ALBERT

Description: The effect of nonequilibrium flow chemistry on the surface temperature distribution over the forebody heat shield on the Aeroassisted Flight Experiment (AFE) vehicle was investigated using a reacting boundary-layer code. Computations were performed by using boundary-layer-edge properties determined from global iterations between the boundary-layer code and flow field solutions from a viscous shock layer (VSL) and a full Navier-Stokes solution. Surface temperature distribution over the AFE heat shield was calculated for two flight conditions during a nominal AFE trajectory. This study indicates that the surface temperature distribution is sensitive to the nonequilibrium chemistry in the shock layer. Heating distributions over the AFE forebody calculated using nonequilibrium edge properties were similar to values calculated using the VSL program.

Description: Phenolic Impregnated Carbon Ablator (PICA) was the forebody heatshield material on the Stardust sample-return capsule and is also a primary candidate material for the Mars Science Lander (MSL), the Orion Crew Module, and the SpaceX Dragon vehicle. As part of the heatshield qualification for Orion, physical and thermal properties of virgin and charred PICA were measured, and an ablation and thermal response model was developed. We validated the model by comparing it with recession and temperature data from stagnation arcjet tests conducted over a wide range of stagnation heat flux of 107 to 1102 W/sq cm. The effect of orthotropic thermal conductivity was evident in the thermal response of the arcjet models. In general, model predictions compared well with the data; however, the uncertainty of the recession prediction was greatest for heat fluxes below 200 W/sq cm. More recent MSL testing focused on the low heat flux regime of 45 to 250 W/sq cm. The new results confirm the recession uncertainty, especially for pressures below 6 kPa. In this work we focus on improving the model predictions for MSL and Orion tests below 250 W/sq cm.

Description: Two sets of finite-rate gas-surface interaction model between air and the carbon surface are studied. The first set is an engineering model with one-way chemical reactions, and the second set is a more detailed model with two-way chemical reactions. These two proposed models intend to cover the carbon surface ablation conditions including the low temperature rate-controlled oxidation, the mid-temperature diffusion-controlled oxidation, and the high temperature sublimation. The prediction of carbon surface recession is achieved by coupling a material thermal response code and a Navier-Stokes flow code. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and Ablation Program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting full Navier-Stokes equations using Data Parallel Line Relaxation method. Recession analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities with heat fluxes ranging from 45 to 1100 wcm2 are performed and compared with data for model validation. The ablating material used in these arc-jet tests is Phenolic Impregnated Carbon Ablator. Additionally, computational predictions of surface recession and shape change are in good agreement with measurement for arc-jet conditions of Small Probe Reentry Investigation for Thermal Protection System Engineering.

Description: The fully implicit ablation and thermal response code simulates pyrolysis and ablation of thermal protection materials and systems. The governing equations, which include energy conservation, a three-component decomposition model, and a surface energy balance, are solved with a moving grid.This work describes new modeling capabilities that are added to a special version of code. These capabilities include a time-dependent pyrolysis gas flow momentum equation with Darcy-Forchheimer terms and pyrolysis gas species conservation equations with finite rate homogeneous chemical reactions. The total energy conservation equation is also enhanced for consistency with these new additions. Two groups of parametric studies of the phenolic impregnated carbon ablator are performed. In the first group, an Orion flight environment for a proposed lunar-return trajectory is considered. In the second group, various test conditions for arcjet models are examined. The central focus of these parametric studies is to understand the effect of pyrolysis gas momentum transfer on material in-depth thermal responses with finite-rate, equilibrium, or frozen homogeneous gas chemistry. Results indicate that the presence of chemical nonequilibrium pyrolysis gas flow does not significantly alter the in-depth thermal response performance predicted using the chemical equilibrium gas model.

Description: In previous work, the authors documented the Multicomponent Ablation Thermochemistry (MAT) and Fully Implicit Ablation and Thermal response (FIAT) programs. In this work, key features from MAT and FIAT were combined to create the new Fully Implicit Ablation, Thermal response, and Chemistry (FIATC) program. FIATC is fully compatible with FIAT (version 2.5) but has expanded capabilities to compute the multispecies surface chemistry and ablation rate as part of the surface energy balance. This new methodology eliminates B' tables, provides blown species fractions as a function of time, and enables calculations that would otherwise be impractical (e.g. 4+ dimensional tables) such as pyrolysis and ablation with kinetic rates or unequal diffusion coefficients. Equations and solution procedures are presented, then representative calculations of equilibrium and finite-rate ablation in flight and ground-test environments are discussed.

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: A test campaign was conducted placing meteorites in the 60 MW plasma Arcjet Interaction Heating Facility at NASA Ames Research Center, with the aim to achieve flight-relevant conditions for asteroid impacts in Earth's atmosphere and to provide insight into how meteoritic materials respond to extreme entry heating environments. The test conditions at heat flux of 4000 W/m2 and 140 kPa stagnation pressure are comparable to those experienced by a 30-meter diameter asteroid moving at 20 km/s velocity at 65 km altitude in the Earth's atmosphere. Test objects were a stony type H5 ordinary chondrite (Tamdakht) and an iron type IAB-MG meteorite (Campo Del Cielo), and included the terrestrial analogs Dense Flood Basalt and Fused Silica. All samples were exposed for only a few seconds in the plasma stream. Significant melt flow and vaporization was observed for both the stony and iron meteorites during exposure. Mass loss from spallation of fragments was also observed. Vapor emitted atomic lines from alkali metals and iron, but did not emit the expected MgO molecular band emissions. The meteoritic melts flowed more rapidly, indicating lower viscosity, than those of Fused Silica. The surface recession was mapped. The effective heat of ablation derived from this showed that ablation under these conditions occurred in the melt-dominated regime. Ablation parameters have an effect on ground damage estimates. A bias in ablation parameters towards the melt-dominated regime would imply that impacting asteroids survive to lower altitude, and therefore could possibly have airbursts with a larger ground damage footprint.

Description: A test campaign was conducted placing meteorites in the 60 MW plasma Arcjet Interaction Heating Facility at NASA Ames Research Center, with the aim to achieve flight-relevant conditions for asteroid impacts in Earth's atmosphere and to provide insight into how meteoritic materials respond to extreme entry heating environments. The test conditions at heat flux of 4000 W/m2 and 140 kPa stagnation pressure are comparable to those experienced by a 30-meter diameter asteroid moving at 20 km/s velocity at 65 km altitude in the Earth's atmosphere. Test objects were a stony type H5 ordinary chondrite (Tamdakht) and an iron type IAB-MG meteorite (Campo Del Cielo), and included the terrestrial analogs Dense Flood Basalt and Fused Silica. All samples were exposed for only a few seconds in the plasma stream. Significant melt flow and vaporization was observed for both the stony and iron meteorites during exposure. Mass loss from spallation of fragments was also observed. Vapor emitted atomic lines from alkali metals and iron, but did not emit the expected MgO molecular band emissions. The meteoritic melts flowed more rapidly, indicating lower viscosity, than those of Fused Silica. The surface recession was mapped. The effective heat of ablation derived from this showed that ablation under these conditions occurred in the melt-dominated regime. Ablation parameters have an effect on ground damage estimates. A bias in ablation parameters towards the melt-dominated regime would imply that impacting asteroids survive to lower altitude, and therefore could possibly have airbursts with a larger ground damage footprint.

Description: A fast code was developed to calculate the forebody heating environment and heat shielding that is required for Jupiter atmospheric entry probes. A carbon phenolic heat shield material was assumed and, since computational efficiency was a major goal, analytic expressions were used, primarily, to calculate the heating, ablation and the required insulation. The code was verified by comparison with flight measurements from the Galileo probe's entry; the calculation required 3.5 sec of CPU time on a work station. The computed surface recessions from ablation were compared with the flight values at six body stations. The average, absolute, predicted difference in the recession was 12.5% too high. The forebody's mass loss was overpredicted by 5.5% and the heat shield mass was calculated to be 15% less than the probe's actual heat shield. However, the calculated heat shield mass did not include contingencies for the various uncertainties that must be considered in the design of probes. Therefore, the agreement with the Galileo probe's values was considered satisfactory, especially in view of the code's fast running time and the methods' approximations.

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