Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Enthalpy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Measurement of Material Recession and Shock Standoff in Plasma Windtunnel using Neural Nets

Arcjets are plasma wind tunnels used to test the performance of heatshield materials for spacecraft atmospheric entry. These facilities present an extremely harsh flow environment with heat fluxes up to 109 W/m2 for up to 30 minutes. The plasma is low-temperature (∼1 eV) but high pressure (> 10 kPa) creating high-enthalpy supersonic flows similar to atmospheric entry conditions. Typically, material samples are measured before and after a test to characterize the total recession. However, this does not capture time-dependent effects such as material expansion and non-linear recession. This work will present new analysis of arcjet test videos which measure both the time-dependent 2D recession of the material samples and the shock standoff distance. The results show non-linear time-dependent effects are present for some conditions. The material and shock edges are extracted from the videos by training and applying a convolutional neural network. Due to the consistent camera settings, the machine learning model achieves high accuracy (± 2 px) relative to manually segmented images with only a small number of training frames (80).

Neural network↗
NASA NTRSDate not supplied

CFD Simulations of the Truncated IHF Semi-Elliptical Nozzle Flow: Convective Calibration Plate

This paper reports computational simulations of a high-enthalpy panel test facility at NASA Ames Research Center. The recently developed facility configuration uses a truncated version of an existing semi-elliptical nozzle in the 60-MW Interaction Heating Facility, and it provides testing capability for panel test articles as large as 43.2 cm x 43.2 cm (17 in x 17 in). A water-cooled calibration plate, which is attached to the semi-elliptical nozzle bottom surface and instrumented with twenty Gardon-gage heat flux calorimeters and twelve pressure gages, is used to characterize the convective heating environment over the panel test articles. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, and flowfield over panel test articles. Both laminar and turbulent simulations are performed, and the computational results are compared with the calibration plate measurements. Comparisons of computed and measured heat fluxes indicate that the boundary layer in the semi-elliptical nozzle is transitional at certain conditions (relatively higher arc-heater pressures or mass flow rates).

Arc-jets↗
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

Computing the effective elasticity of anisotropic porous media from X-ray computed micro-tomography images

Development and optimization of composite materials designed for thermal protection of NASA’s spacecraft requires the understanding of their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the micro-scale on realistic representations of these composites. The digital micro-structures are generated either synthetically or through X-ray micro-tomography reconstructions. One of the main challenges in the prediction of heatshield material structural response is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the micro-scale. These are driven by the fiber layout within the micro-structure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled through the use of a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme. The MPSA, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA), was discussed in a previous presentation at the 15th USCCN and in reference. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Federico Semeraro↗
NASA NTRSDate not supplied

Numerical Investigation of Film Coefficient Engineering Methodology for Dissociated, Chemically Reacting Boundary Layers

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flow field for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for their application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. The sensitivity of the engineering model is assessed based on the impact of the CFD assumptions as an effort to make it available to design engineers across industry.

numerical boundary layer↗
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

CFD Simulations of the Truncated IHF Semi-Elliptical Nozzle Flow: Convective Calibration Plate

This paper reports computational simulations of a high-enthalpy panel test facility at NASA Ames Research Center. The recently developed facility configuration uses a truncated version of an existing semi-elliptical nozzle in the 60-MW Interaction Heating Facility, and it provides testing capability for panel test articles as large as 43.2 cm x 43.2 cm (17 in x 17 in). A water-cooled calibration plate, which is attached to the semi-elliptical nozzle bottom surface and instrumented with twenty Gardon-gage heat flux calorimeters and twelve pressure gages, is used to characterize the convective heating environment over the panel test articles. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, and the flowfields over panel test articles. Both laminar and turbulent simulations are performed, and the computational results are compared with the calibration plate measurements. Comparisons of computed and measured heat fluxes indicate that the boundary layer in the semi-elliptical nozzle is transitional at certain conditions (relatively higher arc-heater pressures or mass flow rates).

Arc-jets↗
NASA NTRSDate not supplied

Using Thermally Imperfect and Ideal Gas Equations for Determining Conditions in Slotted Throat Wind Tunnels

The derivation and application of two back pressure controlling methods for setting the dynamic pressure, Mach number, and Reynolds number of subsonic wind tunnels with a slotted-throat test section are discussed. The first method uses a pressure difference to determine the tunnel dynamic pressure via the compressible Bernoulli equations and is used in the ideal gas path in the CFD code. The second method uses an enthalpy difference in the flow to determine the velocity in the test section and can be applied to multispecies, calorically imperfect and calorically perfect gas flow simulations. For a tunnel operating point close to standard atmosphere conditions, Mach number, pressure, and temperature along the tunnel centerline were very closely matched between the various thermodynamic equation assumption controller methods.

Ground Testing↗
NASA NTRSDate not supplied

Validating Arc-jet Condition Targets Using Inverse 1D Material Response Models

Experiments were conducted using both the Interaction Heating Facility (IHF) and Aerodynamic Heating Facility (AHF), two hypersonic arc-jet facilities at NASA Ames. These tests were run to understand the material response of ablative thermal protection system materials by exposing them to high-enthalpy flows. Using thermocouple data from the materials, test conditions were inversely determined and compared with expected results from a one-dimensional model using facility target conditions; this allowed an approximate environmental scale factor to be determined and applied to more accurately predict material response in additional tests.

material response↗
NASA NTRSDate not supplied

On Characterization of Laminar to Turbulent Transitional Flow in the Truncated IHF Semi-Elliptical Nozzle

This paper reports computational analysis and experimental data to characterize the laminar to turbulent flow transition in a high-enthalpy panel test facility at NASA Ames Research Center. The recently developed facility configuration uses a truncated version of an existing semi-elliptical nozzle in the 60-MW Interaction Heating Facility, and it provides testing capability for panel test articles as large as 43.2 cm x 43.2 cm (17 in x 17 in). A water-cooled calibration plate, which is attached to the semi-elliptical nozzle bottom surface and instrumented with twenty Gardon-gage heat flux calorimeters and twelve pressure gages, is used to characterize the convective heating environment over the panel test articles. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, and the flowfields over panel test articles, and comparisons with the measured test data. Both laminar and turbulent simulations are performed. Constant current studies are performed: the measured data are grouped at the same arc current levels with varying mass flow rates, and they are compared with corresponding computations. For all currents studied, these studies indicate that the semi-elliptical nozzle flow is transitional at relatively higher mass flow rates within the facility operating envelope.

Arc-jets↗
NASA NTRSDate not supplied

Thermal Protection Materials and Systems at NASA Ames Research Center

Thermal Protection Systems (TPS) are critical for enabling NASA missions involving high-speed atmospheric flight where the entries usually include descending into the atmosphere followed by a trajectory that aims to burn off energy and result in a controlled landing. NASA Ames focuses on qualifying and certifying TPS for current missions, sustaining TPS for future missions, and developing new TPS for upcoming missions where a heritage solution is not viable. More recently there is also a focus on advancing and transferring technologies that can benefit both commercial and government space needs. Developing mature thermal protection systems is a lengthy process involving advanced tools, extensive research, and testing. Design and analysis tools are used to predict aerothermal environments, aid the design of test and flight hardware, and support the testing for the thermal/mechanical response of thermal protection systems. More recently, advances in computational methods help reduce the time and cost of technological advances, aid in optimized material architecture design, and improve material properties and performance. While high-enthalpy testing that simulates the conditions of space flight remains essential for the evaluation and development of TPS materials, computational tools are already showing promise in reducing the need for widespread testing and can help fast-track the design cycle. With the exploration of new destinations EDL instrumentation remains an important element of the heatshield and NASA Ames and partners have developed and delivered instrumentation flight hardware in support of recent Mars missions (MSL and M2020) as well as Artemis Orion. Sensors installed on the heatshield and backshell of spacecraft provide coveted information about the aerodynamic and aerothermal environment during entry. Over the years NASA Ames has brought several reusable and ablative TPS materials to a level of readiness to hand off to missions and the branch continues to serve as a TPS steward for the agency. This presentation will cover current developments in the above areas that have enabled recent missions, and look to future TPS needs for missions such as Ice Giants.

thermal protection materials↗
NASA NTRSDate not supplied

On Characterization of Laminar to Turbulent Transitional Flow in the Truncated IHF Semi-Elliptical Nozzle

This paper reports computational analysis and experimental data to characterize the laminar to turbulent flow transition in a high-enthalpy panel test facility at NASA Ames Research Center. The recently developed facility configuration uses a truncated version of an existing semi-elliptical nozzle in the 60-MW Interaction Heating Facility, and it provides testing capability for panel test articles as large as 43.2 cm x 43.2 cm (17 in x 17 in). A water-cooled calibration plate, which is attached to the semi-elliptical nozzle bottom surface and instrumented with twenty Gardon-gage heat flux calorimeters and twelve pressure gages, is used to characterize the convective heating environment over the panel test articles. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, and the flowfields over panel test articles, and comparisons with the measured test data. Both laminar and turbulent simulations are performed. Constant current studies are conducted: the measured data are grouped at the same arc current levels with varying mass flow rates, and they are compared with corresponding computations. For all currents considered, these studies indicate that the semi-elliptical nozzle flow is transitional at relatively higher mass flow rates within the facility operating envelope. An approximate facility map for the transitional flow characterization is developed, showing laminar and transitional to turbulent flow regions.

Arc-jets↗
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

Thermal Protection Materials and Systems at NASA Ames Research Center

Thermal Protection Systems (TPS) are critical for enabling NASA missions involving high-speed atmospheric flight where the entries usually include descending into the atmosphere followed by a trajectory that aims to burn off energy and result in a controlled landing. NASA Ames focuses on qualifying and certifying TPS for current missions, sustaining TPS for future missions, and developing new TPS for upcoming missions where a heritage solution is not viable. More recently there is also a focus on advancing and transferring technologies that can benefit both commercial and government space needs.Developing mature thermal protection systems is a lengthy process involving advanced tools, extensive research, and testing. Design and analysis tools are used to predict aerothermal environments, aid the design of test and flight hardware, and support the testing for the thermal/mechanical response of thermal protection systems. More recently, advances in computational methods help reduce the time and cost of technological advances, aid in optimized material architecture design, and improve material properties and performance. While high-enthalpy testing that simulates the conditions of space flight remains essential for the evaluation and development of TPS materials, computational tools are already showing promise in reducing the need for widespread testing and can help fast-trackthe design cycle.With the exploration of new destinations EDL instrumentation remains an important element of the heatshield and NASA Ames and partners have developed and delivered instrumentation flight hardware in support of recent Mars missions,including Mars Science Lab (MSL) and Mars 2020 (M2020),as well as Artemis Orion. Sensors installed on the heatshield and backshell of spacecraft provide coveted information about the aerodynamic and aerothermal environment during entry.Overthe years NASA Ames has brought several reusable and ablative TPS materials to a level of readiness to hand off to missions and the Thermal Protection Materials branch continues to serve as a TPS steward for the agency.

Thermal protection materials↗
NASA NTRS2023

Numerical Investigation of Film Coefficient Approximation for Chemically Reacting Boundary-Layer Flows

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flowfield for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for his/her application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. Within the subset of tested energy boundary conditions, it is found that the most desirable energy boundary condition is the radiative equilibrium boundary condition, which permits conservative estimates of convective heat flux, but also generates flowfield-dependent spatial thermal distributions along the surface. Other key findings are presented in an effort to make the film coefficient engineering model readily available to design engineers across industry.

Film coefficient↗
NASA NTRSDate not supplied

High-precision predictions of properties of chemically disordered crystals

Multiple scattering theory (MST) combined with density functional theory (DFT) allows to predict properties of chemically disordered materials from the first principles. However, such predictions often suffer from the systematic errors, which depend on crystal geometry. Each computed property of a particular crystal structure typically has a relatively small random error and a larger systematic error. Cancellation of systematic errors allows more accurate predictions. We propose a computational methodology based on the subtraction of the systematic errors in MST. To exemplify it, we apply it to the precipitated alloys. Considering precipitation strengthening in Ni superalloys, we compute the relative enthalpies of the competing Ni3(Al,Ti)1 crystal structures with a chemical disorder on the Al+Ti sublattice. Such predicted composition-structure-property dependencies are useful for the guided design of the next-generation alloys with improved strength. Our predictions are validated by comparison with the results of other DFT methods (having a higher computational cost) and with experiment. We acknowledge funding of the guided design of stronger superalloys for airspace by NASA Aeronautics Research Mission Directorate (ARMD) via Transformational Tools and Technologies (TTT) Project.

Computational↗
NASA NTRS2022

High-precision predictions of properties of chemically disordered crystals

Multiple scattering theory (MST) combined with density functional theory (DFT) allows to predict properties of chemically disordered materials from the first principles. However, such predictions often suffer from the systematic errors, which depend on crystal geometry. Each computed property of a particular crystal structure typically has a relatively small random error and a larger systematic error. Cancellation of systematic errors allows more accurate predictions. We propose a computational methodology based on the subtraction of the systematic errors in MST. To exemplify it, we apply it to the precipitated alloys. Considering precipitation strengthening in Ni superalloys, we compute the relative enthalpies of the competing Ni_3(Al_{1-x}Ti_x)_1 crystal structures with a chemical disorder on the Al+Ti sublattice. Such predicted composition-structure-property dependencies are useful for the guided design of the next-generation alloys with improved strength. Our predictions are validated by comparison with the results of other DFT methods (having a higher computational cost) and with experiment.

density functional theory↗
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) [1,2,3]. 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) [4] 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) [5] 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 [6]. 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 [7]. 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.

Thermal Protection Systems↗