Validating Advanced Thermophysics Models
Aerospace America year in review article for Thermophysics covering the COMARS data from ExoMars, meteor testing in arcj ets and potential energy surface calculations.
PES↗
Engineering topics
Brandis, Aaron M.
Publications and source records attributed to Brandis, Aaron M..
Measurement and Prediction of Radiative Non-Equilibrium for Air Shocks Between 7-9 km/s
The present paper describes a recent characterization of thermochemical non-equilibrium for shock speeds between 7 and 9 km/s in the NASA Ames Electric Arc Shock Tube (EAST) Facility. Data are spectrally resolved from 190-1450 nm and spatially resolved behind the shock front. The data are analyzed in terms of a spectral non-equilibrium metric, defined as the average radiance within +/- 2 cm of the peak. Simulations with DPLR/NEQAIR using different rate chemistries show these conditions to be poorly replicated. The sources of discrepancy are examined, leading to an update to the NEQAIR non-Boltzmann model and DPLR rate chemistry. New parameters for the rate chemistry and non-Boltzmann modeling are reported.
Cruden, Brett A.↗
Experimental and Numerical Investigation of Air Radiation in Superorbital Expanding Flow
To investigate air radiation in expanding flows and provide experimental data for validating associated computational models, experiments were conducted in the X2 expansion tunnel facility at the Centre for Hypersonics of the University of Queensland. A 54 turning angle wedge model was employed to generate steady expanding flows with in flow total enthalpies of 50.7, 63.4 and 75.4 MJkg. VUV spectra from 118 to 180 nm were acquired across the wedge at three equispaced distances away from the top of the model, as well as through its top surface. High speed filtered images were also obtained by coupling a Shimadzu 1 MHz high speed camera to a bandpass filter to obtain calibrated images of the 777 nm oxygen triplet. Both the across-wedge VUV spectra and filtered images of the 777 nm atomic oxygen were compared with NEQAIR simulations, which were performed using flow field data from two-dimensional CFD simulations with two-temperature 11-species air chemistry utilizing the in-house Navier-Stokes flow solver Eilmer3. Data extracted from consecutive frames of the filtered high speed images confirmed up to 8 s of available test time for the flow conditions tested. For the strongly radiating 149 and 174 nm atomic nitrogen lines, large disagreement between experimental data and NEQAIR predictions can be observed from the start of the expansion fan where the electron-ion recombination process commences. The spatial extent, or spans of the radiance profiles of the 149 and 174 nm N lines are significantly under predicted by NEQAIR, and are very close to those of N, N+ and electron number density profiles, which follow that of flow density. The electron-ion recombination process is proposed as the main reason for these discrepancies. The comparisons between NEQAIR simulations and filtered images of the 777 nm oxygen triplet show good agreement in the post-shock compression region and the start of the expansion fan for the 63.4 MJkg condition, but with up to a factor of three over prediction by NEQAIR further downstream, which is attributed to electron-impact excitation. Similar trends are found with the 75.4 MJkg condition, with reduced level of agreement in the compression region, which can be due to uncertainties in inflow condition.
Earth entry↗
Measurement and Prediction of Radiative Non-Equilibrium for Air Shocks Between 7-9 km/s
The present paper describes a recent characterization of thermochemical non-equilibrium for shock speeds between 7 and 9 km/s in the NASA Ames Electric Arc Shock Tube (EAST) Facility. Data are spectrally resolved from 190-1450 nm and spatially resolved behind the shock front. The data are analyzed in terms of a spectral non-equilibrium metric, defined as the average radiance within +/-2 cm of the peak. Simulations with DPLR/NEQAIR using different rate chemistries show these conditions to be poorly replicated. The sources of discrepancy are examined, leading to an update to the NEQAIR non-Boltzmann model and DPLR rate chemistry. New parameters for the rate chemistry and non-Boltzmann modeling are reported.
Cruden, Brett A.↗
Experimental and Numerical Investigation of Air Radiation in Superorbital Expanding Flow
To investigate air radiation in expanding flows and provide experimental data for validating associated computational models, experiments were conducted in the X2 expansion tunnel facility at the Centre for Hypersonics of the University of Queensland. A 54deg turning angle wedge model was employed to generate steady expanding flows with in flow total enthalpies of 50.7, 63.4 and 75.4 MJ/kg. VUV spectra from 118 to 180 nm were acquired across the wedge at three equispaced distances away from the top of the model, as well as through its top surface. High speed filtered images were also obtained by coupling a Shimadzu 1 MHz high speed camera to a bandpass filter to obtain calibrated images of the 777 nm oxygen triplet. Both the across-wedge VUV spectra and filtered images of the 777 nm atomic oxygen were compared with NEQAIR simulations, which were performed using flow field data from two-dimensional CFD simulations with two-temperature 11-species air chemistry utilising the in-house Navier-Stokes flow solver Eilmer3. Data extracted from consecutive frames of the filtered high speed images confirmed up to 8 s of available test time for the flow conditions tested. For the strongly radiating 149 and 174 nm atomic nitrogen lines, large disagreement between experimental data and NEQAIR predictions can be observed from the start of the expansion fan where the electron-ion recombination process commences. The spatial extent, or spans of the radiance profiles of the 149 and 174 nm N lines are significantly underpredicted by NEQAIR, and are very close to those of N, N+ and electron number density profiles, which follow that of flow density. The electron-ion recombination process is proposed as the main reason for these discrepancies. The comparisons between NEQAIR simulations and filtered images of the 777 nm oxygen triplet show good agreement in the post-shock compression region and the start of the expansion fan for the 63.4 MJ/kg condition, but with up to a factor of three overprediction by NEQAIR further downstream, which is attributed to electron-impact excitation. Similar trends are found with the 75.4 MJ/kg condition, with reduced level of agreement in the compression region, which can be due to uncertainties in inflow condition.
expansion tube↗
Measurement of Radiative Non-Equilibrium for Air Shocks Between 7-9 Km/s
This paper describes a recent characterization of non-equilibrium radiation for shock speeds between 7 and 9 km/s in the NASA Ames Electric Arc Shock Tube (EAST) Facility. Data is spectrally resolved from 190- 1450 nm and spatially resolved behind the shock front. Comparisons are made to DPLR/NEQAIR simulations using different modeling options and recommendations for future study are made based on these comparisons.
Radiation↗
Advancements in Afterbody Radiative Heating Simulations for Earth Entry
Four advancements to the simulation of backshell radiative heating for Earth entry are presented. The first of these is the development of a flow field model that treats electronic levels of the dominant backshell radiator, N, as individual species. This is shown to allow improvements in the modeling of electron-ion recombination and two-temperature modeling, which are shown to increase backshell radiative heating by 10 to 40%. By computing the electronic state populations of N within the flow field solver, instead of through the quasi-steady state approximation in the radiation code, the coupling of radiative transition rates to the species continuity equations for the levels of N, including the impact of non-local absorption, becomes feasible. Implementation of this additional level of coupling between the flow field and radiation codes represents the second advancement presented in this work, which is shown to increase the backshell radiation by another 10 to 50%. The impact of radiative transition rates due to non-local absorption indicates the importance of accurate radiation transport in the relatively complex flow geometry of the backshell. This motivates the third advancement, which is the development of a ray-tracing radiation transport approach to compute the radiative transition rates and divergence of the radiative flux at every point for coupling to the flow field, therefore allowing the accuracy of the commonly applied tangent-slab approximation to be assessed for radiative source terms. For the sphere considered at lunar-return conditions, the tangent-slab approximation is shown to provide a sufficient level of accuracy for the radiative source terms, even for backshell cases. This is in contrast to the agreement between the two approaches for computing the radiative flux to the surface, which differ by up to 40%. The final advancement presented is the development of a nonequilibrium model for NO radiation, which provides significant backshell radiation at velocities below 10 km/s. The developed model reduces the nonequilibrium NO radiation by 50% relative to the previous model.
Johnston, Christopher O.↗
DSMC Shock Simulation of Saturn Entry Probe Conditions
This work describes the direct simulation Monte Carlo (DSMC) investigation of Saturn entry probe scenarios and the influence of non-equilibrium phenomena on Saturn entry conditions. The DSMC simulations coincide with rarefied hypersonic shock tube experiments of a hydrogen-helium mixture performed in the Electric Arc Shock Tube (EAST) at the NASA Ames Research Center. The DSMC simulations are post-processed through the NEQAIR line-by-line radiation code to compare directly to the experimental results. Improved collision cross-sections, inelastic collision parameters, and reaction rates are determined for a high temperature DSMC simulation of a 7-species H2-He mixture and an electronic excitation model is implemented in the DSMC code. Simulation results for 27.8 and 27.4 km/s shock waves are obtained at 0.2 and 0.1 Torr, respectively, and compared to measured spectra in the VUV, UV, visible, and IR ranges. These results confirm the persistence of non-equilibrium for several centimeters behind the shock and the diffusion of atomic hydrogen upstream of the shock wave. Although the magnitude of the radiance did not match experiments and an ionization inductance period was not observed in the simulations, the discrepancies indicated where improvements are needed in the DSMC and NEQAIR models.
Higdon, Kyle J.↗
DSMC Shock Simulation of Saturn Entry Probe Conditions
This work describes the direct simulation Monte Carlo (DSMC) investigation of Saturn entry probe scenarios and the influence of non-equilibrium phenomena on Saturn entry conditions. The DSMC simulations coincide with rarefied hypersonic shock tube experiments of a hydrogen-helium mixture performed in the Electric Arc Shock Tube (EAST) at NASA Ames Research Center. To directly compare to the experimental results, the DSMC simulations are post-processed through the NEQAIR line-by-line radiation code. Improved collision cross-sections, inelastic collision parameters, and reaction rates are determined for a high temperature DSMC simulation of a 7-species H2-He mixture and an electronic excitation model is implemented in the DSMC code. Simulation results for 27.8 and 27.4 kms shock waves are obtained at 0.2 and 0.1 Torr respectively and compared to measured spectra in the VUV, UV, visible, and IR ranges. These results confirm the persistence of non-equilibrium for several centimeters behind the shock and the diffusion of atomic hydrogen upstream of the shock wave. Although the magnitude of the radiance did not match experiments and an ionization inductance period was not observed in the simulations, the discrepancies indicated where improvements are needed in the DSMC and NEQAIR models.
Saturn↗
Measurement of Ultraviolet Radiative Heating Augmentation in HIEST Reflected Shock Tunnel
Radiance measurements in air at enthalpies from 8-20 MJkg have been made over a 250mm diameter flat-faced test article in Japan Aerospace Exploration Agency's HIgh-Enthalpy Shock Tunnel (HIEST). Measurements were made in the ultraviolet region (200-400 nm wavelength) in an attempt to resolve the long-standing discrepancy between theoryand measurements of heat flux over a blunt body; this discrepancy is often attributed toradiation. The spectra obtained indicate the presence of atomic iron vapor in the flowfield.At the highest enthalpies, the radiance is at the blackbody limit. An attempt to model theradiance is made by taking a nominal CFD flowfield without any contamination productsand processing it through a line-by-line radiation simulation tool. Iron vapor is introducedinto the shocked gas ahead of the model and radiation computations are repeated; the molefraction of iron vapor is adjusted to match the data. For the higher enthalpy conditions, theradiance was strongly absorbed and it was necessary to adjust the temperature and NOdensity in the freestream to match the signal below 300 nm. Once the observed spectrawere satisfactorily matched, the radiance to the stagnation point was then computed. It isshown that the impurity radiation is sufficiently large to explain the discrepancy.
Shock Tunnel↗
Physics-Based Modeling of Meteor Entry and Breakup
A new research effort at NASA Ames Research Center has been initiated in Planetary Defense, which integrates the disciplines of planetary science, atmospheric entry physics, and physics-based risk assessment. This paper describes work within the new program and is focused on meteor entry and breakup. Over the last six decades significant effort was expended in the US and in Europe to understand meteor entry including ablation, fragmentation and airburst (if any) for various types of meteors ranging from stony to iron spectral types. These efforts have produced primarily empirical mathematical models based on observations. Weaknesses of these models, apart from their empiricism, are reliance on idealized shapes (spheres, cylinders, etc.) and simplified models for thermal response of meteoritic materials to aerodynamic and radiative heating. Furthermore, the fragmentation and energy release of meteors (airburst) is poorly understood. On the other hand, flight of human-made atmospheric entry capsules is well understood. The capsules and their requisite heatshields are designed and margined to survive entry. However, the highest speed Earth entry for capsules is less than 13 km/s (Stardust). Furthermore, Earth entry capsules have never exceeded diameters of 5 m, nor have their peak aerothermal environments exceeded 0.3 atm and 1 kW/cm2. The aims of the current work are: (i) to define the aerothermal environments for objects with entry velocities from 13 to greater than 20 km/s; (ii) to explore various hypotheses of fragmentation and airburst of stony meteors in the near term; (iii) to explore the possibility of performing relevant ground-based tests to verify candidate hypotheses; and (iv) to quantify the energy released in airbursts. The results of the new simulations will be used to anchor said risk assessment analyses. With these aims in mind, state-of-the-art entry capsule design tools are being extended for meteor entries. We describe: (i) applications of current simulation tools to spherical geometries of diameters ranging from 1 to 100 m for an entry velocity of 20 km/s and stagnation pressures ranging from 1 to 100 atm; (ii) the influence of shape and departure of heating environment predictions from those for a simple spherical geometry; (iii) assessment of thermal response models for silica subject to intense radiation; and (iv) results for porosity-driven gross fragmentation of meteors, idealized as a collection of smaller objects. Lessons learned from these simulations will be used to help understand the Chelyabinsk meteor entry up to its first point of fragmentation.
Aerothermodynamics↗
Physics-Based Modeling of Meteor Entry and Breakup
A new research effort at NASA Ames Research Center has been initiated in Planetary Defense, which integrates the disciplines of planetary science, atmospheric entry physics, and physics-based risk assessment. This paper describes work within the new program and is focused on meteor entry and breakup. Over the last six decades significant effort was expended in the US and in Europe to understand meteor entry including ablation, fragmentation and airburst (if any) for various types of meteors ranging from stony to iron spectral types. These efforts have produced primarily empirical mathematical models based on observations. Weaknesses of these models, apart from their empiricism, are reliance on idealized shapes (spheres, cylinders, etc.) and simplified models for thermal response of meteoritic materials to aerodynamic and radiative heating. Furthermore, the fragmentation and energy release of meteors (airburst) is poorly understood. On the other hand, flight of human-made atmospheric entry capsules is well understood. The capsules and their requisite heat shields are designed and margined to survive entry. However, the highest speed Earth entry for capsules is 13 kms (Stardust). Furthermore, Earth entry capsules have never exceeded diameters of 5 m, nor have their peak aerothermal environments exceeded 0.3 atm and 1 kWcm2. The aims of the current work are: (i) to define the aerothermal environments for objects with entry velocities from 13 to 20 kms; (ii) to explore various hypotheses of fragmentation and airburst of stony meteors in the near term; (iii) to explore the possibility of performing relevant ground-based tests to verify candidate hypotheses; and (iv) to quantify the energy released in airbursts. The results of the new simulations will be used to anchor said risk assessment analyses.With these aims in mind, state-of-the-art entry capsule design tools are being extended for meteor entries. We describe: (i) applications of current simulation tools to spherical geometries of diameters ranging from 1 to 100 m for an entry velocity of 20 kms and stagnation pressures ranging from 1 to 100 atm; (ii) the influence of shape and departure of heating environment predictions from those for a simple spherical geometry; (iii) assessment of thermal response models for silica subject to intense radiation; and (iv) results for porosity-driven gross fragmentation of meteors, idealized as a collection of smaller objects. Lessons learned from these simulations will be used to help understand the Chelyabinsk meteor entry up to its first point of fragmentation.
Atmospheric Entry↗
Physics-Based Modeling of Meteor Entry and Breakup
A new research effort at NASA Ames Research Center has been initiated in Planetary Defense, which integrates the disciplines of planetary science, atmospheric entry physics, and physics-based risk assessment. This paper describes work within the new program and is focused on meteor entry and breakup.Over the last six decades significant effort was expended in the US and in Europe to understand meteor entry including ablation, fragmentation and airburst (if any) for various types of meteors ranging from stony to iron spectral types. These efforts have produced primarily empirical mathematical models based on observations. Weaknesses of these models, apart from their empiricism, are reliance on idealized shapes (spheres, cylinders, etc.) and simplified models for thermal response of meteoritic materials to aerodynamic and radiative heating. Furthermore, the fragmentation and energy release of meteors (airburst) is poorly understood.On the other hand, flight of human-made atmospheric entry capsules is well understood. The capsules and their requisite heatshields are designed and margined to survive entry. However, the highest speed Earth entry for capsules is 13 kms (Stardust). Furthermore, Earth entry capsules have never exceeded diameters of 5 m, nor have their peak aerothermal environments exceeded 0.3 atm and 1 kW/sq cm. The aims of the current work are: (i) to define the aerothermal environments for objects with entry velocities from 13 to 20 kms; (ii) to explore various hypotheses of fragmentation and airburst of stony meteors in the near term; (iii) to explore the possibility of performing relevant ground-based tests to verify candidate hypotheses; and (iv) to quantify the energy released in airbursts. The results of the new simulations will be used to anchor said risk assessment analyses. With these aims in mind, state-of-the-art entry capsule design tools are being extended for meteor entries. We describe: (i) applications of current simulation tools to spherical geometries of diameters ranging from 1 to 100 m for an entry velocity of 20 kms and stagnation pressures ranging from 1 to 100 atm; (ii) the influence of shape and departure of heating environment predictions from those for a simple spherical geometry; (iii) assessment of thermal response models for silica subject to intense radiation; and (iv) results for porosity-driven gross fragmentation of meteors, idealized as a collection of smaller objects. Lessons learned from these simulations will be used to help understand the Chelyabinsk meteor entry up to its first point of fragmentation.
Near-Earth Asteroids↗
Status of Shock Layer Radiation Validation Studies
This paper presents a status update for the shock layer radiation validation studies conducted at NASA. A review of the present capability for the simulation and validation of shock layer radiation is presented as well as providing an overview of the data obtained from the Electric Arc Shock Tube (EAST). The paper will include details covering updated convective and radiative heating correlations, provide an overview of the development of new kinetics for Mars entry and detail some recent work calculating after-body radiation. Furthermore, the paper will highlight conditions where there is high confidence in the validation of EAST data (e.g. Earth entry for speeds greater than approximately 10 kms and for many Mars entry conditions) and where further experimental data would be highly beneficial (e.g. lower speed Earth entry around 7.5 to 10 kms and higher speed CO2 entries relevant to Venus). Nominal test conditions for both Earth and Mars are provided for future potential facility-to-facility comparisons.
Entry systems modeling↗
Radiative Heating in MSL Entry: Comparison of Flight Heating Discrepancy to Ground Test and Predictive Models
During the recent entry of the Mars Science Laboratory (MSL), the heat shield was equipped with thermocouple stacks to measure in-depth heating of the thermal protection system (TPS). When only convective heating was considered, the derived heat flux from gauges in the stagnation region was found to be underpredicted by as much as 17 W/sq cm, which is significant compared to the peak heating of 32 W/sq cm. In order to quantify the contribution of radiative heating phenomena to the discrepancy, ground tests and predictive simulations that replicated the MSL entry trajectory were performed. An analysis is carried through to assess the quality of the radiation model and the impact to stagnation line heating. The impact is shown to be significant, but does not fully explain the heating discrepancy.
MSL↗
Shock Layer Radiation Modeling and Uncertainty for Mars Entry
A model for simulating nonequilibrium radiation from Mars entry shock layers is presented. A new chemical kinetic rate model is developed that provides good agreement with recent EAST and X2 shock tube radiation measurements. This model includes a CO dissociation rate that is a factor of 13 larger than the rate used widely in previous models. Uncertainties in the proposed rates are assessed along with uncertainties in translational-vibrational relaxation modeling parameters. The stagnation point radiative flux uncertainty due to these flowfield modeling parameter uncertainties is computed to vary from 50 to 200% for a range of free-stream conditions, with densities ranging from 5e-5 to 5e-4 kg/m3 and velocities ranging from of 6.3 to 7.7 km/s. These conditions cover the range of anticipated peak radiative heating conditions for proposed hypersonic inflatable aerodynamic decelerators (HIADs). Modeling parameters for the radiative spectrum are compiled along with a non-Boltzmann rate model for the dominant radiating molecules, CO, CN, and C2. A method for treating non-local absorption in the non-Boltzmann model is developed, which is shown to result in up to a 50% increase in the radiative flux through absorption by the CO 4th Positive band. The sensitivity of the radiative flux to the radiation modeling parameters is presented and the uncertainty for each parameter is assessed. The stagnation point radiative flux uncertainty due to these radiation modeling parameter uncertainties is computed to vary from 18 to 167% for the considered range of free-stream conditions. The total radiative flux uncertainty is computed as the root sum square of the flowfield and radiation parametric uncertainties, which results in total uncertainties ranging from 50 to 260%. The main contributors to these significant uncertainties are the CO dissociation rate and the CO heavy-particle excitation rates. Applying the baseline flowfield and radiation models developed in this work, the radiative heating for the Mars Pathfinder probe is predicted to be nearly 20 W/cm2. In contrast to previous studies, this value is shown to be significant relative to the convective heating.
Johnston, Christopher O.↗
Assessment of Radiative Heating Uncertainty for Hyperbolic Earth Entry
This paper investigates the shock-layer radiative heating uncertainty for hyperbolic Earth entry, with the main focus being a Mars return. In Part I of this work, a baseline simulation approach involving the LAURA Navier-Stokes code with coupled ablation and radiation is presented, with the HARA radiation code being used for the radiation predictions. Flight cases representative of peak-heating Mars or asteroid return are de ned and the strong influence of coupled ablation and radiation on their aerothermodynamic environments are shown. Structural uncertainties inherent in the baseline simulations are identified, with turbulence modeling, precursor absorption, grid convergence, and radiation transport uncertainties combining for a +34% and ..24% structural uncertainty on the radiative heating. A parametric uncertainty analysis, which assumes interval uncertainties, is presented. This analysis accounts for uncertainties in the radiation models as well as heat of formation uncertainties in the flow field model. Discussions and references are provided to support the uncertainty range chosen for each parameter. A parametric uncertainty of +47.3% and -28.3% is computed for the stagnation-point radiative heating for the 15 km/s Mars-return case. A breakdown of the largest individual uncertainty contributors is presented, which includes C3 Swings cross-section, photoionization edge shift, and Opacity Project atomic lines. Combining the structural and parametric uncertainty components results in a total uncertainty of +81.3% and ..52.3% for the Mars-return case. In Part II, the computational technique and uncertainty analysis presented in Part I are applied to 1960s era shock-tube and constricted-arc experimental cases. It is shown that experiments contain shock layer temperatures and radiative ux values relevant to the Mars-return cases of present interest. Comparisons between the predictions and measurements, accounting for the uncertainty in both, are made for a range of experiments. A measure of comparison quality is de ned, which consists of the percent overlap of the predicted uncertainty bar with the corresponding measurement uncertainty bar. For nearly all cases, this percent overlap is greater than zero, and for most of the higher temperature cases (T >13,000 K) it is greater than 50%. These favorable comparisons provide evidence that the baseline computational technique and uncertainty analysis presented in Part I are adequate for Mars-return simulations. In Part III, the computational technique and uncertainty analysis presented in Part I are applied to EAST shock-tube cases. These experimental cases contain wavelength dependent intensity measurements in a wavelength range that covers 60% of the radiative intensity for the 11 km/s, 5 m radius flight case studied in Part I. Comparisons between the predictions and EAST measurements are made for a range of experiments. The uncertainty analysis presented in Part I is applied to each prediction, and comparisons are made using the metrics defined in Part II. The agreement between predictions and measurements is excellent for velocities greater than 10.5 km/s. Both the wavelength dependent and wavelength integrated intensities agree within 30% for nearly all cases considered. This agreement provides confidence in the computational technique and uncertainty analysis presented in Part I, and provides further evidence that this approach is adequate for Mars-return simulations. Part IV of this paper reviews existing experimental data that include the influence of massive ablation on radiative heating. It is concluded that this existing data is not sufficient for the present uncertainty analysis. Experiments to capture the influence of massive ablation on radiation are suggested as future work, along with further studies of the radiative precursor and improvements in the radiation properties of ablation products.
Johnston, Christopher O.↗