Distributed Roughness-Induced Transition on a Blunt Body at Mach 6: a Numerical Investigation

Este estudio presenta la primera simulación numérica directa de rugosidad distribuida en un cilindro romo a Mach 6, revelando que la disposición de la rugosidad dicta el mecanismo de transición ya sea promoviendo modos de estela sinuosa o ondas de T-S 2D, siendo estas últimas sostenidas de manera única por un bucle de retroalimentación acústica interna que elimina la necesidad de un forzamiento externo.

Lo esencial

Imagina una nave espacial precipitándose a través de la atmósfera a seis veces la velocidad del sonido. Para sobrevivir al intenso calor, su superficie está cubierta por un material especial que se quema lentamente (se ablaciona) para proteger la nave. Sin embargo, a medida que este material se quema, no deja una superficie perfectamente lisa, sino que deja una textura rugosa y con bultos, algo así como papel de lija.

Este artículo es una simulación computacional de alta velocidad que plantea una pregunta simple pero crítica: ¿Cómo es que estos diminutos bultos en la superficie convierten un flujo de aire suave y ordenado en un flujo de aire turbulento y caótico?

Aquí está el desglose de sus hallazgos utilizando analogías de la vida cotidiana:

1. La configuración: El cilindro de "papel de lija"

Los investigadores construyeron un modelo digital de un cilindro romo (como la nariz de un cohete) volando a Mach 6. En lugar de una superficie lisa, cubrieron la superficie con diminutos "bultos" artificiales (rugosidad) para imitar la textura arenosa dejada por el material que se quema.

Probaron tres formas diferentes de organizar estos bultos:

• Alineados: Como soldados parados en filas y columnas perfectas.
• Escalonados: Como una pared de ladrillos, donde los bultos de una fila están desplazados respecto a la fila de atrás.
• Aleatorios: Como guijarros esparcidos en una acera sin ningún patrón.

2. La teoría antigua frente al nuevo descubrimiento

Durante mucho tiempo, los científicos pensaron que la transición a la turbulencia era causada por una "acumulación lenta" de energía, similar a cómo un columpio gana altura si se le empuja de la manera justa a lo largo del tiempo. Esto se llama "crecimiento transitorio".

El hallazgo del artículo:
La simulación mostró que esta teoría de la "acumulación lenta de energía" no explica realmente lo que está sucediendo aquí. Los bultos en la superficie no solo amplificaron la energía lentamente; actuaron como desestabilizadores. Tomaron el flujo de aire e inmediatamente lo volvieron inestable, convirtiéndolo en un tipo específico de onda que crece muy rápido.

Piénselo de esta manera: la vieja teoría pensaba que los bultos estaban empujando suavemente un dominó para que cayera. El nuevo descubrimiento muestra que los bultos en realidad están pateando el dominó, haciendo que choque contra el siguiente de inmediato.

3. Los dos tipos de "ondas"

Dependiendo de cómo se organizaran los bultos, el flujo de aire reaccionaba de dos maneras diferentes:

• La "Serpiente" (Modo Sinuoso): Cuando los bultos estaban alineados (filas perfectas), el flujo de aire comenzó a zigzaguear de lado a lado como una serpiente. Este es un bamboleo muy específico y organizado.
• La "Onda Plana" (Ondas de Tollmien-Schlichting o T-S): Cuando los buntos estaban escalonados o eran aleatorios, el flujo de aire comenzó a ondular hacia arriba y hacia abajo en un patrón de onda 2D plana. Este es un tipo clásico de onda que suele encontrarse en flujos de aire mucho más lentos y de baja velocidad, lo cual fue sorprendente de encontrar en este entorno de alta velocidad.

La idea clave: La disposición de los bumentos dictaba qué "danza" realizaría el aire. La danza de la "serpiente" ocurrió con los bultos alineados, mientras que la danza de la "onda plana" ocurrió con los demás.

4. El final de la "Horquilla"

Una vez que estas ondas crecieron lo suficiente, desencadenaron la etapa final del choque. Las "vetas" constantes de aire creadas por los bultos (que son como largas cintas invisibles de aire lento) de repente se retorcieron y se rompieron en vórtices de horquilla.

Imagine una banda elástica que ha sido estirada con fuerza. De repente, se retuerce y forma un bucle que parece una horquilla. Estos bucles son el nacimiento de la turbulencia. Una vez que estas horquillas se forman, el aire suave se rompe completamente en el caos. Y el calor en la superficie de la nave espacial aumenta drásticamente.

5. La sorpresa de la "Cámara de Eco"

Uno de los descubrimientos más fascinantes fue cómo comenzó la turbulencia en primer lugar para los casos escalonados y aleatorios.

Normalmente, los científicos piensan que se necesita un "empujón" externo (como una ráfaga de viento o una vibración) para iniciar estas ondas. Pero la simulación mostró algo autosustentable:

1. La turbulencia comienza en un punto del cilindro.
2. Debido a que el aire detrás de la onda de choque se mueve más lento que el sonido (subsónico), el ruido de esta turbulencia viaja hacia atrás corriente arriba como un eco.
3. Este "eco" golpea la parte lisa de la superficie por delante de la turbulencia y excita el aire allí, causando que se inicie nueva turbulencia.
4. Esto crea un bucle de retroalimentación: la turbulencia genera ruido, el ruido viaja hacia atrás y el ruido crea más turbulencia.

Es como un micrófono captando su propia salida de altavoz y creando un bucle de retroalimentación chirriante, pero en este caso, el "chirrido" es el aire volviéndose turbulento.

Resumen

Este artículo utilizó una supercomputadora para observar el flujo de aire sobre un cilindro irregular a alta velocidad. Descubrieron que:

• El patrón de los bultos decide exactamente cómo el aire se vuelve turbulento.
• La vieja idea de la "acumulación lenta de energía" no es el principal culpable; en cambio, los bultos desestabilizan directamente ondas específicas.
• Estas ondas crecen hasta que se retuercen en formas de "horquilla", causando que el aire se vuelva caótico.
• En algunos casos, la turbulencia crea su propio "eco" que viaja hacia atrás para iniciar el proceso de nuevo, sin necesidad de ninguna ayuda externa.

Esto ayuda a los ingenieros a comprender que los diminutos y aleatorios bultos dejados por los escudos térmicos que se queman no son solo imperfecciones menores; son los arquitectos principales de cómo y cuándo la superficie de una nave espacial se calienta de manera peligrosa.

Resumen técnico

Resumen Técnico: Transición inducida por rugosidad distribuida en un cuerpo romo a Mach 6

Planteamiento del Problema
La predicción de la transición laminar-turbulenta (LTT) sobre geometrías romas de alta velocidad sigue siendo un desafío crítico para el diseño de ingeniería, particularmente para vehículos de entrada atmosférica donde la ablación superficial crea rugosidad distribuida. Aunque existen correlaciones empíricas (por ejemplo, Hollis 2017, 2019, 2025), estas carecen de un fundamento físico claro. Esfuerzos teóricos previos, como la correlación de Reshotko & Tumin (2004), dependen de mecanismos de crecimiento transitorio pero enfrentan limitaciones: asumen una escala lineal de las perturbaciones con la altura de la rugosidad (válida solo para $k/\delta$ pequeño), requieren constantes de ajuste con datos experimentales y tienen dificultades para explicar la transición en casos donde las alturas de rugosidad exceden el espesor de la capa límite ($k/\delta > 1$). Además, investigaciones numéricas previas a menudo simularon parches de rugosidad en lugar de rugosidad totalmente distribuida, o no observaron la amplificación de la inestabilidad modal, dejando sin aclarar los mecanismos específicos que impulsan la transición en flujos de cuerpo romo hipersónico con rugosidad totalmente distribuida.

Metodología
Este estudio presenta la primera Simulación Numérica Directa (DNS) de un cilindro romo en flujo transversal a Mach 6 con elementos de rugosidad distribuidos sobre toda la superficie. Las simulaciones utilizan el solver CHAMPS para integrar las ecuaciones de Navier-Stokes compresibles (CNSE) para un gas ideal. La geometría es un cilindro de 0.1524 m de diámetro, modelado con base en las condiciones experimentales de Hollis (2017) (Run 38, Test 6975), con un número de Mach de flujo libre de 5.98.

La rugosidad se modela como una distribución sinusoidal con una amplitud fija ($A \approx 35\%$ del espesor de la capa límite) y una longitud de onda correspondiente a la rugosidad nominal de grano de arena. Se investigaron tres configuraciones distintas de desfase en la dirección transversal:

1. Escalonada (Staggered): Las filas están desfasadas $180^\circ$.
2. Alineada (Aligned): Las filas están en fase (desplazamiento de $0^\circ$).
3. Aleatoria (Random): Las filas tienen desfases de fase distribuidos aleatoriamente.

Para analizar el proceso de transición, los autores emplearon:

• DNS: Para capturar la evolución no lineal completa del flujo.
• Ecuaciones de Perturbación Lineal (LDE): Resueltas sobre flujos medios promediados en el tiempo de la DNS para aislar los mecanismos de inestabilidad mediante forzamiento de pulso.
• Teoría de Estabilidad Lineal (LST): Aplicada a perfiles promediados en la dirección transversal para comparación.
• Análisis Espectral: Transformadas Rápidas de Fourier (FFT) de los campos de velocidad y presión para identificar frecuencias y números de onda dominantes.

Resultos Clave

• Flujo Medio y Estrias (Streaks): Todas las configuraciones generaron estrias de velocidad en la dirección de la corriente mediante un mecanismo de levantamiento (lift-up) impulsado por pares de vórtices contra-rotatorios (CVPs). La configuración alineada produjo las estrias estacionarias más fuertes debido a un desfase constante que eleva el fluido de bajo flujo de masa a una mayor altura dentro de la capa límite. Las configuraciones escalonada y aleatoria produjeron patrones de estrias más débiles y complejos. Contrario a las predicciones de la teoría de crecimiento transitorio óptimo, las amplitudes de las estrias exhibieron un crecimiento relativo bajo (menos de un factor de 2) y permanecieron no lineales en todo el dominio, desafiando la utilidad del crecimiento transitorio lineal como herramienta predictiva para este régimen.

• Mecanismos de Inestabilidad: El proceso de transición fue impulsado por inestabilidades convectivas, pero el tipo de modo específico dependió del desfase de la rugosidad:

  • Caso Alineado: Dominado por un modo de estria sinuoso fundamental. Las ondas de Tollmien-Schlichting (T-S) fueron débilmente desestabilizadas.
  • Casos Escalonado y Aleatorio: Dominados por ondas T-S 2D (modo varicoso). Estas configuraciones exhibieron una desestabilización significativa de las ondas T-S, llevando a un rápido crecimiento exponencial.
  • Ruptura (Breakdown): En todos los casos, una vez que las ondas de inestabilidad dominante crecieron hasta una amplitud suficiente, desencadenaron la formación de vórtices de tipo "hairpin" en las estrias estacionarias, llevando a la ruptura y a la turbulencia.

• Localizaciones de la Transición: La DNS predijo localizaciones de transición ($\Theta_{tr}$) significativamente más avanzadas (aproximadamente $11^\circ$–$17^\circ$) que las correlaciones empíricas de vanguardia (que predecían $6^\circ$–$13^\circ$). El caso alineado transicionó más temprano entre los casos de rugosidad debido al forzamiento requerido para iniciar la LTT, mientras que los casos escalonado y aleatorio transicionaron más tarde.

• Mecanismo de Retroalimentación: Un hallazgo novedoso en los casos escalonado y aleatorio fue la identificación de un bucle de retroalimentación. La turbulencia generada en la región subsónica del cilindro dispersó ondas acústicas que se propagaron hacia arriba (aguas arriba). Estas ondas acústicas lentas excitaron las ondas T-S en el lado opuesto del cilindro (a través de la línea de estancamiento) en su punto neutro. Este mecanismo, que no requiere forzamiento externo, sostiene el proceso de transición.

Significancia y Reivindicaciones
Los autores afirman que este trabajo proporciona una nueva comprensión física de la transición inducida por rugosidad en cuerpos romos hipersónicos, desafiando la dependencia prevalente en la teoría de crecimiento transitorio. Las contribuciones clave son:

1. Rechazo del Crecimiento Transitorio como Impulsor Primario: El estudio demuestra que para una rugosidad totalmente distribuida con $k/\delta \approx 0.35$, el desarrollo no lineal de las estrias y las bajas relaciones de amplificación observadas hacen que la teoría de crecimiento transitorio óptimo sea un predictor insuficiente del proceso de transición real.
2. Desestabilización de Inestabilidades Modales: El artículo establece que la rugosidad distribuida puede desestabilizar eficazmente las inestabilidades modales (específicamente las ondas T-S) incluso en cuerpos romos donde el flujo de pared lisa es típicamente estable de modo propio. La configuración de la rugosidad dicta si la transición es impulsada por modos de estria o por ondas T-S.
3. Descubrimiento de un Bucle de Retroalimentación Acústica: La identificación de un mecanismo de retroalimentación autosostenido, donde la acústica generada por la turbulencia excita las ondas T-S en el lado opuesto del cuerpo, ofrece una explicación potencial para la transición en regiones subsónicas de flujo hipersónico, un fenómeno previamente no explicado por mecanismos lineales.
4. Sensibilidad al Desfase: Los resultados resaltan la extrema sensibilidad de la LTT al desfase relativo de los elementos de rugosidad, sugiriendo que topologías simples pueden delimitar el comportamiento de patrones de rugosidad más complejos y aleatorios.

Los autores concluyen que, si bien el mecanismo de retroalimentación acústica propuesto requiere una prueba más rigurosa, la exposición de estos mecanismos de inestabilidad mediante DNS de alta fidelidad sugiere un cambio fundamental en la comprensión de las perturbaciones responsables de la LTT en cuerpos romos hipersónicos.