Solution par éléments spectraux des fonctions noyau dépendant de la profondeur dans la théorie d’intégration en nombre d’ondes de la propagation acoustique sous-marine

Écouter le son sous les vagues

Le son est le principal moyen par lequel nous percevons ce qui se passe sous l’eau, du repérage des sous-marins à l’écoute des baleines. Mais prédire comment le son se propage dans l’océan réel, avec ses températures changeantes et son fond en couches, constitue un défi informatique sérieux. Cet article présente un nouvel outil numérique qui promet des simulations du son sous-marin plus nettes et plus rapides, aidant scientifiques et ingénieurs à concevoir de meilleurs sonars, liaisons de communication et systèmes de surveillance.

Pourquoi le son sous l’eau est difficile à prévoir

Dans l’océan, le son ne voyage pas en lignes droites. Il se courbe, se réfléchit et se diffuse en traversant des couches d’eau aux températures et salinités différentes, et lorsqu’il rencontre la surface ou le fond marin. Pour prévoir comment le son émis par un navire ou un instrument va se propager, les chercheurs utilisent des modèles mathématiques qui résolvent l’équation d’onde, règle fondamentale du comportement acoustique. Une famille puissante de modèles, dite d’intégration en nombre d’ondes, sépare le problème en parties horizontale et verticale. La partie verticale, qui décrit comment le son varie avec la profondeur, est particulièrement délicate et conditionne en grande partie la précision et la rapidité d’une simulation.

Anciennes méthodes et leurs compromis

Deux approches principales ont dominé ce calcul vertical. Les modèles par éléments finis découpent la colonne d’eau en nombreuses couches fines et approchent le champ acoustique dans chaque couche par des fonctions simples. Ils sont efficaces en calcul mais exigent un maillage très fin pour atteindre une grande précision, si bien que leurs erreurs diminuent lentement quand on ajoute du détail. Les modèles spectraux suivent la voie opposée : ils représentent le champ sonore par des formes globales lisses construites à partir de polynômes spéciaux, atteignant une très grande précision avec relativement peu d’inconnues. Cependant, ils génèrent des matrices denses coûteuses à résoudre, ce qui les rend lents pour des problèmes vastes ou détaillés. Jusqu’à présent, il fallait généralement choisir entre rapidité et précision.

Une voie médiane avec les éléments spectraux

Les auteurs présentent SemWI, une nouvelle version du modèle d’intégration en nombre d’ondes qui utilise la méthode des éléments spectraux pour traiter le calcul vertical. L’idée est de diviser la colonne d’eau en éléments, à la manière des éléments finis, mais de représenter le son à l’intérieur de chaque élément par des courbes d’ordre élevé construites à partir de points d’interpolation soigneusement choisis. Ces points se regroupent près des bords d’élément, ce qui améliore la précision là où les caractéristiques acoustiques varient rapidement. Une fois les éléments assemblés, le système d’équations obtenu forme une matrice symétrique par blocs diagonaux beaucoup plus creuse que dans les modèles spectraux standards. Cette structure se résout plus rapidement tout en conservant la réduction d’erreur rapide propre aux techniques spectrales.

Tester la nouvelle méthode

Pour évaluer SemWI, l’équipe a exécuté trois séries d’expériences numériques reproduisant des contextes océaniques courants. D’abord, ils ont examiné une colonne d’eau monolayer simple dont le comportement est connu exactement via les fonctions d’Airy. SemWI a reproduit presque parfaitement à la fois les fonctions « noyau » dépendant de la profondeur et la perte de transmission sonore globale, correspondant à la solution exacte. Ensuite, ils ont modélisé un cas d’eau peu profonde plus réaliste avec un conduit en surface au-dessus d’un fond absorbant. Ici ils ont comparé SemWI à deux programmes établis : SCOOTER, un code d’éléments finis, et WISpec, un code spectral. Les trois ont fourni des champs sonores presque identiques, y compris la fuite d’énergie subtile vers le fond, pour des sources ponctuelles et linéiques. Enfin, ils ont étudié un profil d’océan profond avec un canal acoustique marqué et ont constaté que SemWI capturait les zones de convergence lointaines jusqu’à 100 kilomètres aussi bien que les modèles de référence.

Équilibrer vitesse et précision

Au-delà de la correspondance avec des résultats connus, les auteurs ont exploré le comportement de SemWI en ajustant ses paramètres numériques. En maintenant fixe le nombre d’éléments en profondeur et en augmentant le nombre de points d’interpolation à l’intérieur de chaque élément, ils ont observé une chute d’erreur rapide, quasi exponentielle, reflétant un véritable comportement spectral rivalisant avec le modèle spectral pur WISpec. Lorsqu’ils ont au contraire maintenu constant le nombre de points par élément et augmenté le nombre d’éléments, SemWI a adopté un comportement plus proche d’un code d’éléments finis tout en convergeant plus vite que SCOOTER. Les tests de temps d’exécution ont montré que SemWI est généralement aussi rapide ou plus rapide que SCOOTER dans les cas d’eau peu profonde, et se situe entre SCOOTER et WISpec pour un cas exigeant en eau profonde, tout en offrant une précision notablement supérieure à l’approche par éléments finis.

Ce que cela signifie pour la surveillance océanique

En termes simples, ce travail montre qu’il est possible d’obtenir le meilleur des deux mondes pour simuler le son sous-marin. SemWI offre un réglage flexible entre vitesse et précision en modifiant le nombre d’éléments et de points d’interpolation, et peut même reproduire les modèles d’éléments finis ou spectraux existants comme cas particuliers. Parce que les calculs en profondeur pour différents échantillons en nombre d’ondes peuvent être effectués en parallèle, la méthode est aussi bien adaptée aux ordinateurs multicœurs modernes. Cela fait de SemWI un nouvel outil pratique et puissant pour les scientifiques et ingénieurs qui ont besoin de prédictions fiables sur la propagation du son dans des océans complexes.

Citation: 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

Mots-clés: acoustique sous-marine, propagation du son, modélisation numérique, méthode des éléments spectraux, guide d’ondes océanique