Search NASA⌕ Search

SEARCH · Search NASA

Results for “material response simulation”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Characterizing Porous and Nonporous Phenolic Resins from Molecular Dynamics Simulations

Phenolic resins are an important component of many ablative heat shield materials, which protect spacecrafts from the extreme temperatures reached during atmospheric entry. Examples include the high-density Heritage Carbon Phenolic (HCP) used in the Pioneer-Venus and Galileo missions, as well as the low-density Phenolic Impregnated Carbon Ablator (PICA) used in the Mars Science Laboratory and Mars 2020 missions. Additionally, recent developments within NASA have produced the mid-density Heatshield for Extreme Entry Environment Technology (HEEET) and its derivative 3D Woven Mid-Density Carbon Phenolic (3MDCP). Unlike the nonporous phenolic in HCP, PICA and HEEET/3MDCP are fabricated by infusing preforms with diluted phenolic formulations to obtain a lower density porous matrix. Despite the importance of the phenolic phase to the material response during entry, the variation in properties of porous and nonporous phenolic is not well understood. Here, we present an investigation of porous and nonporous phenolic resins using molecular dynamics (MD) simulations. Resin cure is mimicked in the simulations through the inclusion of representative reaction templates to generate accurate models of the complex crosslinked structures. To create porous models, explicit solvent molecules are included during the cure simulations. We observe nanoscale separation of the phenolic and solvent phases, which results in significant differences in the final structures of porous and nonporous models. In addition to a quantitative assessment of the network structure and porosity, we elucidate the effects of the phenolic formulation on the final material properties. These results are compared with experimental data as appropriate.

phenolic↗
DOE OSTI2023

Hybrid Particle Method for Computational Shock Physics

A long-standing area of research for Eulerian shock wave physics codes has been the treatment of strength and damage for materials. Here we present a method that will aid in the analysis of strength and failure in shock physics applications where excessive diffusion of critical variables can occur and control the solution outcome. Eulerian methods excel for large deformation simulations in general but are inaccurate in capturing structural behavior. Lagrangian methods provide better structural response, but finite element meshes can become tangled. Therefore, a technique for merging Lagrangian and Eulerian treatments of material response, within a single numerical framework, was implemented in the Multiple Component computational shock physics hydrocode. The capability is a Lagrangian/Eulerian Particle Method (LEPM) that uses particles to interface a Lagrangian treatment of material strength with a more traditional Eulerian treatment of the Equation of State (EOS). Lagrangian numerical methods avoid the advection diffusion found in Eulerian methods, which typically strongly affects strength constitutive law internal variables, such as equivalent plastic strain, porosity and/or damage. The Lagrangian capability enhances existing capabilities and permits accurate predictions of high rate, large deformation and/or shock of mechanical structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Structural evolutions of small aromatic mixtures under extreme temperature conditions: Insights from ReaxFF molecular dynamics investigations

Conversion mechanisms of carbonaceous materials during controlled pyrolysis under high temperature currently rely on empirical understanding, rather than a detailed atomistic description of the underlying chemical evolution. Here in this work, a series of molecular dynamics simulations, based on reactive force fields, were performed to probe the structural evolutions of small aromatic (number of aromatic rings ≤4) mixtures under extreme conditions (temperature 3000–4000 K). For any temperature studied and independent of the presence of alkane radicals, three major gas-phase byproducts, H 2 , C 2 H 2 , and C 2 H 4 , were found in bulk treated by rapid heating. Temperatures as well as alkane addition can affect the production/consumption rates and fractionation of these gas species. Compared to alkane radicals, temperatures also play important roles in modulating the initial growth pathways of solid carbon clusters. At 3000 K, aromatic carbon cluster formation is mainly initiated through the direct combination of aryl radicals. Conversely, at 4000K, carbon cluster formation is initiated by the assembly of chained radicals, since substantial polyacetylenic chains were produced. Through cross comparison on different simulation models/methods, as well as with experiments, the effect of processing conditions on carbonaceous materials conversion was revealed, and underlying mechanisms responsible for the observed phenomena are proposed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Finite-element modeling for an explosively loaded ferroelectric generator

A preliminary finite-element model has been developed using the ALEGRA-FE code for explosive driven depoling of a PZT 95/5 ferroelectric generator. The ferroelectric material is characterized using hysteresis-loop and hydrostatic depoling tests. These characteristics are incorporated into ALEGRA-FE simulations that model the explosive drive mechanism and shock environment in the material leading to depoling, as well as the ferroelectric response and the behavior of a coupled circuit. The ferroelectric-to-antiferroelectric phase transition is captured, producing an output voltage pulse that matches experimental data to within 10% in rise time, and to within about 15% for the final voltage. Both experimental and modeled pulse magnitudes are less than the theoretical maximum output of the material. Observations from materials characterization suggest that unmodeled effects such as trapped charge in the stored FEG material may have influenced the experimentally observed output.

42 ENGINEERING↗
DOE OSTI2021

VERA-Grizzly Ex-Core Calculations: Watts Bar Unit 1 Cycles 1-2

The critical structures that comprise light-water reactor (LWR) nuclear power plants are subjected to operating environments that can challenge their integrity. Structures in close proximity to the reactor core, such as the reactor pressure vessel (RPV) and the biological shield wall, are subjected to high levels of radiation emanating from the core, as well as elevated temperatures. As the US fleet of operating LWRs ages, the effects of these operating environments on the integrity of these structures must be considered to ensure their continued safe operation. Extending the lifetime of commercial reactors and maintaining the aging reactor fleet require accurate prediction of the exposure of ex-core components to neutron and photon radiation. In particular, concrete degradation studies must be performed to evaluate the safety and long-term operation of reactors with lifetime extensions. The concrete reactor bioshield is important for providing radiological protection during operation and must last for the entire lifetime of the reactor. Recent interest in lifetime extensions furthers the need to accurately simulate concrete material degradation in the reactor bioshield. As a result of this need, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded this study to couple its tools, Virtual Environment for Reactor Applications (VERA) and Grizzly. VERA allows users to set up models to calculate time-dependent and fully coupled solutions (with thermal feedback) for ex-core quantities of interest such as vessel and coupon fluence and detector responses for multiple statepoints and cycles. Grizzly is a finite-element application based on the Multiphysics Object Oriented Simulation Environment (MOOSE) framework that is used to enable aging materials calculations. This report highlights the work performed to calculate the fluence in the vessel and concrete for Watts Bar Nuclear Plant Unit 1 (WBN1) Cycles 1 and 2. The fluences obtained from VERA were successfully transferred to Grizzly using a Python script. Four simulations were run with Grizzly: (1) the Mazars model with the initial Young’s modulus being the instantaneous modulus, (2) the Mazars model with the initial Young’s modulus being the delayed modulus, (3) the Mazars model with the initial Young’s modulus being the delayed modulus with the addition of the effects of micro-damage caused by irradiation, and (4) the Mazars model with the initial Young’s modulus being the instantaneous modulus, and with the addition of micro-damage and creep. Details regarding the methods used to obtain the fluence and the statistical errors associated with the VERA Monte Carlo Shift calculations are discussed in greater detail in this report. The results obtained from the four Grizzly models are also presented in this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1980

Simulation of charging response of SCATHA (P78-2) satellite

A model of the satellite charging at high altitudes (SCATHA P78-2) satellite was used to simulate the charging response of SCATHA at geosynchronous orbit. The model includes a description of the geometry, currents to exposed surface materials, and electrical connections on the spacecraft. The charging response of the vehicle to that predicted by the NASCAP model for the Day 87, 1979 eclipse charging event, in which the spacecraft charged to several kilovolts negative during a magnetospheric substorm are compared. Double Maxwellian representations of the plasma environment reproduce the charging response observed experimentally.

Schnuelle, G. W.↗
NASA NTRS1982

Nonlinear structural and life analyses of a combustor liner

Three dimensional, nonlinear finite element structural analyses were performed for a simulated combustor liner specimen to assess the capability of nonlinear analyses using classical inelastic material models to represent the thermoplastic creep response of the one half scale component. Results indicate continued cyclic hardening and ratcheting while experimental data suggested a stable stress strain response after only a few loading cycles. The computed stress strain history at the critical location was put into two life prediction methods, strainrange partitioning and a Pratt and Whitney combustor life prediction method to evaluate their ability to predict cyclic crack initiation. It is found that the life prediction analyses over predicted the observed cyclic crack initiation life.

Moreno, V.↗
NASA NTRS1989

The effects of simulated space environmental parameters on six commercially available composite materials

The effects of simulated space environmental parameters on microdamage induced by the environment in a series of commercially available graphite-fiber-reinforced composite materials were determined. Composites with both thermoset and thermoplastic resin systems were studied. Low-Earth-Orbit (LEO) exposures were simulated by thermal cycling; geosynchronous-orbit (GEO) exposures were simulated by electron irradiation plus thermal cycling. The thermal cycling temperature range was -250 F to either 200 F or 150 F. The upper limits of the thermal cycles were different to ensure that an individual composite material was not cycled above its glass transition temperature. Material response was characterized through assessment of the induced microcracking and its influence on mechanical property changes at both room temperature and -250 F. Microdamage was induced in both thermoset and thermoplastic advanced composite materials exposed to the simulated LEO environment. However, a 350 F cure single-phase toughened epoxy composite was not damaged during exposure to the LEO environment. The simuated GEO environment produced microdamage in all materials tested.

Funk, Joan G.↗
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.↗
DOE OSTI2022

Hierarchical Materials from High Information Content Macromolecular Building Blocks: Construction, Dynamic Interventions, and Prediction

Hierarchical materials that exhibit order over multiple length scales are ubiquitous in nature. Because hierarchy gives rise to unique properties and functions, many have sought inspiration from nature when designing and fabricating hierarchical matter. More and more, however, nature’s own high-information content building blocks, proteins, peptides, and peptidomimetics, are being coopted to build hierarchy because the information that determines structure, function, and interfacial interactions can be readily encoded in these versatile macromolecules. Here, we take stock of recent progress in the rational design and characterization of hierarchical materials produced from high-information content blocks with a focus on stimuli-responsive and “smart” architectures. We also review advances in the use of computational simulations and data-driven predictions to shed light on how the side chain chemistry and conformational flexibility of macromolecular blocks drive the emergence of order and the acquisition of hierarchy and also on how ionic, solvent, and surface effects influence the outcomes of assembly. Furthermore, continued progress in the above areas will ultimately usher in an era where an understanding of designed interactions, surface effects, and solution conditions can be harnessed to achieve predictive materials synthesis across scale and drive emergent phenomena in the self-assembly and reconfiguration of high-information content building blocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Modeling brittle fracture due to anisotropic thermal expansion in polycrystalline materials

Here, this work investigated brittle fracture of polycrystalline materials due to thermal stresses arising from anisotropic thermal expansion. We used phase-field fracture simulations with the properties of alpha-uranium (α-U) and assumed a linear elastic mechanical response. Three-dimensional simulations were used to predict fracture for various conditions and crystallographic textures. We found that fracture was more pronounced during cooling than during heating because the anisotropy increased with temperature. We also found that the total crack surface area increased with increasing average misorientation, while the net shape change of the material decreased with increasing misorientation. Two-dimensional simulations in which one crystallographic coefficient of thermal expansion (CTE) was set to zero indicated that the expansion behavior in the crystallographic direction with the smallest CTE was the primary cause of fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

What are the Principles Controlling Biomimetic Heteropolymer Secondary Structure? (Final Technical Report)

The goal of the project was to develop improved theories to understand how nonbiological oligomers could be designed to cooperatively fold into 3D structures. These studies would lay the groundwork for materials made of such molecules, making it possible to create controlled and ordered materials for electron transport, efficient protein-like catalysts that work under extreme conditions, and sensors with highly-specific chemical responsiveness. Two different simulation thrusts were investigated, one focused on programs to identify stable low energy folded structures at a coarse-level of description of oligomers, and another to calculate thermodynamics of such oligomers. We used these theories to answer several specific questions about what properties of oligomers lead to cooperative transitions, and to identify how oligomer knots could serve as secondary structure elements. We also carried out significant collaborative investigation with Dr. Samuel Gellman (UW-Madison, National Academy of Sciences member) on stability for foldamers of interest to them. Only one of the experimentally tested foldamers stably folded, which was indicated by simulations as being the most likely to fold. Finally, we developed new theoretical descriptions of foldamers, showing how cooperativity was determined primarily by the entropy difference between the folded and unfolded state. The research did not answer all questions laid out in the original proposal but laid the groundwork for later efforts to design folded oligomers materials with high switchability.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analytical homogenization techniques applied to the Fickian diffusion: Effective diffusivity coefficient

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. In this work, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing TRISO particles. An extensive survey on effective thermal conductivity modeling was performed in our previous study, considering a random distribution of mono-sized spherical inclusions in a continuous matrix. Using the analogy between heat conduction and the simplified Fickian diffusion (or fission product species conservation), we can use the same analytical homogenization methods to obtain ETC as for the effective diffusivity coefficient (EDC). We performed several numerical experiments at varying conditions to assess the validity of our hypothesis for EDC calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1994

Computational simulation of composite structures with and without damage

A methodology is described which uses finite element analysis of various laminates to computationally simulate the effects of delamination damage initiation and growth on the structural behavior of laminated composite structures. The delamination area is expanded according to a set pattern. As the delamination area increases, how the structural response of the laminate changes with respect to buckling and strain energy release rate are investigated. Rules are presented for laminates of different configurations, materials and thickness. These results demonstrate that computational simulation methods can provide alternate methods to investigate the complex delamination damage mechanisms found in composite structures.

Wilt, Thomas F.↗
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↗
DOE OSTI2025

Establishing a process-structure-property-performance framework for SLS additive manufacturing through integrated multiscale modeling

This study presents a comprehensive suite of high-fidelity computational models that integrate multiscale and multiphysics simulations to capture the full Selective Laser Sintering (SLS) additive manufacturing process—from initial melting and solidification to mechanical response under external loads. Process simulations are linked with mechanical analysis through Representative Volume Elements (RVEs), establishing a process-structure–property-performance framework. The interaction between laser light and polyamide 12 (PA12) powder is modeled, accounting for laser characteristics and the optical, thermal, and geometrical properties of the powder. The heat source is incorporated into a heat transfer model, coupled with crystallization kinetics and densification models to predict material density and crystallinity. The porosity distribution from the densification model and crystallinity interpolated from experimental data are used to construct the RVEs. A multi-mechanism constitutive model is then calibrated using mechanical tests to predict the stress–strain response. Simulation results show good agreement with experimental data in terms of porosity, crystallinity, and mechanical performance when sufficient laser power (62 W or higher) is used. This research supports the inverse design of 3D-printed structures by introducing a high-fidelity framework that combines multiscale and multiphysics modeling with experimental calibration for predictive and performance-driven additive manufacturing.

SLS↗
DOE OSTI2021

Quantifying Temperature Dependence of Electronic Band Gaps and Optical Properties in SnO 2 and SnO via First-Principles Simulations

Tin metal oxides SnO x ( x = 1, 2) have gained interest as gas-sensing materials. For their applications as high-temperature sensors, however, a better understanding of their temperature dependence sensing responses is needed. Here, we comparatively quantify the temperature-dependent electronic band gaps and optical properties of SnO 2 and SnO using first-principles calculations. Without considering the temperature effect, SnO 2 and SnO are predicted to have direct and indirect band gaps of 2.18 and 1.66 eV, respectively, at the PBE + U -GGA level that we employed. The temperature effect on the electronic and optical properties is captured by taking account of the electron–phonon interaction. Band gap renormalization with temperature is calculated via the Allen–Heine–Cardona theory. For both oxides, we find a monotonic decrease in the electronic band gap such that renormalization at zero point (0 K) is ~–0.17 and ~–0.52 eV at 1000 K. These results are also analyzed by employing an analytical equation that helps characterize the band gap shift with temperature. In addition, the optical properties at finite temperatures are simulated using the frozen-phonon method that combines electron–phonon coupling with the momentum matrix. As temperature increases, the optical property spectra are smoothed because of the smearing effect, which diminishes optical constants at shorter wavelengths. Our results are of interest for high-temperature functional materials in applications of optical detection.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗