Spektralno-elementowe rozwiązanie zależnych od głębokości funkcji jądra w teorii całkowania po liczbach falowych dla propagacji akustyki podwodnej

Słuchając dźwięków pod falami

Dźwięk jest głównym sposobem, w jaki wyczuwamy wydarzenia pod wodą — od śledzenia okrętów podwodnych po nasłuch wielorybów. Jednak przewidywanie, jak dźwięk rozchodzi się w rzeczywistym oceanie, z jego zmiennymi temperaturami i warstwowym dnem, stanowi poważne wyzwanie obliczeniowe. Artykuł przedstawia nowe narzędzie numeryczne, które obiecuje jaśniejsze i szybsze symulacje dźwięku podwodnego, pomagając naukowcom i inżynierom projektować lepszy sonar, łącza komunikacyjne i systemy monitoringu.

Dlaczego dźwięk pod wodą trudno przewidzieć

W oceanie dźwięk nie podróżuje po liniach prostych. Załamuje się, odbija i rozchodzi, przechodząc przez warstwy wody o różnej temperaturze i zasoleniu oraz uderzając w powierzchnię czy dno morskie. Aby przewidzieć, jak dźwięk z statku lub instrumentu się rozproszy, badacze używają modeli matematycznych rozwiązujących równanie falowe — podstawową zasadę zachowania dźwięku. Jedna z potężnych rodzin modeli, zwana całkowaniem po liczbach falowych, rozdziela problem na części poziomą i pionową. Część pionowa, opisująca zmiany dźwięku z głębokością, jest szczególnie trudna i w dużej mierze decyduje o dokładności i szybkości symulacji.

Stare metody i ich kompromisy

Dwie główne podejścia dominowały w tej pionowej kalkulacji. Modele elementów skończonych dzielą kolumnę wodną na wiele cienkich warstw i aproksymują dźwięk w każdej z nich prostymi funkcjami. Są efektywne komputerowo, ale wymagają bardzo drobnego podziału, aby osiągnąć wysoką dokładność, więc ich błąd zmniejsza się tylko powoli wraz ze wzrostem szczegółowości. Modele spektralne idą przeciwną drogą: reprezentują pole dźwiękowe za pomocą gładkich, globalnych kształtów zbudowanych ze specjalnych wielomianów, osiągając bardzo wysoką dokładność przy stosunkowo niewielkiej liczbie niewiadomych. Jednak generują gęste macierze kosztowne do rozwiązania, co spowalnia je dla dużych lub szczegółowych problemów. Do tej pory użytkownicy zwykle musieli wybierać między szybkością a precyzją.

Droga pośrednia z wykorzystaniem elementów spektralnych

Autorzy przedstawiają SemWI, nową wersję modelu całkowania po liczbach falowych, która wykorzystuje metodę elementów spektralnych do obliczeń w pionie. Pomysł polega na podziale kolumny wodnej na elementy, podobnie jak w metodzie elementów skończonych, ale na reprezentowaniu dźwięku wewnątrz każdego elementu za pomocą krzywych wysokiego rzędu zbudowanych z ostro dobranych punktów interpolacyjnych. Punkty te skupiają się przy krawędziach elementów, co poprawia dokładność tam, gdzie cechy dźwięku zmieniają się szybko. Po złożeniu wszystkich elementów powstały układ równań tworzy blokowo diagonalną, symetryczną macierz, znacznie rzadszą niż w standardowych modelach spektralnych. Taka struktura może być rozwiązana szybciej, zachowując jednocześnie szybkie zmniejszanie błędu charakterystyczne dla technik spektralnych.

Testowanie nowej metody

Aby ocenić SemWI, zespół przeprowadził trzy zestawy eksperymentów numerycznych naśladujących typowe warunki oceaniczne. Najpierw przeanalizowano prostą jednolitą kolumnę wodną, której zachowanie jest znane dokładnie z funkcji Airy. SemWI niemal idealnie odtworzył zarówno szczegółowe zależne od głębokości „funkcje jądra”, jak i ogólną stratę transmisji dźwięku, zgodnie z rozwiązaniem dokładnym. Następnie zamodelowano bardziej realistyczny przypadek płytkiej wody z kanałem przy powierzchni nad absorbującym dnem. Tutaj porównano SemWI z dwoma ugruntowanymi programami: SCOOTER, kodem elementów skończonych, oraz WISpec, kodem spektralnym. Wszystkie trzy dały niemal identyczne pola dźwiękowe, włącznie z subtelnym wypływem energii do dna, zarówno dla źródeł punktowych, jak i liniowych. Na koniec rozważono profil oceaniczny głębokiej wody z silnym kanałem dźwiękowym i stwierdzono, że SemWI uchwycił strefy konwergencji w odległościach do 100 kilometrów równie dobrze jak modele wzorcowe.

Balans między szybkością a dokładnością

Ponadto autorzy zbadali zachowanie SemWI przy zmianie ustawień numerycznych. Przy stałej liczbie elementów głębokości i zwiększaniu liczby punktów interpolacyjnych w każdym elemencie zaobserwowali szybki, niemal wykładniczy spadek błędu, odzwierciedlający prawdziwe spektralne zachowanie konkurujące z czystym spektralnym modelem WISpec. Gdy natomiast utrzymano stałą liczbę punktów na element i zwiększono liczbę elementów, SemWI zachowywał się bardziej jak kod elementów skończonych, ale wciąż zbiegał szybciej niż SCOOTER. Testy czasowe wykazały, że SemWI jest zwykle tak szybki lub szybszy niż SCOOTER w przypadkach płytkiej wody, a w wymagającym przypadku głębokiej wody plasuje się między SCOOTER a WISpec, dostarczając przy tym zauważalnie wyższą dokładność niż podejście elementów skończonych.

Co to oznacza dla monitoringu oceanicznego

Mówiąc prosto, praca ta pokazuje, że można uzyskać najlepsze cechy obu podejść przy symulacji dźwięku pod wodą. SemWI oferuje elastyczne nastawienie między szybkością a precyzją przez zmianę liczby elementów i punktów interpolacyjnych, a nawet może naśladować istniejące modele elementów skończonych lub spektralne jako przypadki szczególne. Ponieważ obliczenia głębokości dla różnych próbek liczby falowej można wykonywać równolegle, metoda dobrze pasuje do nowoczesnych wielordzeniowych komputerów. To czyni SemWI praktycznym i potężnym nowym narzędziem dla naukowców i inżynierów, którzy potrzebują wiarygodnych przewidywań rozchodzenia się dźwięku w złożonych warunkach oceanicznych.

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

Słowa kluczowe: akustyka podwodna, propagacja dźwięku, modelowanie numeryczne, metoda elementów spektralnych, przewodnik fal oceanicznych