Solución por elementos espectrales de las funciones núcleo dependientes de la profundidad en la teoría de integración por número de onda de la propagación acústica submarina

Escuchar el sonido bajo las olas

El sonido es la vía principal para percibir lo que ocurre bajo el agua, desde rastrear submarinos hasta escuchar ballenas. Pero predecir cómo se desplaza el sonido en el océano real, con sus temperaturas cambiantes y fondos marinos estratificados, es un reto computacional serio. Este artículo presenta una nueva herramienta numérica que promete simulaciones del sonido submarino más claras y rápidas, ayudando a científicos e ingenieros a diseñar mejor sonares, enlaces de comunicación y sistemas de monitoreo.

Por qué es difícil predecir el sonido submarino

En el océano, el sonido no viaja en línea recta. Se dobla, se refleja y se dispersa al pasar por capas de agua con distintas temperaturas y salinidad, y al chocar con la superficie o el lecho marino. Para pronosticar cómo se extenderá el sonido de un barco o un instrumento, los investigadores usan modelos matemáticos que resuelven la ecuación de ondas, una regla fundamental sobre el comportamiento del sonido. Una familia potente de modelos, llamada integración por número de onda, separa el problema en partes horizontales y verticales. La parte vertical, que describe cómo cambia el sonido con la profundidad, es especialmente compleja y determina en gran medida la precisión y la velocidad de una simulación.

Métodos antiguos y sus compensaciones

Dos enfoques principales han dominado este cálculo vertical. Los modelos de elementos finitos dividen la columna de agua en muchas capas delgadas y aproximan el sonido en cada capa con funciones simples. Son eficientes en ordenador pero requieren un mallado muy fino para alcanzar alta precisión, por lo que sus errores disminuyen solo lentamente al aumentar el detalle. Los modelos espectrales siguen la vía opuesta: representan el campo sonoro usando formas globales suaves construidas con polinomios especiales, alcanzando una precisión extremadamente alta con relativamente pocas incógnitas. Sin embargo, generan matrices densas costosas de resolver, lo que los hace lentos para problemas grandes o muy detallados. Hasta ahora, los usuarios generalmente han debido elegir entre rapidez y precisión.

Un camino intermedio con elementos espectrales

Los autores presentan SemWI, una nueva versión del modelo de integración por número de onda que emplea el método de elementos espectrales para resolver la parte de profundidad. La idea es dividir la columna de agua en elementos, como en elementos finitos, pero representar el sonido dentro de cada elemento con curvas de alto orden construidas a partir de puntos de interpolación cuidadosamente elegidos. Estos puntos se agrupan cerca de los bordes de los elementos, lo que mejora la precisión donde las características sonoras cambian rápidamente. Cuando se ensamblan todos los elementos, el sistema resultante de ecuaciones forma una matriz simétrica por bloques y prácticamente diagonal por bloques, mucho más dispersa que en los modelos espectrales estándar. Esta estructura puede resolverse con mayor rapidez manteniendo la rápida reducción del error típica de las técnicas espectrales.

Poniendo a prueba el nuevo método

Para evaluar SemWI, el equipo realizó tres series de experimentos numéricos que imitan escenarios oceánicos comunes. Primero, examinaron una columna de agua simple de una sola capa cuyo comportamiento se conoce exactamente mediante funciones de Airy. SemWI reprodujo casi a la perfección tanto las detalladas funciones núcleo dependientes de la profundidad como la pérdida de transmisión sonora global, igualando la solución exacta. A continuación, modelaron un caso de agua somera más realista con un conducto superficial sobre un lecho absorbente. Aquí compararon SemWI con dos programas consolidados: SCOOTER, un código de elementos finitos, y WISpec, un código espectral. Los tres produjeron campos sonoros casi idénticos, incluida la sutil fuga de energía hacia el lecho, para fuentes puntuales y lineales. Finalmente, abordaron un perfil de océano profundo con un fuerte canal acústico y encontraron que SemWI capturó zonas de convergencia a distancia hasta 100 kilómetros igual de bien que los modelos de referencia.

Equilibrando velocidad y precisión

Más allá de reproducir resultados conocidos, los autores exploraron cómo se comporta SemWI al ajustar sus parámetros numéricos. Manteniendo fijo el número de elementos en profundidad e incrementando el número de puntos de interpolación dentro de cada elemento, observaron una caída rápida, casi exponencial, del error, reflejo de un comportamiento espectral auténtico que rivaliza con el modelo espectral puro WISpec. Cuando, en cambio, fijaron los puntos por elemento y aumentaron el número de elementos, SemWI se comportó más como un código de elementos finitos pero aún así convergió más rápido que SCOOTER. Pruebas de tiempo indicaron que SemWI suele ser tan rápido o más rápido que SCOOTER en casos de aguas someras, y se sitúa entre SCOOTER y WISpec en un caso exigente de aguas profundas, al tiempo que ofrece una precisión claramente superior al enfoque de elementos finitos.

Qué significa esto para la observación oceánica

En términos sencillos, este trabajo muestra que es posible obtener lo mejor de ambos mundos al simular el sonido submarino. SemWI ofrece un dial flexible entre velocidad y precisión cambiando cuántos elementos y puntos de interpolación se usan, y puede incluso imitar modelos existentes de elementos finitos o espectrales como casos particulares. Dado que los cálculos de profundidad para diferentes muestras de número de onda pueden realizarse en paralelo, el método también es adecuado para los ordenadores modernos con múltiples núcleos. Esto convierte a SemWI en una herramienta práctica y potente para científicos e ingenieros que necesitan predicciones fiables sobre cómo viaja el sonido en océanos complejos.

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

Palabras clave: acústica submarina, propagación del sonido, modelado numérico, método de elementos espectrales, guía de ondas oceánica