Assessment of Radiative Heating for Hypersonic Earth Reentry Using Nongray Step Models

: Accurate prediction of the aerothermal environment is of great signiﬁcance to space exploration and return missions. The canonical Fire II trajectory points are simulated to investigate the radiative transfer in the shock layer for Earth reentry at hypervelocity above 10 km/s using a developed radiation–ﬂowﬁeld uncoupling method. The thermochemical nonequilibrium ﬂow is solved by an in-house PHAROS Navier–Stokes code, while the nongray radiation is integrated by the tangent slab approximation, respectively, combined with the two-, ﬁve-, and eight-step models. For the convective heating, the present results agree well with the data of Anderson’s relation. For the radiative heating, the two-step model predicts the closest values with the results of Tauber and Sutton’s relationship, while the ﬁve-and eight-step models predict far greater. The three-step models all present the same order of magnitude of radiative heating of 1 MW/m 2 and show a consistent tendency with the engineering estimation. The Planck-mean absorption coefﬁcient is calculated to show the radiative transfer signiﬁcantly occurs in the shock layer. By performing the steady simulation at each ﬂight trajectory point, the present algorithm using a nongray step model with moderate efﬁciency and reasonable accuracy is promising to solve the real-time problem in engineering for predicting both convective and radiative heating to the atmospheric reentry vehicle in the future.


Introduction
Today, many countries are energetically developing planetary exploration and subsequent return projects [1].The mission spacecraft, usually with a blunt nose, must undergo the harsh thermal environment at hypervelocity even greater than 10 km/s for Earth atmosphere reentry [2].In this situation, the strong bow shock wave around the reentry forebody can stimulate an extremely high temperature with an order of magnitude of 10,000 K and very complicated thermochemical nonequilibrium phenomena including the high excitations of vibrational and electronic energy modes, dissociation, and ionization reactions of molecules and atoms of air species [3][4][5].Except for the convective and chemical diffusive heating, radiative heat transfer could be a considerable or even dominated contributor to the total aerodynamic heating [6].
In fact, the occurrence of air radiation is not only one simple additional mode of heat transfer but also has important impacts on the flow characteristics of the shock layer [7].On one hand, the radiative energy emission and absorption make the flowfield nonadiabatic and lead to the radiation cooling effect [8], which results in a lower temperature, higher density, and subsequently thinner shock layer.On the other hand, the thermochemical reactions and radiation of air species may overlap in some concerned flight regimes.The radiation limits the chemical reaction processes of air species, whereas the air thermodynamic states and composition, in turn, directly affect the level of radiance.The radiation-flowfield coupling should be evaluated [9].Therefore, in engineering, the thermal protection system design of the entry vehicle requires the accurate aerothermal analysis of the hypersonic thermochemical nonequilibrium flowfield including radiation to ensure safe entry [10][11][12].
However, until now, the solution of radiation in the hypersonic shock layer has not been an easy task [13], even for the uncoupled radiation calculation only [14].First, although the scattering could be neglected, the absorption and emission must be accounted for in the air radiation in the hypersonic shock layer.The radiation properties of air, which are emission and absorption coefficients, always vary dramatically with the frequency ranging from zero to infinity [15].The accurate calculation of the air absorption coefficients requires a very large number of frequency points with an order of one million, which is undoubtedly time-consuming and even unaffordable [16].Second, the radiative transfer depends on the multi-dimensions of space, angle, and frequency.The algorithm for solving the radiative transfer equation (RTE) has to simultaneously account for the spatial, angular, and frequency discretizations, which makes the radiation computation difficult and inefficient [17].Facing the foregoing barriers, it is necessary to employ the reduced models to calculate the air radiation property and energy transportation for the aerodynamic heating to reentry vehicle [18].
The step model is an efficient method to calculate the air radiation property with moderate accuracy [19].It divides the whole frequency space into several spectral intervals even including atomic lines in which the radiative absorption coefficient is equivalent to a constant [20].Although the step model appears to be coarser than the line-by-line or multi-band models, it usually can predict the reasonable radiative flux of the hypersonic vehicle [21].For solving radiative transfer, there have been various numerical methods, such as the tangent slab (TS) approximation [22], the spherical harmonics method [23], the ray-tracing method [24], the finite volume method [25], the discrete ordinates method [26], and Monte Carlo method [27].Each method has both advantages and disadvantages either in prediction accuracy or in computation efficiency [28].Among these methods, TS has seemed to be the most frequently used method for both coupled and uncoupled flowfield-radiation simulations for decades [2].TS originates from the one-dimensional analytical integrated solution of RTE for the participating medium between two infinite parallel plates with the radiation variations only along the normal direction of the plate [22].Since the shock layer flow and thermodynamic properties mainly vary in the direction normal to the surface of the blunt-nosed reentry vehicle, TS is a good approximation to model the radiation energytransportation [29].Practically, researchers usually use the body-normal grid lines for convenient TS calculation [2].Hartung et al. used TS to evaluate the characteristics of the shock wave precursor ahead of the hypervelocity entry vehicle including the radiative effects [30].Wright et al. found that TS could over-predict the value of the stagnation radiative heat flux by a minimum of 20% for the Titan aerocapture case [31].Johnston et al. used TS combined with a viscous shock layer flowfield model to account for the radiation-flowfield coupling and showed that the coupled simulation reduced the radiative heating by about 30%, while the convective heating decreased slightly [32].Bauman et al. coupled a reacting flow model and a surface ablation model with TS to develop a two-way loose-coupling procedure for simulating the hypersonic flows with radiation and ablation [33].Johnston and Brandis employed TS to calculate the radiative heating and identified the radiation as a major contributor to afterbody heating for Earth entry at velocities above 10 km/s.They also showed that TS overestimates the afterbody radiative heat transfer by as much as 50% [34].Generally, although TS has some deficiencies and researchers have also proposed a non-tangent-slab procedure, TS is still very worthy of being considered as a first choice to solve RTE coupled or uncoupled with flowfield solver to predict the radiative heating for Earth reentry at hypervelocity, especially to obtain a conservative estimation of the aerodynamic heating in the engineering initial design with moderate accuracy and computation efficiency.
In fact, there have been many engineering empirical relations to quickly predict both the wall convective and radiative heat fluxes of the Earth entry capsule [35].Most methods have been too coarse for the modern design of Earth reentry capsule which can only give the value of heat transfer at the stagnation point but has no ability to present the distribution characteristics of aerodynamic heating and some key physical variables of interest in the shock layer, such as electrons number density and radiation absorption coefficient.However, the accurate radiation-flowfield coupling simulation is time-consuming and even unaffordable for the initial thermal protection system design of reentry vehicles [16].Hence, it is more rational and practical in engineering to choose an uncoupled radiationflowfield procedure, in which the radiation is calculated only once based on the flow simulation results [36].
Gupta et al. used a viscous shock-layer code and aerotherm radiation code with nonequilibrium and equilibrium chemistry to estimate the convective and radiative heating of the Fire II vehicle, but the calculations were limited to the stagnation region [37].Olynick et al. developed a nonequilibrium, axisymmetric flow solver coupled with a radiation GIANTS/NOVAR code to obtain values of the stagnation radiation intensity in the 0.2 and 6.2 eV and the total aerodynamic heating over the entire Fire II vehicle, which provided the better predictions than previously numerical simulations [38].Palmer et al. incorporated the NEQAIR line-by-line radiation code into the DPLR flow solver to investigate the effects of fluid dynamics/radiation coupling by comparing coupled and uncoupled results, and they found that the greatest coupling effect of Fire II occurred at the 1643 s trajectory point [39].Soucasse et al. implemented a hybrid statistical narrow band (HSNB) model with a two-temperature nonequilibrium model to calculate the 1D stagnation line radiative transfer.The HSNB model could reproduce the line-by-line results with an accuracy of better than 5% and a computational time speed up about two orders of magnitude [40].Bonin and Mundt developed a full three-dimensional photon Monte Carlo radiative transport solver to study arbitrary thermal radiation within equilibrium and nonequilibrium hypersonic flows, the code of which was line-by-line accurate but time-consuming [41].Although great progress has been made to predict radiative transfer in hypersonic nonequilibrium flow for atmospheric entry over decades, the too sophisticated methods developed with high accuracy are always of much heavy computational burden and even unaffordable, such as the line-by-line and narrow-band radiation models, especially for multi-dimensional coupling simulation.Therefore, there is still a practical demand for developing the reduced models with moderate time efficiency and reasonable accuracy to predict both radiative heating and radiation characteristics in flowfield for the thermal protection design of Earth reentry vehicles in engineering applications.
For solving the above-mentioned problems, the first objective of this paper is to develop an uncoupled radiation-flowfield algorithm for predicting the aerothermal environment of Earth reentry vehicles at hypervelocities above 10 km/s, which consists of the tangent slab approximation, the nongray step model, and a Navier-Stokes solver including the thermochemical nonequilibrium effects.The detailed physical models and numerical schemes are presented in Section 2. Another objective is to evaluate the performance of the radiation-flowfield uncoupling procedure by analyzing the canonical reentry trajectory cases of Fire II capsule, especially focusing on the aerodynamic heating and the radiation characteristics in the shock layer.The convective and radiative heat fluxes at the stagnation point are also estimated throughout the trajectory, respectively, using the Anderson and Tauber and Sutton relations.A comparison between the present results and those of engineering methods and previous studies is thoroughly presented and discussed in Section 3. The final conclusions are drawn in Section 4.

Flow Governing Equations with Thermochemical Nonequilibrium Models
The hypersonic Earth reentry flow is governed by the Navier-Stokes equations with the two-temperature model including the thermochemical nonequilibrium effects in the conservative forms.It assumes that the air is a multi-gas mixture and for all the composite species, the translational and rotational energy modes are in equilibrium in one translational-rotational temperature, T tr , and the vibrational, electronic, and electron energies are uniformly described by one vibrational-electronic temperature, T ve [42].In this manner, the mass, momentum, and energy conservation equations of the hypersonic nonequilibrium flow can be expressed as follows [3]: ∂ρe ∂t ∂ρe ve ∂t where t is the time, x i is the ordinate variable in the i direction, N s is the total number of the air species, ρ s and ρ are the species density and total density, u i is the flow velocity component in the i direction, p is the pressure, τ ij is the viscous stress tensor, e and e ve are the total energy and vibrational-electronic energy, h is the total enthalpy, q j and q ve,j are the total heat flux and vibrational-electronic heat flux, h s and h ve,s are the enthalpy and vibrational-electronic enthalpy of the species s, J s,j is the mass diffusion flux of the species s in the j direction, ω s is the mass production rate of the species s per unit volume, ω ve is the vibrational-electronic energy source term, and ω r is the radiative source term.The state equation of the air follows: where the subscript "e" represents the electron.
In the present study, the thermodynamic properties of all the species are calculated using analytical relations of the translational, rotational, and electronic excitation energy modes based on the Born-Oppenheimer approximation [43].The transport properties of the air mixture including the dynamic viscosity, thermal conductivity, and species diffusion coefficients are calculated via the extension of Yos' formula with the collision integrals [44].
The mass rate of production of species s is expressed as follows [45]: where N r is the total number of chemical reactions, M s is the molecular mass per mole of species s, ν f r,s and ν b r,s are the forward and backward reaction stoichiometric coefficients of species s of the r-th reaction, k f,r and k b,r are the forward and backward reaction rate coefficients of the r-th reaction, respectively.For each chemical reaction, the forward reaction rate coefficient is calculated using the Arrhenius formula [46], while the backward reaction rate coefficient is obtained from the corresponding forward reaction rate coefficient divided by the equilibrium constant, which is computed using the temperature fitting expression [44].
The vibrational-electronic energy source term ω ve is modeled in the expression proposed by Gnoffo et al. [46], in which the translational-vibrational energy exchange part is calculated by the Landau-Teller model [47], and the relaxation time is calculated by the Millikan-White expression [48] with Park's high-temperature correction [49].For the present uncoupled radiation-flowfield simulation, the radiative source term ω r is neglected.

Flowfield Solver
The Navier-Stokes equations with thermochemical models are solved using an inhouse CFD code PHAROS (Parallel Hypersonic Aerothermodynamics and Radiation Optimized Solver) [50].PHAROS is a parallel multi-block finite volume solver, in which the inviscid fluxes are computed by the modified Steger-Warming flux vector splitting scheme [51] using MUSCL extrapolation [52] with minmod limiters for high order accuracy and stability.The viscous fluxes are discretized in the second central difference.A line relaxation approach is employed for the time marching [53].PHAROS has been used to solve many types of hypersonic thermochemical nonequilibrium flowfields [3,5,6,54].More information on PHAROS can be found in Ref. [50].

Step Models for Radiation Properties
Based on recognizing the variation characteristics of the self-absorption property of air radiation with wavelength, the step model selects some demarcation wavelength points and divides the infinite spectrum into several consecutive spectral regions, such as vacuum ultraviolet and visible bands.In each spectral band, the volatile radiation absorption coefficients are averaged into a constant Planck-mean value, which makes the total radiation a sum of the contributions of all the step regions.The absorption coefficient in the spectral region of each step depends on the temperature, density, and air composition parameters [20].It has been shown that the step model can well account for the important effects of shock-layer nongray self-absorption and radiative cooling on radiative heat transfer [19].Additionally, compared to line-by-line calculation, the step model is a reasonable simplification and can reduce the computational time obviously.Therefore, the step model is selected for the present study.In this paper, the radiation absorption coefficients of the high-temperature gas mixture of air are calculated via the two-step, five-step, and eight-step nongray models [55], respectively, for the purpose of comparative study.The two-step and five-step models are developed based on the hightemperature atomic nitrogen mainly accounting for the air radiance at temperatures above 8000-10,000 K, while the eight-step model includes both the atomic and molecular emission absorption of the high-temperature air [56].The wavelength regions of the three foregoing step models are listed in Table 1.The specific formulas and parameters of the three-step models are given in Appendix A. Neglecting the gas scattering effect, the high-temperature air radiation process in the hypersonic reentry shock layer is described via the radiative transfer equation (RTE) as follows [57]: where x represents the spatial position vector, ν is the radiation frequency, I ν and I bν is the spectral radiative intensity and blackbody radiative intensity at frequency ν, respectively, B j is the j-th component of unit vector B in the transmission direction of I ν , and κ ν is the spectral absorption coefficient at x. TS approximates the hypersonic shock layer around the reentry blunt body as an infinite slab with physical variables only changing in the perpendicular direction to the body surface.Hence, there is an integrated solution for RTE as follows [22]: and where τ ν is the optical thickness at frequency ν perpendicular to the body surface (τ ν = 0 at body surface), the subscript "δ" and "w" represent the outer edge of the shock layer and the wall, respectively; t is a dummy variable of integration, and ε w and T w are the wall emissivity and temperature, respectively.E n is the integro-exponential function of order n as follows: According to step models in Section 2.3, the integration (9) can be written as: where m is the index of the step region.In this paper, the integral with respect to optical thickness in Equation ( 11) is calculated using the trapezoidal method to perform TS procedures.

Radiation-Flowfield Uncoupling Algorithm
For radiation-flowfield coupling simulation, TS approximation with a nongray step model should do the upward and downward integrals of radiative heat flux divergence along the ray line normal to the wall for computing the radiative source term ω r in Equation ( 4) [22].Assuming a two-dimensional problem is solved, the mesh has a total number of grid nodes of N ξ × N η , where N ξ is the discretized number parallel to the wall and N η is the discretized number normal to the wall.Thus, the coupling simulation should need about N ξ × N η × M × N η calculations in one radiation iteration, for which the last product factor N η is due to the numerical integration of radiative heat flux divergence and M is the number of radiation step regions.
A radiation-flowfield uncoupling algorithm is proposed in the present study, in which the hypersonic reentry flow with nonequilibrium chemistry is first solved by PHAROS in Section 2.2, and then the radiation transfer is calculated only once based on the convergent flowfield by the nongray step model in Section 2.3 and TS approximation in Section 2.4.It only needs N ξ × M × N η calculations for radiation solution in total.Therefore, in one iteration, the number of calculations for coupled radiation simulation is far greater with an order of magnitude of O(N η ) than that of the present uncoupled algorithm.Particularly, the present scheme only needs to compute radiation one time in total, while the coupled simulation has to update the radiation in each step over the whole computational process.Even for the loosely coupled manner, it still requires a considerable computational cost, in which the radiation is updated one time after a certain number of flow iterations.Therefore, the present method is more time-efficient.Additionally, the present algorithm can provide more radiation information including the radiative heating to the whole surface of the reentry vehicle and the absorption properties distributed in the flowfield, but the engineering methods certainly cannot make it, that will be seen later in Section 3.4.

Results and Discussion
Fire II was a scaled-down Apollo-shaped capsule launched in 1965 with the calorimeter instrumentation to obtain reentry heating at hyperbolic velocities.In the reentry phase, Fire II jettisoned two nonablating heat shields in sequence at selected trajectory points, thereby with various vehicle nose radiuses.The Fire II flight experiment has become a benchmark for investigating the aerothermal environment for hypersonic Earth entry.Six trajectory points of Fire II have been simulated in the present study with the flight conditions listed in Table 2, where H is the flight altitude, V ∞ and Ma are the flight velocity and Mach number, R N is the vehicle nose radius, ρ ∞ and T ∞ are the freestream density and temperature, and T w is the wall temperature.For each case in Table 2, the time of flow over the vehicle, R N /V ∞ , is of an order of magnitude around 10 −4 s, but the time scale of entry down in the altitude is around 0.5 s.The latter is far greater than the former, which means a steady flowfield establishes very quickly.Therefore, we can use the present algorithm to perform steady simulation for each case to fulfill the real-time prediction throughout the Fire II trajectory.Only the Fire II forebody is considered in this paper, the axisymmetric geometry and grids of which are shown in Figure 1.The computational mesh has the dimensions of 153 × 128 (axial × radial) for every case and the spacing of the first grid layer perpendicular to the wall can ensure the cell Reynolds number with an order of magnitude of one in order to predict the reliable aerodynamic heating [58].The noncatalytic wall condition is used and the wall emissivity is uniformly set to be one.

Convective Heating
where RN is the vehicle nose radius; ρ∞ and V∞ are the freestream density and velocity; Hw and He are the enthalpies at the wall and the outer edge of the boundary layer, which can be evaluated by using the wall temperature Tw and freestream total temperature T0, respectively.T0 can be calculated by the freestream temperature T∞ and Mach number Ma. RN, ρ∞, V∞, Tw, T∞, and Ma are all listed in Table 2.The present results agree well with the other three sets of data both in the magnitude of value at each trajectory point and the total variation tendency with time, which verifies the reliability of the present thermochemical nonequilibrium flow solver PHAROS.Throughout the trajectory from t = 1634 s to 1645 s, the stagnation convective heat transfer is always greater than 1 MW/m 2 and increases continuously up to almost 8 MW/m 2 .Figure 3 shows the convective heating over the whole surface of the Fire II forebody for each trajectory point, which presents that the convective heat transfer still maintains a high level of magnitude greater than 1 MW/m 2 in the region outside the stagnation point for all cases.Figure 4 further compares the forebody convective heating line at t = 1636 s predicted by the present method with those obtained by DPLR and LAURA [59], the good agreement of which shows the high prediction accuracy of the present PHAROS solver again.

Convective Heating
where R N is the vehicle nose radius; ρ ∞ and V ∞ are the freestream density and velocity; H w and H e are the enthalpies at the wall and the outer edge of the boundary layer, which can be evaluated by using the wall temperature T w and freestream total temperature T 0 , respectively.T 0 can be calculated by the freestream temperature T ∞ and Mach number Ma. R N , ρ ∞ , V ∞ , T w , T ∞ , and Ma are all listed in Table 2.The present results agree well with the other three sets of data both in the magnitude of value at each trajectory point and the total variation tendency with time, which verifies the reliability of the present thermochemical nonequilibrium flow solver PHAROS.Throughout the trajectory from t = 1634 s to 1645 s, the stagnation convective heat transfer is always greater than 1 MW/m 2 and increases continuously up to almost 8 MW/m 2 .Figure 3 shows the convective heating over the whole surface of the Fire II forebody for each trajectory point, which presents that the convective heat transfer still maintains a high level of magnitude greater than 1 MW/m 2 in the region outside the stagnation point for all cases.Figure 4 further compares the forebody convective heating line at t = 1636 s predicted by the present method with those obtained by DPLR and LAURA [59], the good agreement of which shows the high prediction accuracy of the present PHAROS solver again.

Thermochemical Nonequilibrium Flowfield
Figure 5 compares the translational-rotational temperature Ttr and vibrational-electronic temperature Tve in the Fire II flowfield throughout the trajectory from t = 1634 s to 1645 s.At t = 1634 s, there is a remarkable difference between Ttr and Tve in a distance closely behind the bow shock, which suggests the translational-rotational and vibrational-electronic energy modes are highly in nonequilibrium in this region.As time goes on with the altitude and velocity down, the thermodynamic nonequilibrium tends to be weakened steadily and the total level of magnitude of temperature also decreases gradually.From t = 1637 s to 1645 s, Ttr and Tve have become consistent in most areas of the shock layer.For all cases, both Ttr and Tve exceed 10 4 K in the shock layer around the Fire II forebody, and the peak of the translational-rotational temperature even reaches up to 44,000 K at t = 1634 s.Such extreme high temperature directly leads to the harsh aerothermal environment for Earth reentry and results in significant radiative heating [19].Figure 6 presents the number densities of the species O, N, O2, N2, NO, and e -along the stagnation line for each case.The high atomic and electron concentrations demonstrate the strong dissociation and ionization reactions of air in the reentry shock layer of Fire II.At t = 1634 s, the number density of each species changes remarkably along the stagnation line, while at the following trajectory points, such variations become more and more unnoticeable and the concentration of each species approaches a constant in most regions of the shock layer.Due to the high ionization, the number density of electrons (that is the sum of the number densities of all positive ions, such as O + , N + , and NO + ) is greater than 10 21 m −3 in most areas of the shock layer.The free electrons and ions constitute a plasma sheath around the reentry vehicle, which absorbs the radio-frequency radiation and causes the communication blackout [60].

Thermochemical Nonequilibrium Flowfield
Figure 5 compares the translational-rotational temperature T tr and vibrational-electronic temperature T ve in the Fire II flowfield throughout the trajectory from t = 1634 s to 1645 s.At t = 1634 s, there is a remarkable difference between T tr and T ve in a distance closely behind the bow shock, which suggests the translational-rotational and vibrational-electronic energy modes are highly in nonequilibrium in this region.As time goes on with the altitude and velocity down, the thermodynamic nonequilibrium tends to be weakened steadily and the total level of magnitude of temperature also decreases gradually.From t = 1637 s to 1645 s, T tr and T ve have become consistent in most areas of the shock layer.For all cases, both T tr and T ve exceed 10 4 K in the shock layer around the Fire II forebody, and the peak of the translational-rotational temperature even reaches up to 44,000 K at t = 1634 s.Such extreme high temperature directly leads to the harsh aerothermal environment for Earth reentry and results in significant radiative heating [19].Figure 6 presents the number densities of the species O, N, O 2 , N 2 , NO, and e − along the stagnation line for each case.The high atomic and electron concentrations demonstrate the strong dissociation and ionization reactions of air in the reentry shock layer of Fire II.At t = 1634 s, the number density of each species changes remarkably along the stagnation line, while at the following trajectory points, such variations become more and more unnoticeable and the concentration of each species approaches a constant in most regions of the shock layer.Due to the high ionization, the number density of electrons (that is the sum of the number densities of all positive ions, such as O + , N + , and NO + ) is greater than 10 21 m −3 in most areas of the shock layer.The free electrons and ions constitute a plasma sheath around the reentry vehicle, which absorbs the radio-frequency radiation and causes the communication blackout [60].

Radiative Heating
The radiative transfer is solved by employing the radiation-flow uncoupling algorithm, in which the radiation is calculated only once based on the convergent flowfield solution using the methods in Section 2. Compared with the radiation-flow coupling approach, although the uncoupling simulation compromises some accuracy in physics, the procedure really can greatly reduce the computing power and time.Figure 7 compares the results of radiative heating at the stagnation point calculated by Tauber and Sutton's (T & S) engineering relationship [61], Gupta et al. [37], Olynick et al. [38], and the present TS method with the two-step, five-step, and eight-step models, respectively.The T & S relation is expressed as follows: where R N is the vehicle nose radius; ρ ∞ and V ∞ are the freestream density and velocity; a and b are empirical exponents; f is a tabulated function of freestream velocity V ∞ .The detailed descriptions of a, b, and f can be found in Ref. [61].
T & S relation is expressed as follows: ( ) where RN is the vehicle nose radius; ρ∞ and V∞ are the freestream density and velocity; a and b are empirical exponents; f is a tabulated function of freestream velocity V∞.The detailed descriptions of a, b, and f can be found in Ref. [61].Gupta et al. developed a RAD code accounting for the molecular band, continuum, and atomic line transitions.They employed a detailed frequency dependence of the absorption coefficients to integrate over the radiation spectrum and used the TS approximation for integrating over physical space, but only limited to the stagnation region.Olynick et al. developed a GIANTS/NOVAR code to obtain values of the stagnation radiation intensity in the 0.2 and 6.2 eV range.They used a "smeared band" approximation instead of a line-by-line approach to account for the air radiation properties of the molecular band systems and reduce the total number of spectral points.Even doing so, Olynick et al. still needed to perform TS integration at each radiation grid point consisting of around 1000 frequency points for calculating absorption and emission coefficients, which is very time-consuming.Gupta et al. and Olynick et al. both implemented the radiation-flow coupling simulation with nonequilibrium chemistry.Figure 7 shows that the radiative heat fluxes calculated by the TS approximation with the step models are mostly greater than those obtained by the T & S relation, Gupta et al., and Olynick et al.In fact, compared with the coupling simulation, the present radiation-flow uncoupling method should predict the higher level of radiative heating [22], which is reasonable and can be seen as the conservative upper limit in the engineering application.should predict the higher level of radiative heating [22], which is reasonable and can be seen as the conservative upper limit in the engineering application.
Although the noticeable differences can be seen among the radiative heat fluxes of different methods, the orders of magnitude of all the data are basically the same and they show a similar tendency, first an increase and then decrease.The five-step model gives the greatest stagnation radiative heat transfer, while the results of the two-step model are the smallest.Unexpectedly, although the two-step model is coarsest in spectral space, it predicts the closest values with those of the T & S relation, Gupta et al. [37], and Olynick et al. [38] compared to the other two models.However, we cannot simply deduce that the two-step model is the best option because it is just only applicable for Fire II trajectory points.The final evaluation needs more experimental data and further high-resolution numerical simulation in the future.

Radiation Field
Except for the more affordable time cost than that of the coupling simulation, the present radiation-flow uncoupling algorithm can provide the more detailed information than those of the engineering methods, not only the radiative heat transfer at the stagnation point but also the radiative heating distribution on the vehicle surface and the radiation characteristics in the whole flowfield.The latter radiation distributions in the flowfield are rarely shown in the previous literature [37][38][39][40][41] but may play some important roles in understanding the mechanism of radiative heating and designing the thermal protection system of atmospheric entry vehicles in the future [6].Figure 8 shows the radiative heating over the whole surface of the Fire II forebody predicted by the TS approximation with the two-step, five-step, and eightstep models, respectively.If the heat transfer of magnitude is not taken into consideration, the tendencies obtained by the three models are basically consistent.The radiative heating varies flatly over the whole vehicle surface for the cases from t = 1634 s to 1637 s, while it becomes more bending for the following cases from t = 1640 s to 1645 s. Figure 9 shows the distributions of the Planck-mean absorption coefficient around the Fire II forebody throughout the trajectory calculated by the two-step, five-step, and eight-step models, respectively.The Planck-mean absorption coefficient is one of the most important average forms that describes the total emission from a fluid element and indicates the level of radiative heat loss from the nearly optically thin flowfield.The Planck-mean absorption coefficient for the present step model is defined as follows: where m is the spectral step index, κ m is the absorption coefficient for the m-th step, I bν is the spectral blackbody intensity at radiation frequency ν, σ is the Stefan-Boltzmann constant, and T is the temperature.In the present hypersonic nonequilibrium flow, κ P is calculated using the vibrational-electronic temperature [62].Figure 9 shows that κ P is remarkable only in the shock layer, particularly nearly behind the bow shock with peak values, while κ P is very small in the freestream.It suggests that the radiation energy transfer mainly occurs in the high-temperature shock layer.The five-step model predicts the greatest κ P , the eight-step model second, and the two-step model the smallest.As time goes on with the altitude and velocity down, the total level of κ P grows gradually.Although the temperature in the shock layer decreases as shown in Figure 5 and makes a negative contribution to radiative transfer as the Fire II flight altitude descends, the air density increase significantly promotes the radiation effect in the flowfield.The results of the three-step models all support this point.Another interesting thing is that κ P predicted by the two-step model shows the Fire II shock layer is close to being optically thin with the optical thickness being an order of magnitude of 10 −2 , while the results of the five-step and eight-step models do not agree with this, which can only be clarified further in the future using the more detailed radiation property model, such as the line-by-line or narrowband calculations [57].Generally, the present radiation-flow uncoupling procedure is a good selection with higher efficiency in time cost than the coupling method to provide rich radiation information on the whole hypersonic nonequilibrium flowfield of Fire II reentry that the engineering methods cannot make.
Aerospace 2022, 9, x 14 of 19 uncoupling procedure is a good selection with higher efficiency in time cost than the coupling method to provide rich radiation information on the whole hypersonic nonequilibrium flowfield of Fire II reentry that the engineering methods cannot make.uncoupling procedure is a good selection with higher efficiency in time cost than the coupling method to provide rich radiation information on the whole hypersonic nonequilibrium flowfield of Fire II reentry that the engineering methods cannot make.

Conclusions
A radiation-flowfield uncoupling procedure is developed to simulate the Fire II trajectory points aimed at studying the radiative heating in the thermochemical nonequilibrium flowfield for Earth reentry at hypervelocity above 10 km/s.The radiative transfer is integrated only once by the TS approximation with the nongray step model based on the flow solution obtained using an in-house N-S solver PHAROS.It is naturally more efficient in computational cost than the coupled scheme, and also provides reasonable and more detailed information on the aerothermal environment than the engineering relations which always can only calculate the aerodynamic heat flux at the stagnation point of reentry vehicles.
The results of Fire II cases throughout the trajectory from t = 1634 s to 1645 s show that the radiative heating grows first and then decreases with the order of magnitude of 1 MW/m 2 , which is comparable to the convective heating and even exceeds the latter.Although there are remarkable differences among the two-, five-, and eight-step models, the three models all show essentially consistent trends in predictions of radiative transfer.The uncoupling calculated radiative heating can be regarded as the upper limit in the engineering application.In the future, more efforts need to be made to clarify the level of the optical thickness for the flowfield of Earth reentry vehicles at hypervelocity above 10 km/s, which the two-step model predicts to be optically thin, while the five-and eightstep models do not agree.The present scheme can also provide more radiation information in the nonequilibrium flowfield than the previous engineering relations.
Due to the fact that the flowfield establishes far more quickly than the vehicle flies down for Earth reentry, it can perform steady flow simulation at each trajectory point.Therefore, the present radiation-flow uncoupling algorithm using a nongray step model with moderate efficiency and reasonable accuracy is promising to solve the real-time problem in engineering for predicting both convective and radiative heating to the atmospheric reentry vehicle in the future.

Nomenclature e
total energy e ve vibrational-electronic energy h total enthalpy h s enthalpy for the species s h ve,s vibrational-electronic enthalpy for the species s k f,r forward reaction rate coefficients of the r-th reaction k b,r backward reaction rate coefficients of the r-th reaction p pressure q total heat flux q ve vibrational-electronic heat flux q cw,stag convective heat flux at stagnation point q rw,stag radiative heat flux at stagnation point

Figure 1 .
Figure 1.Geometry and mesh of Fire II.

Figure 2
Figure 2 compares the present CFD results of convective heat transfer at the stagnation point with those data calculated by Anderson's engineering relationship [35], Gupta et al. [37], and Olynick et al. [38].Anderson's relationship is expressed as follows:

Figure 1 .
Figure 1.Geometry and mesh of Fire II.

Figure 2
Figure 2 compares the present CFD results of convective heat transfer at the stagnation point with those data calculated by Anderson's engineering relationship [35], Gupta et al. [37], and Olynick et al. [38].Anderson's relationship is expressed as follows:

Figure 2 .
Figure 2. Convective heating at the stagnation point of Fire II predicted by the present algorithm (present), Anderson's engineering relation (Anderson Eng.), Olynick et al., and Gupta et al.

Figure 3 .
Figure 3. Forebody convective heating along the Fire II surface throughout the trajectory from t = 1634 s to t = 1645 s.

Figure 2 . 19 Figure 2 .
Figure 2. Convective heating at the stagnation point of Fire II predicted by the present algorithm (present), Anderson's engineering relation (Anderson Eng.), Olynick et al., and Gupta et al.

Figure 3 .
Figure 3. Forebody convective heating along the Fire II surface throughout the trajectory from t = 1634 s to t = 1645 s.

Figure 3 .
Figure 3. Forebody convective heating along the Fire II surface throughout the trajectory from t = 1634 s to t = 1645 s.

Figure 7 .
Figure 7. Radiative heating at the stagnation point of Fire II throughout the trajectory predicted by the present algorithm with 2-step model, 5-step model, 8-step model, T & S engineering relation, Gupta et al., and Olynick et al.

Figure 7 .
Figure 7. Radiative heating at the stagnation point of Fire II throughout the trajectory predicted by the present algorithm with 2-step model, 5-step model, 8-step model, T & S engineering relation, Gupta et al., and Olynick et al.Gupta et al. developed a RAD code accounting for the molecular band, continuum, and atomic line transitions.They employed a detailed frequency dependence of the absorption coefficients to integrate over the radiation spectrum and used the TS approximation for integrating over physical space, but only limited to the stagnation region.Olynick et al. developed a GIANTS/NOVAR code to obtain values of the stagnation radiation intensity in the 0.2 and 6.2 eV range.They used a "smeared band" approximation instead of a line-by-line approach to account for the air radiation properties of the molecular band systems and reduce the total number of spectral points.Even doing so, Olynick et al. still needed to perform TS integration at each radiation grid point consisting of around 1000 frequency points for calculating absorption and emission coefficients, which is very time-consuming.Gupta et al. and Olynick et al. both implemented the radiation-flow coupling simulation with nonequilibrium chemistry.Figure 7 shows that the radiative heat fluxes calculated by the TS approximation with the step models are mostly greater than those obtained by the T & S relation, Gupta et al., and Olynick et al.In fact, compared with the coupling simulation, the present radiation-flow uncoupling methodshould predict the higher level of radiative heating[22], which is reasonable and can be seen as the conservative upper limit in the engineering application.
Figure 7 shows that the radiative heat fluxes calculated by the TS approximation with the step models are mostly greater than those obtained by the T & S relation, Gupta et al., and Olynick et al.In fact, compared with the coupling simulation, the present radiation-flow uncoupling method

Table 1 .
Spectral regions of two-, five-, and eight-step models.

Table 2 .
Fire II flight conditions for the present study.
stoichiometric coefficient of the species s in the r-th backward reaction