Impact of Reservoir Inflow Velocity on Flow Dynamics and Turbine Design in Overtopping Breakwater Energy Conversion Systems

Resumen Técnico: Impacto de la Velocidad de Entrada al Reservorio en la Dinámica de Flujo y el Diseño de Turbinas en Sistemas de Conversión de Energía de Rompeolas con Sobrepase (OBREC)

Planteamiento del Problema
Los sistemas de Rompeolas para Conversión de Energía (OBREC, por sus siglas en inglés) integran la protección costera con la captura de energía de las olas, recolectando el agua que sobrepasa la estructura hacia un reservorio elevado para accionar turbinas de baja caída. Si bien la literatura existente caracteriza extensamente el rendimiento de los OBREC utilizando la carga del reservorio y el caudal de sobrepase, estas métricas son insuficientes para definir las condiciones hidráulicas instantáneas en la entrada de la turbina. La entrada al reservorio ocurre típicamente como chorros altamente transitorios, eventos de salpicadura y pulsos de corta duración. En consecuencia, las evaluaciones basadas únicamente en el caudal integrado pueden pasar por alto características críticas de la entrada dinámica—específicamente, la velocidad de entrada al reservorio ($V_{res}$)—que influyen directamente en la eficiencia, la carga y las estrategias de control de la turbina. Además, la influencia de la geometría de la rampa (lineal vs. convexa) sobre estas condiciones de velocidad transitoria permanece insuficientemente investigada.

Metodología
Este estudio emplea la Hidrodinámica de Partículas Suavizadas (SPH), un método numérico basado en partículas libre de malla, para simular la propagación de las olas y los procesos de sobrepase dentro de un sistema OBREC. Las simulaciones se realizaron utilizando el software de código abierto DualSPHysics con el núcleo Wendland para asegurar la estabilidad y un modelado preciso de la superficie libre.

• Verificación del Modelo: El modelo SPH fue verificado mediante soluciones analíticas derivadas de la teoría de olas de segundo orden. Se seleccionó un espaciamiento de partículas de 0.03 m basándose en el error de Raíz Cuadrática Media (RMS) más bajo en comparación con los perfiles de olas analíticos.
• Configuración de la Simulación: Se modeló un canal de olas numérico (20 m de longitud, 0.8 m de profundidad de agua en calma) con dos geometrías de rampa distintas: Lineal y Convexa. Ambas configuraciones compartieron una pendiente de rampa frontal de 23°, un francobordo de la cresta ($R_r$) de 0.025 m y un ancho de reservorio ($Br$) de 0.135 m.
• Condiciones de Ola: Se generaron olas regulares con alturas de ola significativa ($H_s$) de 0.15 m y 0.18 m, y periodos de ola pico ($T_p$) que oscilaron entre 1.59 s y 2.23 s. Estas condiciones de modelo corresponden a condiciones de prototipo de 3.0 m de altura de ola y 8.05 s de periodo (escala 1:20).
• Marco de Análisis: Se aplicó un análisis no dimensional basado en la gravedad. La velocidad de entrada al reservorio ($V_{res}$) se normalizó por la escala de celeridad de la ola ($\sqrt{gH}$) y se analizó en función de la pendiente de la ola ($H/L_m$), la relación relativa de profundidad-a-altura de ola ($h/H$) y la relación relativa de profundidad-a-longitud de onda ($h/L_m$).

Resultos Clave

1. Influencia Geométrica en el Sobrepase: La geometría de rampa Lineal promovió consistentemente un mayor ascenso de la ola (run-up) y eventos de sobrepase más frecuentes en comparación con la geometría Convexa. La rampa Convexa, caracterizada por un camino de flujo un 15% más largo, indujo una mayor disipación de energía mediante un mayor desarrollo de la capa límite y turbulencia, suprimiendo a menudo el sobrepase por completo bajo condiciones de olas de menor energía (por ejemplo, $H_s = 0.15$ m).
2. Características de la Velocidad de Entrada:
   • Geometría Lineal: Produjo velocidades de entrada medibles que oscilaron aproximadamente entre 0.399 m/s y 0.591 m/s a través de las condiciones probadas.
   • Geometría Convexa: Exhibió cero sobrepase en varios casos (particularmente en alturas de ola menores). Cuando ocurrió el sobrepase (en energías de ola más altas), las velocidades fueron generalmente menores que las de su contraparte Lineal (por ejemplo, de 0.281 m/s a 0.512 m/s).
3. Tendencias No Dimensionales:
   • Pendiente de la Ola ($H/L_m$): La velocidad de entrada normalizada ($V_{res}/\sqrt{gH}$) exhibió una tendencia decreciente a medida que la pendiente de la ola aumentaba. Las olas más empinadas provocaron una rotura temprana y una mayor disipación de energía a lo largo de la rampa, reduciendo la transferencia de momento al reservorio.
   • Profundidad Relativa ($h/H$ y $h/L_m$): Relaciones más altas de profundidad-a-altura de ola ($h/H$) se correlacionaron con velocidades de entrada normalizadas más altas, sugiriendo que una mayor profundidad relativa a la altura de la ola ayuda a preservar el momento. Por el contrario, el aumento de la relación relativa de profundidad-a-longitud de onda ($h/L_m$) resultó en una disminución de la velocidad de entrada normalizada.

Contribuciones Clave

• La Velocidad como Variable Primaria: El estudio trata explícitamente la velocidad de entrada al reservorio ($V_{res}$) como una variable de respuesta hidráulica crítica, yendo más allá de las evaluaciones tradicionales centradas en el caudal para abordar las condiciones de contorno específicas requeridas para el diseño de turbinas.
• Sensibilidad Geométrica: Cuantifica cómo la geometría de la rampa (Lineal vs. Convexa) altera no solo el volumen de agua que sobrepasa, sino también la velocidad y el momento de la entrada, destacando que caudales similares pueden generar condiciones de entrada de turbina sumamente distintas.
• Escalamiento No Dimensional: La investigación establece un marco no dimensional que vincula $V_{res}$ con la pendiente de la ola y los parámetros de profundidad relativa, proporcionando un método escalable para predecir el comportamiento de la entrada en diversos estados de mar.

Significancia y Afirmaciones
Los autores sostienen que tratar $V_{res}$ como un parámetro de diseño explícito es esencial para el diseño hidráulico de sistemas OBREC integrados con turbinas. Mientras que las métricas de caudal estiman la energía potencial disponible, no capturan las condiciones de forzamiento no estacionarias que dictan los umbrales de activación de la turbina, los envolventes de eficiencia y la confiabilidad operativa.

El estudio concluye que la geometría Lineal ofrece un rango operativo más amplio, manteniendo la entrada bajo condiciones de olas de menor energía donde las geometrías Convexas pueden fallar en generar sobrepase. Por el contrario, las geometrías Convexas concentran la operación en estados de mayor energía. Por lo tanto, la velocidad de entrada normalizada ($V_{res}/\sqrt{gH}$) debe considerarse un parámetro hidráulico gobernante junto con la carga y el caudal. Este enfoque cierra la breancia entre la hidrodinámica del sobrepase y el diseño del sistema de toma de potencia (PTO), permitiendo una selección de turbina y una optimización del rendimiento más consistentes desde el punto de vista físico. Los autores señalan que, si bien estos hallazgos aclaran las dependencias sistemáticas, se requiere más investigación para generalizar estos resultados en un espacio de diseño más amplio, incluyendo perfiles de rampa y profundidades de agua variadas.