Fully Turbulent Wakes at Low Reynolds Numbers: the Case of the Thin Flat Plate

Este artículo demuestra, mediante simulación numérica directa y comparación experimental, que el flujo de estela detrás de una placa plana bidimensional delgada se vuelve completamente turbulento a un número de Reynolds relativamente bajo de 400, exhibiendo características estadísticas y espectrales indistinguibles de las estelas turbulentas de números de Reynolds más altos, una vía de transición que difiere significamente de la de los cilindros canónicos circulares o cuadrados.

Lo esencial

Imagina que sostienes un trozo de cartón delgado y plano (como una carta de juego) frente a un viento fuerte. A medida que el viento golpea la carta, crea un rastro de aire desordenado y agitado detrás de ella llamado "estela". Durante mucho tiempo, los científicos creyeron que para que esta estela se volviera verdaderamente caótica o "turbulenta", el viento tendría que soplar muy rápido, o la carta tendría que tener una forma específica como un tubo redondo o un bloque cuadrado.

Este artículo cuenta una historia diferente. Los investigadores descubrieron que, si utilizas una placa delgada y plana, el aire detrás de ella se vuelve completamente caótico y turbulento a una velocidad de viento mucho menor de lo que nadie esperaba. De hecho, sucede a una velocidad en la que, para otras formas, el aire todavía es relativamente tranquilo y ordenado.

Aquí hay un desglose de sus hallazgos utilizando analogías simples:

1. La sorpresa del "Umbral de Turbulencia"

Piensa en la turbulencia como una pista de baile abarrotada.

• La creencia antigua (Cilindros redondos): Si tienes un poste redondo en el viento, el aire detrás de él comienza como un baile tranquilo y rítmico (balanceándose de un lado a otro). Se necesita mucha energía (alta velocidad) antes de que los bailarines comiencen a chocar entre sí, girar salvajemente y crear un caos desordenado (turbulencia). Esta transición ocurre lentamente a través de un amplio rango de velocidades.
• El nuevo descubrimiento (Placa plana delgada): Los investigadores descubrieron que, para una placa delgada y plana, la "pista de baile" pasa de estar tranquila a ser un mosh pit salvaje casi instantáneamente. Incluso a una velocidad de viento relativamente baja (número de Reynolds 400), el aire detrás de la placa ya es completamente caótico. No pasa por las etapas lentas y rítmicas que hacen los postes redondos; salta directamente a la fiesta.

2. Cómo lo demostraron

Para estar seguros de que no estaban imaginando cosas, el equipo actuó como detectives comparando escenas del crimen.

• La simulación (El laboratorio virtual): Utilizaron supercomputadoras para simular el viento golpeando la placa a bajas velocidades (Re 150 y Re 400).
• La prueba del mundo real (El túnel de viento): También observaron experimentos reales donde el viento soplaba mucho más rápido (Re 12,500 y Re 19,700).
• La coincidencia: Cuando compararon la simulación de baja velocidad por computadora (Re 400) con los experimentos del mundo real de alta velocidad, los patrones coincidieron perfectamente. Las "huellas dactilares" de la turbulencia —cómo se movía el aire, cuánta energía tenía y cómo giraba— eran idénticas.
• El grupo de control: Cuando observaron la simulación a una velocidad aún más baja (Re 150), los patrones eran totalmente diferentes. Todavía estaba en la "fase tranquila", no era caótica aún. Esto demostró que la transición al caos ocurre en algún punto entre 150 y 400.

3. La "Huella Dactilar" de la turbulencia

¿Cómo sabes si un flujo es verdaderamente turbulento? El artículo busca "signos de vida" específicos en los datos:

• El espectro de energía (El sonido del viento): En un flujo tranquilo, la energía se concentra en unas pocas notas específicas (como una flauta tocando un solo tono). En un flujo turbulento, suena como ruido blanco o estática, con la energía distribuida en un enorme rango de frecuencias. Los investigadores encontraron que, en Re 400, el "sonido" del viento detrás de la placa ya estaba lleno de esta estática caótica, tal como en los experimentos de alta velocidad.
• La "intermitencia" (El grito ocasional): En un flujo verdaderamente turbulento, el aire no solo gira suavemente; tiene ráfagas repentinas e intensas de velocidad y rotación. Los investigadores encontraron estos "gritos" en los datos en Re 400, pero estaban ausentes en Re 150.

4. ¿Por qué es esto diferente?

El artículo sugiere que la razón de este salto repentino es la forma del objeto.

• Objetos redondos/cuadrados: Cuando el viento golpea un objeto redondo o cuadrado, la parte trasera del objeto actúa como un escudo, estabilizando el flujo de aire detrás de él. Se necesita mucha energía para romper esa estabilidad.
• La placa delgada: Debido a que la placa es tan delgada, no hay una "parte trasera" que proteja el aire. Las fluctuaciones de presión (el empuje y la tracción del aire) están conectadas directamente con los vórtices giratorios desde el principio. Es como intentar equilibrar un lápiz sobre su punta frente a equilibrar una bola de boliche; el lápiz (la placa delgada) es inherentemente inestable y se vuelca hacia el caos mucho más rápido.

La conclusión fundamental

Este artículo cambia nuestra comprensión de cómo fluye el aire alrededor de objetos planos. Demuestra que las placas planas delgadas crean estelas completamente turbulentas a velocidades sorprendentemente bajas, mucho menores que los objetos redondos o cuadrados. La transición al caos no es un proceso lento y gradual para estas formas; es un cambio repentino y fundamental que ocurre muy temprano en el rango de velocidad.

Los investigadores no discutieron cómo esto se aplica a la construcción de puentes, el diseño de autos o dispositivos médicos. Se centraron estrictamente en demostrar que este fenómeno existe y cómo difiere la física del flujo de aire de lo que pensábamos anteriormente.

Resumen técnico

Resumen Técnico: Estelas Completamente Turbulentas a Bajos Números de Reynolds: El Caso de la Placa Plana Delgada

Planteamiento del Problema
La simulación numérica del flujo nominalmente bidimensional (2D) alrededor de una placa plana delgada normal a una corriente uniforme ha resultado históricamente desafiante, con estudios de dinámica de fluidos computacional (CFD) que arrojan resultados contradictorios respecto a la dinámica de la estela a bajos números de Reynolds ($Re < 1000$). Mientras que los datos experimentales sugieren que los parámetros integrales (como los coeficientes de arrastre y los números de Strouhal) dependen débilmente de $Re$ en el régimen 3D, las simulaciones han reportado comportamientos de estela notablemente distintos. Una cuestión central permanece: ¿transiciona la estela de una placa plana delgada a un estado completamente turbulento inmediatamente después del inicio de las inestabilidades 3D, de forma similar a las estelas de cilindros canónicos, o sigue una trayectoria de transición de múltiples etapas distinta? Los autores investigan si la estela tras una placa plana delgada ya es completamente turbulenta a un número de Reynolds relativamente bajo de $Re = 400$, un régimen donde las estelas de cilindros circulares y cuadrados canónicos suelen exhibir una fuerte dependencia de $Re$ y características de fluctuación no turbulentas.

Metodología
El estudio emplea un enfoque multifacético que combina Simulaciones Numéricas Directas (DNS) y validación experimental:

1. Simulaciones Numéricas: Se utilizaron dos códigos de DNS independientes para asegurar la robustez:

   • Método de Elementos Espectrales/hp (SEM): Implementado en Nektar++, utilizado para la DNS principal de alta precisión a $Re = 400$ (DNS400) y $Re = 150$ (DNS150).
   • Método de Volúmenes Finitos (FVM): Implementado en OpenFOAM, utilizado para estudios de convergencia de malla y validación cruzada (DNS400V).
   • El dominio computacional fue truncado y sometido a condiciones de no deslizamiento en la placa, condiciones de entrada de flujo libre y condiciones de salida específicas propuestas por Dong et al. Se impuso la periodicidad en la dirección de la envergadura.
   • La convergencia de la malla fue probada rigurosamente, asegurando la resolución de la escala de longitud de Kolmogorov ($\eta$) con valores de $\Delta x/\eta$ comparables o superiores a la literatura existente para geometrías similares.

2. Datos Experimentales: Se adquirieron dos conjuntos de datos experimentales utilizando PIV estereoscópica de tiempo resuelto en un túnel de viento en la Universidad de Calgary:

   • EX12K: $Re = 12,500$
   • EX20K: $Re = 19,700$
   • Estos conjuntos de datos representan regímenes completamente turbulentos y sirven como referencia para las simulaciones de bajo $Re$.

3. Análisis: El estudio compara cantidades integrales globales, campos de velocidad media, tensiones de Reynolds, presupuestos de energía cinética turbulenta (TKE), espectros de energía y propiedades estadísticas de los gradientes de velocidad (específicamente la no-gaussianidad) a través de los diferentes casos de $Re$.

Contribuciones Clave y Resultados

• Inicio Temprano de la Turbulencia: El hallazgo principal es que la estela de una placa plana delgada se vuelve completamente turbulenta a $Re = 400$. Esto es significativamente menor que la transición observada en las estelas de cilindros circulares ($Re \approx 10^5$) o cuadrados ($Re \approx 6000$) canónicos.
• Acuerdo Cuantitativo con Experimentos de Alto $Re$: Los resultados de la DNS a $Re = 400$ (DNS400) muestran un acuerdo notable con los datos experimentales a $Re = 12,500$y$19,700$ (EX12K, EX20K) en cuanto a:
  • Cantidades Globales: El coeficiente de arrastre medio ($C_D$) y el número de Strouhal ($St$) son funciones de variación lenta de $Re$ para $Re \ge 400$.
  • Estructura del Flujo: Las líneas de corriente medias, los perfiles de velocidad y las distribuciones de tensión de Reynolds ($u'^2, v'^2, w'^2$) son estadísticamente indistinguibles entre el caso $Re=400$ y los experimentos de alto $Re$ dentro de la incertidumbre de la medición.
  • Presupuestos de Energía: La distribución espacial de los términos de producción y advección de TKE en la ecuación de transporte de $k$ es consistente entre el caso $Re=400$ y los casos de alto $Re$.
• Firmas Turbulentas a Bajo $Re$: El caso $Re=400$ exhibe características definitivas de turbulencia ausentes en el caso $Re=150$:
  • Espectros de Energía: Los componentes de velocidad fluctuante muestran espectros de banda ancha con un escalamiento de rango inercial aproximadamente de $f^{-5/3}$ (ley de Kolmogorov), indicativo de un cascada de energía. En contraste, los espectros de $Re=150$ muestran transferencias de energía entre solo unos pocos modos.
  • No-Gaussianidad: Las Funciones de Densidad de Probabilidad (PDF) de las tasas de deformación ($\Sigma$) y rotación ($\Omega$) fluctuantes exhiben colas pesadas, lo que significa intermitencia y eventos de alta intensidad poco comunes, característicos de la turbulencia. El caso $Re=150$ muestra distribuciones gaussianas.
• Trayectoria de Transición Distinta: La transición a $Re=150$ involucra inestabilidades 3D pero carece de la estructura coherente y plenamente turbulenta vista en $Re=400$. La estela de $Re=150$ presenta un flujo de retorno más débil, una recuperación de la estela más lenta y una evolución no monotónica del ancho de la estela, diferenciándose fundamentalmente del comportamiento de los cilindros canónicos donde las inestabilidades 3D (Modos A, B, C) evolucionan a través de un amplio rango de $Re$ antes de alcanzar la turbulencia total.

Significancia y Reivindicaciones
El artículo argumenta que el camino hacia la transición a la turbulencia para una placa plana delgada es fundamentalmente diferente al de los cilindros canónicos. Mientras que las estelas de los cilindros pasan por una secuencia prolongada de inestabilidades (Modo A, B, C) durante un amplio rango de números de Reynolds antes de alcanzar un estado totalmente turbulento, la estela de la placa plana delgada transiciona a un estado completamente turbulento casi inmediatamente después del inicio de las inestabilidades 3D (entre $Re=150$y$Re=400$).

Los autores sugieren que esta diferencia puede atribuirse a la falta de un cuerpo posterior (afterbody) en la dirección de la corriente de la geometría. Para los cilindros, el cuerpo posterior estabiliza el flujo y desacopla las fluctuaciones de presión del punto de separación. En contraste, para las placas planas delgadas, las fluctuaciones de presión están directamente acopladas a la región de formación de vórtices, lo que podría impulsar una transición a la turbulencia más temprana y rápida. Estos hallazgos ayudan a reconciliar las discrepancias en la literatura previa sobre la sensibilidad de las estelas de placas delgadas a las relaciones de aspecto y números de Reynolds, estableciendo que el régimen "totalmente turbulento" para esta geometría comienza a un $Re$ mucho más bajo de lo asumido anteriormente para cuerpos de choque canónicos.