Спектрально-элементное решение глубинно-зависимых керновых функций в теории интегрирования по волновому числу для подводного акустического распространения

Слушая звук под волнами

Звук — главный способ узнавать, что происходит под водой: от обнаружения подлодок до прослушивания китов. Но предсказать, как звук распространяется в реальном океане с его изменяющейся температурой и слоистым дном, — серьёзная вычислительная задача. В этой статье представлен новый численный инструмент, который обещает более чёткие и быстрые моделирования подводного звука, помогая учёным и инженерам разрабатывать лучшие сонары, каналы связи и системы мониторинга.

Почему предсказать подводный звук сложно

В океане звук не идёт по прямым линиям. Он преломляется, отражается и рассеивается при прохождении через слои воды с разной температурой и солёностью, а также при взаимодействии с поверхностью и дном. Чтобы предсказать, как звук от корабля или источника распространится, исследователи используют математические модели, решающие волновое уравнение — базовое правило поведения звука. Одна мощная группа моделей, называемая интегрированием по волновому числу, разделяет задачу на горизонтальную и вертикальную части. Вертикальная часть, описывающая изменение звука с глубиной, особенно сложна и во многом определяет точность и скорость моделирования.

Старые методы и их компромиссы

Два основных подхода доминировали в этих глубинных расчётах. Методы конечных элементов разбивают водяной столб на множество тонких слоёв и аппроксимируют звук в каждом слое простыми функциями. Они вычислительно эффективны, но требуют очень тонкой разметки для высокой точности, поэтому погрешность уменьшается медленно по мере увеличения детализации. Спектральные модели идут в противоположную сторону: они представляют звуковое поле гладкими глобальными формами, построенными из специальных полиномов, достигая чрезвычайно высокой точности при относительно небольшом числе неизвестных. Однако они дают плотные матрицы, которые дорого решать, что делает их медленными для больших или детализированных задач. До сих пор пользователям обычно приходилось выбирать между скоростью и точностью.

Промежуточный путь со спектральными элементами

Авторы предлагают SemWI — новую версию модели интегрирования по волновому числу, использующую спектрально-элементный метод для глубинных расчётов. Идея в том, чтобы разделить водяной столб на элементы, как в методе конечных элементов, но представлять звук внутри каждого элемента с помощью высокопорядковых кривых, построенных на специально подобранных интерполяционных точках. Эти точки скапливаются у краёв элементов, что повышает точность там, где звуковые характеристики быстро меняются. После сборки всех элементов система уравнений образует блочно-диагональную симметричную матрицу, гораздо более разреженную, чем в стандартных спектральных моделях. Такая структура решается быстрее, при этом сохраняя быстрое уменьшение ошибки, характерное для спектральных методов.

Тестирование нового метода

Чтобы оценить SemWI, команда провела три набора численных экспериментов, имитирующих типичные океанские условия. Сначала они рассмотрели простой однородный водяной столб, поведение которого известно точно из функций Эйри. SemWI почти идеально воспроизвёл как детальные глубинно-зависимые «керновые» функции, так и общие потери передачи звука, совпав с точным решением. Далее они смоделировали более реалистичный мелководный случай с поверхностным каналом над поглощающим дном. Здесь сравнивали SemWI с двумя устоявшимися программами: SCOOTER (код на методе конечных элементов) и WISpec (спектральный код). Все три инструмента дали почти идентичные звуковые поля, включая тонкое утекание энергии в донные слои, как для точечных, так и для линейных источников. Наконец, они изучили глубоководный профиль с сильным звуковым каналом и показали, что SemWI захватывает зоны схождения на расстояниях до 100 км так же хорошо, как эталонные модели.

Баланс скорости и точности

Помимо соответствия известным результатам, авторы исследовали поведение SemWI при изменении численных параметров. При фиксированном числе глубинных элементов и увеличении числа интерполяционных точек внутри каждого элемента они наблюдали быстрое, почти экспоненциальное падение ошибки, что отражает истинное спектральное поведение и сопоставимо с чисто спектральной моделью WISpec. Когда же фиксировали число точек на элемент и увеличивали число элементов, SemWI вёл себя ближе к методу конечных элементов, но по-прежнему сходился быстрее, чем SCOOTER. Тесты по времени показали, что в мелководных случаях SemWI обычно работает так же быстро или быстрее SCOOTER, а в требовательном глубоководном случае занимает промежуточное положение между SCOOTER и WISpec, обеспечивая при этом заметно большую точность по сравнению с подходом конечных элементов.

Что это значит для океанской разведки

Проще говоря, работа показывает, что при моделировании подводного звука можно получить лучшее из двух миров. SemWI предлагает гибкую регулировку между скоростью и точностью путем изменения числа элементов и интерполяционных точек, и даже может воспроизводить существующие методы конечных или спектральных элементов как частные случаи. Поскольку глубинные расчёты для разных выборок волнового числа можно выполнять параллельно, метод хорошо подходит для современных многопроцессорных компьютеров. Это делает SemWI практичным и мощным новым инструментом для учёных и инженеров, которым нужны надёжные прогнозы распространения звука в сложных океанских условиях.

Цитирование: 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

Ключевые слова: подводная акустика, распространение звука, численное моделирование, спектрально-элементный метод, океанский волновод