HYMOR: An open-source package for global modal, non-modal, and receptivity analysis in high-enthalpy hypersonic vehicles

El artículo presenta HYMOR, un paquete de código abierto en MATLAB y Julia que realiza análisis de estabilidad lineal global (modal, no modal y de receptividad) para flujos hipersónicos de alta entalpía, incorporando efectos de gas real y una formulación de ajuste de choque para capturar interacciones físicas inaccesibles mediante métodos locales tradicionales.

Lo esencial

¡Hola! Imagina que estás diseñando un cohete o una nave espacial que viaja tan rápido que el aire a su alrededor se calienta tanto como la superficie del Sol. Este es el mundo de la hipersónica.

El problema principal es que, a esas velocidades, el aire no se comporta como un gas normal; se vuelve un "plasma" químico, se ioniza y cambia de forma. Lo más peligroso es que, si la capa de aire pegada al cohete (la capa límite) pasa de ser suave (laminar) a caótica (turbulenta), el calor aumenta drásticamente y podría derretir la nave.

Aquí es donde entra HYMOR, la herramienta que presentan en este artículo. Vamos a explicarlo con analogías sencillas:

1. ¿Qué es HYMOR?

Imagina que HYMOR es un simulador de vuelo ultra-preciso para ingenieros, pero en lugar de pilotar, tú eres un "detective de ondas".

• Su misión: Predecir exactamente cuándo y dónde el aire suave alrededor de la nave se volverá turbulento.
• Su superpoder: Es de código abierto (como Wikipedia o Linux), lo que significa que cualquier científico en el mundo puede usarlo, modificarlo y mejorar el diseño de naves espaciales sin pagar licencias costosas. Está escrito en dos lenguajes de programación populares (MATLAB y Julia), como si tuviera dos idiomas para hablar con diferentes tipos de usuarios.

2. El problema de la "Marea" (La onda de choque)

Cuando una nave viaja más rápido que el sonido, empuja el aire creando una onda de choque (una barrera invisible de aire comprimido) justo frente a su nariz. Es como el pico de agua que se forma frente a un barco que va muy rápido.

• El problema de los métodos antiguos: La mayoría de los programas de computadora tratan esta onda de choque como una "niebla" o un borrón difuso. Imagina intentar medir la velocidad de una pelota de tenis golpeando una pared, pero la cámara de video tiene el enfoque borroso. Eso genera errores.
• La solución de HYMOR (Shock-fitting): HYMOR trata la onda de choque como una cortina de acero perfectamente afilada. No es un borrón; es una línea exacta. Esto permite ver con precisión milimétrica cómo las pequeñas perturbaciones (como una ráfaga de viento en el espacio) rebotan y se transforman al cruzar esa cortina. Es como cambiar de una cámara borrosa a una de alta definición 8K.

3. Los tres modos de detección de HYMOR

HYMOR tiene tres formas de investigar por qué el aire se vuelve turbulento, como si tuviera tres lentes diferentes:

• A. El Lente de "Modos" (Modal Analysis):
  Imagina que el aire es una cuerda de guitarra. Si la tocas, vibra en una nota específica. Este modo busca esas "notas" naturales que hacen que la vibración crezca por sí sola. Si encuentra una nota que crece, ¡alerta! La turbulencia llegará pronto.

• B. El Lente de "Crecimiento Transitorio" (Non-modal):
  A veces, aunque no haya una "nota" que crezca sola, si empujas la cuerda de la forma correcta (con la fuerza y el ángulo adecuados), puede vibrar mucho más fuerte por un momento antes de calmarse. HYMOR busca esos empujones perfectos. Es como encontrar la forma exacta de soplar una vela para que se apague, aunque el viento sea suave.

• C. El Lente de "Receptividad" (Receptivity):
  Esto es crucial. Imagina que el aire exterior tiene pequeñas "suciedades" o ráfagas (disturbios). ¿Cómo entra esa suciedad dentro del sistema de la nave? HYMOR rastrea cómo esas pequeñas ráfagas del espacio exterior cruzan la "cortina de acero" (la onda de choque) y se convierten en grandes problemas dentro de la capa de aire de la nave. Es como ver cómo una pequeña gota de lluvia entra en un barco a través de una escotilla y termina inundando la bodega.

4. ¿Por qué es tan importante esto?

En el pasado, los ingenieros tenían que hacer naves muy pesadas y con escudos térmicos gigantes "por si acaso", porque no podían predecir con exactitud cuándo se calentaría la nave.

Con HYMOR, podemos:

1. Ahorrar peso: Diseñar escudos térmicos más delgados y eficientes.
2. Viajar más lejos: Entender mejor cómo funcionan las naves en la atmósfera de Marte o en viajes a velocidades extremas.
3. Acelerar la ciencia: Al ser gratuito y abierto, todo el mundo puede usarlo para probar nuevas ideas sin tener que construir un programa desde cero.

En resumen

HYMOR es como un laboratorio virtual de alta tecnología que permite a los científicos ver el invisible. Usa matemáticas avanzadas para tratar las ondas de choque como líneas perfectas (no como borrones) y descubre cómo las pequeñas perturbaciones del espacio pueden convertirse en grandes tormentas de calor alrededor de una nave.

Es una herramienta que promete hacer que el viaje al espacio sea más seguro, más ligero y más eficiente para todos. ¡Y lo mejor es que cualquiera puede descargarla y empezar a usarla hoy mismo!

Resumen técnico

Resumen Técnico: HYMOR - Un paquete de código abierto para análisis de estabilidad modal, no modal y receptividad en vehículos hipersónicos de alta entalpía

1. Planteamiento del Problema

La transición de flujo laminar a turbulento en flujos hipersónicos sigue siendo uno de los problemas abiertos más críticos en la aerodinámica de alta velocidad. La capacidad de predecir cuándo y dónde ocurre esta transición es esencial para dimensionar los sistemas de protección térmica y establecer márgenes de diseño.

• Limitaciones de los métodos actuales: Las herramientas tradicionales de estabilidad lineal, como la Teoría de Estabilidad Lineal (LST) y las Ecuaciones de Estabilidad Parabolizadas (PSE), asumen flujos paralelos o débilmente no paralelos. Estas aproximaciones no pueden capturar interacciones entre mecanismos físicos separados espacialmente ni fenómenos globales.
• El desafío de las ondas de choque: En flujos de alta entalpía, las ondas de choque fuertes juegan un papel crucial en la receptividad (cómo las perturbaciones del flujo libre excitan la capa límite). Los métodos numéricos convencionales de "captura de choque" (shock-capturing) introducen un espesor artificial del choque que corrompe la respuesta de las perturbaciones infinitesimales, invalidando el análisis de interacción lineal (LIA).
• Falta de herramientas integradas: No existía un paquete de código abierto que integrara análisis modal global, crecimiento transitorio (no modal) y receptividad del flujo libre para flujos de alta entalpía con efectos de gas real.

2. Metodología y Formulación

El paquete HYMOR (HYpersonic MOdal/non-modal, and Receptivity) es un marco computacional desarrollado en MATLAB y Julia que aborda estas limitaciones mediante las siguientes estrategias:

• Formulación Global: A diferencia de los métodos locales, HYMOR resuelve la dependencia espacial completa del campo de perturbaciones, permitiendo capturar interacciones entre regiones del flujo separadas.
• Enfoque de Ajuste de Choque (Shock-Fitting):
  • El choque de proa se trata como una discontinuidad nítida definida por splines cúbicos, en lugar de ser capturado por la malla.
  • Se imponen las condiciones de salto de Rankine-Hugoniot exactas en la interfaz choque-flujo.
  • Esto elimina artefactos numéricos (como el fenómeno "carbuncle") y asegura que la interacción entre perturbaciones infinitesimales y el choque reproduzca exactamente la respuesta predicha por el Análisis de Interacción Lineal (LIA).
• Modelos Termodinámicos y Químicos:
  • Incluye múltiples modelos para efectos de gas real en regímenes de alta entalpía: Gas perfecto calórico (CPG), química congelada, equilibrio químico y térmico (RTV, RTVE) y no equilibrio químico (NonEq-RTVE).
  • Utiliza propiedades termodinámicas efectivas ($\gamma^*$, $c_v^*$) derivadas de ajustes de funciones de base radial (RBF) para acelerar los cálculos en comparación con resolver ecuaciones de equilibrio en cada celda.
• Análisis de Estabilidad Lineal:
  • Análisis Modal Global: Resuelve el problema de autovalores para identificar modos inestables. Utiliza una transformación espectral exponencial ($\hat{L} \approx \exp(L\tau)$) en lugar del método "shift-and-invert" tradicional, lo que preserva la dispersidad de la matriz y reduce drásticamente los requisitos de memoria, permitiendo el uso de GPUs.
  • Crecimiento Transitorio (No Modal): Calcula la ganancia de energía óptima utilizando la norma de energía de Chu, adaptada para incluir química de equilibrio. Identifica perturbaciones iniciales que maximizan la amplificación de energía en un tiempo dado.
  • Receptividad del Flujo Libre: Modela la interacción entre perturbaciones del flujo libre y el flujo post-choque como un sistema entrada-salida, considerando la transmisión a través del choque y la amplificación posterior.
• Implementación Numérica:
  • Discretización de volumen finito de segundo orden en mallas curvilíneas.
  • Integración temporal explícita (RK4) o implícita (Euler hacia atrás) para estados base.
  • Linearización automática de los operadores discretos mediante diferencias finitas.

3. Contribuciones Clave

1. Integración Unificada: Es el primer paquete de código abierto (licencia MIT) que combina análisis modal, no modal y de receptividad en un solo marco global para flujos hipersónicos de alta entalpía.
2. Tratamiento Exacto de Choques: La implementación de choque ajustado (shock-fitting) garantiza la precisión en la interacción choque-perturbación, algo crítico para la receptividad que los métodos de captura de choque no pueden ofrecer.
3. Eficiencia Computacional y Escalabilidad:
   • Uso de transformaciones espectrales que evitan el llenado (fill-in) de matrices densas, permitiendo análisis en mallas grandes.
   • Soporte nativo para GPU (NVIDIA CUDA) en los solvers de autovalores y crecimiento transitorio, logrando aceleraciones significativas.
   • Implementación dual en MATLAB y Julia con funcionalidad idéntica.
4. Modelado de Gas Real: Inclusión de modelos termodinámicos avanzados (equilibrio y no equilibrio) necesarios para simulaciones de reentrada atmosférica y vuelos a Marte.

4. Resultados y Validación

El código ha sido verificado y validado contra una amplia colección de casos de referencia:

• Geometrías 2D y 3D: Comparación favorable con soluciones de Navier-Stokes inviscidos para choques sobre cilindros y esferas (Mach 2.5 a 8.06), mostrando excelente acuerdo en perfiles de presión y densidad.
• Efectos de Gas Real: Validación de los modelos de no equilibrio químico (NonEq-RTVE) contra el código Eilmer, reproduciendo correctamente los picos de temperatura translacional-rotacional y la relajación química detrás del choque.
• Condiciones de Salto: Verificación de las condiciones de Rankine-Hugoniot con efectos termodinámicos contra la Shock and Detonation Toolbox.
• Estabilidad de Flujos Clásicos: Recuperación precisa de autovalores inestables en flujos de canal incompresible y bifurcaciones de Hopf en cavidades de arrastre de tapa (lid-driven cavity), tanto incompresibles como compresibles.
• Interacción Choque-Perturbación: Confirmación de que el solver reproduce la conversión de modos (entropía a acústica) predicha analíticamente por la LIA a través de un choque 1D.
• Caso de Aplicación: Análisis de un cono romo hipersónico (Mach 12, Re=100,000).
  • El flujo resultó ser modalmente estable (todos los autovalores con parte real negativa).
  • El análisis de crecimiento transitorio mostró una amplificación de energía de ~40, dominada por contribuciones entrópicas y cinéticas, asociada al mecanismo Orr en la capa límite.
  • El análisis de receptividad reveló una amplificación masiva (~30,000) debido a la combinación de la amplificación a través del choque (escala $M_\infty^2$) y la interacción con el gradiente de presión en la región de estancamiento.

5. Significado e Impacto

HYMOR llena un vacío crítico en la comunidad de investigación de hipersónicos al proporcionar una herramienta accesible y rigurosa para estudiar la transición de flujo.

• Reproducibilidad: Al ser de código abierto, fomenta la reproducibilidad en la investigación de transición hipersónica.
• Precisión Física: Su enfoque de choque ajustado es fundamental para estudios de receptividad donde la física de la interacción choque-perturbación es dominante.
• Escalabilidad: La eficiencia en GPU y la capacidad de manejar mallas grandes permiten realizar análisis globales en geometrías realistas que antes eran computacionalmente prohibitivas.
• Base para Futuros Desarrollos: Sirve como plataforma base para que la comunidad desarrolle nuevos modelos físicos y algoritmos de control de flujo.

En resumen, HYMOR representa un avance significativo en la capacidad de predecir y entender los mecanismos de transición en vehículos hipersónicos de alta entalpía, superando las limitaciones de los métodos locales y ofreciendo una precisión sin precedentes en la modelización de interacciones de choque.