Leveraging unstructured grids for direct numerical simulations of wall turbulence

Este trabajo presenta la {\eta}-malla, un marco de generación de mallas no estructuradas para simulaciones numéricas directas de turbulencia de pared que escala los tamaños de malla con la escala local de Kolmogorov, logrando una precisión comparable a la de las mallas cartesianas convencionales mientras reduce significativamente el costo computacional, particularmente a altos números de Reynolds y sobre geometrías complejas de riblets.

Lo esencial

Imagina intentar simular cómo fluye el aire sobre un automóvil o cómo se mueve el agua alrededor del casco de un barco. Para lograrlo con precisión en una computadora, los científicos utilizan una técnica llamada Simulación Numérica Directa (DNS). Piensa en la DNS como crear un microscopio digital gigante en 3D que divide el fluido (aire o agua) en millones de cubos diminutos e invisibles (una cuadrícula). Luego, la computadora calcula cómo se mueve e interactúa cada cubo individual con sus vecinos.

El problema es que el flujo de fluido cerca de una superficie (como el costado de un barco) es increíblemente caótico y detallado. Para obtener una imagen clara, necesitas una cantidad masiva de estos cubos diminutos justo al lado de la superficie. Sin embargo, a medida que te alejas de la superficie, el caos se suaviza y ya no necesitas cubos tan pequeños.

La Vieja Forma: El "Muro de Ladrillos Rígidos"

Tradicionalmente, los científicos utilizaban una cuadrícula cartesiana. Imagina construir un muro con ladrillos idénticos y rígidos.

• El Problema: Para ver los detalles diminutos cerca de la superficie, tienes que hacer los ladrillos de la parte inferior muy pequeños. Pero como estos ladrillos son rígidos y están conectados en una línea recta, se te obliga a usar esos mismos ladrillos diminutos hasta la parte superior de tu muro, incluso donde los detalles no son importantes.
• El Resultado: Terminas con un muro hecho de miles de millones de ladrillos diminutos, la mayoría de los cuales son innecesarios. Esto hace que la simulación por computadora sea increíblemente lenta y costosa, como intentar contar cada grano de arena en una playa solo para medir la marea.

La Nueva Solución: La "Red Inteligente y Elástica"

Este artículo introduce un nuevo método llamado $\eta$-grid (red eta). En lugar de ladrillos rígidos, imagina una red de pesca inteligente y elástica.

• Cómo funciona: Los autores diseñaron un sistema donde el tamaño de los agujeros de la red cambia automáticamente según la cantidad de detalle necesaria.
  • Cerca de la superficie (La "Capa Interna"): La red tiene agujeros muy pequeños y ajustados para atrapar los remolinos caóticos y diminutos del fluido.
  • Más lejos (La "Capa Externa"): A medida que el fluido se calma, la red se estira automáticamente, haciendo los agujeros mucho más grandes.
• El Ingrediente Secreto: El tamaño de estos agujeros se basa en algo llamado escala de Kolmogorov (denotada como $\eta$). Piensa en esto como el "remolino más pequeño posible" que puede existir en el fluido a cualquier altura dada. La nueva cuadrícula simplemente dice: "Haz el tamaño del agujero lo suficientemente grande como para atrapar el remolino más pequeño a esta altura específica, y no más grande".

Por Qué Esto es Importante

Los autores probaron esta "red inteligente" en dos tipos diferentes de códigos informáticos (uno similar al método de elementos espectrales, el otro similar al método de volúmenes finitos) y lo compararon con el antiguo método de "ladrillos rígidos".

1. Precisión: Los resultados fueron casi idénticos. La "red inteligente" capturó la física tan bien como los "ladrillos rígidos", con menos de un 1% de diferencia en mediciones clave como la fricción y la velocidad.
2. Ahorros Masivos: Aquí es donde ocurre la magia.
   • Para superficies lisas (como una pared plana), la nueva cuadrícula redujo la cantidad de "ladrillos" necesarios (puntos de la cuadrícula) en aproximadamente un 90% a altas velocidades.
   • Para superficies rugosas (como una pared con pequeñas ranuras llamadas "riblets" diseñadas para reducir la resistencia), los ahorros fueron aún más dramáticos: hasta un 97% menos de puntos de cuadrícula.

La Analogía de la "Pared Surcada"

Para entender la parte de los riblets, imagina una pared cubierta de pequeñas ranuras paralelas (como la textura de una pelota de golf o la piel de un tiburón).

• La Vieja Forma: Para simular esto, el método de ladrillos rígidos tenía que mantener los ladrillos diminutos en todas partes porque las ranuras obligaban a que la cuadrícula fuera fina hasta arriba. Era como intentar contar cada hilo de un suéter, incluso las partes lejos de la tela.
• La Nueva Forma: La "red inteligente" sabe que una vez que estás a unos centímetros por encima de esas pequeñas ranuras, el flujo se vuelve suave nuevamente. Estira los agujeros inmediatamente por encima de las ranuras, ignorando los detalles diminutos que ya no importan.

La Conclusión

Los autores han creado un marco que actúa como un zoom inteligente para las simulaciones de fluidos. Concentra su poder computacional exactamente donde se necesita (cerca de la pared) y se relaja donde no es necesario.

• Para Paredes Lisas: Escala el esfuerzo mucho más lentamente a medida que la simulación se hace más grande.
• Para Paredes Rugosas: Escala aún mejor, haciendo factibles en computadoras que anteriormente no podían manejarlas las simulaciones de superficies complejas que reducen la resistencia.

En resumen, encontraron una manera de realizar el mismo trabajo de alta calidad con una fracción de la potencia de computación, convirtiendo una tarea que podría llevar a un superordenador un mes en una que podría realizarse en unos pocos días.

Resumen técnico

Resumen Técnico: Aprovechamiento de Mallas No Estructuradas para Simulaciones Numéricas Directas de Turbulencia de Pared

Planteamiento del Problema
La Simulación Numérica Directa (DNS) de la turbulencia de pared es computacionalmente costosa, principalmente debido a la necesidad de resolver la escala de longitud de Kolmogorov ($\eta$) en todo el dominio. Los solucionadores de mallas cartesianas convencionales, ampliamente utilizados para superficies lisas y no lisas (por ejemplo, mediante Métodos de Frontera Inmersa), suelen imponer una resolución de malla fija o de variación lenta en las direcciones streamwise ($\Delta x^+$) y spanwise ($\Delta z^+$) a lo largo de todo el dominio. Este enfoque conduce a una sobre-resolución significativa en las regiones externas donde la escala de Kolmogorov crece, creando un cuello de botella computacional. Esto es particularmente agudo en simulaciones que involucran superficies complejas como riblets, donde la malla debe ser extremadamente fina cerca de la superficie para resolver micro-grooves, sin embargo, esta resolución fina se mantiene innecesariamente a lo largo de todo el espesor de la capa límite en marcos cartesianos.

Metodología
Los autores proponen un nuevo marco de generación de mallas no estructuradas denominado $\eta$-grid. El principio fundamental consiste en establecer los tamaños de malla en la dirección normal a la pared ($\Delta y^+$) y spanwise ($\Delta z^+$) proporcionales a la escala local de Kolmogorov ($\eta^+$), permitiendo el ensanchamiento de la malla en las regiones de flujo externo donde las escalas de turbulencia son mayores.

1. Formulación de la Malla:

   • Capa Interna: Una capa delgada ($\sim 50$ unidades viscosas) adyacente a la pared utiliza tamaños de malla escalados viscosamente similares a la DNS convencional ($0.3 \lesssim \Delta y^+ \lesssim 4$, $\Delta z^+ \simeq 5$). Para superficies no lisas (por ejemplo, riblets), la malla resuelve la longitud de onda superficial más pequeña ($\ell^+$) con $\ell^+/30 \lesssim \Delta y^+, \Delta z^+ \lesssim 4$.
   • Región Externa: Por encima de la capa interna, los tamaños de malla aumentan proporcionalmente a la escala local de Kolmogorov ($\Delta y^+ \simeq \Delta z^+ \simeq 2\eta^+$).
   • Leyes de Escala: Los autores refinan ajustes semi-empíricos para $\eta^+$ tanto en la región logarítmica como en la región de estela para el flujo turbulento en canal y las capas límite turbulentas (TBL) con gradiente de presión cero (ZPG). Proponen ajustes de ley de potencia por tramos para $\eta^+$ que transicionan en $y^+ \approx \delta^+/2$.
   • Adaptación a Riblets: Para geometrías de riblets, el marco introduce una "subcapa de riblet" ($\delta_r^+$) llena de celdas cuadradas para resolver flujos secundarios y rodillos de Kelvin-Helmholtz, transicionando suavemente hacia la escala basada en $\eta$ por encima.

2. Implementación:

   • El marco utiliza el paquete de código abierto Gmsh para generar mallas no estructuradas multi-bloque.
   • Las mallas se prueban utilizando dos solucionadores distintos: SOD2D (Método de Elementos Espectrales) y OpenFOAM (Método de Volúmenes Finitos).
   • Las simulaciones abarcan flujo turbulento en canal y TBLs ZPG sobre paredes lisas y diversas geometrías de riblets (surcos triangulares) hasta un número de Reynolds de fricción $\delta_0^+ = 1000$.

Contribuciones Clave

• Desarrollo del Marco: Una estrategia robusta de generación de mallas no estructuradas multi-bloque que adapta dinámicamente el espaciado de la malla a la escala local de Kolmogorov, aplicable tanto a superficies lisas como a superficies rugosas complejas.
• Leyes de Escala Refinadas: Ajustes semi-empíricos actualizados para la escala de Kolmogorov en la región de estela de TBLs y flujos en canal, mejorando la precisión del dimensionamiento de la malla en la capa externa.
• Agnosticismo del Solucionador: Demostración de que el marco $\eta$-grid produce resultados consistentes y de alta fidelidad a través de diferentes esquemas de discretización numérica (SEM y FVM).
• Formulación Específica para Riblets: Una estrategia de malla adaptada para riblets que resuelve la física de la subcapa (flujos secundarios, rodillos KH) mientras permite un ensanchamiento rápido en el flujo externo, superando las limitaciones de los enfoques cartesianos/IBM.

Resultados

• Precisión: El $\eta$-grid produce resultados con una diferencia menor al 1% en comparación con mallas cartesianas convencionales y datos de referencia de DNS/experimentales. Las métricas incluyen coeficientes de fricción superficial ($C_f$), perfiles de velocidad media ($U^+$), tensiones turbulentas y espectrogramas.
• Eficiencia de la Malla:
  • Paredes Lisas: El número de puntos de malla ($N_\eta$) escala como $\propto \delta_0^{+2.5}$, en comparación con $\propto \delta_0^{+3.0}$ para mallas cartesianas con estiramiento de tangente hiperbólica. En $\delta_0^+ = 6000$, el $\eta$-grid requiere solo $\sim 10\%$ de los puntos de malla de una malla cartesiana estándar ($N_\eta/N_{Tanh} \simeq 0.1$).
  • Riblets: La ganancia en eficiencia es aún más pronunciada. Para riblets triangulares que reducen la resistencia, $N_\eta$ escala como $\propto \delta_0^{+2.0}$, mientras que las mallas cartesianas escalan como $\propto \delta_0^{+3.0}$. En $\delta_0^+ = 6000$, el $\eta$-grid requiere solo $\sim 3\%$ de los puntos de malla de una malla cartesiana ($N_\eta/N_{Tanh} \simeq 0.03$).
• Física de Riblets: El marco captura con éxito mecanismos de flujo complejos sobre riblets, incluyendo el origen virtual, flujos secundarios y rodillos de Kelvin-Helmholtz, logrando la convergencia de la malla utilizando parámetros de malla moderados ($C_y=2.0, C_z=2.5, y_{in}^+=50$).
• Validación de TBL: La DNS de TBLs ZPG sobre riblets mostró un buen acuerdo con datos experimentales (Baron & Quadrio; Choi & Orchard) respecto a los porcentajes de reducción de arrastre y perfiles de velocidad, a pesar de las discrepancias inherentes entre experimento-DNS en las regiones de estela.

Significado
El artículo afirma que el marco $\eta$-grid ofrece un "enorme ahorro de malla", haciendo computacionalmente factible la DNS de turbulencia de pared a altos números de Reynolds sobre superficies complejas. Al alinear la resolución de la malla con las escalas físicas locales (escala de Kolmogorov) en lugar de mantener resoluciones fijas, el método reduce drásticamente el costo computacional (puntos de malla y restricciones de paso de tiempo) sin sacrificar la precisión. Los autores enfatizan que este enfoque permite la viabilidad económica de realizar DNS de TBLs sobre superficies complejas a números de Reynolds que coinciden con las condiciones experimentales, una tarea previamente obstaculizada por el costo prohibitivo de los métodos basados en cartesián. El trabajo valida la eficacia de las mallas no estructuradas para simulaciones de turbulencia de alta fidelidad, superando las restricciones de las mallas cartesianas estructuradas.