🔬 physics

Potential flow reconstruction of plunging breaking wave measured with PIV

Este estudio combina mediciones experimentales de Velocimetría de Imágenes de Partículas (PIV) con el modelado de flujo potencial totalmente no lineal (FNPF) para investigar olas de rotura de tipo plunging, revelando que el criterio de inicio de rotura cuasi universal de $B=0.85$ se cumple solo después de que la ola ya ha comenzado a volcarse.

Autores originales: Jeffrey Harris, Lili Kimmoun, Sunil Mohanlal, Jensen McTighe, Jason Dahl, Stephan Grilli

Publicado 2026-06-29

📖 7 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Jeffrey Harris, Lili Kimmoun, Sunil Mohanlal, Jensen McTighe, Jason Dahl, Stephan Grilli

Artículo original bajo licencia CC BY 4.0 (https://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

La visión general: Capturar una ola antes de que rompa

Imagine que está tratando de predecir exactamente cuándo romperá una ola gigante en el océano. Para los ingenieros que construyen plataformas petrolíferas en alta mar o defensas costeras, saber cuándo y con qué fuerza golpea una ola es crucial. Si se equivoca, la estructura podría sufrir daños.

Los científicos han intentado durante mucho tiempo crear modelos informáticos para simular estas olas. Sin embargo, una vez que una ola comienza a romper y convertirse en espuma blanca, la física se vuelve desordenada y caótica. La mayoría de los modelos informáticos tienen dificultades con esta parte "desordenada".

Este artículo trata sobre un equipo de investigadores que quería probar si su modelo informático súper preciso podía predecir el momento exacto en que una ola decide romper, y si esa predicción coincide con lo que ven en un tanque de olas real.

El experimento: Una "máquina de olas" en un laboratorio

Los investigadores construyeron una piscina de agua larga y estrecha (un canal de olas) en una universidad en Francia. Piense en ello como una bañera gigante y de alta tecnología.

El generador de olas: En un extremo, utilizaron un enorme remo que oscila de un lado a otro. Al mover el remo con un ritmo muy específico y complejo (como un DJ mezclando pistas), pudieron forzar al agua a crear una única ola "enfocada" y masiva.
El objetivo: Querían que esta ola creciera más y más alta hasta que alcanzara un punto específico en el medio de la piscina y luego rompiera (se desplomara) hacia adelante, tal como una ola para surfistas.
La medición: No se limitaron a observar; utilizaron dos herramientas:
1. Sondas de oleaje: Como reglas que sobresalen del agua para medir la altura.
2. PIV (Velocimetría de Imágenes de Partículas): Esta es la parte genial. Añadieron partículas diminutas e inofensivas al agua y proyectaron una lámina de láser a través de ella. Una cámara de alta velocidad tomó miles de fotos por segundo. Al rastrear cómo se movían las partículas, pudieron mapear la velocidad y la dirección del agua dentro de la ola, llegando justo hasta la parte superior de la cresta.

El modelo informático: El simulador "perfecto"

En el lado informático, utilizaron un modelo de "Flujo Potencial Total No Lineal" (FNPF, por sus siglas en inglés).

La analogía: Imagine el agua como una hoja de seda perfectamente lisa y sin fricción. En el mundo real, el agua es pegajosa y desordenada (viscosa). En este modelo informático, fingen que el agua es una seda perfecta.
El desafío: Las olas reales rompen debido a esa viscosidad y caos. Dado que el modelo informático ignora la viscosidad, no puede simular naturalmente el "choque" y la espuma. Si simplemente dejaran correr el modelo, la ola seguiría curvándose eternamente sin llegar a salpicar nunca.
La solución: Para hacer que el modelo actúe como una ola real que rompe, los investigadores añadieron un "freno". Cuando el ordenador detectaba que una ola estaba a punto de romper, aplicaba una "presión de absorción" especial (como una esponja) a la parte superior de la ola para drenar su energía, imitando la pérdida de energía de un choque real.

El gran descubrimiento: El "número mágico"

Durante años, los científicos han debatido una regla específica para predecir cuándo romperá una ola. La regla se basa en una relación llamada B.

La relación: Compara qué tan rápido se mueve el agua en la parte superior de la ola ( $u_c$ ) frente a qué tan rápido se desplaza la forma de la ola en sí misma ( $c_c$ ).
La teoría: Si el agua en la parte superior se mueve más rápido que el 85% de la velocidad de la ola ( $B = 0.85$ ), la ola está destinada a romper.

Lo que los investigadores encontraron:

El tiempo es ajustado: Descubrieron que la ola realmente comienza a volcarse (overturn) casi en el mismo momento en que se alcanza este ratio de 0.85. Es como un dominó cayendo; en el momento en que el dominó alcanza la marca de 0.85, ya se está volcando.
El inconveniente: Debido a que la ola comienza a volcarse tan rápidamente después de alcanzar el 0.85, es extremadamente difícil para los modelos informáticos capturar este momento exacto. Si el modelo se retrasa aunque sea un poco, la ola ya ha roto y el modelo de "seda perfecta" deja de funcionar.
La coincidencia de datos: A pesar de la dificultad, su modelo informático coincidió increíblemente bien con las mediciones de la cámara y el láser de la vida real. La velocidad del agua dentro de la ola, medida por el láser, coincidió casi perfectamente con la predicción del ordenador, incluso cerca de la parte superior de la ola.

La pérdida de energía: ¿Con qué fuerza golpea?

El artículo también analizó cuánta energía se pierde cuando la ola rompe.

La analogía: Piense en un accidente de coche. La "fuerza de rotura" es como la fuerza con la que el coche golpea la pared.
El hallazgo: Intentaron usar una fórmula que otros científicos habían creado para predecir esta fuerza de choque. Sin embargo, debido a que su ola era tan intensa y enfocada (como un choque de coche a alta velocidad), la fórmula antigua no funcionó del todo bien. Tuvieron que ajustar el "frenado" en su modelo informático para que coincidiera con la pérdida de energía del mundo real. Descubrieron que, para estas olas específicas e intensas, la disipación de energía seguía un patrón diferente al observado en olas más suaves.

El "truco de magia" de la simplificación

Una de las partes más ingeniosas del artículo es cómo analizaron los datos.

El problema: El modelo informático genera millones de puntos de datos, lo cual es difícil de comparar con las mediciones del láser.
La solución: Los investigadores se dieron cuenta de que el movimiento complejo del agua cerca de la cresta podía describirse mediante solo tres "polos" invisibles (puntos matemáticos) flotando en el aire sobre el agua.
La analogía: Imagine intentar describir la forma de una nube. En lugar de mapear cada gota de agua, podría simplemente decir: "Hay tres imanes invisibles tirando de la nube para darle esta forma". Descubrieron que usar solo estos tres "imanes" matemáticos era suficiente para recrear perfectamente la velocidad y la forma de la ola al romper. Esta es una simplificación enorme que hace que entender la ola sea mucho más fácil.

Resumen

Los investigadores demostraron con éxito que:

Su modelo informático es lo suficientemente preciso como para predecir el momento exacto en que una ola comienza a romper.
La "regla del 0.85" (donde la velocidad del agua alcanza el 85% de la velocidad de la ola) es un indicador muy fiable de que una ola está a punto de romper.
Aunque el agua se vuelve desordenada cuando rompe, el flujo justo antes del choque sigue reglas simples y predecibles que pueden describirse con muy pocos "polos" matemáticos.

Han puesto todos sus datos y mediciones a disposición de otros científicos para que sirvan como un estándar de referencia ("gold standard") para probar sus propios modelos de olas.

Resumen Técnico: Reconstrucción del Flujo Potencial de una Ola de Rotura de Tipo Plunging Medida con PIV

Planteamiento del Problema
El modelado preciso de las olas de rotura es crítico para el diseño de estructuras marinas y la ingeniería costera, pero sigue siendo un desafío debido a la compleja física de la disipación de energía y el vuelco (overturning). Si bien los modelos de Dinámica de Fluidos Computacional (CFD) pueden simular la rotura, son computacionalmente costosos y a menudo están restringidos a escalas pequeñas. En consecuencia, se utilizan ampliamente modelos de orden reducido, como los modelos de Flujo Potencial Completamente No Lineal (FNPF), pero estos dependen de parametrizaciones aproximadas para la disipación de energía y el inicio de la rotura. Se ha propuesto un criterio cinemático específico, $B = u_c/c_c \approx 0.85$ (donde $u_c$ es la velocidad del fluido en la cresta y $c_c$ es la celeridad de la cresta), como un precursor universal de la rotura. Sin embargo, existe una falta de evidencia experimental de alta resolución que valide este criterio y la parametrización de la disipación de energía asociada, particularmente respecto a la cinemática precisa cerca de la cresta justo antes y durante el vuelco. Además, medir estas velocidades experimentalmente es difícil debido a la complejidad de la superficie libre y las limitaciones de la Velocimetría de Imágenes de Partículas (PIV) cerca de la cresta.

Metodología
El estudio emplea un enfoque dual que combina experimentos de laboratorio de precisión con simulaciones numéricas de alta fidelidad:

Configuración Experimental: Los experimentos se llevaron a cabo en el canal de olas de precisión de l'École Centrale de Marseille (ECM). Un generador de olas de tipo paleta generó olas de rotura de tipo plunging focalizadas utilizando un espectro de densidad de amplitud de Ricker modificado para lograr el enfoque de frecuencia. La profundidad del agua se fijó en $h = 0.667$ m.
- Mediciones: La elevación de la superficie fue registrada por nueve sensores de oleaje (wave gauges). La cinemática del flujo en la región de la cresta se midió mediante una cámara de video de alta velocidad y PIV con partículas trazadoras.
- Repetibilidad: Se realizaron múltiples pruebas para asegurar la repetibilidad, mostrando las mediciones de los sensores de oleaje diferencias menores a 0.5 mm en el punto de rotura.
Modelo Numérico: Se utilizó un modelo de Flujo Potencial Completamente No Lineal (FNPF) 2D basado en el Método de Elementos de Frontera (BEM) para simular los experimentos.
- Ecuaciones de Gobierno: El modelo resuelve la ecuación de Laplace para el potencial de velocidad con condiciones de contorno cinemáticas y dinámicas completamente no lineales en la superficie libre, utilizando un marco de trabajo de Eulerian-Lagrangiano mixto.
- Inicio de la Rotura y Disipación: El modelo incorpora un mecanismo de detección de rotura basado en el criterio $B = u_c/c_c$ . Una vez alcanzado el umbral $B_t = 0.85$ , se aplica una presión superficial "absorbente" para simular la disipación de energía. La fuerza de disipación fue parametrizada utilizando un parámetro de fuerza de rotura no dimensional constante $b$ , calibrado con datos experimentales.
- Reconstrucción de la Velocidad: Para comparar los resultados numéricos con los datos de PIV (que a menudo presentan vacíos o ruido cerca de la cresta), los autores desarrollaron una representación analítica de orden reducido del campo de velocidad cercano a la cresta. Este método aproxima el flujo mediante una superposición de polos complejos discretos (singularidades) ubicados fuera del dominio del fluido, ajustados a las velocidades de la superficie del BEM.

Contribuciones Clave

Validación de Alta Resolución: El estudio proporciona una comparación directa entre los resultados del modelo FNPF y las mediciones de precisión de PIV de una ola de rotura de tipo plunging, centrándose específicamente en la cinemática cerca de la cresta en el inicio de la rotura.
Validación del Criterio $B$ : La investigación investiga el tiempo del criterio $B = 0.85$ en relación con el vuelco físico de la ola.
Modelado de Velocidad de Orden Reducido: Se introduce un nuevo enfoque analítico utilizando polos complejos para reconstruir el campo de velocidad interno cerca de la cresta de vuelco. Este método permite una comparación precisa con los datos experimentales sin el costo computacional de evaluar integrales singulares en cada punto y maneja la geometría compleja de la superficie de vuelco.
Parametrización de la Disipación: El estudio evalúa diferentes formulaciones para el parámetro de fuerza de rotura ( $\gamma$ y $b$ ), encontrando que para estas olas de tipo plunging focalizadas, una parametrización basada en Romero et al. (2022) proporcionó un mejor acuerdo con la pérdida de energía experimental que la formulación de Derakhti et al. (2019).

Resultados

Precisión del Modelo: El modelo FNPF-BEM reprodujo con precisión la evolución de la superficie libre y la fase de las olas hasta el punto de vuelco. Los resultados numéricos mostraron un excelente acuerdo con las series temporales de los sensores de oleaje y los campos de velocidad derivados de PIV, tanto antes como después del evento de rotura (cuando la disipación estaba activa).
Tiempo de Inicio de la Rotura: El análisis reveló que el criterio $B = 0.85$ se satisface casi simultáneamente con el momento en que el frente de la ola se vuelve vertical y comienza el vuelco. La diferencia de tiempo entre alcanzar $B=0.85$ y que el frente se vuelva vertical fue mínima (aprox. 0.02 s), lo que sugiere que para olas de rotura de tipo plunging rápidas, el criterio puede alcanzarse demasiado tarde para servir como un precursor predictivo en algunos esquemas numéricos sin una disipación inmediata.
Comparación de Velocidad: El modelo de polo analítico reconstruyó con éxito el campo de velocidad. Las comparaciones entre el modelo de polo, la solución completa del BEM y las mediciones de PIV mostraron una fuerte concordancia para las velocidades horizontales y verticales en ubicaciones ligeramente por debajo de la superficie libre. Se observaron discrepancias muy cerca de la cresta, consistentes con las limitaciones de medición de PIV (ruido, reflexión).
Calibración de la Disipación: El análisis de sensibilidad indicó que un parámetro de fuerza de rotura constante $b \approx 0.0431$ (derivado de la relación empírica de Romero et al.) produjo el mejor ajuste con las alturas de ola experimentales aguas abajo del rompedor. El estudio señaló que la alta intensidad de la ola focalizada resultó en un parámetro de fuerza de rotura $\gamma$ que estaba fuera del rango de los conjuntos de datos previos utilizados para derivar leyes de disipación alternativas.

Significancia
El artículo concluye que el conjunto de datos generado a partir de estos experimentos, caracterizado por una alta repetibilidad y mediciones cinemáticas detalladas, sirve como un referente experimental robusto para validar criterios de inicio de rotura y modelos de disipación de energía en simulaciones de olas de orden reducido. El estudio confirma que los modelos FNPF, cuando se acoplan con mecanismos de disipación adecuados y el criterio $B$ , pueden simular con precisión las etapas de pre-rotura y el inicio del vuelco de las olas de tipo plunging. La representación analítica de polos propuesta ofrece una herramienta prometedora para el modelado de orden reducido, capaz de capturar el comportamiento localizado y casi singular de las olas de rotura con un pequeño número de parámetros. Los autores señalan modestamente que, si bien la representación simplificada captura bien la cinemática, aún no proporciona una evaluación detallada de los criterios de rotura en sí, dejando esto para trabajos futuros. El conjunto de datos se pone a disposición como material suplementario para facilitar la investigación futura en este dominio.