Soluzione tramite elementi spettrali delle funzioni kernel dipendenti dalla profondità nella teoria d’integrazione per numero d’onda della propagazione acustica subacquea

Ascoltare il suono sotto le onde

Il suono è il principale mezzo con cui percepiamo ciò che accade sott’acqua, dal tracciamento di sottomarini all’ascolto delle balene. Ma prevedere come il suono si propaga nell’oceano reale, con temperature variabili e fondali stratificati, è una sfida computazionale importante. Questo articolo introduce un nuovo strumento numerico che promette simulazioni del suono subacqueo più chiare e rapide, aiutando scienziati e ingegneri a progettare sonar, collegamenti di comunicazione e sistemi di monitoraggio migliori.

Perché il suono subacqueo è difficile da prevedere

In mare il suono non percorre linee rette. Si piega, si riflette e si disperde mentre attraversa strati d’acqua con temperature e salinità diverse, e quando incontra la superficie o il fondale. Per prevedere come il suono emesso da una nave o da uno strumento si diffonderà, i ricercatori usano modelli matematici che risolvono l’equazione d’onda, la regola fondamentale del comportamento acustico. Una famiglia potente di modelli, detta integrazione per numero d’onda, separa il problema nelle componenti orizzontale e verticale. La parte verticale, che descrive come il suono varia con la profondità, è particolarmente complessa e determina in gran parte l’accuratezza e la velocità della simulazione.

Metodi tradizionali e i loro compromessi

Due approcci principali hanno dominato questo calcolo verticale. I modelli agli elementi finiti suddividono la colonna d’acqua in molti strati sottili e approssimano il campo sonoro in ciascuno con funzioni semplici. Sono efficienti dal punto di vista computazionale ma richiedono una stratificazione molto fine per raggiungere alta accuratezza, perciò gli errori diminuiscono lentamente con l’aumentare della risoluzione. I modelli spettrali seguono la via opposta: rappresentano il campo sonoro usando forme globali lisce costruite con polinomi speciali, raggiungendo elevatissima accuratezza con relativamente poche incognite. Tuttavia generano matrici dense costose da risolvere, rendendoli lenti per problemi grandi o dettagliati. Finora gli utenti dovevano scegliere tra velocità e precisione.

Una via di mezzo con elementi spettrali

Gli autori presentano SemWI, una nuova versione del modello di integrazione per numero d’onda che utilizza il metodo degli elementi spettrali per il calcolo in profondità. L’idea è dividere la colonna d’acqua in elementi, come negli elementi finiti, ma rappresentare il suono all’interno di ciascun elemento con curve di alto ordine costruite a partire da punti di interpolazione accuratamente scelti. Questi punti si concentrano vicino ai bordi degli elementi, migliorando l’accuratezza dove le caratteristiche del suono cambiano rapidamente. Quando tutti gli elementi sono assemblati, il sistema risultante di equazioni forma una matrice simmetrica a blocchi diagonali molto più sparsa rispetto ai modelli spettrali standard. Questa struttura può essere risolta più rapidamente pur mantenendo la rapida riduzione dell’errore tipica delle tecniche spettrali.

Mettere alla prova il nuovo metodo

Per valutare SemWI, il team ha eseguito tre serie di esperimenti numerici che riproducono scenari oceanici comuni. Prima hanno esaminato una semplice colonna d’acqua a strato singolo il cui comportamento è noto esattamente dalle funzioni di Airy. SemWI ha riprodotto sia le dettagliate funzioni “kernel” dipendenti dalla profondità sia la perdita di trasmissione complessiva quasi perfettamente, corrispondendo alla soluzione esatta. Poi hanno modellato un caso di acque costiere più realistico con un duct superficiale sopra un fondale assorbente. Qui hanno confrontato SemWI con due programmi affermati: SCOOTER, un codice a elementi finiti, e WISpec, un codice spettrale. Tutti e tre hanno prodotto campi sonori quasi identici, compresa la sottile perdita di energia nel fondale, sia per sorgenti puntiformi sia per sorgenti lineari. Infine, su un profilo oceanico profondo con un canale acustico marcato, SemWI ha catturato le zone di convergenza a grande distanza fino a 100 chilometri tanto bene quanto i modelli di riferimento.

Bilanciare velocità e accuratezza

Oltre a riprodurre risultati noti, gli autori hanno esplorato il comportamento di SemWI variando i parametri numerici. Mantenendo fisso il numero di elementi in profondità e aumentando il numero di punti di interpolazione all’interno di ciascun elemento, hanno osservato una rapida, quasi esponenziale riduzione dell’errore, rispecchiando il comportamento spettrale che rivaleggia con il modello spettrale puro WISpec. Quando invece hanno fissato i punti per elemento e aumentato il numero di elementi, SemWI si è comportato più come un codice ad elementi finiti ma ha comunque convinto più rapidamente di SCOOTER. I test temporali hanno mostrato che SemWI è di norma veloce quanto o più veloce di SCOOTER nei casi di acque poco profonde, e si colloca tra SCOOTER e WISpec in un caso impegnativo di acque profonde, offrendo al contempo un’accuratezza sensibilmente superiore rispetto all’approccio agli elementi finiti.

Cosa significa per il monitoraggio oceanico

In termini pratici, questo lavoro mostra che è possibile ottenere il meglio di entrambi i mondi nella simulazione del suono subacqueo. SemWI offre una manopola flessibile tra velocità e precisione cambiando il numero di elementi e di punti di interpolazione utilizzati, e può anche riprodurre come casi particolari i modelli a elementi finiti o spettrali esistenti. Poiché i calcoli in profondità per diversi campioni di numero d’onda possono essere eseguiti in parallelo, il metodo è inoltre ben adatto ai computer multicore moderni. Ciò rende SemWI uno strumento pratico e potente per scienziati e ingegneri che necessitano di previsioni affidabili su come il suono si propaga in oceani complessi.

Citazione: Tu, H., Wang, Y., Wang, Y. et al. Spectral element solution of the depth-dependent kernel functions in wavenumber integration theory of underwater acoustic propagation. npj Acoust. 2, 19 (2026). https://doi.org/10.1038/s44384-026-00055-8

Parole chiave: acustica subacquea, propagazione del suono, modellazione numerica, metodo degli elementi spettrali, guidavia oceanica