🔬 physics

Steady, Quasi-Periodic, and Chaotic Regimes of Natural Convection Around a Heated Cylinder in a Square Cavity

Este estudio numérico de la convección natural alrededor de un cilindro calentado en una cavidad cuadrada revela que el flujo transiciona de un estado estacionario al caos a través de la ruta de Ruelle-Takens-Newhouse a medida que el número de Rayleigh aumenta de $10^4$ a $5 \times 10^7$ , estableciendo al mismo tiempo una correlación específica de ley de potencia para el número de Nusselt del cilindro y destacando el papel dominante de la pared superior en la transferencia de calor debido a la estratificación térmica.

Autores originales: Mannu Yadav, Anil Lal Sadasivan

Publicado 2026-07-07

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Mannu Yadav, Anil Lal Sadasivan

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

Imagina una habitación cuadrada con paredes frías y un calentador brillante y caliente situado justo en el medio. Esta es la "cavidad cuadrada" en el estudio. El aire dentro no está simplemente quieto; debido a que el calentador está caliente y las paredes están frías, el aire quiere moverse. El aire caliente sube, el aire frío baja, creando una circulación natural. Esto es convección natural.

Los investigadores querían ver qué le sucede a este baile del aire a medida que aumentan el calor (representado por algo llamado número de Rayleigh). No se limitaron a suponer; construyeron un modelo digital superpreciso utilizando un programa informático para observar el movimiento del aire en tiempo real, sin forzarlo a comportarse de ninguna manera específica.

Aquí está la historia de lo que encontraron, desglosada en pasos sencillos:

1. El comienzo tranquilo (Estado estacionario)

Cuando el calor es bajo o moderado, el aire se comporta como un ballet bien ensayado.

El patrón: Dos bucles gigantes de aire, imagen especular, se forman a ambos lados del calentador. El aire sube por los lados del calentador, golpea el techo frío, se desliza por las paredes y regresa al fondo.
La pluma: Una única columna constante de aire caliente (una "pluma") sube recta desde la parte superior del calentador hacia el techo.
El resultado: Todo es predecible y simétrico. Si dibujaras una línea por el medio de la habitación, el lado izquierdo se vería exactamente igual al derecho. Esto ocurre hasta cierto nivel de calor (unos 5 millones de unidades de calor).

2. La primera señal de problemas (El punto de inflexión)

A medida que aumentaban el calor solo un poco más, algo interesante sucedió. La simetría perfecta empezó a agrietarse, aunque aún no se podía ver a simple vista.

El cambio: Los bucles de aire de la izquierda y la derecha dejaron de ser perfectamente iguales en fuerza. Uno empezó a volverse ligeramente más fuerte que el otro.
El descubrimiento: Los investigadores encontraron que el "punto de inflexión" donde el flujo tranquilo y constante se rompe hacia el caos ocurre en algún lugar entre 5 millones y 7,5 millones de unidades de calor. Antes de esto, el aire estaba tranquilo; después, empezó a oscilar.

3. El baile se complica (Cuasi-periódico)

Una vez que el calor superó ese punto de inflexión, el aire dejó de ser un ballet tranquilo y empezó a realizar un baile rítmico complejo.

El ciclo límite: Al principio, el aire empezó a balancearse de un lado a otro en un único ritmo perfecto (como un péndulo).
El toroide doble: A medida que aumentaban el calor más allá de eso, el baile se volvió más complicado. Ya no era solo balancearse de un lado a otro; empezó a tambalearse mientras se balanceaba. Imagina un trompo que también se tambalea de lado a lado. El aire seguía ahora dos ritmos diferentes a la vez que no terminaban de encajar. Esto se llama "cuasi-periódico".

4. La tormenta caótica (Caos)

En los niveles de calor más altos (50 millones de unidades), el baile ordenado se desmoronó por completo.

El atractor extraño: El aire se volvió caótico. No era ruido aleatorio, sino que era impredecible. La pluma caliente no solo subía recta; se inclinaba a la izquierda, luego a la derecha, y luego se desprendía y giraba en pequeños grumos.
La propagación del caos: Curiosamente, el caos no ocurrió en todas partes al mismo tiempo. Comenzó en la parte superior, justo encima del calentador (la plema). El aire en la parte inferior de la habitación se mantuvo tranquilo y estratificado (en capas) durante mucho tiempo. Solo cuando el calor llegó a ser extremadamente alto, el remolino caótico finalmente alcanzó el fondo de la habitación.

5. La simetría "fantasma"

Aquí hay una parte truculenta que el artículo destaca:

Instante vs. Promedio: Si tomaras una instantánea del aire en cualquier momento único durante la fase caótica, se vería desordenado y asimétrico (desequilibrado). Sin embargo, si tomaras una foto de "lapso de tiempo" (promediando todo durante un largo periodo), la confusión se cancelaría y la imagen se vería perfectamente simétrica de nuevo, igual que al principio tranquilo.
La metáfora: Imagina una multitud de personas corriendo en un desorden caótico y revuelto. Si tomas una foto de un segundo, es un borrón. Si tomas una foto durante una hora, el borrón se promedia y parece un círculo perfecto. El promedio parece tranquilo, pero la realidad es caótica.

6. La transferencia de calor (¿A dónde va el calor?)

Los investigadores también rastrearon cuánto calor absorbían las paredes.

El techo es el ganador: Independientemente de lo caótico que se volviera el aire, la pared superior siempre absorbió aproximadamente el 55% del calor. Las paredes laterales compartieron el resto (alrededor del 22% cada una).
El suelo es ignorado: La pared inferior absorbió casi nada (menos del 2%).
¿Por qué?: A medida que el aire se calentaba, el aire frío en la parte inferior se quedó "atascado" en una capa. El aire caliente que subía del calentador y el aire frío que bajaba por los lados crearon una barrera que mantuvo aislado el suelo. Es como si el fondo de la habitación se hubiera convertido en una "zona muerta" donde el calor no podía llegar, a pesar de que el calentador trabajaba más que nunca.

Resumen

El artículo cuenta la historia del aire en una caja con un calentador.

Calor bajo: Bucles tranquilos y simétricos.
Calor medio: Los bucles empiezan a tambalearse y a romper la simetría.
Calor alto: El aire entra en un baile rítmico complejo con dos frecuencias.
Calor extremo: El aire entra en un caos total, comenzando en la parte superior y extendiéndose lentamente hacia abajo.
La lección: Incluso cuando el aire es caótico, la distribución de calor promedio parece tranquila, y la parte inferior de la habitación permanece sorprendentemente fría porque el calor queda atrapado en las capas superiores.

Los investigadores utilizaron un método computacional especial (Método de Elementos Finitos) para asegurar que su modelo digital fuera preciso, verificándolo con puntos de referencia conocidos antes de realizar sus nuevos experimentos. Descubrieron que estudios anteriores podrían haber pasado por alto el inicio de este caos porque forzaban al aire a mantenerse simétrico, lo que impedía que el "tambaleo" ocurriera de forma natural.

Resumen Técnico: Regímenes Estacionarios, Cuasi-Periódicos y Caóticos de la Convección Natural alrededor de un Cilindro Caliente en una Cavidad Cuadrada

Planteamiento del Problema
Este estudio investiga la convección natural no estacionaria del aire ($Pr = 0.71$) dentro de una cavidad cuadrada fría que contiene un cilindro circular calentado centralmente con un radio $R = 0.2L$ . Si bien la cavidad cuadrada diferencialmente calentada es un estándar de referencia, la adición de un cilindro interno calentado introduce una reorganización compleja del flujo. La literatura previa se ha centrado principalmente en soluciones de estado estacionario o rangos no estacionarios limitados, imponiendo a menudo restricciones de simetría artificiales (por ejemplo, fijar la función de corriente a cero en el cilindro) que suprimen los modos de flujo asimétricos. En consecuencia, el número de Rayleigh crítico ( $Ra_c$ ) para el inicio de la inestabilidad, la ruta específica hacia el caos y la secuencia completa de regímenes de flujo para esta configuración permanecieron sin resolver. Este trabajo tiene como objetivo cerrar estas brechas mediante la simulación de la transición desde el flujo laminar estacionario hacia el caos sin imponer simetría, permitiendo que los modos asimétricos crezcan de forma natural.

Metodología
Los autores emplearon un método de elementos finitos (FEM) basado en presión sobre una malla conforme al cuerpo. Las características metodológicas clave incluyen:

Ecuaciones de Gobierno: Se resolvieron las ecuaciones de Navier-Stokes 2D no estacionarias y de energía bajo la aproximación de Boussinesq en variables primitivas. No se impusieron condiciones de simetría en la línea de centro vertical.
Generación de Malla: Se generó una malla conforme al cuerpo a partir de una rejilla de fondo cartesiana mediante una proyección de ray-cast seguida de un suavizado armónico. Este enfoque resuelve la frontera curva del cilindro sin necesidad de transformaciones de coordenadas o interpolación de frontera inmersa.
Esquema Numérico: Se utilizó un algoritmo de paso fraccionado explícito (proyección) con un tratamiento de Adams-Bashforth de segundo orden para los términos de convección-difusión. La ecuación de Poisson de la presión se resolvió mediante un método de Gradiente Conjugado con inicio en caliente (warm-started) utilizando productos de matriz libre de Elemento por Elemento (EbE) para manejar mallas no estructuradas de manera eficiente.
Validación e Independencia de la Malla: El solver fue validado contra el benchmark de De Vahl Davis ( $Ra = 10^6$ ). Un estudio del Índice de Convergencia de la Malla (GCI) confirmó la independencia de la malla con una incertidumbre de malla fina del 0.84% en el número de Nusselt del cilindro.
Alcance de la Simulación: Se simularon ocho números de Rayleigh, que oscilan entre $10^4$ y $5 \times 10^7$ . Se registraron las historias temporales del número de Nusselt promedio del cilindro ( $Nu_{cyl}$ ), y se generaron densidades espectrales de energía (ESD) y mapas de Poincaré utilizando datos de once puntos de sonda distribuidos en la cavidad.

Resultajes Clave
El estudio identifica tres regímenes de flujo distintos y caracteriza la transición entre ellos:

Régimen Estacionario ( $10^4 \le Ra \le 5 \times 10^6$ ):
- El flujo permanece estacionario y simétrico respecto al eje vertical, caracterizado por dos celdas de recirculación de imagen especular y una única pluma térmica que asciende desde el cilindro hacia la pared superior.
- A $Ra = 5 \times 10^6$ , aparece la primera firma física de la ruptura de simetría: la magnitud del núcleo de la función de corriente negativa ( $|\psi_{min}|$ ) excede al núcleo positivo ( $\psi_{max}$ ) en aproximadamente un 3.6%, aunque el flujo permanece visualmente simétrico.
- La transferencia de calor está dominada por la pared superior (absorbiendo ~55% del calor total de la pared), mientras que la pared inferior representa menos del 2% debido a la estratificación térmica en la cavidad inferior.
Inicio de la Inestabilidad ( $5 \times 10^6 < Ra_c < 7.5 \times 10^6$ ):
- El número de Rayleigh crítico para el inicio de la inestabilidad se encuentra entre $5 \times 10^6$ y $7.5 \times 10^6$ . Esto se determina por un salto brusco (de más de tres órdenes de magnitud) en la desviación estándar de $Nu_{cyl}$ .
- Este valor crítico es mayor que el implícito en estudios previos de solo estado estacionario que imponían simetría, confirmando que la supresión de los modos asimétricos retrasa artificialmente la transición.
Régimen No Estacionario ( $Ra \ge 7.5 \times 10^6$ ):
- Ciclo Límite ( $Ra = 7.5 \times 10^6$ ): El flujo exhibe una oscilación periódica de frecuencia única (ciclo límite). El desorden se origina en la pluma sobre el cilindro.
- Cuasi-Periódico ( $Ra = 10^7$ ): Emerge una segunda frecuencia inconmensurable, situando el flujo en un toro de dos dimensiones. Las secciones de Poincaré muestran una curva cerrada, y el estanque estratificado inferior comienza a registrar la frecuencia secundaria.
- Caos ( $Ra = 5 \times 10^7$ ): El flujo transiciona hacia un atractor extraño de banda ancha. Las secciones de Poincaré revelan una nube dispersa de puntos, y la densidad espectral de energía se vuelve continua. El desorden se propaga desde la pluma hacia las paredes laterales, alcanzando finalmente el estanque estratificado inferior, el cual se vuelve activo a través de una amplia gama de escalas temporales.
- Simetría en los Promedios: A pesar de que las instantáneas son fuertemente asimétricas, los campos promediados en el tiempo recuperan la simetría bilateral. Este es un efecto estadístico donde el sistema visita estados de imagen especular con la misma frecuencia a lo largo del tiempo.

Correlaciones de Transferencia de Calor

Una única ley de potencia, $Nu_{cyl} = 0.764 Ra^{0.226}$ , se ajusta a los números de Nusselt del cilindro promediados en el tiempo a través de todo el rango ( $10^4 \le Ra \le 5 \times 10^7$ ) con una desviación de menos del 6%.
El exponente para la rama estacionaria es $n \approx 0.224$ , y para la rama no estacionaria es $n \approx 0.243$ . Ambos valores se sitúan entre la predicción clásica de la capa límite laminar ( $n=0.25$ ) y el escalamiento de turbulencia suave de Castaing ( $n \approx 0.286$ ), lo que sugiere que el flujo no es totalmente turbulento pero está influenciado por campos fluctuantes.
El número de Nusselt de la pared inferior ( $Nu_{bot}$ ) disminuye a medida que $Ra$ aumenta, una tendencia contraintuitiva atribuida al fortalecimiento de la estratificación térmica y al atrapamiento de chorros fríos en los vórtices de las esquinas inferiores.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera caracterización numérica integral de la secuencia completa de transición desde la convección natural estacionaria hasta la caótica en una configuración de cilindro calentado centralmente sin restricciones de simetría artificiales. Las contribuciones clave son:

Número de Rayleigh Crítico: Establecer un rango más preciso para $Ra_c$ al permitir el desarrollo de modos asimétricos, corrigiendo sobreestimaciones previas causadas por las restricciones de simetría.
Ruta al Caos: Confirmar que la transición sigue la ruta de Ruelle-Takens-Newhouse (Ciclo Límite $\to$ Toro de dos dimensiones $\to$ Atractor Extraño) en lugar de una cascada de duplicación de periodo.
Propagación Espacial del Desorden: Demostrar que el caos se inicia en la pluma térmica y se propaga hacia abajo, manteniendo el estanque estratificado inferior relativamente quieto hasta los números de Rayleigh más altos.
Datos de Referencia: Proporcionar datos de referencia fiables para el número de Rayleigh crítico, las transiciones de régimen y las correlaciones globales de transferencia de calor para configuraciones donde la geometría interna y la dinámica no estacionaria interactúan.

Los autores señalan que el estudio está limitado a 2D, un único radio de cilindro ( $R=0.2L$ ) y $Pr=0.71$, sugiriendo extensiones a 3D, posiciones descentradas y variaciones en los números de Prandtl como direcciones futuras.