Search NASA⌕ Search

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.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Improving dynamic collision frequencies: Impacts on dynamic structure factors and stopping powers in warm dense matter

Simulations and diagnostics of high-energy-density plasmas and warm dense matter rely on models of material response properties, both static and dynamic (frequency-dependent). Here, in this work, we systematically investigate variations in dynamic electron–ion collision frequencies ν(ω) in warm dense matter using data from a self-consistent-field average-atom model. We show that including the full quantum density of states, strong collisions, and inelastic collisions lead to significant changes in ν(ω). These changes result in red shifts and broadening of the plasmon peak in the dynamic structure factor, an effect observable in x-ray Thomson scattering spectra, and modify stopping powers around the Bragg peak. These changes improve the agreement of computationally efficient average-atom models with first-principles time-dependent density functional theory in warm dense aluminum, carbon, and deuterium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Analysis of MSL/MEDLI Entry Data with Coupled CFD and Material Response

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performances for protecting probes during planetary entry. 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]. The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on OpenFOAM (PATO) [4,5,6] using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA [7], is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Thermal Protection Systems↗
NASA NTRSDate not supplied

ANALYSIS OF THE MSL/MEDLI ENTRY DATA WITH COUPLED CFD AND MATERIAL RESPONSE.

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heat-shield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material. PICA is a lightweight carbon fiber/polymeric resin material that offers out-standing performances for protecting probes during planetary entry. 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). The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on Open-FOAM (PATO) using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA, is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Mars Science Laboratory↗
NASA NTRSDate not supplied

Overview of Ablative TPS Modeling at NASA Ames

Over the past decade, NASA has invested in efforts to build predictive thermal protection system (TPS) material models from the micro-scale to the macro-scale. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. Traditional material response and ablation modeling tools, such as the heritage code FIAT, and its multidimensional siblings, TITAN and 3dFIAT, are being complemented with newly developed software such as Icarus and PATO. Both of these programs are three-dimensional, finite-volume solvers that use unstructured meshes and 21st century programming paradigms to allow for efficient parallel simulations. FIAT and Icarus are also used for TPS sizing purposes. Today, these traditional tools are being supplemented with computational materials models at the atomistic level. The scales of interest range from computational chemistry (Density Functional Theory [DFT]), to atomistic simulations (Molecular Dynamics [MD]), to the microscale with the Porous Microstructure Analysis (PuMA) software that was recently awarded the 2022 NASA Software of the Year award. Finally, thermo-structural modeling is also of interest to the TPS Materials branch and done using commercial tools such as MSC MARC, MENTAT, NASTRAN and PATRAN. The present talk will also link the use of these computational tools to current NASA missions and projects associated with challenging and complex vehicles entries/reentries.

materials modeling↗
NASA NTRS2012

Influence of Catalysis and Oxidation on Slug Calorimeter Measurements in Arc Jets

Arc jet tests play a critical role in the characterization and certification of thermal protection materials and systems (TPS). The results from these arc jet tests feed directly into computational models of material response and aerothermodynamics to predict the performance of the TPS in flight. Thus the precise knowledge of the plasma environment to which the test material is subjected, is invaluable. As one of the environmental parameters, the heat flux is commonly measured. The measured heat flux is used to determine the plasma enthalpy through analytical or computational models. At NASA Ames Research Center (ARC), slug calorimeters of a geometrically similar body to the test article are routinely used to determine the heat flux. A slug calorimeter is a thermal capacitance-type calorimeter that uses the temperature rise in a thermally insulated slug to determine the heat transfer rate, see Figure 1(left). Current best practices for measuring the heat flux with a slug calorimeter are described in ASTM E457 - 96. Both the calorimeter body and slug are made of Oxygen Free High Conductivity Copper, and are cleaned before each run.

Nawaz, Anuscheh↗
NASA NTRS2017

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗
NASA NTRSDate not supplied

Overview of Additively Manufactured TPS Proposed Flight Test and Earth Re-Entry Capsule Design

A flight mechanics overview is presented of an Earth flight test designed to investigate a novel, 3D printed thermal protection system (TPS) that is currently in development at NASA as part of the Additive Manufacturing of Thermal Protection Systems project. The project is pioneering a method to print a thermal protection system onto an entry vehicle forebody one layer at a time. This method reduces labor and complexity as compared to traditional manufacturing methods while increasing mission-dependent customization of through-depth materials properties. The flight test has three objectives. First, subject the forebody stagnation point of a capsule equipped with additively manufactured TPS (AMTPS) material to peak heat fluxes in excess of 100 $W/cm^2$. Second, capture in-flight data to enable flight reconstruction and AMTPS material thermal response model improvement. Third, recover the capsule with data storage and forebody AMTPS intact to enable post-flight inspection and analysis of AMTPS performance. The flight test trajectory is designed to achieve a peak stagnation point, cold-wall, entry heat flux of 135~$W/cm^2$. Flight mechanics simulations are performed using the Program to Optimize Simulated Trajectories II (POST2) and Monte-Carlo analysis yields statistical percentiles on vehicle performance at key points along the trajectory. Based on the flight mechanics analysis presented in this paper, a prototype capsule was designed, partially fabricated, and underwent preliminary component stress testing in preparation for fabrication of the flight unit capsule. The capsule outer mold line is a modified version of the heritage Mars Microprobe geometry. The capsule has a 0.356~m diameter, a 30~kg mass, and a hypersonic ballistic coefficient of 300~$kg/m^2$. Sensor selection is guided by flight dynamics simulations with the goal of resolving the re-entry heating pulse. On-board instrumentation include forebody and aftbody pressure sensors and thermocouples, a 9-axis IMU, a GPS receiver, and an Iridium satellite modem, all of which collect and store data throughout flight via on-board avionics systems. A two-stage parachute system is designed to decelerate the capsule to touchdown velocities that will not result in significant fracture or deformation of the charred AMTPS material at ground impact.

Flight Mechanics↗
NASA NTRS2001

Optimized shielding for space radiation protection

Future deep space mission and International Space Station exposures will be dominated by the high-charge and -energy (HZE) ions of the Galactic Cosmic Rays (GCR). A few mammalian systems have been extensively tested over a broad range of ion types and energies. For example, C3H10T1/2 cells, V79 cells, and Harderian gland tumors have been described by various track-structure dependent response models. The attenuation of GCR induced biological effects depends strongly on the biological endpoint, response model used, and material composition. Optimization of space shielding is then driven by the nature of the response model and the transmission characteristics of the given material.

NASA Discipline Radiation Health↗
DOE OSTI2021

Systems Engineering Approach for Design and Implementation of a Gas Breech for Actinide Experiments at the JASPER Facility

The Joint Actinide Shock Physics Experimental Research (JASPER) Facility is located approximately 65 miles north of Las Vegas, Nevada at the Nevada National Security Site (NNSS). The primary mission is to conduct shock physics research on actinide materials in support of NNSA’s Stockpile Stewardship Program. JASPER experiment uses a two-stage light gas gun to accelerate projectiles into targets at velocities up to 8 km/s (17,000 mph). The first stage uses an ignited propellent to drive a piston to compress gas in the pump tube. At the second stage, the gas compression exceeds a specified pressure and a rupture valve at the end of the pump tube opens. This launches a projectile to impact the target. The target is housed inside the target assembly and the debris field is contained inside the primary target chamber (PTC). The PTC is placed inside the secondary confinement chamber (SCC) as an added protection against possible contamination. These major components are illustrated in Figure 1.JASPER is capable of generating and measuring data on the properties of radioactive chemical elements at high shock pressures, temperatures, and strain rates approximating the conditions in nuclear weapons by using a two-stage gas gun. The data is used to determine material equations-of-state and validate computer models of material response. The work advances predictive capability, thus ensuring confidence in the nuclear stockpile.

42 ENGINEERING↗
NASA NTRS2011

Comparative Measurements of Earth and Martian Entry Environments in the NASA Langley HYMETS Facility

Arc-jet facilities play a major role in the development of heat shield materials for entry vehicles because they are capable of producing representative high-enthalpy flow environments. Arc-jet test data is used to certify material performance for a particular mission and to validate or calibrate models of material response during atmospheric entry. Materials used on missions entering Earth s atmosphere are certified in an arc-jet using a simulated air entry environment. Materials used on missions entering the Martian atmosphere should be certified in an arc-jet using a simulated Martian atmosphere entry environment, which requires the use of carbon dioxide. Carbon dioxide has not been used as a test gas in a United States arc-jet facility since the early 1970 s during the certification of materials for the Viking Missions. Materials certified for the Viking missions have been used on every entry mission to Mars since that time. The use of carbon dioxide as a test gas in an arc-jet is again of interest to the thermal protection system community for certification of new heat shield materials that can increase the landed mass capability for Mars bound missions beyond that of Viking and Pathfinder. This paper describes the modification, operation, and performance of the Hypersonic Materials Environmental Test System (HYMETS) arc-jet facility with carbon dioxide as a test gas. A basic comparison of heat fluxes, various bulk properties, and performance characteristics for various Earth and Martian entry environments in HYMETS is provided. The Earth and Martian entry environments consist of a standard Earth atmosphere, an oxygen-rich Earth atmosphere, and a simulated Martian atmosphere. Finally, a preliminary comparison of the HYMETS arc-jet facility to several European plasma facilities is made to place the HYMETS facility in a more global context of arc-jet testing capability.

Splinter, Scott C.↗
NASA NTRS2013

Thermal and Structural Performance of Woven Carbon Cloth For Adaptive Deployable Entry and Placement Technology

Arcjet testing and analysis of a three-dimensional (3D) woven carbon fabric has shown that it can be used as a thermal protection system and as a load bearing structural component for a low ballistic coefficient hypersonic decelerator called ADEPT (Adaptive Deployable Entry and Placement Technology). Results of arcjet tests proved that the 3D woven carbon fabric can withstand flight-like heating while under flight-like biaxial mechanical loads representative of those encountered during shallow entry flight path angles into the atmosphere of Venus. Importantly, the arcjet test results have been used to extend a preliminary material thermal response model based on previous testing of the same 3D woven carbon fabric under uni-axial mechanical loading.

Carbon fabric↗
NASA NTRS2000

Optimized Shielding for Space Radiation Protection

Abstract. Future deep space mission and International Space Station exposures will be dominated by the high-charge and -energy (HZE) ions of the Galactic Cosmic Rays (GCR). A few mammalian systems have been extensively tested over a broad range of ion types and energies. For example, C3H10T1/2 cells, V79 cells, and Harderian gland tumors have been described by various track-structure dependent response models. The attenuation of GCR induced biological effects depends strongly on the biological endpoint, response model used, and material composition. Optimization of space shielding is then driven by the nature of the response model and the transmission characteristics of the given material.

Wilson, J. W.↗
DOE OSTI2020

Multiscale-Informed Modeling of High Temperature Component Response with Uncertainty Quantification

This report summarizes a joint effort between Argonne National Laboratory, Idaho National Laboratory, and Los Alamos National Laboratory to develop and deploy constitutive models targeted at predicting the life of Grade 91 alloy components subjected to high temperature environments typical of those that structural components in advanced nuclear reactors would experience. Two distinct, but complementary constitutive modeling approaches have been taken here. The first employs a phenomenological viscoplastic model for which parameters have been calibrated based on experimental data for a wide range of Grade 91 alloy that has undergone a variety of processing. A Bayesian approach was used to derive distributions of uncertain parameters for this model based on this data set. The second approach is a reduced order model suitable for engineering-scale analysis that is based on the results of a large set of mesoscale simulations. Mesoscale models allow for the microstructure and composition of a particular alloy to be directly taken into account in the computation of the viscoplastic response, but are computationally expensive, which makes it impractical to directly call those models for the material constitutive response in an engineering-scale simulation. The reduced-order representation of the response of the underlying model used here allows for an engineering-scale model to take into account the characteristics of the underlying microstructure, while only incurring a reasonable computational expense. Both of these approaches have different strengths, and are applicable for different parts of the design/analysis process. The phenomenological models can be readily parameterized based on a set of experimental data for a given class of materials and used for scoping calculations. Once a specific material is chosen and adequately characterized, the reduced order models can accurately predict the response of that specific alloy, and because the models are based on predictive models of the underlying microstructure, they can be used to more confidently predict the response under conditions in regions where there is limited experimental data. Both of these models have been integrated in the Grizzly code, which is used here to perform proof-of-concept uncertainty quantification analyses of a simple component under prototypical conditions. The built- in stochastic analysis capabilities in the MOOSE framework that Grizzly is built on are used here to run large sets of simulations for this uncertainty quantification analysis. As would be expected, because the reduced order models are developed for a much more tightly defined alloy, they predict tighter distributions of the time to failure than the phenomenological models, which are calibrated to a broader set of data. Also important is that these simulations demonstrate that a reduced order modeling approach can be successfully deployed to propagate uncertainties from the material scale to practical engineering-scale component simulations.

42 ENGINEERING↗
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
NASA NTRS2020

Quasi-Static and Dynamic Analysis of Composite Panels: Challenges in Engineering-Level Commercial Software

In engineering applications, where analysis time and efficiency are essential, analytical models and commercially available finite element (FE) material models are the two mainstream methods used to analyze composite crush scenarios. A reliable numerical model addressing the crushing response of composite materials requires a nonlinear progressive damage analysis which includes appropriate failure mechanisms such as matrix cracking, fragmentation, delamination, and fiber breakage. One might be able to address all these interacting failure responses through implementation/validation of a user defined material model or by utilization of a commercially available model, the latter being the most efficient engineering approach. In this paper, the quasi-static and dynamic crushing response of composite specimens, consisting of a C-channel and a corrugated geometry manufactured with carbon fiber/epoxy matrix fabric, are modeled using the progressive composite damage model within ABAQUS. The novel part of this work, which has not been previously addressed in the literature, is the utilization of a detailed parametric study to optimize composite damage model parameters for a given quasi-static crush simulation in ABAQUS. C-channel experiments were used to calibrate the model and then refine/optimize a set of material parameters that best fit the experimental results. These identical parameters were then used to simulate the crush response of the corrugated channel specimens, thereby making it a pure prediction and not a double recalibration. Next, the dynamic response of both panels were predicted, further confirming the predictive capability of the developed numerical model. The effect of parameter variation on the crush response is fully documented for future use in engineering applications utilizing this specific ABAQUS damage model.

Paria Naghipour↗