RAPRAL v1.0: RAdiation Prediction using RAy tracing and Line-by-line methods for hypersonic air flows
Pith reviewed 2026-05-10 16:00 UTC · model grok-4.3
The pith
RAPRAL combines line-by-line spectra with ray tracing to predict radiative heat flux in hypersonic air flows.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
RAPRAL integrates detailed line-by-line spectral modeling with a ray-tracing solution of the radiative transfer equation. When applied to predict afterbody radiative heating in the Fire II flight experiment based on a two-temperature, 11-species air flowfield, the approach yields reliable radiative heat flux predictions and captures the dominant radiation mechanisms.
What carries the argument
Line-by-line computation of atomic and molecular spectral coefficients combined with ray tracing to integrate the radiative transfer equation over space and spectrum.
If this is right
- Bulk spectral coefficients for air species can be computed accurately over a wide range of temperatures and pressures.
- Dominant radiation mechanisms in afterbody regions of hypersonic flows are identified and quantified.
- The solver serves as a reliable tool for radiative heating predictions in Earth-atmosphere hypersonic entries.
- Planned extensions will incorporate additional species for planetary entry simulations.
Where Pith is reading between the lines
- The same framework could be tested on other flight data sets to check whether spectral agreement continues to produce correct total heat loads.
- Coupling RAPRAL directly to a flow solver instead of using a precomputed field might reveal sensitivity of heating predictions to flowfield details.
- The ray-tracing approach may allow efficient parallelization for three-dimensional vehicle geometries.
Load-bearing premise
The supplied two-temperature, 11-species air flowfield is accurate, and matching isolated spectral coefficients guarantees correct integrated radiative heating.
What would settle it
Measurement of afterbody radiative heat flux on a different reentry vehicle where an independent, high-fidelity flowfield solution is available for direct comparison.
Figures
read the original abstract
A new radiation solver, RAPRAL (RAdiation Prediction based on RAy tracing and Line-by-line) implemented in C++, is developed for simulating high-temperature thermochemical nonequilibrium radiative processes. RAPRAL integrates detailed line-by-line spectral modeling with a ray-tracing solution of the radiative transfer equation, enabling accurate resolution of both spectral features and spatial radiation transport. The adopted methods and their implementation are described in detail. To assess the overall capability and accuracy of RAPRAL, we first focus on the computation of atomic and molecular bulk spectral coefficients. Through comparison with the established code in the literature, RAPRAL demonstrates its ability to accurately capture key spectral features across a wide range of conditions. Moreover, RAPRAL is applied to predict afterbody radiative heating in the Fire II flight experiment, based on a two-temperature, 11-species air flowfield. The results demonstrate that the present approach provides reliable predictions of radiative heat flux and effectively captures the dominant radiation mechanisms. Overall, the presented results demonstrate that RAPRAL is a robust tool for simulating radiative processes in hypersonic air flows, and future versions will extend its capabilities to include species relevant to planetary atmospheres.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces RAPRAL v1.0, a C++ radiation solver combining line-by-line spectral modeling with ray-tracing solution of the radiative transfer equation for thermochemical nonequilibrium hypersonic air flows. It reports qualitative agreement of atomic and molecular bulk spectral coefficients with established literature codes across a range of conditions, then applies the solver to afterbody radiative heating predictions for the Fire II flight experiment using an externally supplied two-temperature, 11-species air flowfield, concluding that the approach yields reliable radiative heat flux predictions and captures dominant mechanisms.
Significance. If the integrated flux predictions can be shown to be reliable, RAPRAL would provide a useful new open-source tool for detailed radiation modeling in hypersonic flows, with its explicit C++ implementation and combined spectral-spatial treatment offering a clear alternative to existing solvers. The detailed description of methods is a strength, but the current evidence base does not yet establish the reliability of the integrated quantities.
major comments (2)
- [Abstract and Fire II results section] Abstract and Fire II results section: The claim that 'the present approach provides reliable predictions of radiative heat flux' for Fire II afterbody heating is not supported by the presented results. Only computed flux values are reported; no quantitative comparison to flight data or to benchmark codes (e.g., NEQAIR) on identical input flowfields is given. Because integrated flux depends on ray-tracing through the full spatial domain, agreement on isolated spectral coefficients does not establish accuracy of the integrated heating.
- [Spectral validation section] Spectral validation section: While qualitative matches to literature codes on spectral features are shown, the manuscript provides no quantitative error metrics (e.g., RMS or percentage differences) or sensitivity studies to uncertainties in the supplied 2T/11-species flowfield. This is load-bearing for the reliability claim, as the radiation module operates on externally supplied flow data whose accuracy directly determines the predicted fluxes.
minor comments (2)
- [Abstract] The abstract could more precisely state the specific conditions and spectral features used in the code-to-code comparisons.
- [Implementation and numerical methods] A table summarizing numerical parameters (e.g., ray discretization, spectral resolution) used in the ray-tracing and line-by-line calculations would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive review of our manuscript on RAPRAL v1.0. We address the major comments point by point below, indicating where revisions will be made to improve the presentation and clarify the scope of the validation.
read point-by-point responses
-
Referee: [Abstract and Fire II results section] Abstract and Fire II results section: The claim that 'the present approach provides reliable predictions of radiative heat flux' for Fire II afterbody heating is not supported by the presented results. Only computed flux values are reported; no quantitative comparison to flight data or to benchmark codes (e.g., NEQAIR) on identical input flowfields is given. Because integrated flux depends on ray-tracing through the full spatial domain, agreement on isolated spectral coefficients does not establish accuracy of the integrated heating.
Authors: We agree that the phrasing 'provides reliable predictions of radiative heat flux' in the abstract and Fire II section overstates what the presented results demonstrate. The Fire II case is included as an application example using an externally supplied flowfield to illustrate the solver's use on a realistic problem and to identify dominant mechanisms, rather than as a comprehensive validation of integrated quantities. We will revise the abstract and results section to remove or qualify this claim, emphasizing instead the demonstration of the method and its ability to capture key features. We will also add a limitations paragraph noting that quantitative benchmarking of integrated fluxes against flight data or other codes on identical flowfields is not performed here. revision: partial
-
Referee: [Spectral validation section] Spectral validation section: While qualitative matches to literature codes on spectral features are shown, the manuscript provides no quantitative error metrics (e.g., RMS or percentage differences) or sensitivity studies to uncertainties in the supplied 2T/11-species flowfield. This is load-bearing for the reliability claim, as the radiation module operates on externally supplied flow data whose accuracy directly determines the predicted fluxes.
Authors: We concur that quantitative error metrics are missing from the spectral validation section and that their addition would strengthen the assessment. We will compute and report RMS differences and/or average percentage deviations for the atomic and molecular spectral coefficient comparisons in the revised manuscript. For sensitivity to uncertainties in the supplied 2T/11-species flowfield, we will include a brief discussion of how variations in temperature and species concentrations could propagate to the radiative quantities. A full sensitivity study, however, lies beyond the scope of the present work, which centers on the radiation solver implementation and its basic validation against established spectral data. revision: partial
- Direct quantitative comparison of the integrated afterbody radiative heat fluxes to Fire II flight data or to benchmark codes such as NEQAIR on identical input flowfields cannot be added in the current revision, as this would require new simulations and access to matching external flowfield data and code implementations.
Circularity Check
No circularity: radiation module computes outputs from externally supplied flowfield inputs
full rationale
The paper describes a radiation solver that takes a two-temperature, 11-species air flowfield as input and applies line-by-line spectral modeling plus ray-tracing RTE solution to produce radiative heat flux predictions. Spectral coefficient comparisons to literature codes are presented as validation of the module, but the integrated Fire II afterbody results are computed values, not fitted or redefined from those inputs. No self-definitional equations, fitted parameters renamed as predictions, or load-bearing self-citations appear in the provided text. The derivation chain is self-contained as a forward computation tool.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Two-temperature model for thermochemical nonequilibrium in air
- domain assumption Line-by-line spectral modeling accurately captures atomic and molecular transitions in high-temperature air
Reference graph
Works this paper leans on
-
[1]
James A Fay and Frederick R Riddell. Theory of stagnation point heat transfer in dissociated air.Journal of the Aerospace Sciences, 25(2):73–85, 1958
work page 1958
-
[2]
Critical review of stagnation point heat transfer theory
H Hoshizaki, YS Chou, NG Kulgein, and JW Meyer. Critical review of stagnation point heat transfer theory. Technical report, NASA, 1975
work page 1975
-
[3]
Leonidas Gkimisis, Bruno Dias, James B Scoggins, Thierry Magin, Miguel A Mendez, and Alessandro Turchi. Data-driven modeling of hypersonic reentry flow with heat and mass transfer.AIAA Journal, 61(8):3269–3286, 2023
work page 2023
-
[4]
Odelma Teixeira and José Páscoa. Catalytic wall effects for hypersonic nozzle flow in thermochemical non- equilibrium.Acta Astronautica, 203:48–59, 2023. 33 APREPRINT- APRIL15, 2026
work page 2023
-
[5]
Zhenyuan Tang, Haonan Xu, Xueying Li, and Jing Ren. Aerodynamic heating in hypersonic shock wave and turbulent boundary layer interaction.Journal of Fluid Mechanics, 999:A66, 2024
work page 2024
-
[6]
Christopher O Johnston and Aaron M Brandis. Features of afterbody radiative heating for earth entry.Journal of Spacecraft and Rockets, 52(1):105–119, 2015
work page 2015
-
[7]
Refinements to afterbody radiative heating simulations for earth entry
Christopher O Johnston, Aaron M Brandis, and Marco Panesi. Refinements to afterbody radiative heating simulations for earth entry. In46th AIAA Thermophysics Conference, page 3693, 2016
work page 2016
-
[8]
Englewood Cliffs, NJ: Prentice-Hall, 1966
John J Martin.Atmospheric reentry: an introduction to its science and engineering. Englewood Cliffs, NJ: Prentice-Hall, 1966
work page 1966
-
[9]
Characterization of stagnation-point heat flux for earth entry
Aaron M Brandis and Christopher O Johnston. Characterization of stagnation-point heat flux for earth entry. In 45th AIAA Plasmadynamics and Lasers Conference, page 2374, 2014
work page 2014
-
[10]
Kenneth Sutton. Air radiation revisited. In19th Thermophysics Conference, page 1733, 1984
work page 1984
-
[11]
The aerothermodynamic environment of the apollo command module during superorbital entry
Dorothy B Lee and Winston D Goodrich. The aerothermodynamic environment of the apollo command module during superorbital entry. Technical report, NASA, 1972
work page 1972
-
[12]
Investigation of nonequilibrium radiation for mars entry
Aaron Brandis, Christopher Johnston, Marco Panesi, Brett Cruden, Dinesh Prabhu, and Deepak Bose. Investigation of nonequilibrium radiation for mars entry. In51st AIAA Aerospace Sciences Meeting including the New Horizons F orum and Aerospace Exposition, page 1055, 2013
work page 2013
-
[13]
Stardust: Discovery’s interstellar dust and cometary sample return mission
Kenneth L Atkins, Donald E Brownlee, Tom Duxbury, C-W Yen, Peter Tsou, and Joseph M Vellinga. Stardust: Discovery’s interstellar dust and cometary sample return mission. In1997 IEEE Aerospace Conference, volume 4, pages 229–245. IEEE, 1997
work page 1997
-
[14]
Nonequilibrium air radiation (nequair) program: User’s manual
Chul Park. Nonequilibrium air radiation (nequair) program: User’s manual. Technical report, NASA, 1985
work page 1985
-
[15]
Neqair96, nonequilibrium and equilibrium radiative transport and spectra program: User’s manual
Ellis E Whiting, Chul Park, Yen Liu, James O Arnold, and John A Paterson. Neqair96, nonequilibrium and equilibrium radiative transport and spectra program: User’s manual. Technical report, NASA, 1996
work page 1996
-
[16]
PhD thesis, The University of Queensland, 2011
Daniel Potter.Modelling of radiating shock layers for atmospheric entry at Earth and Mars. PhD thesis, The University of Queensland, 2011
work page 2011
-
[17]
Spark: a software package for aerodynamics, radiation and kinetics
Bruno Lopez and Mario Lino Da Silva. Spark: a software package for aerodynamics, radiation and kinetics. In 46th AIAA thermophysics conference, page 4025, 2016
work page 2016
-
[18]
PhD thesis, Virginia Polytechnic Institute and State University, 2006
Christopher Owen Johnston.Nonequilibrium shock-layer radiative heating for Earth and Titan entry. PhD thesis, Virginia Polytechnic Institute and State University, 2006
work page 2006
-
[19]
Christophe O Laux, TG Spence, CH Kruger, and RN Zare. Optical diagnostics of atmospheric pressure air plasmas.Plasma Sources Science and Technology, 12(2):125–138, 2003
work page 2003
-
[20]
Kazuhisa Fujita and Takashi Abe.SPRADIAN, structured package for radiation analysis: Theory and application. Number 669. Institute of Space and Astronautical Science, 1997
work page 1997
-
[21]
Erwan Pannier and Christophe O Laux. Radis: A nonequilibrium line-by-line radiative code for co2 and hitran-like database species.Journal of Quantitative Spectroscopy and Radiative Transfer, 222:12–25, 2019
work page 2019
-
[22]
PhD thesis, School of Engineering, The University of Queensland, 2008
Rowan Gollan.Computational Modelling of High-temperature Gas Effects with Application to Hypersonic Flows. PhD thesis, School of Engineering, The University of Queensland, 2008
work page 2008
-
[23]
Rowan J Gollan, Peter A Jacobs, Nicholas N Gibbons, and Kyle A Damm. The eilmer user guide. Technical report, The University of Queensland, 2026
work page 2026
-
[24]
University of Queensland. Gas dynamics toolkit (gdtk). https://gdtk.uqcloud.net/. Accessed: 2026-04-05
work page 2026
-
[25]
Yuzhe Zhang, Qizhen Hong, Xiaoyong Wang, Chao Yang, and Quanhua Sun. Uncertainty quantification in electronic state-specific modeling of high-enthalpy nitrogen flows.AIAA Journal, 63(10):4479–4486, 2025
work page 2025
-
[26]
Gerhard Herzberg.Molecular Spectra and molecular structure-V ol I, volume 1. Read Books Ltd, 2013
work page 2013
-
[27]
Chul Park.Nonequilibrium hypersonic aerothermodynamics. Wiley, New York, 1989
work page 1989
-
[28]
Walter G Vincenti, Charles H Kruger Jr, and T Teichmann.Introduction to physical gas dynamics. Wiley, New York, 1966
work page 1966
- [29]
-
[30]
Sophie Chauveau, Christine Deron, Marie-Yvonne Perrin, Philippe Rivière, and Anouar Soufiani. Radiative transfer in lte air plasmas for temperatures up to 15,000 k.Journal of Quantitative Spectroscopy and Radiative Transfer, 77(2):113–130, 2003. 34 APREPRINT- APRIL15, 2026
work page 2003
-
[31]
A. Kramida, Yu. Ralchenko, J. Reader, and and NIST ASD Team. NIST Atomic Spectra Database (ver. 5.12), [Online]. Available: https://physics.nist.gov/asd [2025, November 6]. National Institute of Standards and Technology, Gaithersburg, MD., 2024
work page 2025
-
[32]
Springer Science & Business Media, 2007
Gordon WF Drake.Springer handbook of atomic, molecular , and optical physics. Springer Science & Business Media, 2007
work page 2007
-
[33]
Line-by-line transport calculations for jupiter entry probes
JO Arnold, DM Cooper, Chul Park, and SG Prakash. Line-by-line transport calculations for jupiter entry probes. In American Institute of Aeronautics and Astronautics, Thermophysics Conference, number AIAA PAPER 79-1082, page 58, 1979
work page 1979
-
[34]
The Opacity Project Team.https://cds.unistra.fr//topbase/topbase.html. On-line Resources, 2009
work page 2009
-
[35]
G Peach. Continuous absorption coefficients for non-hydrogenic atoms.Memoirs of the Royal Astronomical Society, V ol. 73, p. 1, 73:1–123, 1970
work page 1970
-
[36]
Forrest R Gilmore, Russ R Laher, and Patrick J Espy. Franck–condon factors, r-centroids, electronic transition moments, and einstein coefficients for many nitrogen and oxygen band systems.Journal of physical and chemical reference data, 21(5):1005–1107, 1992
work page 1992
-
[37]
Rihui Liang, Yanming Liu, and Fangyan Li. Radiative transition probabilities of diatomic N2, O2, N + 2 and NO in air plasmas.Physica Scripta, 96(12):125402, 2021
work page 2021
-
[38]
The NIST Team.https://webbook.nist.gov/chemistry/form-ser/. On-line Resources, 2025
work page 2025
-
[39]
JO Arnold, EE Whiting, and GC Lyle. Line by line calculation of spectra from diatomic molecules and atoms assuming a voigt line profile.Journal of Quantitative Spectroscopy and Radiative Transfer, 9(6):775–798, 1969
work page 1969
-
[40]
István Kovács.Rotational structure in the spectra of diatomic molecules. Academiai Kiado, 1969
work page 1969
-
[41]
Christopher O Johnston and Alireza Mazaheri. Impact of non-tangent-slab radiative transport on flowfield– radiation coupling.Journal of Spacecraft and Rockets, 55(4):899–913, 2018
work page 2018
-
[42]
Álvaro González. Measurement of areas on a sphere using fibonacci and latitude–longitude lattices.Mathematical Geosciences, 42(1):49–64, 2010
work page 2010
-
[43]
Sergey F Gimelshein and Ingrid J Wysong. Validation of high-temperature air reaction and relaxation models using emission data.Journal of Thermophysics and Heat Transfer, 33(3):606–616, 2019
work page 2019
-
[44]
Modeling of no emissions from a hypersonic ground- based facility
Ozgur Tumuklu, Matthew Jouffray, and Deborah A Levin. Modeling of no emissions from a hypersonic ground- based facility. InAIAA Scitech 2021 F orum, page 0105, 2021
work page 2021
-
[45]
Radiation measurements of shockwaves in synthetic air and pure nitrogen
Alex B Glenn, Orin Varley, Peter Collen, and Matthew McGilvray. Radiation measurements of shockwaves in synthetic air and pure nitrogen. InAIAA SCITECH 2024 F orum, page 0226, 2024
work page 2024
-
[46]
Measurement of radiative non-equilibrium for air shocks between 7-9 km/s
Brett A Cruden and Aaron M Brandis. Measurement of radiative non-equilibrium for air shocks between 7-9 km/s. In47th AIAA Thermophysics Conference, page 4535, 2017
work page 2017
-
[47]
Analysis of shockwave radiation data in nitrogen
Brett A Cruden and Aaron M Brandis. Analysis of shockwave radiation data in nitrogen. InAIAA Aviation 2019 F orum, page 3359, 2019
work page 2019
-
[48]
S Chauveau, M-Y Perrin, Ph Riviere, and A Soufiani. Contributions of diatomic molecular electronic systems to heated air radiation.Journal of Quantitative Spectroscopy and Radiative Transfer, 72(4):503–530, 2002
work page 2002
-
[49]
Project FIRE integrated post flight evaluation report, flight no
MA Minovitch. Project FIRE integrated post flight evaluation report, flight no. II.NASA Contractor Report, page 05936, 1965
work page 1965
-
[50]
Flight parameters and vehicle performance for project Fire flight II, launched may 22, 1965
John H Lewis Jr and William I Scallion. Flight parameters and vehicle performance for project Fire flight II, launched may 22, 1965. Technical report, NASA, 1966
work page 1965
-
[51]
ES Cornette and Forebody Temperatures. Calorimeter heating rates measured during project FIRE II reentry at 11.35 kilometers per second,".NASA TM X-1305, National Air and Space Administration, 1966
work page 1966
-
[52]
T Slocumb. Project Fire flight 2 afterbody temperatures and pressures at 11.35 kilometers per second (37,200 feet per second). Technical report, NASA, 1966
work page 1966
-
[53]
Dona L Cauchon. Radiative heating results from the FIRE II flight experiment at a reentry velocity of 11.4 kilometers per second. Technical report, NASA, 1967
work page 1967
-
[54]
David Hash, Joseph Olejniczak, Michael Wright, Dinesh Prabhu, Maria Pulsonetti, Brian Hollis, Peter Gnoffo, Michael Barnhardt, Ioannis Nompelis, and Graham Candler. FIRE ii calculations for hypersonic nonequilibrium aerothermodynamics code verification: Dplr, laura, and us3d. In45th AIAA Aerospace Sciences Meeting and Exhibit, page 605, 2007. 35 APREPRINT...
work page 2007
-
[55]
Xiaoyong Wang, Qizhen Hong, Yuan Hu, and Quanhua Sun. On the accuracy of two-temperature models for hypersonic nonequilibrium flow.Acta Mechanica Sinica, 39(2):122193, 2023
work page 2023
-
[56]
Chul Park, Richard L Jaffe, and Harry Partridge. Chemical-kinetic parameters of hyperbolic earth entry.Journal of Thermophysics and Heat transfer, 15(1):76–90, 2001
work page 2001
-
[57]
Elden S Cornette.F orebody temperatures and calorimeter heating rates measured during project FIRE 2 reentry at 11.35 kilometers per second. AIAA, 1966
work page 1966
-
[58]
Modification of chemical-kinetic parameters for 11-air species in re-entry flows
Jae Gang Kim and Sung Min Jo. Modification of chemical-kinetic parameters for 11-air species in re-entry flows. International Journal of Heat and Mass Transfer, 169:120950, 2021
work page 2021
-
[59]
Bianca R Capra and Richard G Morgan. Total heat transfer measurements on a flight investigation of reentry environment model.Journal of Spacecraft and Rockets, 50(3):494–503, 2013
work page 2013
-
[60]
Features of afterbody radiative heating for titan entry
Christopher O Johnston, Thomas K West, and Aaron M Brandis. Features of afterbody radiative heating for titan entry. InAIAA Aviation 2019 F orum, page 3010, 2019
work page 2019
-
[61]
Features of afterbody radiative heating for jupiter entry
Aaron J Erb, Thomas K West, and Christopher O Johnston. Features of afterbody radiative heating for jupiter entry. InAIAA SciTech F orum and Exposition, 2021
work page 2021
-
[62]
Convective and radiative heat transfer analysis for the fire ii forebody
Robert B Greendyke and Lin C Hartung. Convective and radiative heat transfer analysis for the fire ii forebody. Journal of Spacecraft and Rockets, 31(6):986–992, 1994. 36
work page 1994
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.