RAPRAL v1.0: RAdiation Prediction using RAy tracing and Line-by-line methods for hypersonic air flows

El artículo presenta RAPRAL v1.0, un nuevo solver de radiación en C++ que combina el modelado espectral línea por línea con trazado de rayos para simular con precisión procesos radiativos en flujos hipersónicos, validando su capacidad mediante la comparación con códigos existentes y la predicción del calentamiento radiativo en el experimento de vuelo Fire II.

Lo esencial

¡Hola! Imagina que estás diseñando un vehículo espacial que viaja tan rápido que el aire a su alrededor se convierte en fuego líquido. Cuando una nave regresa a la Tierra desde el espacio (como las misiones a Marte o trayendo muestras de asteroides), viaja a velocidades increíbles (más de 10 km/s). En ese momento, el aire no solo se calienta por fricción, sino que se "rompe" en sus partes más pequeñas (átomos y electrones) y emite una luz brillante y peligrosa.

Este documento presenta RAPRAL, un nuevo programa informático creado por científicos chinos para predecir exactamente cuánto calor radiante recibirá la nave.

Aquí tienes la explicación sencilla, usando analogías de la vida diaria:

1. El Problema: El "Fuego" Invisible

Cuando la nave entra a la atmósfera, el aire se comprime y se calienta tanto que se convierte en un plasma (gas ionizado). Este plasma emite luz y calor, como una estufa encendida.

• La analogía: Imagina que la nave es un coche que entra en un túnel de viento súper caliente. No solo el viento quema (calor por fricción), sino que el aire mismo se vuelve una bombilla gigante que irradia calor hacia el coche. Si no calculamos bien cuánto calor emite esa "bombilla", el escudo térmico de la nave podría derretirse.

2. La Solución: RAPRAL (El "Detective de la Luz")

Los autores crearon un software llamado RAPRAL. Su nombre es un acrónimo que significa "Predicción de Radiación usando Trazado de Rayos y Método Línea por Línea".

Imagina que RAPRAL es un detective muy detallista que tiene dos herramientas principales:

A. La herramienta "Línea por Línea" (El Microscopio)

El aire caliente no emite luz de un solo color uniforme. Emite millones de colores específicos, como si fuera una orquesta tocando millones de notas diferentes al mismo tiempo.

• La analogía: Si usas un prisma para descomponer la luz del sol, ves un arcoíris. El aire caliente es como un arcoíris supercomplejo con millones de líneas finas.
• Qué hace RAPRAL: En lugar de mirar el arcoíris de lejos (como hacen otros programas más simples), RAPRAL toma un microscopio y mira cada una de esas líneas individuales. Cuenta cuánta energía hay en cada nota musical. Esto le permite ver detalles que otros programas pierden, como si distinguiera entre un violín y una viola en medio de una sinfonía.

B. La herramienta "Trazado de Rayos" (El Mapa de Luces)

Una vez que sabe qué colores emite el aire, necesita saber cómo viaja esa luz hasta la nave.

• La analogía: Imagina que estás en una habitación llena de niebla brillante y quieres saber cuánta luz llega a tu cara. No basta con saber qué tan brillante es la niebla; tienes que trazar líneas imaginarias desde tu cara hacia todas las direcciones posibles para ver qué trozos de niebla te están "mirando" y enviando luz.
• Qué hace RAPRAL: Dispara miles de "rayos láser" imaginarios desde la superficie de la nave hacia el aire caliente. Sigue cada rayo, calcula cuánto calor absorbe o emite el aire en cada punto, y suma todo el calor que llega a la nave.

3. ¿Cómo funciona por dentro? (El Motor del Coche)

El programa simula dos cosas principales:

1. Los átomos y moléculas: Cuando el aire se calienta, los átomos de nitrógeno y oxígeno se excitan (saltan a niveles de energía más altos) y luego caen, soltando luz. RAPRAL calcula esto con una precisión matemática extrema, considerando si el aire está en equilibrio o si está "desordenado" (fuera de equilibrio), lo cual es muy común en velocidades supersónicas.
2. La geometría: Usa un método llamado "Esfera de Fibonacci" para disparar sus rayos de luz de manera uniforme en todas direcciones, como si estuviera llenando una esfera con puntos perfectamente distribuidos para no dejar ningún ángulo sin revisar.

4. La Prueba de Fuego: El Experimento "Fire II"

Para ver si RAPRAL funciona de verdad, los científicos lo pusieron a prueba contra un experimento real de los años 60 llamado Fire II.

• La analogía: Fue como tomar un mapa nuevo y compararlo con un viaje real que ya se hizo.
• El resultado: El programa logró predecir con gran precisión cuánto calor radiante golpeaba la parte trasera de la nave (el "afterbody"). En algunos momentos, sus predicciones coincidieron casi perfectamente con los sensores reales de la misión.
• El detalle: En un momento específico, el programa predijo un poco menos de calor del que se midió. Los autores explican que esto podría ser porque el programa no contó con los gases que salen del escudo térmico de la nave (como si el coche estuviera soltando humo que también brilla), algo que se corregirá en futuras versiones.

5. ¿Por qué es importante?

Hoy en día, queremos ir a Marte y traer muestras de asteroides. Estas misiones viajan más rápido que las misiones anteriores a la Luna.

• La conclusión: Si no calculamos bien el calor de la "bombilla" de plasma, la nave se quemará. RAPRAL es una herramienta nueva, más precisa y flexible (escrita en un lenguaje moderno llamado C++) que ayuda a los ingenieros a diseñar escudos térmicos más seguros y ligeros.

En resumen:
RAPRAL es como un simulador de realidad virtual ultra-realista para el calor. En lugar de adivinar cuánto calor recibe una nave, lo "ve" píxel por píxel y rayo por rayo, asegurando que cuando enviemos humanos o robots al espacio profundo, vuelvan sanos y salvos.

Resumen técnico

Aquí presento un resumen técnico detallado del artículo "RAPRAL V1.0: RADIATION PREDICTION USING RAY TRACING AND LINE-BY-LINE METHODS FOR HYPERSONIC AIR FLOWS", estructurado según los puntos solicitados.

1. El Problema

En los flujos hipersónicos, particularmente durante la reentrada atmosférica a altas velocidades (superiores a 9-10 km/s), el calentamiento radiativo se convierte en un componente crítico, y en algunos casos dominante, de la carga térmica total sobre las naves espaciales. A diferencia del calentamiento convectivo, la radiación depende fuertemente de la velocidad y la escala de longitud, y su predicción requiere resolver procesos termoquímicos fuera de equilibrio en condiciones de alta entalpía.

Los desafíos principales incluyen:

• Complejidad Espectral: La necesidad de modelar transiciones atómicas y moleculares (enlazadas, libres-enlazadas y libres-libres) con alta precisión en el rango espectral VUV/UV/Vis.
• No Equilibrio Termodinámico: En la estela de la nave (región de expansión), las poblaciones de estados electrónicos, vibracionales y rotacionales no siguen distribuciones de Boltzmann locales, lo que invalida modelos simplificados.
• Limitaciones de Herramientas Existentes: Muchos códigos de código abierto existentes (como RADIS o SPARK) tienen limitaciones en la portabilidad, la eficiencia computacional para análisis independientes, o la capacidad de manejar poblaciones no equilibradas de manera integrada en flujos complejos.

2. Metodología

Los autores desarrollaron RAPRAL (RAdiation Prediction based on RAy tracing and Line-by-line), un nuevo solucionador de radiación implementado en C++ con utilidades en Python. La metodología se basa en tres pilares fundamentales:

A. Módulo de Colisiones y Poblaciones de Estados

• Modelado de Estados: Se utilizan distribuciones de Boltzmann para condiciones de equilibrio local (LTE). Para condiciones de fuerte no equilibrio, se emplea un modelo Colisional-Radiativo (CR) bajo la aproximación de estado cuasi-estacionario (QSS).
• Mecanismos: El modelo CR resuelve la ecuación maestra considerando procesos elementales como excitación/impacto de partículas pesadas y electrones, disociación, ionización, recombinación y transiciones radiativas.
• Especies: Se modelan especies atómicas (N, O) y diatómicas (N2, O2, NO, N2+, etc.). Para moléculas, se consideran modos electrónicos, vibracionales y rotacionales, utilizando distribuciones de Boltzmann caracterizadas por temperaturas vibracionales ($T_v$) y rotacionales ($T_r$) distintas de la temperatura traslacional ($T_t$).

B. Módulo Radiativo (Espectroscopía)

• Método Línea por Línea (LBL): Se adopta el método LBL para calcular coeficientes de emisión y absorción con alta precisión, evitando las aproximaciones de modelos de bandas.
• Transiciones Atómicas: Se incluyen transiciones enlazadas-enlazadas (B-B), libres-enlazadas (B-F) y libres-libres (F-F). Se utilizan perfiles de línea Voigt para combinar el ensanchamiento Doppler y de presión (Stark, resonancia, Van der Waals).
• Transiciones Moleculares: Se calculan coeficientes para sistemas de bandas diatómicas (ej. N2 First/Second Positive, N2+ First Negative, NO, O2). Se utilizan factores de Hönl-London y constantes de acoplamiento de Hund (casos a, b y transiciones a→b) para determinar la estructura rotacional.
• Datos: Se utilizan bases de datos actualizadas (NIST, TOPbase) y momentos de transición electrónica recientes (Liang et al., Chauveau et al.).

C. Módulo de Transferencia Radiativa (Ray Tracing)

• Ecuación de Transferencia Radiativa (RTE): Se resuelve la RTE en flujos estructurados utilizando un enfoque de trazado de rayos.
• Discretización Angular: Se utiliza el método de la esfera de Fibonacci para generar direcciones de rayos (líneas de visión, LOS) uniformemente distribuidas en el espacio 3D, evitando la agrupación no uniforme típica de la discretización esférica estándar.
• Integración: La RTE se integra a lo largo de cada LOS desde el campo de flujo hasta un punto objetivo en la pared. Se asume que la emisividad varía linealmente y la absorbancia es constante entre nodos adyacentes para equilibrar precisión y costo computacional.
• Interfaz: El código es independiente del solucionador de dinámica de fluidos computacional (CFD), permitiendo su uso como una herramienta de post-procesamiento.

3. Contribuciones Clave

1. Desarrollo de RAPRAL V1.0: Un solucionador robusto en C++ que integra modelado espectral detallado (LBL) con un esquema de trazado de rayos eficiente para flujos hipersónicos.
2. Implementación de No Equilibrio: Capacidad para manejar poblaciones de estados electrónicos no equilibradas mediante modelos colisionales-radiativos (CR) y distribuciones de Boltzmann no locales para modos vibracionales/rotacionales.
3. Optimización Computacional: Uso de paralelización OpenMP para la evaluación de propiedades espectrales y cálculos de transferencia radiativa, junto con algoritmos eficientes de búsqueda de rejilla e interpolación.
4. Validación Exhaustiva: Comparación directa con el código de código abierto SPARK y datos experimentales del vuelo Fire II, demostrando la capacidad del código para capturar características espectrales finas y flujos de calor radiativos complejos.

4. Resultados

• Comparación Espectral (RAPRAL vs. SPARK):
  • Se observaron diferencias en la intensidad de pico y la forma de las bandas rovibracionales.
  • Las discrepancias se atribuyeron principalmente al uso de una rejilla espectral más fina, modelos de ensanchamiento de línea actualizados (Stark) y diferentes tratamientos del número cuántico rotacional máximo.
  • RAPRAL mostró una mejor resolución de las características espectrales finas.
• Predicción del Experimento Fire II:
  • Se simuló el calentamiento en la parte trasera (afterbody) de la cápsula Fire II en dos instantes de tiempo (1637.5 s y 1640.5 s).
  • A 1640.5 s: El código predijo con precisión el flujo de calor radiativo, coincidiendo bien con las mediciones experimentales.
  • A 1637.5 s: Se subestimó significativamente el flujo de calor radiativo en la región de la parte trasera.
  • Análisis de la discrepancia: Los autores atribuyen la subestimación a dos factores principales:
    1. La ruptura de la aproximación de estado cuasi-estacionario (QSS) en la región de expansión, donde las poblaciones de estados excitados están "congeladas" desde la región frontal de alta temperatura y no evolucionan localmente.
    2. La omisión de los productos de ablación de la pared, que pueden contribuir con hasta un 25% adicional al flujo de calor radiativo.
  • Se identificó que la radiación atómica (N, O) en el rango VUV domina la contribución en la estela, y que la radiación de continuo (B-F y F-F) contribuye con más del 50% en ciertas líneas de visión.

5. Significancia

El trabajo presenta una herramienta avanzada y flexible para la predicción de calentamiento radiativo en misiones de reentrada hipersónica y exploración planetaria.

• Precisión: Al utilizar el método LBL y trazado de rayos, RAPRAL supera las limitaciones de los modelos de bandas aproximados, ofreciendo una resolución espectral y espacial superior.
• Aplicabilidad: Su diseño modular y su independencia del solucionador CFD lo hacen ideal para el análisis de sensibilidad y la validación de modelos de física en diferentes entornos de simulación.
• Futuro: El código sienta las bases para futuras extensiones que incluirán radiación de continuo molecular (B-F y F-F) y especies relevantes para atmósferas planetarias (como CO2 y CH4), lo cual es crucial para el diseño de sistemas de protección térmica (TPS) en misiones de retorno de muestras de Marte o asteroides.

En resumen, RAPRAL demuestra ser una herramienta robusta para simular procesos radiativos en flujos de aire hipersónicos, proporcionando predicciones confiables que son esenciales para la seguridad y el diseño de vehículos de reentrada de alta velocidad.