🔢 matemáticas

SIREN: A nonlinear Rankine-source boundary element solver for the stationary motion of a marine vehicle in calm water

SIREN es un resolvedor de elementos de frontera de fuente de Rankine no lineal que mejora los métodos tradicionales mediante el uso de un problema de Neumann actualizado iterativamente con una superficie libre de tapa rígida y un marcado dinámico de paneles húmedos/secos para calcular con precisión la resistencia de onda estacionaria de vehículos marinos en aguas tranquilas.

Autores originales: Daniele Peri

Publicado 2026-10-14

📖 7 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Daniele Peri

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

Cada barco que atraviesa aguas tranquilas deja una estela, una ondulación permanente de energía por la que el navío debe pagar en forma de resistencia por arrastre. Esta resistencia por formación de olas es uno de los enigmas más antiguos y persistentes en la ingeniería naval. Durante más de un siglo, los científicos han dependido de modelos matemáticos para predecir cuánta potencia necesita un barco para superar estas olas, con la esperanza de diseñar buques que se deslicen de manera más eficiente. El desafío central reside en el agua misma: a medida que un barco se desplaza, empuja la superficie hacia arriba y hacia abajo, creando una frontera compleja y cambiante entre el casco sólido y el aire fluido superior. Los métodos tradicionales a menudo han simplificado esta interacción, tratando la superficie del agua como un plano plano e invariable o asumiendo que el área mojada del barco termina exactamente donde el agua está tranquila. Si bien estos atajos hicieron posibles los cálculos, introdujeron errores, particularmente cerca de la proa, donde las olas son más altas y las fuerzas de presión son más críticas.

Una nueva herramienta computacional llamada SIREN, desarrollada por el investigador Daniele Peri, aborda este problema abandonando esas simplificaciones. En lugar de obligar al agua a permanecer plana o al barco a detenerse en una línea fija, este software permite que la superficie del agua suba y baje naturalmente alrededor de la forma real del casco. Trata a todo el barco, desde la quilla hasta la cubierta, como un único objeto sólido al que el agua no puede penetrar, independientemente de si una parte específica del casco está actualmente sumergida o expuesta al aire. Al resolver la física del flujo sobre esta geometría real y móvil, el método captura las interacciones sutiles que los modelos anteriores pasan por alto. El resultado es una predicción más precisa de cuánta resistencia enfrentará un barco, ofreciendo un camino más claro hacia la comprensión del verdadero costo de desplazarse a través del agua.

El corazón de este nuevo enfoque es un cambio en la forma en que la computadora visualiza el agua alrededor del barco. En los métodos antiguos, el software calculaba el flujo de agua alrededor de un "doble cuerpo"—una forma teórica creada al reflejar el barco a través de una línea de flotación plana. Esto funcionaba lo suficientemente bien para casos simples, pero fallaba cuando la proa del barco se elevaba fuera del agua o cuando las olas subían lo suficiente como para tocar la cubierta. SIREN evita esto resolviendo el problema sobre la forma real y actual del barco y la superficie del agua tal como existen en ese momento. Imagine la superficie del agua como una tapa rígida que se sitúa exactamente donde las olas la han empujado; la computadora calcula cómo fluye el agua bajo esta tapa y alrededor del casco, tratando a todo el barco como una barrera sólida. Esto elimina la necesidad de trucos matemáticos complejos para manejar la frontera entre el agua y el aire, permitiendo que la simulación gestione las partes húmedas y secas del casco sin problemas.

Para encontrar la respuesta final, el software no solo calcula una vez; itera, o repite, el proceso para refinar el resultado. Comienza con una suposición inicial de la altura de la ola, calcula el flujo y luego mueve la superficie del agua virtual hacia arriba o hacia abajo para coincidir con la nueva altura de la ola que acaba de calcular. Luego, recalcula el flujo alrededor de esta nueva forma. Este ciclo continúa, con la computadora ajustando lentamente el perfil de la ola hasta que la solución se estabiliza. Una innovación clave en este proceso es cómo el software maneja las partes del barco que están temporalmente secas. En versiones anteriores, la computadora eliminaba por completo los paneles secos de sus cálculos, lo que causaba inestabilidad y errores cuando esos paneles volvían a mojarse. El nuevo método mantiene cada parte del barco en el cálculo en todo momento, marcándolas como húmedas o secas solo para decidir qué partes contribuyen a la medición de la fuerza final. Este simple cambio evita los errores matemáticos que causaban que las simulaciones anteriores fallaran o produjeran resultados poco fiables.

Los investigadores probaron este nuevo solver en un casco de prueba estándar conocido como casco Wigley, una forma simple y simétrica que se utiliza a menudo para evaluar teorías hidrodinámicas. Realizaron simulaciones a través de una gama de velocidades, representadas por un número llamado número de Froude, que compara la velocidad del barco con el tamaño de las olas que genera. Las pruebas cubrieron desde velocidades muy lentas hasta bastante rápidas, utilizando tres niveles diferentes de detalle en la rejilla computacional, desde una gruesa hasta una muy fina. Los resultados mostraron que el nuevo método era notablemente estable. Completó con éxito los complejos cálculos repetitivos para casi todas las velocidades y tamaños de rejilla probados, una mejora significativa respecto a intentos anteriores que a menudo no lograban converger en ciertas velocidades.

Cuando los investigadores compararon los nuevos resultados con datos físicos de tanques de remolque, donde se arrastran modelos de barcos reales a través del agua, el nuevo método mostró mejoras claras. En el rango de velocidad más crítico, donde la resistencia por formación de olas alcanza su punto máximo, la versión no lineal del solver —que tiene en cuenta la superficie de agua en movimiento— predijo valores de resistencia más cercanos a las mediciones del mundo real que los métodos lineales más antiguos. El nuevo enfoque predijo correctamente que la resistencia aumentaría a medida que las olas crecieran, un comportamiento físico que los modelos antiguos a veces erraban. Sin embargo, el estudio también destacó dónde el método aún tiene límites. A las velocidades más altas probadas, los resultados comenzaron a variar dependiendo de qué tan detallada fuera la rejilla computacional, lo que sugiere que incluso se podrían necesitar rejillas más finas para capturar los comportamientos de las olas más extremos.

El estudio también analizó qué tan bien coincidía la resistencia por formación de olas calculada con verificaciones independientes basadas en la energía de las olas que quedan muy atrás del barco. Para la mayoría de las velocidades probadas, la presión calculada en el casco del barco concordó bien con la energía medida en la estela, confirmando que la simulación es físicamente sólida. Sin embargo, a las velocidades más altas, este acuerdo se volvió menos consistente, indicando que la resolución de la rejilla actual podría no ser suficiente para capturar la complejidad total de las olas en esas condiciones. Los investigadores señalaron que, si bien el método es un gran paso adelante, aún no es una solución final para todos los tipos de barcos. Aún no ha sido probado en embarcaciones complejas de doble casco o en barcos con transomas planos en la popa, y el software actualmente asume que el ángulo del barco en el agua es fijo en lugar de permitir que se hunda o se incline de forma natural.

A pesar de estas interrogantes abiertas, el trabajo demuestra una nueva forma robusta de simular la hidrodinámica de los barcos. Al tratar la superficie del agua como una frontera dinámica y móvil que interactúa con todo el casco, el solver SIREN proporciona una representación más fiel de la realidad que los modelos anteriores. Resuelve con éxito los problemas de inestabilidad que plagaban las simulaciones no lineales anteriores y ofrece una herramienta más precisa para predecir la resistencia por formación de olas. Aunque se requiere un mayor refinamiento para manejar las velocidades más extremas y las geometrías de barcos más complejas, este enfoque proporciona una base sólida para futuros desarrollos en el diseño de barcos, prometiendo una comprensión más clara de las fuerzas invisibles que dan forma a cómo los navíos se desplazan por el mar.

Resumen Técnico: SIREN – Un Solucionador de Elementos de Contorno de Fuente de Rankine No Lineal

Planteamiento del Problema
El artículo aborda el cálculo de la resistencia de formación de olas en régimen estacionario para un vehículo marino que avanza en aguas tranquilas. Mientras que el marco clásico establecido por Michell y Havelock se basa en la teoría de naves delgadas, y los métodos de paneles subsiguientes (Dawson, 1977) mejoraron la satisfacción de las condiciones de contorno, tres aproximaciones específicas en el esquema clásico de fuentes de Rankine limitan la precisión:

El flujo base se linealiza sobre un flujo de "doble cuerpo" asumiendo una superficie libre plana en $z=0$ .
Los paneles de la superficie libre están fijados en $z=0$ .
Se asume que la superficie de la carena termina exactamente en el plano de agua quieta ( $z=0$ ).

Esta tercera aproximación es particularmente problemática porque la integral de presión para la resistencia de olas es altamente sensible a la región de proa, donde las crestas de las olas suelen elevarse por encima de $z=0$ . Además, el flujo de doble cuerpo clásico requiere un sistema de imágenes que carece de sentido físico cuando el casco se extiende por encima de $z=0$ , lo que requiere paneles desingularizados para evitar singularidades numéricas.

Metodología
El artículo presenta SIREN (Ship Inviscid Resistance and Elevation eNgine), un solucionador de paneles de fuentes de Rankine que relaja las tres aproximaciones clásicas mediante un esquema de iteración no lineal específico.

Formulación del Flujo Base: En lugar de un flujo de doble cuerpo, el flujo base se resuelve como un problema de Neumann sobre la geometría actual, real. Las fuentes se distribuyen sobre todo el casco (tanto paneles húmedos como secos) y los paneles de la superficie libre. Se impone la impermeabilidad en todo el casco y en la superficie libre (actuando como una tapa rígida en su posición actual). Esto elimina la necesidad de sistemas de imágenes y de paneles desingularizados, permitiendo que la carena se extienda por encima de $z=0$ sin solaparse con su propia imagen.
Iteración No Lineal: La condición de la superficie libre se linealiza sobre este flujo base. La iteración procede de la siguiente manera:
1. Se computa el flujo base.
2. Los paneles de la superficie libre se mueven rígidamente en la dirección vertical hacia la elevación de la ola computada en sus centroides.
3. Los paneles del casco se marcan como "húmedos" o "secos" basándose en el nuevo nivel de la superficie. Crucialmente, la impermeabilidad se impone en todos los paneles del casco independientemente de esta marca; la marca se utiliza únicamente para seleccionar qué paneles contribuyen a la integral de fuerza.
4. Se vuelve a computar el flujo base sobre la geometría movida.
5. Se aplica la sub-relajación solo en las pasadas donde la corrección no lineal se acumula, evitando la inestabilidad observada en configuraciones anteriores.
Implementación Numérica: El solucionador utiliza una distribución de fuerza de fuente constante en paneles cuadriláteros. La condición de radiación se impone mediante un desplazamiento hacia arriba (upstream) de los puntos de colocación (escalonamiento). Los sistemas lineales se resuelven mediante una eliminación gaussiana incompleta o una factorización LU de OpenBLAS con refinamiento iterativo.
Evaluación de Fuerza: La resistencia de olas se calcula mediante la integración de presión en campo cercano y, de forma independiente, mediante un análisis de corte de onda transversal en campo lejano (método de Eggers/Sharma). La resistencia por fricción se estima utilizando ITTC-1957 y correlaciones de placa plana panel por panel.

Contribuciones Clave

Flujo Base Dependiente de la Geometría: La principal contribución es un flujo base definido sobre la geometría real y actual (incluyendo secciones secas del casco y paneles de superficie libre movidos) en lugar de una tapa plana fija. Esto elimina la necesidad de sistemas de imágenes y desingularización.
Manejo Robusto de Húmedo/Seco: Al imponer la impermeabilidad en todo el casco (húmedo y seco) en cada iteración y utilizar la marca húmedo/seco solo para la integración de la fuerza, el método evita defectos estructurales en el sistema lineal que anteriormente causaban picos de presión espurios en los paneles de proa recién mojados.
Iteración No Lineal Estable: La restricción específica de la sub-relajación a las pasadas de acumulación, combinada con la impermeabilidad de todo el casco, estabiliza la iteración no lineal en un amplio rango de números de Froude ( $Fn = 0.150–0.450$ ), un régimen donde las configuraciones previas del proyecto divergían.
Marco de Verificación: El artículo implementa un riguroso estudio de convergencia de malla utilizando el procedimiento de Richardson de tres mallas y compara los resultados contra la teoría de naves delgadas de Michell y datos de tanques de remolque.

Resultados
El solucionador fue probado en un casco tipo Wigley con tres densidades de malla de superficie libre ( $Fn = 0.150–0.450$ ).

Estabilidad: La iteración no lineal fue estable en todas las mallas para todos los números de Froude, excepto para el más bajo ( $Fn=0.150$ ) en la malla más fina.
Convergencia de Malla: La convergencia de la malla fue estrecha tanto para las soluciones lineales como para las no lineales en la mayoría de las velocidades. Sin embargo, en los dos números de Froude más altos ( $Fn=0.408, 0.450$ ) en la malla más fina, los resultados no lineales fueron inconsistentes con la malla, mostrando una desviación significativa respecto a las mallas más gruesas.
Comparación con Experimentos: La corrección no lineal generalmente aumentó el coeficiente de resistencia de olas ( $C_w$ ), lo cual es el signo físicamente esperado. En la principal joroba de resistencia ( $Fn \approx 0.316–0.354$ ), la solución no lineal mejoró la concordancia con los datos de los tanques de remolque en comparación con la solución lineal. Sin embargo, en los extremos del rango, la corrección no lineal no mejoró la concordancia de manera uniforme.
Perfiles de Ola: La solución no lineal siguió mejor los valles de las olas y los ascensos de popa medidos que la solución lineal, pero subestimó consistentemente la altura de la cresta de proa, probablemente debido a las limitaciones de resolución de la malla de la superficie libre en la proa.
Consistencia de Campo Lejano: La estimación del flujo de energía (campo lejano) concordó bien con los resultados de la integración de presión (campo cercano) en el rango de velocidad de diseño, validando la solución de campo cercano. El método de corte de onda transversal mostró un desfase multiplicativo consistente, y el corte de onda longitudinal resultó ser una herramienta de diagnóstico para errores de implementación más que un estimador de resistencia fiable en su forma actual.

Significado y Limitaciones
El artículo afirma que la formulación propuesta logra estabilizar la iteración de la fuente de Rankine no lineal y proporciona un tratamiento físicamente consistente de la superficie mojada sin sistemas de imágenes. Los resultados demuestran que la corrección no lineal mejora las predicciones de la resistencia en la región crítica de la joroba de resistencia.

Sin embargo, el autor establece explícitamente que estos resultados son preliminares. Se señalan varias limitaciones y puntos abiertos:

La tolerancia iterativa no se alcanzó completamente dentro del número de pasadas utilizadas.
La inconsistencia de la malla persiste en los números de Froude más altos para la solución no lineal.
La base de consistencia entre campo cercano y campo lejano (relación cercana a la unidad) establecida en formulaciones anteriores no se ha restablecido con el nuevo tratamiento de húmedo/seco.
El solucionador aún no ha sido probado en catamaranes, cascos de popa de transón, o con hundimiento y trim gratuito (aunque la infraestructura de código para esto existe).
Los paneles de la superficie libre satisfacen la condición en sus centroides movidos rígidamente, no en la verdadera superficie inclinada.

El trabajo se presenta como un paso fundacional hacia un solucionador de Rankine-source no lineal más robusto, identificando como siguientes pasos necesarios las extensiones de hundimiento/trim libre y formas de casco complejas (transón, casco doble) antes de que los hallazgos puedan generalizarse más allá de los monocascos de actitud fija.

Más como este