Results for “NEQAIR” — Page 4

NASA NTRS2017

Applications of Radiative Heating for Space Exploration

Vehicles entering planetary atmospheres at high speeds (6 - 12 kms) experience intense heating by flows with temperatures of the order 10 000K. The flow around the vehicle experiences significant dissociation and ionization and is characterized by thermal and chemical non-equilibrium near the shock front, relaxing toward equilibrium. Emission from the plasma is intense enough to impart a significant heat flux on the entering spacecraft, making it necessary to predict the magnitude of radiative heating. Shock tubes represent a unique method capable of characterizing these processes in a flight-similar environment. The Electric Arc Shock tube (EAST) facility is one of the only facilities in its class, able to produce hypersonic flows at speeds up to Mach 50. This talk will review the characterization of radiation measured in EAST with simulations by the codes DPLR and NEQAIR, and in particular, focus on the impact these analyses have on recent missions to explore the solar system.

simulations↗

NASA NTRS2018

Overview of Recent EAST Testing, Modeling & Analysis

Shock Layer Radiation at NASA Ames - Background: Complex aerothermal and thermochemical phenomena of planetary entry define convective and radiative heating. A spacecraft's TPS (Thermal Protective System) mitigates heat transfer to substructure. Successful TPS design relies on verifiable characterization of these phenomena in the anticipated flight environment. - Approach: EAST (Electric Arc Shock Tube) [facility at Ames] simulates high-enthalpy, real-gas phenomena encountered by hypersonic vehicles entering planetary atmospheres by spectrally imaging the flow behind a moving shock wave. - Goal: Validate aerothermal models (DPLR (Data Parallel Line Relaxation Code) and NEQAIR (Nonequilibrium Radiative Transport and Spectra Program)), inform model improvements, reduce uncertainty and quantify design uncertainties. - Recent Relevant Projects: MSL (Mars Science Laboratory) & Mars 2020, InSight, OSIRIS-REx (Origins Spectral Interpretation Resource Identification Security Regolith Explorer), Orion EFT-1 (Exploration Flight-1) and EM-1 (Exploration Mission-1) and New Frontiers.

Brandis, Aaron↗

NASA NTRS2018

Aerothermodynamics for Dragonfly's Titan Entry

Dragonfly is a proposed spacecraft and mission that would send a mobile robotic rotorcraft lander to Titan, the largest moon of Saturn, in order to study prebiotic chemistry and extraterrestrial habitability at various locations. Titan is unique in having an abundant, complex, and diverse carbon-rich chemistry on the surface of a water-ice-dominated world with an interior water ocean, making it a high priority target for astrobiology and origin of life studies. The mission was initially proposed in April 2017 to NASA's New Frontiers program by the Johns Hopkins Applied Physics Laboratory. In December 2017, it was selected as one of two finalists (out of twelve proposals) to further refine the mission's concept. NASA Ames Research Center and NASA Langley Research Center are partnering as the leads for the Dragonfly's entry system to provide the completed EDL Assembly. The aerothermal analysis for Dragonfly utilizes four simulation tools from NASA Ames Research Center. Traj for calculating the trajectory, DPLR 4.04.0 for calculating the flowfield around the vehicle and convective heating, NEQAIR V15.0 for calculating the radiative heating, and FIAT for calculating the material response and thermal protection system (TPS) sizing for the heatshield. The entry conditions are relatively benign and can readily be accommodated with a tiled PICA heatshield similar to MSL and a number of flight proven materials for the backshell. This work will demonstrate that the aerothermal entry environments can be readily solved using heritage materials and techniques.

New Frontiers↗

NASA NTRS2018

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.

Sharma, Maitreyee↗

NASA NTRS2018

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for a range of high speed Earth entry experiments conducted in the Electric Arc Shock Tube at NASA Ames. The nonequilibrium radiative transport equations are solved in NEQAIR using flow field variables from the full facility CFD simulations of the EAST shock tube performed by US3D ow solver. These physics-based flow calculations lead to a significantly different post-shock gas state and associated radiation field as compared to that based on a simplified but computationally inexpensive calculation for flow over a blunt-body with appropriate initial conditions. The radiation spectra and radiance profiles are computed for an extensive range of wavelengths, from deep VUV to IR, which are pertinent to the emission characteristics of high enthalpy shock waves in air. The radiation properties of the shocked gas are calculated both in the nonequilibrium region at the shock, and in the equilibrium region behind the shock. Numerical predictions are found to be consistent with the experimental observations.

Simulations↗

NASA NTRS2019

InSight's Reconstructed Aerothermal Environments

The InSight Mars Lander successfully landed on the surface on November 26, 2018. This poster will describe the methodologies and margins used in developing the aerothermal environments for design of the thermal protection systems (TPS), as well as a prediction of as-flown environments based on the best estimated trajectory. The InSight mission spacecraft design approach included the effects of radiant heat flux to the aft body from the wake for the first time on a US Mars Mission, due to overwhelming evidence in ground testing for the European ExoMars mission (2009/2010) [1] and 2010 tests in the Electric Arc Shock Tube (EAST) facility [2]. The radiant energy on an aftbody was also recently confirmed via measurement on the Schiaparelli mission [3]. In addition, the InSight mission expected to enter the Mars atmosphere during the dust storm season, so the heatshield TPS was designed to accommodate the extra recession due to the potential dust impact. This poster will compare the predicted aerothermal environments using the reconstructed best estimated trajectory to the design environments. Design Approach: The InSight spacecraft was planned to be a near-design-to-print copy of the Phoenix spacecraft. The determination of the heatshield TPS requirements was approached as if it was a new design due to the new requirement of flying through a dust storm. The baseline for aftbody was build-to-print, and all analyses focused on ensuring adequate margin. This proved to be a challenge because the Phoenix aftbody was designed to withstand only convective heating and the InSight aftbody was evaluated for both convective and radiative heating. Aerothermal environments were predicted using the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) and the Data Parallel Line Relaxation (DPLR) CFD codes, and the Nonequilibrium Radiative Transport and Spectra Program (NEQAIR) utilizing bounding design trajectories derived from Monte Carlo analyses from the Program to Optimize Simulated Trajectories II (POST2). In all cases, super-catalytic flowfields were assigned to ensure the most conservative heating results. Two trajectories were evaluated: 1) the trajectory with the maximum heat flux was utilized to determine the flowfield characteristics and the viability of the selection of TPS materials; and 2) the trajectory with the maximum heat load was used to determine the required thicknesses of the TPS materials. Evaluation of the MEDLI data [4], along with ground test data [5] led to the determination of whether or not the flow would transition from laminar to turbulent on the heatshield, which also determined the TPS sizing location for the heatshield. Aerothermal margins were added for the convective heating and developed for the radiative heating. TPS material sizing was determined with the Reaction Kinetic Ablation Program (REKAP) and the Fully Implicit Ablation and Thermal Analysis program (FIAT) using a three-branched approach to account for aerothermal, material response, and material properties uncertainties. In addition, the heatshield recession was augmented by an analysis of the effect of entry through a potential dusty atmosphere using a methodology developed in References [6] and [7]. These analyses resulted in an increase to the Phoenix heatshield TPS thickness. Reconstruction Efforts: Once the best estimated trajectory is reconstructed by the team, the LAURA/HARA (High-Temperature Aerothermo-dynamic Radiation model) and DPLR/NEQAIR code pairs will be used to predict the as-flown aerothermal conditions. In these runs, fully-catalytic flowfields will be assigned because it is a more physically accurate description of the chemistry in the flow. Once again, determination of the onset of turbulence on the heatshield will be evaluated. The as-flown aerothermal environments will then be compared to the design environments.

Beck, R. A.↗

NASA NTRS2018

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.

Sharma, Maitreyee↗

NASA NTRS2019

Radiation Modeling for the Reentry of the Hayabusa Sample Return Capsule

Predicted shock-layer emission signatures during the reentry of the Japanese Hayabusa capsule are presented and compared with flight measurements conducted during an airborne observation mission in NASA's DC-8 Airborne Laboratory. For selected altitudes at 11 points along the flight trajectory of the capsule, lines of sight were extracted from flow field solutions computed using the in-house high-fidelity CFD code, DPLR. These lines of sight were used as inputs for the line-by-line radiation code NEQAIR, and emission spectra of the air plasma were computed in the wavelength range from 300 nm to 1600 nm, a range which covers all of the different experiments onboard the DC-8. In addition, the computed flow field solutions were post-processed with the material thermal response code FIAT, and the resulting surface temperatures of the heat shield were used to generate thermal emission spectra based on Planck radiation. Both spectra were summed and integrated over the flow field. The resulting emission at each trajectory point was propagated to the DC-8 position and transformed into incident irradiance to be finally compared with experimental data.

shock-layer radiation↗

NASA NTRS2019

Reconstruction of Schiaparelli and Comars Flight Data

ESA recently flew an entry, descent, and landing demonstrator module called Schiaparelli that entered the atmosphere of Mars on the 19th of October, 2016. The instrumentation suite included heatshield and backshell pressure transducers and thermocouples (known as AMELIA) and backshell radiation and direct heatflux-sensing sensors (known as COMARS and ICOTOM). Due to the failed landing of Schiaparelli, only a subset of the flight data was transmitted before and after plasma black-out. The goal of this paper is to present comparisons of the flight data with calculations from NASA simulation tools, DPLR/NEQAIR and LAURA/HARA. DPLR and LAURA are used to calculate the flowfield around the vehicle and surface properties, such as pressure and convective heating. The flowfield data are passed to NEQAIR and HARA to calculate the radiative heat flux. Comparisons will be made to the COMARS total heat flux, radiative heat flux and pressure measurements. Results will also be shown against the reconstructed heat flux which was calculated from an inverse analysis of the AMELIA thermocouple data performed by Astrium. Preliminary calculations are presented in this abstract. The aerodynamics of the vehicle and certain as yet unexplained features of the inverse analysis and forebody data will be investigated.

Schiaparelli↗

NASA NTRS2019

Current Status of Shock Layer Radiation Studies for Planetary Probes

Radiative heating by the shock layer gases surrounding an entry probe is an important heating mechanism for many planetary entries. The magnitude of radiative heating for planetary atmospheres may be characterized with flight similarity in a shock tube capable of obtaining flight relevant velocities and atmospheric density/composition. Over the last decade, studies performed in the Electric Arc Shock Tube (EAST) at NASA Ames have characterized the radiative heating magnitudes for most relevant planetary destinations: Mars, Venus, Titan, Saturn, Uranus. Predictions of heating for a given entry trajectory are simulated by a combination of hypersonic flowfield modeling and the radiation code, NEQAIR, and then compared to the EAST test data. This paper will discuss the measurements, summarize the current simulation capability and provide indications for destination specific radiative heating uncertainties.

Cruden, Brett A.↗

NASA NTRS2020

InSight Aerothermal Environment Assessment

The InSight spacecraft was proposed to be a build-to-print copy of the Phoenix vehicle due to the knowledge that the lander payload would be similar and the trajectory would be similar. However, the InSight aerothermal analysts, based on tests performed in CO2 during the Mars Science Laboratory mission (MSL) and completion of Russian databases, considered radiative heat flux to the aftbody from the wake for the first time for a US Mars mission. The combined convective and radiative heat flux was used to determine if the as-flown Phoenix thermal protection system (TPS) design would be sufficient for InSight. All analyses showed that the design would be adequate. Once the InSight lander was successfully delivered to Mars on November 26, 2018, work began to reconstruct the atmosphere and trajectory in order to evaluate the aerothermal environments that were actually encountered by the spacecraft and to compare them to the design environments.The best estimated trajectory (BET) reconstructed for the InSight atmospheric entry fell between the two trajectories considered for the design, when looking at the velocity versus altitude values. The maximum heat rate design trajectory (MHR) flew at a higher velocity and the maximum heat load design trajectory (MHL) flew at a lower velocity than the BET. For TPS sizing, the MHL trajectory drove the design. Reconstruction has shown that the BET flew for a shorter time than either of the design environments, hence total heat load on the vehicle should have been less than used in design. Utilizing the BET, both DPLR and LAURA were first run to analyze the convective heating on the vehicle with no angle of attack. Both codes were run with axisymmetric, laminar flow in radiative equilibrium and vibrational non-equilibrium with a surface emissivity of 0.8. Eight species Mitcheltree chemistry was assumed with CO2, CO, N2, O2, NO, C, N, and O. Both codes agreed within 1% on the forebody and had the expected differences on the aftbody. The NEQAIR and HARA codes were used to analyze the radiative heating on the vehicle using full spherical ray-tracing. The codes agreed within 5% on most aftbody points of interest.The LAURA code was then used to evaluate the conditions at angle of attack at the peak heating and peak pressure times. Boundary layer properties were investigated to confirm that the flow over the forebody was laminar for the flight.Comparisons of the aerothermal heating determined for the reconstructed trajectory to the design trajectories showed that the as-flown conditions were less severe than design

Beck, Robin↗

NASA NTRS2019

InSight's Reconstructed Aerothermal Environments

The InSight spacecraft was proposed to be a build-to-print copy of the Phoenix vehicle due to the knowledge that the lander payload would be similar and the trajectory would be similar. However, the InSight aerothermal analysts, based on tests performed in CO2 during the Mars Science Laboratory mission (MSL) and completion of Russian databases, considered radiative heat flux to the aftbody from the wake for the first time for a US Mars mission. The combined convective and radiative heat flux was used to determine if the as-flown Phoenix thermal protection system (TPS) design would be sufficient for InSight. All analyses showed that the design would be adequate. Once the InSight lander was successfully delivered to Mars on November 26, 2018, work began to reconstruct the atmosphere and trajectory in order to evaluate the aerothermal environments that were actually encountered by the spacecraft and to compare them to the design environments.The best estimated trajectory (BET) reconstructed for the InSight atmospheric entry fell between the two trajectories considered for the design, when looking at the velocity versus altitude values. The maximum heat rate design trajectory (MHR) flew at a higher velocity and the maximum heat load design trajectory (MHL) flew at a lower velocity than the BET. For TPS sizing, the MHL trajectory drove the design. Reconstruction has shown that the BET flew for a shorter time than either of the design environments, hence total heat load on the vehicle should have been less than used in design. Utilizing the BET, both DPLR and LAURA were first run to analyze the convective heating on the vehicle with no angle of attack. Both codes were run with axisymmetric, laminar flow in radiative equilibrium and vibrational non-equilibrium with a surface emissivity of 0.8. Eight species Mitcheltree chemistry was assumed with CO2, CO, N2, O2, NO, C, N, and O. Both codes agreed within 1% on the forebody and had the expected differences on the aftbody. The NEQAIR and HARA codes were used to analyze the radiative heating on the vehicle using full spherical ray-tracing. The codes agreed within 5% on most aftbody points of interest.The LAURA code was then used to evaluate the conditions at angle of attack at the peak heating and peak pressure times. Boundary layer properties were investigated to confirm that the flow over the forebody was laminar for the flight.Comparisons of the aerothermal heating determined for the reconstructed trajectory to the design trajectories showed that the as-flown conditions were less severe than design

radiation↗

NASA NTRSDate not supplied

Predictive Modeling of Carbon Ablators

Efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. To reduce the need for extensive testing, accelerate the design cycle process, and reduce uncertainty margins applied to final designs, NASA is developing simulation and modeling tools that enable characterization of material properties and response to high-enthalpy environments. The Porous Microstructure Analysis (PuMA) code has been developed for computing macroscale (volume averaged) properties of porous materials using microscale images from micro-computed tomography (micro-CT). Microscale modeling requires a realistic representation of a material microstructure; these are obtained either synthetically during the design of the material or through X-ray micro-CT. Volume averaged properties are then used to inform macroscale material response models, such as those implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO) software, also actively developed by NASA. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale assuming local thermal equilibrium. These tools were developed to efficiently interface with other pre-existing codes such as SPARTA (direct simulation Monte Carlo), DPLR (hypersonic CFD), NEQAIR (radiative transport) and DAKOTA (uncertainty quantification and optimization). Detailed flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI]) is critical for validating these computational tools for NASA applications. Examples of modeling ablative material response using these codes will be presented including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrate the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains, through the use of massively parallel computations.

Thermal Protection Systems↗

NASA NTRSDate not supplied

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performance for protecting probes during planetary entry. Data from the Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. A space-grade silicone-based coating commercially known as NuSil CV-1144-0 [3] was applied to the entire MSL heatshield, including the MEDLI plugs, to mitigate the spread of dust from PICA. Modeling the thermal response of PICA-NuSil (PICA-N) system is still an open challenge. Ground testing of PICA-N models exhibited surface temperature jumps of the order of 150 K due to oxide scale formation and sub-sequent NuSil burn-off. It is therefore critical to include a validated model for the material response of the coating in engineering codes. A test campaign has been conducted at the NASA’s Langley HyMETS [4] facility to screen the response of PICA-N and gather detailed data on its behavior [5]. A first model of PICA-N thermal response has been developed using the Hy-METS experiments [6]. The objective of this work is to analyze the material response of the latest PICA-N model compared to the engineering model used to simulate the entry of MSL. The environment and material response around the MSL aeroshell during Mars atmospheric entry is simulated using a collection of tools. The Direct Simulation Monte Carlo SPARTA code [7] is used in the rarefied regime, the Data Parallel Line Relaxation (DPLR) code [8] is used in the continuum regime and radiative heating conditions are provided by the Nonequilibrium air radiation (NEQAIR) code [9] to estimate the environmental conditions. The thermal response inside the material is computed using the Porous material Analysis Toolbox based on Open-FOAM (PATO) [10,11,12]. Thermodynamic and chemistry properties are estimated using the Mutation++ library [13]. The approach implemented in PATO as a first cut PICA-N thermal response model is outlined in Figure 1. While the recession is less than the coating thickness, the Surface mass and energy balance Boundary Condition (SBC) uses the NuSil B’ tables. Once the recession removes the coating, the usual PICA B’ tables are used for the SBC. The B’ tables are computed using an equilibrium solver implemented in Mutation++, given the temperature, pressure, blowing rate, composition of the pyrolysis and environment gases, and the condensed species at the surface. Preliminary results of the 3D material response of the MSL heat-shield at the peak heating (80 sec after Entry Interface) are shown in Figure 2. Current NASA’s mission to Mars, Mars 2020, used the spare heatshield of MSL for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the PATO high-fidelity material response capability was benchmarked against flight data from MEDLI. This effort represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes. This bench-marking is awaiting the final release of the MEDLI-2 data.

Aerospace↗

NASA NTRSDate not supplied

Predictive Modeling of Carbon Ablators Using Micro and Macro-Scale Modeling

Efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. To reduce the need for extensive testing, accelerate the design cycle process, and reduce uncertainty margins applied to final designs, NASA is developing simulation and modeling tools that enable characterization of material properties and response to high-enthalpy environments. The Porous Microstructure Analysis (PuMA) code has been developed for computing macroscale (volume averaged) properties of porous materials using microscale images from micro-computed tomography (micro-CT). Microscale modeling requires a realistic representation of a material microstructure; these are obtained either synthetically during the design of the material or through X-ray micro-CT. Volume averaged properties are then used to inform macroscale material response models, such as those implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO) software, also actively developed by NASA. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale assuming local thermal equilibrium. These tools were developed to efficiently interface with other pre-existing codes such as SPARTA (direct simulation Monte Carlo), DPLR (hypersonic CFD), NEQAIR (radiative transport) and DAKOTA (uncertainty quantification and optimization). Detailed flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI]) is critical for validating these computational tools for NASA applications. Examples of modeling ablative material response using these codes will be presented including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrate the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains, through the use of massively parallel computations.

Thermal Protection Systems↗

NASA NTRSDate not supplied

Analysis of Nonequilibrium Molecular Nitrogen Ultraviolet Radiation in Pure N2 Shockwaves

This paper presents an analysis of data previously reported from a test series in the Electric Arc Shock Tube Facility for incident shocks composed of pure Nitrogen. The present work focuses on the lowest velocity measured in the test (approximately 7 km/s) where molecular radiation is significant. Analysis of the spectral data with the NEQAIR radiative transport code obtains post-shock trends in temperature and number densities of the states of N2. The results suggest that chemical, and possibly thermal, equilibrium has not been obtained within the distance/time measured (about 5 cm/120 s). In the nonequilibrium region of the shock, differences between rotational temperatures of different species are observed, as well as differences between rotational and vibrational temperatures. Population of rotational states above the dissociation limit is observed to occur over a distance of several cm. Finally, these data suggest that predissociation strongly affects the C^3 Π_u state of N2, leading to a lower measured radiation than predicted. An effective average pre-dissociation rate is shown to produce reasonable agreement under quasi-steady state approximation.

Augustin Tibere Inglesse↗

NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

NASA NTRSDate not supplied

Analysis of nonequilibrium molecular nitrogen ultraviolet radiation in pure N2 shockwaves

An analysis of data previously reported from a test series in the Electric Arc Shock Tube Facility for incident shocks composed of pure Nitrogen was performed. The present work focuses on the lowest velocity measured in the test (approximately 7 km/s) where molecular radiation is significant. Analysis of the spectral data with the NEQAIR radiative transport code obtains post-shock trends in temperature and number densities of the states of N2. The results suggest that chemical, and possibly thermal, equilibrium has not been obtained within the distance/time measured (about 5 cm/120 ms). In the nonequilibrium region of the shock, differences between rotational temperatures of different species are observed, as well as differences between rotational and vibrational temperatures. Population of rotational states above the dissociation limit is observed to occur over a distance of several cm. Finally, these data suggest that predissociation strongly affects the C^3 Π_u state of N2, leading to a lower measured radiation than predicted. An effective average pre-dissociation rate is shown to produce reasonable agreement under quasi-steady state approximation.

Augustin Tibere-Inglesse↗