ALBERT

Description: Spacecraft entering a planetary atmosphere dissipate a great deal of energy into the surrounding gas. In the frame of reference of the vehicle, the atmospheric gas suddenly decelerates from hypersonic (Mach ~5-50) to subsonic velocities. The kinetic energy of the gas is rapidly converted to thermal and chemical energy, forming a bow shock behind which a plasma with energies on the order of one electron volt (eV) is produced. The resulting shock layer relaxes from strong thermal non-equilibrium that is translationally hot but internally cold and un-ionized toward a thermochemically equilibrated plasma over a distance of a few centimeters. Composition is dependent upon the planetary atmosphere Air for Earth, CO2/N2 for Mars and Venus, N2/CH4 for Titan and H2/He/CH4 for Saturn, Neptune and Jupiter. Typical velocities of entry may range from 3-7 km/s (4-25 MJ/kg) for Titan/Mars, 8-14 km/s (30-100 MJ/kg) for Earth/Venus, and 25-40 km/s (300-800 MJ/kg) for outer planets. The equilibrium plasmas produced from these conditions are highly dissociated (up to and above 99%) and ionized (0.1- 15%), with temperatures from 7,000-15,000K and pressures from 0.1-1.0 bar. Understanding the behavior of these plasmas the way in which they approach equilibrium, how they radiate, and how they interact with materials is an active area of research necessitated by requirements to predict and test the performance of thermal protection systems (TPS) that enable spacecraft to deliver scientific instruments, and people, to foreign worlds and back to Earth. The endeavor is a multi-physics problem, with key processes highlighted in Fig. 1. This white paper describes the current state of the art in simulating shock layer plasmas both computationally and in ground test facilities. Gaps requiring further research and development are identified.

Description: Data from a pure nitrogen test series in the Electric Arc Shock Tube Facility were previously reported for velocities spanning 6-12 km/s at a free-stream pressure of 0.2 Torr. This test series provides validation data for a range of physical phenomena to investigate, including vibrational relaxation, molecular radiation, nitrogen dissociation and ionization, and atomic radiation and ionization. This paper details analysis of data obtained at a nominal velocity of 10.3 km/s. The spectra are analyzed to extract temperatures and the densities of excited states as a function of position behind the shock. The effect of different methods for calculating state populations and ionization processes is assessed, as is a rigorous assessment of the atomic line lists, with both missing and extra lines identified.

Description: 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.

Description: 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.

Description: The NEQAIR code is one of the original heritage solvers for radiative heating prediction in aerothermal environments, and is still used today for mission design purposes. This paper discusses the implementation of the first major revision to the NEQAIR code in the last five years, NEQAIR v14.0. The most notable features of NEQAIR v14.0 are the parallelization of the radiation computation, reducing runtimes by about 30, and the inclusion of mid-wave CO2 infrared radiation.

Description: 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.

Description: This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA's exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA's Entry Systems Modeling (ESM) Project, which serves as the Agency's only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specific models for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures.

Description: 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.

Description: 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.

Description: Time accurate simulation of non-equilibrium flows inside shock tube facilities presents several challenges from both physical and mathematical aspects. Furthermore, the large computational cost makes it impractical to support a real-time experimental test campaign. In this work, we explore other methods for modeling the shock tube problem with the main focus on the post-shock region and the absolute radiation emanating from it. The proposed alternative approach is several orders of magnitude less computationally expensive while still accurate enough with regards to the quantities of interest. Excellent agreement is found with the established stagnation-line approach. Comparison with time-accurate simulations shows good agreement close to the peak values and disagreement of the temperatures relaxation and radiance profiles toward equilibrium.