Spectral-Element-Lösung der tiefenabhängigen Kernfunktionen in der Wellenzahl-Integrationstheorie der Unterwasser-Akustik

Dem Schall unter den Wellen lauschen

Schall ist die wichtigste Quelle, um unter Wasser zu erkennen, was geschieht – von der Ortung von U-Booten bis zum Lauschen auf Wale. Doch vorherzusagen, wie sich Schall durch den realen Ozean bewegt, mit seinen wechselnden Temperaturen und geschichteten Meeresböden, ist eine anspruchsvolle Rechenaufgabe. Dieser Beitrag stellt ein neues numerisches Werkzeug vor, das klarere und schnellere Simulationen der Unterwasser-Schallausbreitung verspricht und Wissenschaftlern sowie Ingenieuren hilft, bessere Sonare, Kommunikationsverbindungen und Überwachungssysteme zu entwickeln.

Warum Unterwasser-Schall schwer vorherzusagen ist

Im Ozean bewegt sich Schall nicht geradlinig. Er wird gebogen, reflektiert und breitet sich aus, wenn er durch Wasser mit unterschiedlichen Temperaturen und Salzgehalten läuft oder auf die Oberfläche bzw. den Meeresboden trifft. Um zu prognostizieren, wie sich Schall von einem Schiff oder einer Messquelle ausbreitet, nutzen Forscher mathematische Modelle, die die Wellengleichung lösen – die grundlegende Regel für das Verhalten von Schall. Eine leistungsfähige Modellfamilie, die Wellenzahl-Integration, zerlegt das Problem in horizontale und vertikale Anteile. Der vertikale Anteil, der beschreibt, wie sich der Schall mit der Tiefe ändert, ist besonders knifflig und bestimmt weitgehend die Genauigkeit und Geschwindigkeit einer Simulation.

Alte Methoden und ihre Kompromisse

Zwei Hauptansätze haben diese vertikale Berechnung dominiert. Finite-Elemente-Modelle teilen die Wassersäule in viele dünne Schichten und approximieren den Schall in jeder Schicht mit einfachen Funktionen. Sie sind rechnerisch effizient, benötigen aber sehr feine Schichtungen, um hohe Genauigkeit zu erreichen, so dass ihre Fehler nur langsam schrumpfen, wenn mehr Details hinzukommen. Spektrale Modelle verfolgen den umgekehrten Ansatz: Sie stellen das Schallfeld mit glatten globalen Formen dar, die aus speziellen Polynomen aufgebaut sind, und erreichen sehr hohe Genauigkeit mit vergleichsweise wenigen Unbekannten. Allerdings erzeugen sie dichte Matrizen, deren Lösung teuer ist, was sie für große oder detaillierte Probleme langsam macht. Bislang mussten Anwender in der Regel zwischen Geschwindigkeit und Präzision wählen.

Ein Mittelweg mit Spektralelementen

Die Autoren stellen SemWI vor, eine neue Variante des Wellenzahl-Integrationsmodells, das die Spektralelementmethode für die Tiefenberechnung nutzt. Die Idee ist, die Wassersäule in Elemente zu teilen, ähnlich wie bei Finite-Elementen, aber den Schall innerhalb jedes Elements mit hochordnigen Kurven darzustellen, die aus sorgfältig gewählten Interpolationspunkten aufgebaut sind. Diese Punkte häufen sich an den Elementrändern, was die Genauigkeit dort verbessert, wo sich Schallmerkmale schnell ändern. Wenn alle Elemente zusammengesetzt sind, ergibt das System eine blockdiagonale, symmetrische Matrix, die deutlich dünner besetzt ist als bei Standard-Spektralmodellen. Diese Struktur lässt sich schneller lösen und bietet dennoch die schnelle Fehlerreduktion, die für spektrale Techniken typisch ist.

Prüfung der neuen Methode

Um SemWI zu bewerten, führte das Team drei Reihen numerischer Experimente durch, die typische Ozeanbedingungen nachbilden. Zuerst untersuchten sie eine einfache Einzelschicht-Wassersäule, deren Verhalten sich genau durch Airy-Funktionen beschreiben lässt. SemWI reproduzierte sowohl die detaillierten tiefenabhängigen „Kern“-Funktionen als auch den gesamten Schalltransmissionsverlust nahezu perfekt und stimmte mit der exakten Lösung überein. Anschließend modellierten sie einen realistischeren Flachwasserfall mit einem Oberflächendukt über einem absorbierenden Meeresboden. Hier verglichen sie SemWI mit zwei etablierten Programmen: SCOOTER, einem Finite-Elemente-Code, und WISpec, einem Spektralcode. Alle drei lieferten nahezu identische Schallfelder, einschließlich subtiler Energieverluste in den Meeresboden, sowohl für Punkt- als auch Linienquellen. Schließlich betrachteten sie ein Tiefseeprofil mit einem starken Schallkanal und fanden, dass SemWI entfernte Konvergenzzonen bis zu 100 Kilometern ebenso gut erfasste wie die Referenzmodelle.

Geschwindigkeit und Genauigkeit ausbalancieren

Über das Abgleichen bekannter Ergebnisse hinaus untersuchten die Autoren, wie sich SemWI beim Anpassen numerischer Einstellungen verhält. Hielt man die Anzahl der Tiefenelemente konstant und erhöhte die Anzahl der Interpolationspunkte pro Element, beobachtete man einen raschen, nahezu exponentiellen Fehlerabfall – ein echtes spektrales Verhalten, das dem reinen Spektralmodell WISpec gleicht. Wenn stattdessen die Punkte pro Element konstant blieben und die Anzahl der Elemente erhöht wurde, ähnelte SemWI eher einem Finite-Elemente-Code, konvergierte jedoch immer noch schneller als SCOOTER. Zeitmessungen zeigten, dass SemWI in Flachwasserfällen meist so schnell oder schneller als SCOOTER ist und in einem anspruchsvollen Tiefwasserfall zwischen SCOOTER und WISpec liegt, dabei jedoch deutlich höhere Genauigkeit als der Finite-Elemente-Ansatz liefert.

Was das für die Ozeanbeobachtung bedeutet

Kurz gesagt zeigt diese Arbeit, dass es möglich ist, das Beste aus beiden Welten bei der Simulation von Unterwasser-Schall zu erreichen. SemWI bietet einen flexiblen Regelbereich zwischen Geschwindigkeit und Präzision durch die Wahl der Anzahl der Elemente und Interpolationspunkte und kann sogar bestehende Finite-Elemente- oder Spektralmodelle als Spezialfälle nachbilden. Da die Tiefenberechnungen für verschiedene Wellenzahlproben parallel ausgeführt werden können, eignet sich die Methode auch gut für moderne Mehrkernrechner. Das macht SemWI zu einem praktischen und leistungsfähigen neuen Werkzeug für Wissenschaftler und Ingenieure, die zuverlässige Vorhersagen darüber benötigen, wie sich Schall in komplexen Ozeanen ausbreitet.

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

Schlüsselwörter: Unterwasserakustik, Schallausbreitung, numerische Modellierung, Spektralelementmethode, Ozean-Waveguide