🔢 mathematics

Cavitation Acoustic Perturbation Equations: A Computational Framework for Source-Resolved Multiphase Hydroacoustics

Este artículo presenta las Ecuaciones de Perturbación Acústica de Cavitación (CAPE), un marco computacional unificado que integra directamente la física de la cavitación en el modelado acústico para predecir, localizar y analizar con precisión las fuentes hidroacústicas y su propagación en flujos multifásicos con cavitación.

Autores originales: Zhi Cheng, Rajeev K. Jaiman

Publicado 2026-07-23

📖 7 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Zhi Cheng, Rajeev K. Jaiman

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Imagina el océano como una gigantesca orquesta invisible. A veces, la música es un zumbido suave, pero otras veces, es un rugido ensordecedor que puede escucharse a kilómetros de distancia. No se trata solo de las canciones de las ballenas; se trata de las máquinas que construimos para movernos a través del agua, como las hélices de los barcos y las turbinas submarinas. Cuando estas máquinas giran demasiado rápido o golpean el ángulo incorrecto, crean un extraño gemido de alta frecuencia conocido como "canto". Esto sucede debido a un fenómeno escurridizo llamado cavitación. Piensa en la cavitación como agua hirviendo, pero sin el calor. Cuando el agua se mueve lo suficientemente rápido, la presión baja tanto que literalmente hierve a temperatura ambiente, convirtiéndose en diminutas burbujas de vapor. Estas burbujas son inestables; crecen durante una fracción de segundo y luego colapsan violentamente, estallando como micro-petardos. Cada estallido crea una pequeña onda de choque de sonido. Cuando miles de millones de estas burbujas estallan en un ritmo sincronizado, crean un tono fuerte y molesto que puede dañar las piezas de los barcos, revelar la ubicación de un submarino ante los enemigos o dañar la vida marina.

Durante décadas, los científicos han intentado predecir y comprender este ruido utilizando matemáticas complejas. Los métodos tradicionales eran como intentar escuchar un susurro en medio de un huracán escuchando únicamente el viento golpeando las paredes de una habitación; eran excelentes para escuchar el ruido causado por la forma del objeto (como la pala de una hélice), pero pasaban por alto el ruido proveniente de las propias burbujas estallando dentro del agua. Este artículo presenta una nueva forma más inteligente de escuchar. Los investigadores construyeron un marco de trabajo "consciente de la cavitación" llamado CAPE (Ecuaciones de Perturbación Acústica de Cavitación). En lugar de ignorar las burbujas o tratarlas como una simple nota al pie, esta nueva matemática trata la creación y el estallido de las burbujas de vapor como una fuente de sonido directa y fuerte justo donde ocurre. Es como actualizar un micrófono que solo capta el viento a uno que también puede escuchar los diminutos petardos estallando en medio de la habitación.

El Nuevo Detective de Sonidos

El logro central de este trabajo es la creación de una herramienta computacional que puede simular exactamente cómo la cavitación genera ruido, llegando incluso a los tonos específicos de "canto" que atormentan a los ingenieros marinos. Los autores, Zhi Cheng y Rajeev K. Jaiman de la Universidad de British Columbia, desarrollaron un sistema que acopla la física del flujo de agua con la física de las ondas sonoras en una única ecuación unificada.

En el pasado, los científicos tenían que elegir entre dos caminos difíciles. Podían ejecutar una simulación extremadamente costosa y lenta que intentara calcular cada molécula de aire y de agua moviéndose y estallando (lo cual es como intentar contar cada grano de arena en una playa para predecir una ola). O bien, podían usar un método de atajo que asumía que el agua era simplemente una forma sólida y calculaba el ruido basándose en cómo el agua empujaba contra el objeto (como adivinar el sonido de un tambor mirando únicamente el baquet). El atajo omitía por completo el "ruido de las burbujas".

Este nuevo marco CAPE cierra esa brecha. Toma el flujo del agua (el "estado base") y añade una capa de "perturbaciones acústicas" (las ondas sonoras) encima de ello. Crucialmente, añade un término especial a la matemática que da cuenta de la transferencia de masa. En lenguaje sencillo, esto significa que la matemática rastrea explícitamente cuándo el agua se convierte en vapor y cuándo el vapor vuelve a convertirse en agua. Debido a que convertir un líquido en un gas cambia el volumen de ese espacio instantáneamente, actúa como un altavoz gigante e invisible que empuja el aire (o el agua) hacia afuera. El artículo muestra que este "cambio de volumen" es la razón principal por la cual la cavitación crea un sonido de "monopolo" —un tipo de ruido que se irradia de manera igual en todas las direcciones, como un globo estallando—. Esto es diferente al sonido de "dipolo" de un objeto sin cavitación, que es más direccional, como un altavoz apuntando en una dirección.

Los Experimentos Virtuales

Para demostrar que su nueva matemática funciona, el equipo no solo escribió ecuaciones; realizaron una serie de experimentos virtuales.

Primero, probaron el sistema en una línea unidimensional simple. Enviaron una onda sonora a lo largo de esta línea y comprobaron si su computadora podía manejarla sin que la onda rebotara de forma extraña o perdiera su forma. Descubrieron que su técnica de "Capa de Coincidencia Perfecta" (PML, por sus siglas en inglés) funcionaba como una espuma absorbente de sonido, tragándose las ondas en el borde de la simulación para que no rebotaran y arruinaran los resultados. También comprobaron si el sonido se extinguía a la velocidad correcta, coincidiendo con una regla clásica de la física llamada ley de Stokes, que predice cómo el sonido se desvanece en un fluido debido a la fricción. La simulación coincidió perfectamente con esta regla, demostrando que la matemática era físicamente sólida.

Luego, pasaron a una prueba más compleja: el flujo de agua pasando por un cilindro circular (un poste redondo).

Sin Cavitación: Cuando el agua fluía suavemente, el patrón de ruido era un "dipolo". Era fuerte en los lados del poste y silencioso al frente y detrás, causado por los vórtices giratorios (remolinos) que se desprendían de la parte posterior. Esto coincidía con lo que los científicos ya sabían.
Con Cavitación: Cuando activaron la configuración de "cavitación", el patrón de ruido cambió drásticamente. La simulación mostró que el ruido se convirtió en un "monopolo", irradiándose casi por igual en todas las direcciones. El tono de "canto" apareció, y los investigadores pudieron señalar exactamente de dónde provenía.

Descubrieron que el "canto" no era causado por los grandes remolinos giratorios en la estela, como alguien podría suponer. En cambio, era causado por eventos de colapso localizados. Imagina una nube de burbujas formándose en la superficie del poste, creciendo y luego implosionando repentinamente. Esta implosión crea un pico de presión masivo y rápido (un aumento de la "tasa de presión"). El artículo muestra que estos estallidos repentinos y violentos son los principales motores del ruido tonal fuerte. La matemática logró separar el "respiro" de toda la nube de burbujas (un ritmo lento) del "estallido" violento de las burbujas individuales (un ritmo rápido y fuerte).

Finalmente, probaron el sistema en un hidroala NACA 6412 (una forma similar a la pala de una hélice de barco). Aquí, los resultados fueron aún más interesantes. Debido a que las burbujas solo se formaban en el "lado de succión" (la curva trasera) de la pala, el sonido no se irradiaba uniformemente. El hidroala actuaba como un escudo, bloqueando parte del sonido en una dirección y permitiéndolo pasar en la otra. La simulación capturó esta asimetría, mostrando que el ruido era más fuerte cerca de la zona de colapso de las burbujas y más silencioso en el lado opuesto. Esto confirmó que el nuevo marco podía manejar no solo formas simples, sino geometrías complejas del mundo real donde la fuente del sonido está oculta y la forma del objeto importa.

Lo Que Esto Significa

El artículo concluye que este nuevo marco CAPE es una herramienta poderosa y eficiente para comprender el ruido submarino. No requiere la potencia de cómputo imposible de simular cada molécula, pero tampoco ignora las burbujas. Al tratar el cambio de fase (agua a vapor y viceversa) como una fuente directa de sonido, permite a los ingenieros ver por qué una hélice está cantando y de dónde proviene el ruido.

Los autores advierten cuidadosamente que estos resultados se basan en simulaciones de escenarios específicos y controlados (como un cilindro y un hidroala en un flujo recto). Si bien la matemática es robusta y ha sido verificada contra leyes físicas conocidas, el artículo sugiere que el trabajo futuro deberá probar esto en tres dimensiones con hélices rotatorias y turbulencia del mundo real. Sin embargo, por ahora, este marco ofrece una forma clara y computacionalmente eficiente de resolver el misterio del "canto" inducido por la cavitación, convirtiendo un complejo problema acústico en algo que puede ser visualizado, medido y, eventualmente, tal vez, silenciado.

Resumen Técnico: Ecuaciones de Perturbación Acústica por Cavitación (CAPE)

Planteamiento del Problema

El ruido hidrodinámico proveniente de flujos turbulentos y cavitantes plantea desafíos significativos para los sistemas marinos e hidráulicos, particularmente en lo que respecta a la integridad estructural, el impacto ambiental y el sigilo acústico. Un fenómeno específico y problemático es el "canto" (singing): un ruido tonal de banda estrecha y fuerte generado por la cavitación que puede propagarse a largas distancias. Aunque la Computación de Ruido Directo (resolver las ecuaciones de Navier-Stokes totalmente compresibles) es teóricamente completa, resulta computacionalmente prohibitiva para flujos cavitantes de bajo número de Mach debido a la vasta disparidad entre las escalas hidrodinámicas y acústicas.

Los métodos híbridos existentes enfrentan limitaciones:

Analogías de Ffowcs Williams-Hawkings (FW-H): Aunque son eficientes para la predicción en campo lejano, las formulaciones clásicas tienen dificultades para representar con precisión las fuentes monopolares inducidas por la cavitación que surgen del cambio de fase volumétrico (generación y colapso de vapor) dentro del campo de flujo.
Ecuaciones de Perturbación Acústica (APE): Estas permiten la resolución de fuentes en el dominio, pero históricamente se han desarrollado para flujos de una sola fase o débilmente compresibles. Típicamente, omiten la transferencia de masa entre fases y las variaciones de compresibilidad de la mezcla, fallando en capturar matemáticamente las contribuciones acústicas dominantes de la vaporización rápida, la condensación y el colapso de la cavidad.

Consecuentemente, existe una brecha entre las simulaciones de flujo cavitante de alta fidelidad y los modelos acústicos capaces de acoplar consistentemente la dinámica de la cavitación con la propagación de ondas acústicas para explicar el fenómeno del canto tonal.

Metodología

Los autores proponen un marco de Ecuaciones de Perturbación Acústica por Cavitación (CAPE), el cual extiende el enfoque convencional de APE a mezclas homogéneas de líquido-vapor.

Formulación Matemática

Descomposición: El flujo se descompone en un estado base hidrodinámico incompresible de variación lenta (resuelto mediante URANS con un modelo de cavitación) y pequeñas perturbaciones acústicas compresibles.
Ecuaciones de Gobierno: El sistema CAPE se deriva insertando las variables descompuestas en las ecuaciones de Navier-Stokes compresibles. Crucialmente, la formulación incorpora explícitamente:
- Compresibilidad de la Mezcla: La velocidad del sonido en la mezcla se determina mediante una relación de tipo Wood dependiente de la fracción de volumen de vapor local, lo que permite una reducción significativa de la velocidad del sonido en las regiones cavitantes.
- Términos Fuente de Transferencia de Masa: Las ecuaciones de continuidad y cantidad de movimiento incluyen términos fuente ( $h$ ) derivados del modelo de cavitación de Schnerr-Sauer, que representan la tasa neta de transferencia de masa ( $\dot{m}$ ) entre las fases líquida y de vapor.
- Amplificación de la Tasa de Presión: La formulación captura la derivada temporal de la presión ( $\partial P/\partial t$ ) como un motor para la generación acústica.
Balance de Energía Acústica: Un análisis de energía revela que la transferencia de masa inducida por el cambio de fase actúa como una fuente monopolar volumétrica, distinta de las fuentes dipolares asociadas con la carga superficial inestable.

Implementación Numérica

Solucionador: Se emplea un esquema de predictor-corrector de volumen finito. El flujo base incompresible se resuelve utilizando el algoritmo PIMPLE, mientras que el sistema CAPE utiliza un enfoque segregado con subciclos acústicos.
Estrategia de Acoplamiento: Se utiliza un acoplamiento unidireccional donde la solución hidrodinámica impulsa el campo acústico a través de términos fuente, pero se desprecia la retroalimentación acústica hacia el flujo (justificado por las pequeñas amplitudes de presión acústica en relación con las escalas hidrodinámicas).
Paso de Tiempo: Se utiliza un enfoque de paso de tiempo múltiple donde el solucionador acústico avanza con un paso de tiempo $\Delta t = \Delta t_f / m$ (donde $m=1000$ ) relativo al paso de tiempo hidrodinámico $\Delta t_f$ . Los términos fuente se interpolan linealmente entre los pasos de tiempo hidrodinámicos.
Condiciones de Contorno: Se implementa una Capa de Coincidencia Perfecta (PML) para absorber las ondas acústicas salientes y suprimir las reflexiones de frontera, complementada con condiciones estándar no reflectantes.

Contribuciones Clave

Derivación de CAPE: Una derivación rigurosa de las ecuaciones de perturbación acústica para flujos cavitantes que integra directamente la transferencia de masa por cambio de fase y la compresibilidad de la mezcla en las ecuaciones de control acústico.
Identificación de la Fuente Monopolar: La formulación proporciona un vínculo matemático directo que muestra que la transferencia de masa inducida por la cavitación actúa como una fuente monopolar volumétrica, aclarando el mecanismo físico detrás del ruido de cavitación.
Marco Computacional: Implementación de un solucionador de volumen finito estable y eficiente con subciclos acústicos y tratamiento de PML, capaz de resolver los componentes tonales inducidos por la cavitación dentro del dominio computacional.

Resultados y Verificación

Verificación Numérica

Propagación de Ondas 1D: El solucionador fue verificado contra problemas de propagación de ondas 1D. Los resultados demostraron convergencia de malla y de paso de tiempo, supresión efectiva de reflexiones de frontera mediante PML y preservación precisa de las frecuencias dominantes.
Validación de la Ley de Stokes: El solucionador reprodujo con éxito la atenuación acústica dependiente de la frecuencia predicha por la ley de Stokes ( $\alpha \propto f^2$ ), confirmando la implementación correcta de los términos de disipación viscosa.

Aplicación a Flujos Canónicos

El marco se aplicó a dos configuraciones: flujo alrededor de un cilindro circular y flujo alrededor de un perfil hidrodinámico NACA 6412.

1. Cilindro Circular (Re = 1.11 × 10⁵):

Caso No Cavitante: Exhibió un patrón de radiación de tipo dipolar asociado con la desprendimiento de vórtices inestable, consistente con la literatura establecida.
Caso Cavitante:
- Patrón de Radiación: Cambió a un patrón de radiación casi axisimétrico de tipo monopolar, impulsado por efectos de cambio de fase volumétricos.
- Análisis Espectral: Identificó tres frecuencias dominantes: un modo de "respiración" de baja frecuencia de la cavidad, la frecuencia de desprendimiento de vórtices y una frecuencia más alta asociada con el desprendimiento/colapso de la cavidad.
- Mecanismo de Canto: El estudio vinculó el "canto" (ruido tonal) con eventos de colapso de cavidad localizados cerca de la pared. Estos eventos generan fluctuaciones intensas y localizadas de la tasa de presión ( $\partial P/\partial t$ ) y transferencia de masa, que actúan como fuertes impulsores acústicos. Los picos tonales en el espectro acústico se alinearon con la frecuencia de desprendimiento de la cavidad ( $f_c$ ) y sus armónicos, en lugar de la frecuencia de desprendimiento de vórtices.
- Localización de la Fuente: A diferencia de los flujos no cavitantes donde las fuentes de ruido están en la estela, las fuentes de ruido de cavitación se encontraron localizadas cerca de la superficie del cilindro durante los eventos de colapso.

2. Perfil Hidrodinámico NACA:

Asimetría: Debido a la geometría, el campo acústico exhibió asimetría espacial. Los niveles de presión sonora (SPL) más altos se concentraron en el lado de succión donde ocurre la cavitación.
Apantallamiento: La geometría del perfil bloqueó y reflejó parcialmente las ondas acústicas, creando un patrón de radiación sesgado distinto de la simetría monopolar del cilindro.
Consistencia del Mecanismo: El caso del perfil confirmó que los eventos de colapso localizados en el lado de succión generan la misma amplificación de la tasa de presión y respuesta acústica tonal observada en el caso del cilindro.

Significado y Reivindicaciones

El artículo afirma que el marco CAPE propuesto extiende con éxito los métodos de perturbación acústica a sistemas multifásicos cavitantes. Su principal significado radica en:

Resolución Unificada: Permite la resolución directa de las fuentes de ruido inducidas por la cavitación (tanto el monopolo volumétrico como el dipolo de carga) dentro del dominio computacional, evitando las limitaciones de los métodos de integral de superficie como FW-H para flujos cavitantes.
Información sobre el Mecanismo: Proporciona una herramienta cuantitativa para analizar el vínculo físico entre la dinámica de la cavitación (específicamente la amplificación de la tasa de presión por colapso) y la radiación acústica, ofreciendo una explicación clara del surgimiento del "canto" tonal.
Eficiencia Computacional: Al utilizar un acoplamiento unidireccional con flujos base incompresibles, ofrece una alternativa computacionalmente eficiente para la predicción hidroacústica en aplicaciones marinas e hidráulicas.

Los autores concluyen que esta metodología proporciona un marco numérico práctico para la predicción hidroacústica, la localización de fuentes y el análisis de mecanismos, con planes futuros para extender la formulación a flujos turbulentos tridimensionales, propulsores rotatorios e interacciones fluido-estructura.