SEARCH · Search NASA

Results for “Material Response Modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Thermomechanical deformation testing and modeling in the presence of metallurgical instabilities

A number of viscoplastic constitutive models were developed to describe deformation behavior under complex combinations of thermal and mechanical loading. Questions remain, however, regarding the validity of procedures used to characterize these models for specific structural alloys. One area of concern is that the majority of experimental data available for this purpose are determined under isothermal conditions. This experimental study is aimed at determining whether viscoplastic constitutive theories characterized using an isothermal data base can adequately model material response under the complex thermomechanical loading conditions typical of power generation service. The approach adopted was to conduct a series of carefully controlled thermomechanical experiments on a nickel-based superalloy, Hastelloy Alloy X. Previous investigations had shown that this material experiences metallurgical instabilities leading to complex hardening behavior, termed dynamic strain aging. Investigating this phenomenon under full thermomechanical conditions leads to a number of challenging experimental difficulties which up to the present work were unresolved. To correct this situation, a number of advances were made in thermomechanical testing techniques. Advanced methods for dynamic temperature gradient control, phasing control and thermal strain compensation were developed and incorporated into real time test control software. These advances allowed the thermomechanical data to be analyzed with minimal experimental uncertainty. The thermomechanical results were evaluated on both a phenomenological and microstructural basis. Phenomenological results revealed that the thermomechanical hardening trends were not bounded by those displayed under isothermal conditions. For the case of Hastelloy Alloy X (and similar dynamic strain aging materials), this strongly suggests that some form of thermomechanical testing is necessary when characterizing a thermoviscoplastic deformation model. Transmission electron microscopy was used to study the microstructural physics, and analyze the unique phenomenological behavior.

Castelli, Michael G.↗

NASA NTRSDate not supplied

Simulating Ordinary Chondrite (Tamdakht) Ablation at the Hypersonic Materials Environmental Test System (HyMETS) Facility.

NASA’s Science Mission Directorate (SMD) created the Asteroid Threat Assessment Project (ATAP) to inform decision-makers of risks associated with Potentially Hazardous Objects (PHOs), which may pose an existential threat to human civilization. ATAP assesses risk, in part, by developing analytical physics-based damage models which draw upon data collected during material property characterization efforts, entry simulations, hazard simulations, and ground-based arc-jet testing. Therefore, a pathfinder test campaign was con-ducted at the Hypersonic Materials Environmental Test System (HyMETS) at the NASA Langley Re-search Center to investigate the ablation mechanisms of an ordinary chondrite (Tamdakht H5) and a terrestrial analog (basalt). The HyMETS facility is a 400 kW constricted arc heater used to screen material performance under simulated aerothermal conditions of hypersonic flow.1 Facility conditions were chosen to simulate Earth entry conditions be-tween the upper-mesosphere and the lower thermo-sphere. The ablation mechanisms of Tamdakht and bas-alt are illustrated with still images collected from high-speed video cameras (fig. 1, Tamdakht). Tamdakht produces a relatively stable melt flow, as evidenced by the perpetual growth of a flange at the sidewall and the exiguous detachment of molten material into the flow. Furthermore, an inspection of high-speed video and chemical analysis of the post-test melt layer suggests that the dominant mode of mass loss for Tamdakht, at these test conditions, is expressed through the vaporization of volatiles (e.g., iron, potassium, sodium, phosphorus, and sulfur). The vaporization rates of Tamdakht are greater than basalt under the most extreme test conditions, which leads to an enhanced blowing layer, in-creased thermal shielding, and reduced recession. The ablation mechanisms of basalt are markedly different from Tamdakht and are attributed to the presence of hydrated minerals dispersed at irregular intervals throughout the silicate matrix. Rapid de-composition of secondary minerals and the subsequent formation of water lead to the ejection of the subsurface and overlying melt layer. Furthermore, the presence of water is suspected of lowering the viscosity of the melt layer resulting in increased mass loss near the edge of the test article where drag forces overcome rheological properties (fig. 1, bas-alt). Finally, the results will be discussed in the broader context of meteorite-based material response model development and the potential impact on the Asteroid Threat Assessment Project (ATAP).

Brody K. Bessire↗

NASA NTRS2018

Overview of Heatshield for Extreme Entry Environment Technology (HEEET)

The Heatshield for Extreme Entry Environment Technology (HEEET) projects objective is to mature a 3-D Woven Thermal Protection System (TPS) to Technical Readiness Level (TRL) 6 to support future NASA missions to destinations such as Venus and Saturn. The scope of the project, status of which will be discussed, encompasses development of manufacturing and integration processes, fabrication of a prototype 1m diameter engineering test unit (ETU) that will undergo a series of structural tests, characterizing material aerothermal performance including development of a material response model, and structural testing and analysis to develop tools to support design and establish system capability.

Driver, David M.↗

NASA NTRS2014

AVCOAT Density Characterization for Orion Multi-Purpose Crew Vehicle

The Orion Multi-Purpose Crew Vehicle (MPCV) will transport four crew members to and from lunar-class orbital destinations. The first orbital Exploration Flight Test (EFT-1) is scheduled for December 2014 and will provide valuable data on several systems, including the heat shield. The heat shield material is AvcoatTM, a mid-density ablator. Specifically, the heat shield contains a fiberglass-phenolic honeycomb structure filled with an ablative epoxy novolac resin. Post-flight characterization of the EFT-1 thermal protection system will be conducted in order to study the material response and measure the char front of the material. Avcoat density profiling will be conducted to study three critical material zones: char, pyrolysis, and virgin states as a function of material depth. As part of a ground test campaign, Avcoat coupons are tested at the NASA Ames Research Center Arc Jet Complex. The test campaign is set forth to study material response based on environments, perform Avcoat material density characterization, and compare the char depths from the HEAT sensor and density profiling studies. These investigations will then be compared to flight data in order to improve upon TPS material response models.

AVCOAT↗

NASA NTRS2019

Laboratory-Based Thermal Shock Investigation of Heat Flux Sensors for the Mars 2020 Backshell

Miller, R. A.↗

NASA NTRS2019

Laboratory-Based Thermal Shock Investigation of Heat Flux Sensors for the Mars 2020 Backshell

In 2012 during the entry, descent, and landing of the Mars Science Laboratory (MSL), the MSL Entry, Descent, and Landing Instrumentation (MEDLI) sensor suite was collecting in-flight heatshield pressure and temperature data. The data collected by the MEDLI instruments has since been used for reconstruction of vehicle aerodynamics, atmospheric conditions, aerothermal heating, and Thermal Protection System (TPS) performance as well as material response model validation and refinement. The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite for the Mars 2020 heatshield and backshell is being designed to expand on the measurements and knowledge gained from MEDLI. Similar to MEDLI, MEDLI2 will measure the pressure and temperature of the heatshield. MEDLI2 will additionally measure the temperature, pressure, total heat flux, and radiative heat flux on the backshell.Since the backshell instrumentation is new to MEDLI2, Do No Harm (DNH) testing was conducted on instrumented backshell TPS (SLA-561V) panels. The panels consisted of four pressure port holes, one Mars Entry Atmospheric Data System (MEADS) pressure port plug, one MEDLI2 Integrated Sensor Plug (MISP) thermal plug, and one heat flux sensor. DNH testing was conducted to ensure the performance of the TPS was not degraded due to sensor integration and to characterize any TPS performance changes. The testing consisted of environmental testingâ€” vibration, shock, thermal vacuum (TVAC) cyclingâ€” and bounding aerothermal (arc jet) testing. During arc jet testing, the heat flux sensors embedded in the SLA-561V panels exhibited an unexpected temporary reduction in the heat flux sensor temperature and response. After review of the test results, it was determined that this unexpected response was confined to the two heat flux sensors that experienced the greatest thermal shock condition. This condition consisted of a liquid nitrogen (LN2) bath that induced temperatures of approximately -190Â°C, and then a transition (thermal shock) to an arc jet test at a heat rate of approximately 21 W/cm2. Both heat flux sensors that were exposed to this thermal shock experienced a blister in the thermal coating during the arc jet test.Two heat flux sensor thermal shock test series were performed to investigate the cause of the blistering and subsequent energy release. In these tests, the heat flux sensor was first cold soaked in either a dry ice or LN2 bath to induce temperatures of approximately -78Â°C or -190Â°C, respectively. Then the sensors were thermally shocked using two propane torches with a heat rate of either approximately 8 W/cm2 or 21 W/cm2. The key findings indicated that there is a correlation between thermal shock and the blistering observed in the DNH test series, and that the cause appeared to be rooted in the heat flux sensor epoxy that encapsulates the sensor thermopile.Since the heat flux sensors are required to measure heat fluxes up to 15 W/cm2 during the Mars 2020 entry, a third test series was designed to determine if blistering is an issue at this maximum expected flight heat flux. Results from all three thermal shock test series and a discussion about whether or not blistering of the heat flux sensor thermal coating could be an issue for the Mars 2020 mission will be presented.

Miller, R. A.↗

NASA NTRSDate not supplied

MEDLI2: MISP Inferred Aerothermal Environment and Flow Transition Assessment

The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite on the Mars2020 mission contained multiple sensors on the aeroshell to measure the aerothermal environment during entry into the Martian atmosphere. These sensors performed superbly and successfully returned forebody and aftbody heating measurements. Analysis of MEDLI2 data indicated flow transitioning from a laminar to turbulent state on the heatshield. No evidence of flow transition was observed on the backshell. Two methods were used to estimate flow transition times on the heatshield: (1) temperature gradient of near-surface thermocouple data and (2) heat flux gradient from an inverse reconstruction approach using thermocouple data and material response modeling. Both methods produced similar transition times with an estimated accuracy of ±1 s. To assess various transition criteria, transition parameters were evaluated at each sensor location using flow field solutions from computational fluid dynamics (CFD) simulations. The idea was to use conservative values inferred from MEDLI2 data as transition criteria for other Mars missions. To test this hypothesis, MEDLI data from the Mars Science Laboratory (MSL) mission was used to compare predicted vs. actual flow transition times. The comparisons suggest smooth wall transition criteria are not well-suited in modeling the rapid progression of a turbulent transition front. Transition criteria containing a roughness element parameter agreed better with the flight data. In summary, critical transition values derived from MEDLI2 data may be used as a starting point in constructing a flow transition model for future Mars missions.

Chun Y Tang↗

NASA NTRSDate not supplied

Simulating Meteor Ablation at the Hypersonic Materials Environmental Test System

The ablation mechanisms of stony meteorite material and a terrestrial analog were investigated at the Hypersonic Materials Environmental Test System (HyMETS) facility at NASA’s Langley Research Center. The primary analyte (Tamdakht-H5) and arc-jet test conditions were chosen to simulate the entry of an ordinary chondrite through the Earth’s mesosphere. Basalt was selected to compare the meteorite response against material with enhanced polymerization and an elevated concentration of volatile species. The initial ablation process for each material commences with the vaporization of volatile species and germination of a viscous melt layer on the stagnation surface. However, each material exhibited distinct mass loss mechanisms, with basalt recessing at twice the rate of Tamdakht. Inspection of high-speed video and chemical analysis of the post-test melt layer suggests the dominant mode of mass loss of Tamdakht is expressed through the vaporization of volatile compounds. Vaporization rates of Tamdakht are 18x higher than basalt under the most severe test conditions which provides enhanced thermal shielding. Furthermore, Tamdakht produces a relatively stable melt flow, as evidenced by perpetual growth at the sidewall and the exiguous detachment of molten material into the flow. The ablation mechanism of basalt is markedly different from Tamdakht and is governed by volatile species in the virgin material. The dominant mode of mass loss for basalt is attributed to the detachment of large sections of the stagnation surface near the edge radius. An additional but subordinate mode of mass loss is assigned to the removal of melt flow as it passes over the edge radius, where shear forces overcome the rheological properties of the molten matter. Finally, mechanistic details will be discussed in terms of designing future test campaigns, developing meteorite-based materials response models, and impacting the broader scope of the Asteroid Threat Assessment Project (ATAP).

Meteor Ablation↗

NASA NTRSDate not supplied

In Situ Synchrotron Micro-CT of Microballoon-Filled RTV Silicone During Thermal Decomposition

This work investigates the thermal decomposition of microballoon-filled room-temperature-vulcanizing (RTV) silicone using time-resolved synchrotron micro-CT at the Advanced Light Source. In situ imaging of constrained and unconstrained samples captures the evolution of the internal microstructure during heating, including material expansion, microballoon and void growth, tearing, pyrolysis, shrinkage, and increasing porosity. The measurements provide new insight into the mechanisms governing RTV degradation and generate three-dimensional data for future quantitative analysis and material-response model development.

TPS↗

DOE OSTI2024

Peridynamic neural operators: A data-driven nonlocal constitutive model for complex material responses

Neural operators, which can act as implicit solution operators of hidden governing equations, have recently become popular tools for learning the responses of complex real-world physical systems. Nevertheless, most neural operator applications have thus far been data-driven and neglect the intrinsic preservation of fundamental physical laws in data. Here, in this work, we introduce a novel integral neural operator architecture called the Peridynamic Neural Operator (PNO) that learns a nonlocal constitutive law from data. This neural operator provides a forward model in the form of state-based peridynamics, with objectivity and momentum balance laws automatically guaranteed. As applications, we demonstrate the expressivity and efficacy of our model in learning complex material behaviors from both synthetic and experimental data sets. We also compare the performances with baseline models that use predefined constitutive laws. We show that, owing to its ability to capture complex responses, our learned neural operator achieves improved accuracy and efficiency. Moreover, by preserving the essential physical laws within the neural network architecture, the PNO is robust in treating noisy data. The method shows generalizability to different domain configurations, external loadings, and discretizations.

42 ENGINEERING↗

NASA NTRSDate not supplied

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗

NASA NTRS2015

Trajectory Simulation of Meteors Assuming Mass Loss and Fragmentation

Program used to simulate atmospheric flight trajectories of entry capsules [1] Includes models of atmospheres of different planetary destinations - Earth, Mars, Venus, Jupiter, Saturn, Uranus, Titan, ... Solves 3-‐degrees of freedom (3DoF) equations for a single body treated as a point mass. Also supports 6-DoF trajectory simula4on and Monte Carlo analyses. Uses Fehlberg-‐Runge-‐Kuna (4th-5th order) time integraion with automaic step size control. Includes rotating spheroidal planet with gravitational field having a J2 harmonic. Includes a variety of engineering aerodynamic and heat flux models. Capable of specifying events - heatshield jettison, parachute deployment, etc. - at predefined altitudes or Mach number. Has material thermal response models of typical aerospace materials integrated.

Atmospheric Entry↗

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↗

DOE OSTI2022

Sensitivity Study of Multiscale and Phenomenological Elasto-Viscoplastic Grade 91 Material Models for Component-Scale Response

Many advanced nuclear reactor concepts currently being developed are targeting higher operating temperatures relative to the current fleet of light water nuclear reactors, for efficiency gains and other operational considerations. The design of high temperature structural components with reliable long-term operational performance will depend on material models that accurately capture the inelastic deformation mechanisms active in these environments. In this work, we perform a detailed parameter sensitivity analysis of two unified elasto-viscoplastic Grade 91 material models capable of capturing long term high temperature creep deformation. The first model is a phenomelogical material model from the Nuclear Engineering Material Library (NEML) developed at Argonne National Lab. The NEML model parameters and their uncertainty were fit to a range of Grade 91 experimental data using Bayesian Markov Chain Monte Carlo analysis. The second model is a LAROMance data-driven surrogate material model developed at Los Alamos National Lab. The LAROMance model is fit to a large database of responses produced by a mechanistic crystal plasticity based polycrystal model. Parameters for the LAROMance surrogate material model reflect the pedigree of the Grade 91 microstructure. Both material models have been integrated into the Grizzly code, based on the open-source MOOSE multiphysics simulation framework, to simulate both the progression of aging mechanisms and the effects of that aging on nuclear power plant structures. Grizzly is used analyze a three-dimensional Grade 91 piping system to compare the long-term inelastic response predicted by these two fundamentally different models and assess the sensitivity of the material model input parameters on this quantity of interest.

42 ENGINEERING↗

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

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↗

NASA NTRS1984

Constitutive modeling of inelastic anisotropic material response

A constitutive equation was developed to predict the inelastic thermomechanical response of single crystal turbine blades. These equations are essential for developing accurate finite element models of hot section components and contribute significantly to the understanding and prediction of crack initiation and propagation. The method used was limited to unified state variable constitutive equations. Two approaches to developing an anisotropic constitutive equation were reviewed. One approach was to apply the Stouffer-Bodner representation for deformation induced anisotropy to materials with an initial anisotropy such as single crystals. The second approach was to determine the global inelastic strain rate from the contribution of the slip in each of the possible crystallographic slip systems. A three dimensional finite element is being developed with a variable constitutive equation link that can be used for constitutive equation development and to predict the response of an experiment using the actual specimen geometry and loading conditions.

Stouffer, D. C.↗

NASA NTRSDate not supplied

Three Dimensionally Woven Mid-Density Carbon Phenolic (3MDCP) is a Novel Single Piece Ablative TPS for Extreme Entry Environments

Three Dimensionally Woven Mid-Density Carbon Phenolic (3MDCP) is robust, single piece, carbon phenolic ablative thermal protection system under development at NASA Ames Research Center initially for the Mars Sample Return Earth Entry System (MSR EES). The MSR EES requirements drove the need for TPS with no seams, capable of surviving the highest entry conditions for any NASA Earth return capsule with heat fluxes >2000 W/cm2 and pressures >1.5 atmospheres. To produce 3MDCP required development of new weaving infrastructure to enable weaving of preforms large enough to form into a single piece heatshield. It required development of forming techniques to transform a flat woven panel into a sphere cone shape and enhanced infusion processes to support larger scale infusion of resin into the formed preforms. A rigorous performance testing campaign was conducted to develop and validate the materials thermal response model used to determine the required material thickness and to demonstrate the material can survive the extreme entry conditions. The end of the development effort (end of FY26) will result in a TPS at Technical Readiness Level (TRL) 6 and Manufacturing Readiness Level (MRL) 6+ for the MSR EES mission and a mature system ready to support other missions. This poster will provide a snapshot of where 3MDCP is in its development phase.

Ablator↗