Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials simulations”

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 1,117 records · Page 62

Improved Axisymmetric and High Temperature Material Structural Modeling in MOOSE and NEML

This report describes improvements made to the solid mechanics formulation in the MOOSE open source finite element simulation environment and the open source Nuclear Engineering Material model Library (NEML) for mechanical constitutive models. The focus of these improvements is to improve the usability and performance of simulations involving one or both pieces of software. Specifically, this work completes a new system for solid mechanics simulations in the MOOSE ecosystem providing exact linearizations and optimal (quadratic) convergence, for a variety of coordinate systems and material types, including large deformation simulations. This work then provides users a framework to build highly efficient mechanical simulations of structures or materials or to couple in additional MOOSE physics modules to build complex, scalable multiphysics simulations.

36 MATERIALS SCIENCE↗

Thermochemical Ablation Analysis of the Orion Heatshield

The Orion Multi-Purpose Crew Vehicle will one day carry astronauts to the Moon and beyond, and Orion's heatshield is a critical component in ensuring their safe return to Earth. The Orion heatshield is the structural component responsible for absorbing the intense heating environment caused by re-entry to Earth's atmosphere. The heatshield is primarily composed of Avcoat, an ablative material that is consumed during the re-entry process. Ablation is primarily characterized by two processes: pyrolysis and recession. The decomposition of in-depth virgin material is known as pyrolysis. Recession occurs when the exposed surface of the heatshield reacts with the surrounding flow. The Orion heatshield design was changed from an individually filled Avcoat honeycomb to a molded block Avcoat design. The molded block Avcoat heatshield relies on an adhesive bond to keep it attached to the capsule. In some locations on the heatshield, the integrity of the adhesive bond cannot be verified. For these locations, a mechanical retention device was proposed. Avcoat ablation was modelled in CHAR and the in-depth virgin material temperatures were used in a Thermal Desktop model of the mechanical retention device. The retention device was analyzed and shown to cause a large increase in the maximum bondline temperature. In order to study the impact of individual ablation modelling parameters on the heatshield sizing process, a Monte Carlo simulation of the sizing process was proposed. The simulation will give the sensitivity of the ablation model to each of its input parameters. As part of the Monte Carlo simulation, statistical uncertainties on material properties were required for Avcoat. Several properties were difficult to acquire uncertainties for: the pyrolysis gas enthalpy, non-dimensional mass loss rate (B´c), and Arrhenius equation parameters. Variability in the elemental composition of Avcoat was used as the basis for determining the statistical uncertainty in pyrolysis gas enthalpy and B´c. A MATLAB program was developed to allow for faster, more accurate and automated computation of Arrhenius reaction parameters. These parameters are required for a material model to be used in the CHAR ablation analysis program. This MATLAB program, along with thermogravimetric analysis (TGA) data, was used to generate uncertainties on the Arrhenius parameters for Avcoat. In addition, the TGA fitting program was developed to provide Arrhenius parameters for the ablation model of the gap filler material, RTV silicone.

Sixel, William↗

Friction and wear of carbon-graphite materials for high-energy brakes

Caliper type brake simulation experiments were conducted on seven different carbon graphite materials formulations against a steel disk material and against a carbon graphite disk material. The effects of binder level, boron carbide (B4C) additions, SiC additions, graphite fiber additions, and graphite cloth reinforcement on friction and wear behavior were investigated. Reductions in binder level, additions of B4C, and additions of SiC each resulted in increased wear. The wear rate was not affected by the addition of graphite fibers. Transition to severe wear and high friction was observed in the case of graphite-cloth-reinforced carbon sliding against a disk of similar composition. The transition was related to the disruption of a continuous graphite shear film that must form on the sliding surfaces if low wear is to occur.

Bill, R. C.↗

Tunable Nanoscale Evolution and Topological Phase Transitions of a Polar Vortex Supercrystal

Here, understanding the phase transitions and domain evolutions of mesoscale topological structures in ferroic materials is critical to realizing their potential applications in next-generation high-performance storage devices. Here, the behaviors of a mesoscale supercrystal are studied with 3D nanoscale periodicity and rotational topology phases in a PbTiO 3 /SrTiO 3 (PTO/STO) superlattice under thermal and electrical stimuli using a combination of phase-field simulations and X-ray diffraction experiments. A phase diagram of temperature versus polar state is constructed, showing the formation of the supercrystal from a mixed vortex and α-twin state and a temperature-dependent erasing process of a supercrystal returning to a classical α-twin structure. Under an in-plane electric field bias at room temperature, the vortex topology of the supercrystal irreversibly transforms to a new type of stripe-like supercrystal. Under an out-of-plane electric field, the vortices inside the supercrystal undergo a topological phase transition to polar skyrmions. These results demonstrate the potential for the on-demand manipulation of polar topology and transformations in supercrystals using electric fields. The findings provide a theoretical understanding that may be utilized to guide the design and control of mesoscale polar structures and to explore novel polar structures in other systems and their topological nature.

P=phase-field simulations↗

Ingrained: An Automated Framework for Fusing Atomic-Scale Image Simulations into Experiments

To fully leverage the power of image simulation to corroborate and explain patterns and structures in atomic resolution microscopy, an initial correspondence between the simulation and experimental image must be established at the outset of further high accuracy simulations or calculations. Furthermore, if simulation is to be used in context of highly automated processes or high-throughput optimization, the process of finding this correspondence itself must be automated. In this work, "ingrained," an open-source automation framework which solves for this correspondence and fuses atomic resolution image simulations into the experimental images to which they correspond, is introduced. Here, the overall "ingrained" workflow, focusing on its application to interface structure approximations, and the development of an experimentally rationalized forward model for scanning tunneling microscopy simulation are described.

36 MATERIALS SCIENCE↗

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.

Bioinorganic chemistry↗

Quantitative measurements of dislocations in metals for advancing predictive simulations

LLNL applications require scientists to predict how materials evolve under various thermomechanical conditions. While this is achieved through physics-based simulations, uncertainty in the predictions of mechanical properties remains a serious challenge that limits the predictive capabilities of models because we lack methods to compare predictions of atomic-scale defects (dislocations) with experimental measurements. High energy X-ray diffraction (HEXRD) is the most relevant technique that can provide the necessary statistical information on dislocations. However, this technique is not yet quantitative because we lack a precise understanding of the relationship between X-ray diffraction patterns and the underlying material dislocation content and arrangements. To address this need, we used our novel computational X-ray diffraction method to simulate the effect of dislocations on the diffraction patterns. We compared virtual and experimental diffraction patterns. Results allowed us to clearly establish the relationship between X-ray diffraction patterns and the underlying dislocation structures, proving that it is feasible to quantitatively measure dislocation statistics with HEXRD. This project delivered a method that can provide the missing piece to LLNL’s mechanical property simulations in advanced metals by obtaining experimentally long-needed quantitative dislocation data, which could fully enable predictive capabilities.

36 MATERIALS SCIENCE↗

2026 Annual Molten Salt Reactor Campaign Review

This report documents the 2026 Annual Molten Salt Reactor (MSR) Campaign Review held in Albuquerque, New Mexico, from April 21–24, 2026. The review was organized by Dr. Patricia Paviet, National Technical Director of the Advanced Reactor Technology (ART) MSR program for the U.S. Department of Energy Office of Nuclear Energy (DOE-NE), and included participation from principal investigators, federal managers, developers, regulators, and members of the broader MSR community. The 2026 review expanded beyond the ART-MSR campaign to include related DOE-NE programs supporting molten salt reactor advancement, including the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the Advanced Materials and Manufacturing Technologies (AMMT) program, Advanced Reactor Safeguards and Security (ARSS), Material Protection, Accounting and Control Technologies (MPACT), and relevant advanced fuels activities. Three days were dedicated to technical presentations and panel discussions covering thermal properties, off-gas management, modeling and simulation, safety and licensing, safeguards and security, materials and corrosion, and irradiation activities. A fourth day was dedicated to technical tours of Sandia National Laboratories and Kairos Power facilities. Attendance was strong and comparable to the prior annual review, with approximately 75 in-person participants per day, 75–80 virtual attendees per day, and approximately 30 participants in the tour day. The review fostered significant technical exchange across national laboratories, universities, industry, regulators, and international participants. This report summarizes the review structure, technical themes, participation, tours, feedback, and conclusions relevant to future planning for the MSR campaign.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Determining reference standard strength for neutron-irradiated reduced activation ferritic/martensitic steel F82H by Bayesian method

The deterministic approach widely adopted in the design of structural components relies on systematically defined design limits using empirically determined safety factors. However, this approach is not always appropriate because structures are subjected to a variety of loads in the practical environment, which may result in excessively conservative design limits. In recent years, a more rigorous probabilistic approach that incorporates material strength distributions has become an important solution. In the probabilistic approach, the probability density functions of material strength properties underpin the design criteria. Here, the objective of this study is to identify the density distribution functions that best describe tensile properties of irradiated F82H to define a reference strength for DEMO design. Due to the limited number of existing data, this study specifically employs a Bayesian prediction method based on Monte Carlo simulations to determine a material reference value with statistical reliability and to investigate its effectiveness. For example, the dependence of tensile properties of 300 °C irradiated materials on irradiation damage and the range predicted by 95% Bayesian estimation was evaluated. As a statistical model for the dose dependence of statistical parameters, the normal distribution exhibited a better fit for 0.2% proof strength and tensile strength, whereas the distribution of total elongation data gave comparable reference values for both the normal and Weibull distribution models. Both models gave comparable criteria for the distribution of total elongation data. The Weibull model also gave better results for uniform elongation. The function best describing the model was a logarithmic law for both 0.2% proof strength and tensile strength, while a power law for both total and uniform elongation, which allowed for more comprehensive data prediction of irradiation data with statistical accuracy for DEMO reactor design.

36 MATERIALS SCIENCE↗

Lost Circulation Materials in Geothermal Drilling: Thermal Degradation, Viscosity, Compression, and Flow Loop Tests

This dataset provides a set of experimental data on the performance of lost circulation materials (LCMs) used in geothermal drilling operations. It includes results from four distinct tests: thermal degradation, viscosity measurements, compression tests, and high-temperature flow loop experiments. The materials tested, such as cedar fiber, cotton seed hulls, magma fiber, microcellulose, and sawdust, were analyzed under a range of conditions to assess their effectiveness in mitigating fluid losses in geothermal wells. Data on thermal degradation explores mass loss, water loss, and gas release for each material when subjected to temperatures between 90C and 250C. Viscosity measurements capture the behavior of LCM mixtures under different temperatures, while compression tests investigate the structural integrity of both undegraded and thermally degraded materials. The high-temperature flow loop experiments simulate geothermal well conditions, analyzing how different materials and mixtures affect fluid flow and sealing performance in a fractured reservoir environment. In addition to the experimental data, the dataset includes photos documenting the materials, schematics of test setups, and several published manuscripts that provide detailed insights into the experimental methods and findings.

15 GEOTHERMAL ENERGY↗

Development of New Lunar Highland Regolith Simulant, NUW-LHT-5M

Introduction: NASA has a need for large quantities of lunar simulants that closely match future manned lunar missions at landing sites near the Lunar South Pole. The current simulant was designed to approximate NU-LHT-2M, and -4M, except using natural minerals and a fully synthetic, non-basaltic, high calcium glass. Washington Mills in Niagara Falls, NY was contracted to produce this new simulant due to their electric arc furnace (EAF) fusion technology and capabilities for crushing and sizing ceramics. The primary objectives of this work were: (1) to create a commercially produced simulant that contains a high fidelity, high calcium (An90+), low Mg/Fe glass without basaltic constituents; (2) to create a simulant that uses existing, terrestrially mined, and readily available minerals blended with a synthetic glass, and (3) accomplish objectives (1) and (2) using conventional processing methods capable of production levels that can meet substantial current and future demands of NASA projects. Another consideration was that highland simulants closest to meeting these requirements, NU-LHT-2M and -4M types, are out of production and largely out of stock. Design: The design chosen for development and production for this work is NU-LHT-2M, de-signed by Doug Stoeser and Douglas Rickman, which uses Stillwater anorthosite and norite rocks and synthetic glass. The composition targeted the average composition, glass content and particle size distribution of Apollo 16 samples [1]. The design utilizes two Stillwater minerals, anorthosite (37.7 wt.%), norite (17.6 wt.%), a commercially sourced olivine (4.7 wt.%), and a synthetic glass (40 wt.%). Processing: Stillwater rocks were hand collected and initially crushed by the USGS. Further crushing and milling was performed at Washing-ton Mills. For the glass, oxides were blended in the ap-propriate ratios, blended in a v-blender, and placed in a graphite lined, water cooled pot. The mixture was fused in Washington Mills’ pilot scale, 500 kW EAF. The molten glass was then poured into water and quenched to inhibit crystallization. Several pour/quench cycles were required to produce the total quantity of glass needed. Anorthosite, norite and olivine were milled together to the target particle size in a ball mill. The glass was ball milled separately. The complete simulant was produced by blending the milled materials. Particle size distribution, phase and chemical analyses of the materials was were per-formed at Washington Mills. Results: The glass composition was found to match the Apollo 16 oxide content average [1] quite well, with the exception of SiO2 and Fe2O3. The deviation of SiO2 and Fe2O3 from batched content was believed to be due to limitations inherent in the use of the arc furnace technique for these types of materials. A particle size distribution of the full simulant formulation closely matching the average for the Apollo 16 samples was also achieved. Fig. 1 shows the morphology of the final simulant particles. The size distribution is compared to the target in Fig. 2. Conclusions: A new lunar highlands type simulant containing a high-quality synthetic glass, closely matching an average composition of the Apollo 16 regolith was produced. Processing utilized easily scalable mineral and glass making methods.

simulant↗

Crystallization of the P 3 Sn 4 Phase upon Cooling P 2 Sn 5 Liquid by Molecular Dynamics Simulation Using a Machine Learning Interatomic Potential

We performed molecular dynamics simulations to study the crystallization of the P3Sn4 phase from P 2 Sn 5 liquid using a machine learning (ML) interatomic potential with desirable efficiency and accuracy. Our results capture the liquid properties of P 2 Sn 5 at 1300 K, which is well above the melting temperature. The phase separation and crystallization are observed when P 2 Sn 5 liquid is cooled down below 832 and 505 K, respectively. The simulation results are in good agreement with the experimentally observed phase transformation behaviors and provide useful insights into the complex nucleation and crystallization process at the details of atomistic scale. Our work also demonstrated that ML interatomic potentials based on neural network deep learning are robust and capable of accurately describing the energetics and kinetics of complex materials through molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Processing of Lunar Soil Simulant for Space Exploration Applications

NASA's long-term vision for space exploration includes developing human habitats and conducting scientific investigations on planetary bodies, especially on Moon and Mars. To reduce the level of up-mass processing and utilization of planetary in-situ resources is recognized as an important element of this vision. Within this scope and context, we have undertaken a general effort aimed primarily at extracting and refining metals, developing glass, glass-ceramic, or traditional ceramic type materials using lunar soil simulants. In this paper we will present preliminary results on our effort on carbothermal reduction of oxides for elemental extraction and zone refining for obtaining high purity metals. In additions we will demonstrate the possibility of developing glasses from lunar soil simulant for fixing nuclear waste from potential nuclear power generators on planetary bodies. Compositional analysis, x-ray diffraction patterns and differential thermal analysis of processed samples will be presented.

Sen, Subhayu↗

Simulation of Creep Deformation and Failure in Graded AM Microstructures

This report describes modeling tools and techniques developed to simulate the long-term material performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF). A physics-based Crystal Plasticity Finite Element model is used to simulate creep in microstructures and to study the roles of grain morphology, porosity, and texture. We describe our modeling methodology, including an orientation-mapping technique to capture the spatially varying crystallographic orientation that results from the build conditions. Our study of microstructural features shows that AM microstructures produced by LPBF tend to creep faster in the build direction, while texture and grain boundaries strengthen the transverse directions. However, when grain-boundary porosity and the consequent cavity growth are included in the model, the transverse directions begin to creep faster. In examining texture, the results indicate that spatially varying orientation arising from the build conditions increases anisotropy in the material, making it critical to account for orientation gradients in the material to accurately model its mechanical behavior. We also describe a material-model calibration campaign in which we calibrated the constitutive model specifically for LPBF 316H stainless steel at 725℃ for both solution-annealed and as-built conditions. Finally, these tools and techniques are used to model creep in microstructures representing different regions of an LPBF material with graded microstructure, owing to intentional variation in processing conditions. The creep simulation results show good agreement with experimental data across all three microstructures, with future work planned to study rupture in the material.

36 MATERIALS SCIENCE↗

A plane strain analysis of the blunted crack tip using small strain deformation plasticity theory

A deformation plasticity analysis of the tip region of a blunted crack in plane strain is presented. The power hardening material is incompressible both elastically and plastically, in order to simulate behavior of a stress freezing material above critical temperature. Stress and displacement fields surrounding the crack tip are presented. The results indicate that the maximum stress seen at the crack tip is indeed limited and is determined by the tensile properties; however, the scale over which the stresses act is dependent on the loading. Comparisons are good between the forward crack tip displacement and micro-fractographic measurments of stretch zones observed in plane strain fracture toughness tests.

Mcgowan, J. J.↗

Friction-Testing Machine

Testing machine evaluates wear and ignition characteristics of materials in rubbing contact. Offers advantages over other laboratory methods of measuring wear because it simulates operating conditions under which material will actually be used. Machine used to determine wear characteristics, rank and select materials for service with such active oxidizers as oxygen, halogens, and oxides of nitrogen, measure wear characteristics, and determine coefficients of friction.

Benz, F. J.↗

Current Developments in Future Planetary Probe Sensors for TPS

In-situ Thermal Protection System (TPS) sensors are required to provide traceability of TPS performance and sizing tools. Traceability will lead to higher fidelity design tools, which in turn will lead to lower design safety margins, and decreased heatshield mass. Decreasing TPS mass will enable certain missions that are not otherwise feasible, and directly increase science payload. NASA Ames is currently developing two flight measurements as essential to advancing the state of TPS traceability for material modeling and aerothermal simulation: heat flux and surface recession (for ablators). The heat flux gage is applicable to both ablators and non-ablators and is therefore the more generalized sensor concept of the two with wider applicability to mission scenarios. This paper describes the development of a microsensor capable of surface and in-depth temperature and heat flux measurements for TPS materials appropriate to Titan, Neptune, and Mars aerocapture, and direct entry. The thermal sensor will be monolithic solid state devices composed of thick film platinum RTD on an alumina substrate. Choice of materials and critical dimensions are used to tailor gage response, determined during calibration activities, to specific (forebody vs. aftbody) heating environments. Current design has maximum operating temperature of 1500 K, and allowable constant heat flux of q=28.7 watts per square centimeter, and time constants between 0.05 and 0.2 seconds. The catalytic and radiative response of these heat flux gages can also be changed through the use of appropriate coatings. By using several co-located gages with various surface coatings, data can be obtained to isolate surface heat flux components due to radiation, catalycity and convection. Selectivity to radiative heat flux is a useful feature even for an in-depth gage, as radiative transport may be a significant heat transport mechanism for porous TPS materials in Titan aerocapture. This paper also reports on progress to adapt a previously flown surface recession sensor, based on the Jupiter probe Galileo Analog Resistance Ablation Detector (ARAD), to appropriate aerocapture conditions.

Martinez, Ed↗

Development of Solid State Thermal Sensors for Aeroshell TPS Flight Applications

In-situ Thermal Protection System (TPS) sensors are required to provide verification by traceability of TPS performance and sizing tools. Traceability will lead to higher fidelity design tools, which in turn will lead to lower design safety margins, and decreased heatshield mass. Decreasing TPS mass will enable certain missions that are not otherwise feasible, and directly increase science payload. NASA Ames is currently developing two flight measurements as essential to advancing the state of TPS traceability for material modeling and aerothermal simulation: heat flux and surface recession (for ablators). The heat flux gage is applicable to both ablators and non-ablators and is therefore the more generalized sensor concept of the two with wider applicability to mission scenarios. This paper describes the continuing development of a thermal microsensor capable of surface and in-depth temperature and heat flux measurements for TPS materials appropriate to Titan, Neptune, and Mars aerocapture, and direct entry. The thermal sensor is a monolithic solid state device composed of thick film platinum RTD on an alumina substrate. Choice of materials and critical dimensions are used to tailor gage response, determined during calibration activities, to specific (forebody vs. aftbody) heating environments. Current design has maximum operating temperature of 1500K, and allowable constant heat flux of q=28.7 W/cm(sup 2), and time constants between 0.05 and 0.2 seconds. The catalytic and radiative response of these heat flux gages can also be changed through the use of appropriate coatings. By using several co-located gages with various surface coatings, data can be obtained to isolate surface heat flux components due to radiation, catalycity and convection. Selectivity to radiative heat flux is a useful feature even for an in-depth gage, as radiative transport may be a significant heat transport mechanism for porous TPS materials in Titan aerocapture.

Martinez, Ed↗