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

Imagine you are building a spacecraft designed to dive back into Earth's atmosphere from deep space. As it plummets at speeds faster than 10 kilometers per second (that's about 22,000 miles per hour!), the air in front of it doesn't just get hot; it gets so hot that it turns into a glowing, super-charged soup of particles called plasma.

This isn't just a hot wind; it's a blindingly bright fireball. This fireball emits massive amounts of radiation (light and heat energy), which can cook the spacecraft from the outside in, potentially melting it before it even touches the ground.

The paper you provided introduces a new computer program called RAPRAL. Think of RAPRAL as a highly sophisticated "heat-prediction crystal ball" for engineers. Its job is to calculate exactly how much of this blinding radiation will hit the back of the spacecraft, so designers can build shields thick enough to survive.

Here is a breakdown of how RAPRAL works, using simple analogies:

1. The Problem: It's Too Complicated to Guess

When a spacecraft re-enters the atmosphere, the air molecules smash together so hard that they break apart (dissociate) and even rip electrons off atoms (ionize). This creates a chaotic mix of different particles (like Nitrogen, Oxygen, and their ionized versions) moving at different speeds.

To predict the heat, you can't just use a simple thermometer. You have to understand:

The "Music" of the Atoms: Every time an electron jumps between energy levels in an atom, it emits a specific color of light (a spectral line). It's like every atom is a tiny instrument playing a specific note.
The "Traffic" of Light: This light doesn't just travel in a straight line; it gets absorbed and re-emitted by other particles as it travels through the hot gas cloud.

2. The Solution: Two Superpowers Combined

RAPRAL combines two powerful methods to solve this puzzle:

A. The "Line-by-Line" Method (The Microscope)

Imagine trying to listen to a choir where every singer is singing a slightly different note. A simple method might just say, "It's loud." But RAPRAL acts like a super-microscope. It looks at every single note (spectral line) that every single atom and molecule can sing.

It calculates the exact "volume" (intensity) of each note.
It accounts for how the "notes" get blurry due to heat and collisions (like how a siren sounds different when a car zooms past you).
Why this matters: If you miss even a few notes, you might underestimate the total heat by a huge amount.

B. The "Ray Tracing" Method (The Flashlight)

Once RAPRAL knows what notes are being sung, it needs to figure out how that sound travels to the back of the spaceship.

Imagine standing in a foggy room with a flashlight. You shine the light in thousands of different directions.
RAPRAL shoots thousands of "virtual light beams" (rays) from the hot gas cloud toward the back of the spacecraft.
As these rays travel, the program checks: "Did this ray hit a particle that absorbed it? Did it get brighter because a particle emitted light?"
By adding up all these rays, it calculates the total heat hitting a specific spot on the ship.

3. The "Traffic Jam" Analogy (Non-Equilibrium)

One of the hardest parts of this physics is that the gas isn't "calm."

The Scenario: Imagine a highway where cars (atoms) are suddenly forced to speed up by a shockwave. The engines (electronic energy) rev up instantly, but the tires (vibrational energy) take time to heat up.
The Result: For a while, the engines are screaming hot, but the tires are still cool. This is called Thermochemical Nonequilibrium.
RAPRAL's Job: It doesn't assume everything is at the same temperature. It tracks the "hot engines" and "cool tires" separately, which is crucial for getting the heat prediction right. If you assume they are the same, your prediction will be wrong.

4. The Test Drive: The Fire II Experiment

To prove RAPRAL works, the authors tested it against a real-life historical event: the Fire II mission from the 1960s.

The Setup: A scaled-down Apollo capsule was dropped from a rocket to simulate a lunar return. It had sensors to measure how much heat hit its back.
The Result: RAPRAL took the data from the flight and ran its simulation.
- Success: It predicted the heat flux (heat flow) very accurately for one part of the flight.
- The Glitch: For another part, it underestimated the heat. The authors realized this was likely because the simulation didn't account for the "sweat" of the spacecraft (ablation products). When the heat shield burns away, it releases gas that actually makes the radiation hotter. RAPRAL didn't know about this "sweat" yet, so it predicted less heat than actually happened.

5. Why This Matters

Current tools are either too slow (trying to calculate every single note takes forever) or too simple (guessing the notes). RAPRAL is a new, fast, and accurate tool written in modern code (C++) that can handle these complex calculations efficiently.

The Bottom Line:
RAPRAL is a new, high-tech calculator for space engineers. It uses a "microscope" to see every tiny detail of the light emitted by super-hot air and a "flashlight" to trace exactly where that heat goes. This helps ensure that when we send astronauts to Mars or bring back samples from asteroids, our heat shields are strong enough to keep them safe from the fiery sky.

Future Plans:
The creators plan to upgrade RAPRAL to handle even more complex atmospheres (like Mars, which is mostly Carbon Dioxide) and to include the "sweat" (ablation) effects so it never misses a beat again.

1. Problem Statement

In hypersonic atmospheric reentry, radiative heating becomes a dominant component of the total thermal load, particularly at velocities exceeding 10 km/s (e.g., lunar or Mars return missions). Accurate prediction of radiative heat flux is critical for the design of Thermal Protection Systems (TPS). However, existing radiation prediction tools often face limitations:

Spectral Accuracy vs. Efficiency: Many codes rely on band models which sacrifice spectral detail for speed, while high-fidelity Line-by-Line (LBL) methods are computationally expensive.
Nonequilibrium Modeling: Existing open-source tools often lack robust capabilities for handling complex thermochemical nonequilibrium (non-Boltzmann populations) in the Vacuum Ultraviolet (VUV), UV, and Visible spectral ranges.
Geometric Complexity: Traditional ray-tracing or tangent-slab approximations may introduce errors in complex flowfields, such as the expansion and recirculation regions behind a vehicle's shoulder (afterbody).

The authors developed RAPRAL (RAdiation Prediction based on RAy tracing and Line-by-line) to address these gaps, providing a high-fidelity, portable solver for simulating radiative processes in high-temperature air flows.

2. Methodology

RAPRAL is a C++-based solver that integrates detailed spectral modeling with a ray-tracing solution of the Radiative Transfer Equation (RTE). The methodology is structured into three main layers:

A. Collisional-Radiative (CR) Module (Population Distributions)

To determine the population of internal energy states (essential for calculating emission/absorption), RAPRAL employs:

Equilibrium: Uses Boltzmann distributions for electronic states and Saha-Boltzmann distributions for ionization equilibrium.
Nonequilibrium: Implements a Collisional-Radiative (CR) model under the Quasi-Steady-State (QSS) approximation. This solves a master equation balancing source and sink terms for various mechanisms:
- Heavy-particle and electron impact excitation/dissociation.
- Radiative transitions (bound-bound).
- Electron impact ionization and charge exchange.
Molecular States: Accounts for electronic, vibrational, and rotational modes. It utilizes Hund's coupling cases (a, b, and transitions a→b) to calculate rotational energy levels and Hönl-London factors for diatomic molecules (e.g., $N_2$ , $NO$, $N_2^+$ ).

B. Radiative Module (Spectral Properties)

The solver calculates emission ( $j_\lambda$ ) and absorption ( $\kappa_\lambda$ ) coefficients using a Line-by-Line (LBL) approach:

Atomic Radiation: Includes Bound-Bound (B-B), Bound-Free (B-F), and Free-Free (F-F) transitions.
- Line Broadening: Uses a Voigt profile combining Doppler (thermal) and Lorentz (pressure/Stark/Resonance/Van der Waals) broadening. Stark broadening coefficients are updated based on recent literature.
- Continuum: Includes photo-ionization (B-F) and Bremsstrahlung (F-F) transitions, scaling hydrogenic data for other atomic species.
Molecular Radiation: Focuses on B-B transitions for major air species ( $N_2$ $N_{2}$ , $N_2^+$ $N_{2}^{+}$ , $NO$, $O_2$ $O_{2}$ ).
- Uses updated electronic transition moment functions (ETMF) and Franck-Condon factors from recent ab initio calculations (Liang et al.).
- Explicitly handles line alternation factors due to nuclear spin statistics.

C. Radiative Transfer Solver (Ray Tracing)

Algorithm: Solves the RTE along Lines of Sight (LOS) using a ray-tracing approach.
Discretization: Employs the Fibonacci sphere method to generate a uniform distribution of LOS directions in 3D space, avoiding the clustering issues of standard spherical coordinate discretization.
Integration: Interpolates flowfield variables (temperature, density) from structured CFD grids onto LOS nodes. Integrates the RTE using a linear variation assumption for emission and constant absorption between nodes, with special handling for optically thin limits to ensure numerical stability.
Portability: Features a unified interface allowing it to operate independently of specific CFD solvers.

3. Key Contributions

Development of RAPRAL V1.0: A new, open-architecture radiation solver implemented in C++ with Python utilities for database generation, specifically designed for hypersonic air flows.
High-Fidelity Spectral Modeling: Integration of updated line-broadening models (Stark broadening) and high-resolution LBL calculations for both atomic and molecular species, including complex nonequilibrium population kinetics via the CR model.
Robust Ray-Tracing Framework: Implementation of a 3D ray-tracing solver using the Fibonacci sphere method to accurately capture radiative transport in complex flowfields (e.g., vehicle afterbodies) where 1D approximations fail.
Comprehensive Validation: Rigorous testing against the open-source Spark code for spectral coefficients and against experimental data from the Fire II flight experiment for radiative heating.

4. Results

Spectral Coefficients Validation:
- RAPRAL was compared against the Spark code for atomic Nitrogen and major molecular bands ( $N_2$ , $NO$, $O_2$ ).
- RAPRAL showed higher peak intensities and sharper spectral features, attributed to the use of updated line-shape models (Johnston's Stark broadening) and a finer spectral grid.
- Discrepancies in integrated intensities were linked to differences in partition function evaluations and the selection of maximum rotational quantum numbers.
Fire II Flight Experiment Prediction:
- The solver was applied to the Fire II reentry trajectory (1637.5 s and 1640.5 s) using a 2-temperature, 11-species flowfield.
- Flowfield Analysis: Identified significant thermochemical nonequilibrium in the afterbody expansion region, where vibrational-electronic temperatures ( $T_{vee}$ ) remained high while translational temperatures ( $T_{tr}$ ) dropped rapidly.
- Heat Flux Prediction:
  - At t = 1640.5 s, RAPRAL provided reliable predictions of radiative heat flux that agreed well with experimental measurements.
  - At t = 1637.5 s, the solver underpredicted the radiative heat flux. The authors attribute this to the breakdown of the QSS assumption in the expansion region and the neglect of wall ablation products (which can contribute ~25% to radiative flux).
- Dominant Mechanisms: Analysis confirmed that atomic emission (VUV range) dominates the leeside radiation, while continuum radiation (B-F and F-F) contributes significantly (>50%) in highly ionized regions.

5. Significance

Reliability for TPS Design: RAPRAL demonstrates robust capability in simulating radiative processes in hypersonic flows, offering a reliable tool for thermal protection system design, particularly for high-enthalpy reentry scenarios.
Advancement in Nonequilibrium Modeling: By successfully coupling detailed CR kinetics with ray tracing, the code provides deeper insight into the spatial distribution of radiation in complex flowfields, moving beyond simplified 1D assumptions.
Foundation for Future Work: The modular architecture allows for future extensions to include planetary atmospheres (e.g., $CO_2$ , $CH_4$ ) and molecular continuum radiation, making it a versatile platform for deep-space exploration missions.

In conclusion, RAPRAL represents a significant step forward in high-fidelity radiative heat transfer prediction, balancing computational efficiency with the physical accuracy required for next-generation hypersonic vehicle design.