Fast Prediction of Aero-Thermal Performance During Hypersonic Vehicle Design

Resumen Técnico: Predicción Rápida del Rendimiento Aero-Térmico Durante el Diseño de Vehículos Hipersónicos

Planteamiento del Problema
El diseño de vehículos hipersónicos requiere una estimación rápida de las fuerzas aerodinámicas y los flujos de calor superficial durante las etapas tempranas de desarrollo. Si bien las simulaciones de Dinámica de Fluidos Computacional (CFD) de alta fidelidad ofrecen una precisión detallada, resultan computacionalmente prohibitivas para las extensas iteraciones de diseño y optimizaciones de trayectoria requeridas en las fases preliminares. Los métodos tradicionales de nivel de ingeniería, como las técnicas de inclinación de superficie local (cuña tangente, Newton modificado, Prandtl-Meyer), proporcionan velocidad pero a menudo carecen de fidelidad al modelar los efectos viscosos. Específicamente, los enfoques estándar suelen asumir que la capa límite es completamente laminar o completamente turbulenta. La suposición de un flujo totalmente turbulento, aunque es conservadora, frecuentemente conduce a una sobrepredicción del arrastre por fricción superficial y del flujo de calor superficial. Esto resulta en Sistemas de Protección Térmica (TPS) sobredimensionados, mayor masa del vehículo y mayor consumo de combustible. Además, pocas herramientas de ingeniería rápidas incorporan actualmente los efectos de transición de la capa límite dentro de los marcos de inclinación local, a pesar del impacto significativo de los flos transicionales en las cargas aerotérmicas.

Metodología
El artículo presenta HipeX, una herramienta de análisis aerotermodinámico ligera diseñada para cerrar la brecha entre los métodos analíticos clásicos y la CFD de alta fidelidad. HipeX utiliza un enfoque híbrido de nivel de ingeniería basado en métodos de inclinación de superficie local para estimar las propiedades del flujo de entrada y emplea un modelo de capa límite mezclado para abordar los efectos viscosos y la transición.

• Modelado de Flujo Inviscido: El solver procesa mallas de superficie no estructuradas (formato STL) y selecciona el modelo de flujo apropiado basándose en la inclinación de la superficie local con respecto al vector de velocidad incidente:

  • Teoría de Newton Modificada: Aplicada para choques desprendidos (cuerpos romos) y ángulos de deflexión grandes.
  • Relaciones de Choque Oblicuo: Utilizadas para choques adheridos y regiones de compresión.
  • Expansión de Prandtl-Meyer: Aplicada para regiones de expansión de flujo, con una salvaguarda que limita su aplicación al 50% del ángulo de deflexión máximo teórico para asegurar la estabilidad numérica.
  • Fórmula de Gaubeaud: Utilizada como un respaldo empírico cuando las soluciones de Prandtl-Meyer fallan.
  • Las propiedades atmosféricas se determinan mediante modelos estándar basados en la altitud de vuelo.

• Modelado Viscoso y Térmico:

  • Compresibilidad: Corregida mediante el método de temperatura de referencia de Eckert, aplicable tanto a regímenes laminares como turbulentos.
  • Modelado de Transición: Se calcula empíricamente un número de Reynolds crítico ($Re_{cr}$) como función del número de Mach para determinar el inicio de la transición. El modelo emplea un factor de intermitencia ($\gamma_{it}$) que oscila entre 0 (laminar) y 1 (turbulento) para mezclar los coeficientes de fricción superficial de los flujos laminares y turbulentos sobre una región de transición definida.
  • Cálculo de Flujo de Calor: El flujo de calor superficial se deriva del coeficiente de fricción superficial utilizando la analogía de Reynolds y el número de Stanton. El modelo tiene en cuenta tanto la transferencia de calor por convección como la transferencia de calor por radiación (asumiendo una emisividad de 0.4), calculando el flujo de calor neto y la carga térmica total mediante la integración sobre el área mojada.

• Flujo de Trabajo: La herramienta acepta geometría 3D y condiciones de vuelo (Mach, altitud, Ángulo de Ataque, Ángo de Deriva). Calcula las propiedades inviscidas, aplica correcciones viscosas, integra las fuerzas para determinar los coeficientes aerodinámicos y exporta los resultados en formato VTK. El código admite la ejecución en paralelo para el análisis de trayectoria.

Contribuciones Clave

1. Modelo de Transición Mezclado: La principal contribución es la integración de un modelo de capa límite transicional en un marco de inclinación local rápido. Esto permite al solver estimar el cambio de un flujo laminar a uno turbulento basándose en los números de Reynolds locales y los números de Mach, en lugar de depender de suposiciones binarias de laminar/turbulento.
2. Eficiencia Computacional: HipeX está diseñado para una ejecución rápida, capaz de computar el rendimiento de una trayectoria completa (aprox. 300 puntos) en una estación de trabajo estándar en aproximadamente un minuto.
3. Validación contra Datos de Referencia: La herramienta es validada mediante datos de CFD y referencias experimentales para la configuración HEXAFLY-INT (un diseño de cuerpo de sustentación mezclada para un crucero de Mach 7–8).

Resultados
La validación se realizó utilizando la configuración HEXAFLY-INT a través de varios números de Mach (M=2 a M=7.5) y trayectorias de vuelo.

• Coeficientes Aerodinámicos: Para números de Mach altos (M=6 y M=7), las predicciones de HipeX para la sustentación ($C_L$), el arrastre ($C_D$) y el centro de presión ($X_{cp}$) se alinean estrechamente con los datos de CFD. La solución transicional mezclada generalmente proporciona resultados más precisos que las suposiciones puramente laminares o totalmente turbulentas. A números de Mach más bajos (M=4 y M=2), las discrepancias aumentan, particularmente en los valores absolutos de $C_L$ y $C_D$, aunque la relación sustentación/arrastre ($L/D$) se mantiene razonablemente precisa. La herramienta identifica correctamente los límites de los métodos de inclinación local cuando el flujo se aleja de las condiciones fuertemente hipersónicas.
• Cargas Térmicas: En casos donde el flujo es predominantemente laminar (alta altitud), la solución mezclada produce resultados similares al caso totalmente laminar. Sin embargo, en regímenes transicionales, el modelo mezclado captura niveles de flujo de calor intermedios que se encuentran entre los extremos laminar y turbulento.
• Implicaciones del TPS: El estudio demuestra que asumir un flujo totalmente turbulento conduce a una sobrepredicción significativa del calentamiento, particularmente cerca del borde de ataque. Al capturar la naturaleza laminar o parcialmente turbulenta de la capa límite sobre gran parte de la trayectoria, HipeX proporciona estimaciones de flujo de calor más bajas y realistas. Esto sugiere el potencial para reducir la masa del TPS sin comprometer los márgenes de seguridad.
• Inicio de la Transición: Las comparaciones con el código BLITZ muestran tendencias similares en la ubicación de la transición, aunque HipeX tiende a predecir un inicio de turbulencia más temprano. El modelo actualmente tiene dificultades para evitar la predicción de relaminarización en ciertas regiones de expansión.

Significancia y Reivindicaciones
El artículo posiciona a HipeX como una base robusta para la evaluación aerotermodinámica y la optimización de diseño en etapas tempranas. Su significancia radica en ofrecer un "equilibrio balanceado" entre el conservadurismo del diseño y la eficiencia de la misión. Al alejarse de la suposición excesivamente conservadora de un flujo totalmente turbulento, la herramienta permite un dimensionamiento más preciso de los sistemas de protección térmica, reduciendo potencialmente la masa del vehículo y el consumo de combustible.

Los autores reconocen que HipeX hereda las limitaciones intrínsecas de los métodos de inclinación de superficie local, específicamente la incapacidad de predecir interacciones de choque, zonas de separación y efectos tridimensionales fuertes. Además, la determinación precisa del inicio de la transición en condiciones reales de vuelo sigue siendo un desafío debido a las dependencias de la rugosidad de la superficie y las perturbaciones del flujo libre. Sin embargo, la herramienta logra llenar el vacío entre los métodos analíticos clásicos y la CFD de alta fidelidad, permitiendo la evaluación rápida de trayectorias completas y grandes conjuntos de configuraciones geométricas esenciales para el diseño preliminar.