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 235 records · Page 13

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗

Space Environmental Effects on Additively Manufactured Materials

Space simulations including thermal vacuum, atomic oxygen, and ultraviolet radiation were performed to study the durability of various additively manufactured materials. In addition to ground simulations, additively manufactured materials were selected for a one-year flight on the Materials on International Space Station Experiment (MISSE) Flight Facility. The space environment is composed of atomic oxygen, ultraviolet radiation, protons, electrons, meteoroid/space debris impacts, thermal cycling, and hard vacuum. An improved UV sensor is also discussed.

Finckenor, M. M.↗

Application of First-Principles-Based Artificial Neural Network Potentials to Multiscale-Shock Dynamics Simulations on Solid Materials

The use of artificial neural network (ANN) potentials trained with first-principles calculations has emerged as a promising approach for molecular dynamics (MD) simulations encompassing large space and time scales while retaining first-principles accuracy. To date, however, the application of ANN-MD has been limited to near-equilibrium processes. Here we combine first-principles-trained ANN-MD with multiscale shock theory (MSST) to successfully describe far-from-equilibrium shock phenomena. Overall, our ANN-MSST-MD approach describes shock-wave propagation in solids with first-principles accuracy but a 5000 times shorter computing time. Accordingly, ANN-MD-MSST was able to resolve fine, long-time elastic deformation at low shock speed, which was impossible with first-principles MD because of the high computational cost. This work thus lays a foundation of ANN-MD simulation to study a wide range of far-from-equilibrium processes.

36 MATERIALS SCIENCE↗

Compressing Hamiltonians with ab initio downfolding for simulating strongly-correlated materials on quantum computers

The accurate first-principles description of strongly correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving compressed many-body Hamiltonians that maintain the essential physics of strongly correlated materials. The solution of these material-specific models is still exponentially difficult to generate on classical computers, but quantum algorithms allow for a significant speed-up in obtaining the ground states of these compressed Hamiltonians. Here, we demonstrate that using quantum algorithms to obtain the properties of downfolded Hamiltonians can indeed yield high-fidelity solutions. By combining ab initio downfolding and variational quantum eigensolvers, we correctly predict the antiferromagnetic state of one-dimensional cuprate Ca 2 Cu O 3 , the excitonic ground state of monolayer W Te 2 , and the charge-ordered state of correlated metal Sr VO 3 . Numerical simulations using a classical tensor network implementation of variational quantum eigensolvers allow us to simulate large models with up to 54 qubits and encompassing up to four bands in the correlated subspace, which is indicative of the complexity that our framework can address. Through these methods we demonstrate the potential of classical preoptimization and downfolding techniques for enabling efficient materials simulation using quantum algorithms.

Alvertis, Antonios M. [NASA, Ames; LBNL, Berkeley]↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Space simulation test for thermal control materials

Tests were run in TRW's Combined Environment Facility to examine the degradation of thermal control materials in a simulated space environment. Thermal control materials selected for the test were those presently being used on spacecraft or predicted to be used within the next few years. The geosynchronous orbit environment was selected as the most interesting. One of the goals was to match degradation of those materials with available flight data. Another aim was to determine if degradation can adequately be determined with accelerated or short term ground tests.

Hardgrove, W. R.↗

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗

Aladyn - Adaptive Neural Network Molecular Dynamics Simulation Code: Computational Materials Mini-Application

This report provides an overview and commands description of the Computational Materials mini-application, Aladyn. Aladyn is a simple molecular dynamics code written in FORTRAN 2008, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of quantum mechanics-based energy calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. The Aladyn code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing many core architectures based on graphic processing units (GPUs). The effort is related to the High Performance Computing Incubator (HPCI) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗

Reactions and Surface Transformations of a Bone-Bioactive Material in a Simulated Microgravity Environment

A comprehensive program to investigate the expeditious in vitro formation of three-dimensional bone-like tissue is currently underway at the University of Pennsylvania. The study reported here forms a part of that program. Three-dimensional bone-like tissue structures may be grown under the simulated microgravity conditions of NASA designed Rotating Wall Bioreactor Vessels (RWV's). Such tissue growth will have wide clinical applications. In addition, an understanding of the fundamental changes that occur to bone cells under simulated microgravity would yield important information that will help in preventing or minimizing astronaut bone loss, a major health issue with travel or stay in space over long periods of time. The growth of three-dimensional bone-like tissue structures in RWV's is facilitated by the use of microcarriers which provide structural support. If the microcarrier material additionally promotes bone cell growth, then it is particularly advantageous to employ such microcarriers. We have found that reactive, bone-bioactive glass (BBG) is an attractive candidate for use as microcarrier material. Specifically, it has been found that BBG containing Ca- and P- oxides upregulates osteoprogenitor cells to osteoblasts. This effect on cells is preceded by BBG reactions in solution which result in the formation of a Ca-P surface layer. This surface further transforms to a bone-like mineral (i.e., carbonated crystalline hydroxyapatite (c-HA)). At normal gravity, time-dependent, immersion-induced BBG reactions and transformations are greatly affected both by variations in the composition of the milieu in which the glass is immersed and on the immersion conditions. However, the nature of BBG reactions and phase transformations under the simulated microgravity conditions of RWV's are unknown, and must be understood in order to successfully use BBG as microcarrier material in RWV'S. In this paper, we report some of our recent findings in this regard using experimental and numerical methods. BBG composition 45S5, the most reactive among known bone-bioactive glasses, was chosen for the study. BBG 45S5 behavior in physiological solutions was tested in simulated microgravity and compared with that at normal gravity. On the basis of our numerical study, we have chosen the BBG granule size to be in the range 40-70 microns, and a RWV rotational speed of 10 rpm. Our numerical study has shown that these parameters enable the microcarrier to remain suspended in the medium without experiencing collisions with the wall of the vessel. Immersion-induced changes in the solution composition and the material surface were analyzed after immersion.

Radin, S.↗

Impact Cratering Experiments Using Lunar Highland Simulant as Target Material

The properties of a target material have long been known to influence the cratering process for relatively small impact events, where the material “strength” is more dominant than the local gravitational force[1]. Strength in granular materials, however, is dependent on the interactions among packing fraction, grain size, and cohesion[2], and can thus vary for the same bulk material at different porosities. Experimentally, this relationship has been demonstrated for the most commonly used target material in impact cratering studies, unconsolidated sand, where [3] observed a notable decrease in depth-diameter ratios(d/D)for more densely packed (higher packing fraction, or lower porosity) targets. In order to extend the applicability, or properly scale, those findings to relatively small lunar craters, for example, a more realistic analog material(in terms of physical properties) should be interrogated for use as a target medium. To this end, we have conducted a series of impact cratering experiments using a lunar-highlands regolith simulant to investigate the relationship between packing fraction and crater morphology/morphometry.

impact cratering↗

Evaluation of the Mechanical Properties of Innovative Spacecraft Materials under a LEO Simulated Atomic Oxygen Environment

External spacecraft materials are crucial for protecting satellites from the harsh space environment. This study evaluates the radiation effects of low Earth orbit (LEO) space weather, focusing on high-energy electron and atomic oxygen (AO) exposure, on several modern spacecraft materials. The materials under investigation are “flight duplicates" of samples scheduled for launch as part of the Materials International Space Station Experiment Flight Facility (MISSE-FF) mission in 2025. To observe changes in stress distribution within weathered polymers, the experimental setup based on the photo elasticity technique was developed

Yuliya Kuznetsova↗

Organics on Fe-Silicate Grains: Potential Mimicry of Meteoritic Processes?

Currently, it is unknown what exact process or combination of processes produced organics that are found in meteorites or are detected in comets and nebulas. One particular process that forms organics involves Fischer-Tropsch type (FTT) reactions. Fischer-Tropsch type reactions produce hydrocarbons by hydrogenating carbon monoxide via catalytic reactions. The overall reaction is shown. The products of these reactions have been studied using natural catalysts and calculations of the efficiency of FTT synthesis in the Solar Nebula suggest that these types of reactions could make significant contributions to the composition of material near three AU. We use FTT synthesis to coat amorphous Fesilicate grains with organic material to simulate the chemistry in the early Solar Nebula. We used lab-synthesized amorphous Fe-silicate grains for the catalyst because they might better simulate the starting materials found in protostellar nebulas. A brief description of the synthesis of these grains is given in Experiments. This work is different from previous studies because we focus here on the carbonaceous material deposited on the grains. In our experiments, we roughly simulate a model of the nebular environment where grains are successively transported from hot to cold regions of the nebula. In other words, the starting cold regions of the nebula. In other words, the starting gases and FTT products are continuously circulated through the grains at high temperature with intervals of cooling. Overall, organics generated in this manner could represent the carbonaceous material incorporated into comets and meteorites. We present the analyses of the organics produced using pyrolysis gas chromatography mass spectrometry (GCMS) and compare the results with those organics found in the Murchison meteorite.

Johnson, N. M.↗

Tight-Binding-Based Method for Simulations of Complex Materials with Accompanying Database

The NRL-Tight-Binding (NRL-TB) method has been applied successfully to many monatomic, binary, and a few ternary systems. The method consists of a fit to Density Functional Theory (DFT) results and has transferability beyond the fitted DFT database including performing tight-binding molecular dynamics simulations. We are working on a new website with a detailed guide of how to use the method in performing the fit to DFT results (Fitting Code) and obtaining non-orthogonal TB parameters expressed in a polynomial form as environment-dependent Slater-Koster parameters (on-site terms), and neighbor distance-dependent Hamiltonian and Overlap terms. These parameters are then used in a Static Code to verify or make predictions for large systems often beyond the computational capabilities of standard DFT codes. Our starting point was a website developed at NRL some years ago which NRL no longer supports. We have now been augmenting that website with new examples of expanding the original capabilities of the method such as incorporating in the codes f-electron states. Our website (esd.cos.gmu.edu/tb) contains two separate branches: One named NRL-TB with the full documentation of the method including source codes and examples. The other branch, named ESD-dB, is a library of electronic structure calculations for all the monatomic elements (a subset of the elements also found in a recent book by the author of this proposal. See reference 22) and an approximately 700 hundred compounds in specific crystal structures. We discuss below the progress made so far in improving and expanding the capabilities of our TB method and the plans for further extensions of our approach.

36 MATERIALS SCIENCE↗

Initial framework for engineering-scale statistical creep-fatigue modeling

This report describes the integration of new solid and interface-cohesive mechanics systems into MOOSE. The purpose of these new systems is to support the ability of MOOSE to run full-field crystal plasticity finite element method simulations of key material processes in high temperature metallic materials. These simulations could be used to help accurately predict the performance of key high temperature structural materials in future advanced nuclear reactor components. Previous work implemented preliminary versions of many of these systems in MOOSE Apps. The current work reports on their integration into the main MOOSE tensor mechanics module along with associated improvements to the basic formulations and numerical implementations. Finally, the report provides an example of the full-field crystal plasticity simulations now possible in MOOSE, including examples of realistic geometries requiring millions of degrees of freedom to resolve the microstructural features and macroscale geometry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Micromorphic Length-Scale Coupling Framework for the Determination of Higher-Order Constitutive Models and the Multi-Scale Simulation of Heterogeneous Materials [Thesis]

Heterogeneous materials and materials with complex microstructures pose a unique challenge in the development of accurate models of their response to external stimuli. These difficulties principally arise due to the difficulty in characterizing and modeling the constituents and their interactions. It is usually possible, though non-trivial, to construct an explicit representation of the microstructure (a direct numeric simulation or DNS) but the method by which the complex modes of deformation and other processes can be homogenized to a reduced-order approximation is frequently unclear. Many homogenization approaches are ad-hoc and lack a strong justification beyond ease of computation. Furthermore, the homogenization approach can, in some cases, not utilize the full breadth of information available in the computation of the macro-scale stresses and deformation measures. It is also noteworthy that homogenization, by its very nature, will tend to obfuscate details of processes occurring at the lower length-scale. It is therefore of interest to include as much information as possible in the construction of the reduced-order model so as to be predictive in a variety of loading environments. We here present a length-scale bridging technique based upon the micromorphic continuum mechanics of Eringen which incorporates volume and surface area averages as a part of its construction. This approach enforces the balance equations at the micro-scale and then studies the effect of the spatially varying nature on the macro-scale. This leads, naturally, to further balance equations which are solved at the macro-scale. This work details a homogenization framework which arises naturally from the micromorphic construction of Eringen attempting to introduce no definitions beyond which are justifiable from micro-structural considerations. One of the results of this effort is the so-called “micromorphic filter” which has been developed and applied to several DNS to demonstrate its effectively. In order to determine the macroscopic degrees of freedom we utilize the special case of an overlap coupling technique where the macro-scale is fully constrained to the micro-scale. This enables us to study the resulting material properties as expected but also allows us to further study the boundary conditions on the additional degrees of freedom at the macro-scale.

42 ENGINEERING↗