Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “microstructure 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Simultaneous Development and Robust Optimization of a Microstructure Dependent Material

Recent microstructure characterization techniques combined with Symbolic Regression(SR)analysis has been proven to generate white box plasticity models well suited for incorporation into FEA software.The current work builds upon those efforts and demonstrates the applicability of Sequential Monte-Carlo (SMC) methods within SR analysis to condense model development and robust optimization into a single, co-dependent process. In this project, SMC methods provide a mechanism through which the observed microstructure features and associated variability can be incorporated into the discovery phase of model development and simultaneously recover approximate parameter distributions through SR analysis. The demonstration utilized a data set consisting of tensile test results from a limited number of sample specimens with corresponding EBSD data from which microstructure features were characterized.The maximum threshold stress model in the Visco-Plastic Self-Consistent (VPSC) code developed by Los Alamos National Laboratories was calibrated using mechanical test data.Synthetic volume elements with statistically equivalent microstructure were generated with DREAM3Dbased on the observed EBSD data. VPSC was used to simulate the corresponding tensile test response for each of the synthetic volume elements. The simulated microstructure and tensile test data was used astraining datafor SMC-SR algorithm and the resulting model was validated with data from the original empirical data set.

Karl Garbrecht↗
DOE OSTI2022

Understanding Phase and Interfacial Effects of Spall Fracture in Additively Manufactured Ti-5Al-5V-5Mo-3Cr

Additive manufactured Ti-5Al-5V-5Mo-3Cr (Ti-5553) is being considered as an AM repair material for engineering applications because of its superior strength properties compared to other titanium alloys. Here, we describe the failure mechanisms observed through computed tomography, electron backscatter diffraction (EBSD), and scanning electron microscopy (SEM) of spall damage as a result of tensile failure in as-built and annealed Ti-5553. We also investigate the phase stability in native powder, as-built and annealed Ti-5553 through diamond anvil cell (DAC) and ramp compression experiments. We then explore the effect of tensile loading on a sample containing an interface between a Ti-6Al-V4 (Ti-64) baseplate and additively manufactured Ti-5553 layer. Post-mortem materials characterization showed spallation occurred in regions of initial porosity and the interface provides a nucleation site for spall damage below the spall strength of Ti-5553. Preliminary peridynamics modeling of the dynamic experiments is described. Finally, we discuss further development of Stochastic Parallel PARticle Kinteic Simulator (SPPARKS) Monte Carlo (MC) capabilities to include the integration of alpha (α)-phase and microstructural simulations for this multiphase titanium alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Acceleration of Thermochemistry Solves in MOOSE and Pronghorn

This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

Visualizing Uranium Crystallization from Melt: Experiment-Informed Phase Field Modeling and Machine Learning

The focus of this project was to observe and simulate the solidification of uranium metal at the crystallographic level from its molten state. Melting experiments were conducted at two different scales to observe microstructural evolution using either a laboratory-scale induction furnace (hundreds of grams of metal) or a microscope heating stage (hundreds of milligrams of metal), respectively. Experimental parameters and characterization data were then used to inform a phase field model of gamma-U crystal growth as dendrites with or without secondary phase impurities in the form of uranium carbide particles. Finally, training datasets were generated by the phase field model as inputs to a neural network, developed with the aim of providing a faster, cheaper surrogate model for microstructural simulations within a given parameter space. Progress is reported herein for each of these task areas. Ultimately, 1) an optical microscope heating stage capability has been stood-up for uranium metal solidification studies, 2) a phase field model was advanced to simulate multiple uranium grains growing in the presence of carbide impurity particles and 3) a neural network was constructed and optimized to predict the microstructure features of individually growing uranium crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Neural Phase Simulation

The Neural Phase Simulation (NPS) is a package of codes for simulating microstructure evolution and accelerated molecular dynamics with deep neural-networks based surrogate models. NPS is designed to offer quantitatively accurate and computationally efficient simulation capabilities by leveraging modern machine-learning techniques. The primary intended use cases of NPS are training neural network surrogate models, though performing simulations on a single node is also supported. The NPS surrogate models can be trained from ground truth simulation methods, which are supposed to be accurate but expensive, such as molecular dynamics, phase field methods, kinetic Monte Carlo and discrete dislocation dynamics.

Zhou, Fei↗
DOE OSTI2021

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Microstructure Modeling of 3rd Generation Disk Alloys

The objective of this program is to model, validate, and predict the precipitation microstructure evolution, using PrecipiCalc (QuesTek Innovations LLC) software, for 3rd generation Ni-based gas turbine disc superalloys during processing and service, with a set of logical and consistent experiments and characterizations. Furthermore, within this program, the originally research-oriented microstructure simulation tool will be further improved and implemented to be a useful and user-friendly engineering tool. In this report, the key accomplishment achieved during the second year (2008) of the program is summarized. The activities of this year include final selection of multicomponent thermodynamics and mobility databases, precipitate surface energy determination from nucleation experiment, multiscale comparison of predicted versus measured intragrain precipitation microstructure in quench samples showing good agreement, isothermal coarsening experiment and interaction of grain boundary and intergrain precipitates, primary microstructure of subsolvus treatment, and finally the software implementation plan for the third year of the project. In the following year, the calibrated models and simulation tools will be validated against an independently developed experimental data set, with actual disc heat treatment process conditions. Furthermore, software integration and implementation will be developed to provide material engineers valuable information in order to optimize the processing of the 3rd generation gas turbine disc alloys.

Jou, Herng-Jeng↗
NASA NTRS2026

Hybrid Modeling Study on Grain Evolution in the Metal Welding Process and Its Potential Lunar Application

Metal is most commonly used structural material in a wide range of spacecraft, and welding is the principal method for joining metal components into functional systems. However, conducting welding experiments under extreme environments—such as microgravity or vacuum conditions in space—is prohibitively expensive and experimentally challenging. To overcome these limitations, multi-physics computational welding models provide a cost-effective and versatile alternative. In this work, the authors have developed a coupled thermal (fluid) microstructure simulation framework to model metal welding under varying gravity conditions. The framework integrates a mixed-mode heat transfer formulation (conduction, convection, and radiation) with molten pool fluid dynamics, enabling accurate prediction of temperature fields and weld-pool geometry. A grain growth model is further incorporated to capture the spatial and temporal evolution of microstructure, including grain size distribution and morphological transitions during solidification. This approach provides detailed insight into molten pool evolution and grain-level microstructure development throughout the welding process. By explicitly parameterizing environmental conditions, the model supports extrapolation to off-Earth manufacturing scenarios such as welding on the lunar surface. Tantalum—chosen in this study due to its high melting point, oxidation resistance, and mechanical stability at elevated temperatures—serves as the material system for model demonstration. Beyond Tantalum, the integrated multi-physics framework offers broad applicability for predictive welding simulations of various structural and refractory metals or alloys used in extreme terrestrial or extraterrestrial environments.

kinetic Monte Carlo (SPPARKS)↗
DOE OSTI2022

Understanding Uncertainty in Microstructure Evolution and Constitutive Properties in Additive Process Modeling

Coupled process–microstructure–property modeling, and understanding the sources of uncertainty and their propagation toward error in part property prediction, are key steps toward full utilization of additive manufacturing (AM) for predictable quality part development. The OpenFOAM model for process conditions, the ExaCA model for as-solidified grain structure, and the ExaConstit model for constitutive mechanical properties are used as part of the ExaAM modeling framework to examine a few of the various sources of uncertainty in the modeling workflow. In addition to “random” uncertainty (due to random number generation in the orientations and locations of grains present), the heterogeneous nucleation density N 0 and the mean substrate grain spacing S 0 are varied to examine their impact of grain area development as a function of build height in the simulated microstructure. While mean grain area after 1 mm of build is found to be sensitive to N 0 and S 0 , particularly at small N 0 and large S 0 (despite some convergence toward similar values), the resulting grain shapes and overall textures develop in a reasonably similar manner. As a result of these similar textures, ExaConstit simulation using ExaCA representative volume elements (RVEs) from various permutations of N 0 , S 0 , and location within the build resulted in similar yield stress, stress–strain curve shape, and stress triaxiality distributions. It is concluded that for this particular material and scan pattern, 15 layers is sufficient for ExaCA texture and ExaConstit predicted properties to become relatively independent of additional layer simulation, provided that reasonable estimates for N 0 and S 0 are used. However, additional layers of ExaCA will need to be run to obtain mean grain areas independent of build height and baseplate structure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Experimental and Numerical Characterization of High Temperature Deformation Behavior of 347H Stainless Steel

This study investigates how temperature affects the plasticity and thermal creep behavior of 347H stainless steel under uniaxial tension. The research combined experimental testing with advanced computational modeling. Two types of experiments were conducted: uniaxial tensile tests at temperatures from 100°C to 750°C using strain rates of ~10⁻⁴ s⁻¹, and creep tests at temperatures between 600°C and 750°C under various stress levels. These experimental results were used to develop and validate a new integrated mechanistic model that can predict material behavior under any loading condition while accounting for both stress and temperature effects. The model was implemented using a polycrystalline microstructure simulation framework based on elasto-viscoplastic Fast Fourier Transform (EVPFFT). It incorporates three key deformation mechanisms: thermally activated dislocation glide, dislocation climb, and vacancy diffusional creep. The model accounts for internal stress distribution within single crystals and considers how precipitates and solute atoms (both interstitial and substitutional) affect dislocation movement. After validation against experimental data, the model was used to generate Ashby-Weertman deformation mechanism maps for 347H steel, providing new insights into how microstructure influences the activation of different creep mechanisms.

Joy, Jobin K.↗
DOE OSTI2022

An algorithm for physics informed scan path optimization in additive manufacturing

Site specific microstructure control is a critical research area within the field of additive manufacturing due to its potential to revolutionize part performance. One way to achieve site specific microstructure control is through control of the solidification conditions via the construction of intricate scan paths; however, the search space for such a problem is large. Previous attempts only considered the solidification conditions at the top surface while also requiring either lots of manual-fine tuning or large amounts of computational resources. This paper introduces a general method for scan path optimization which considers the solidification conditions in the bulk of the material without an increase in computational expense. This method consists of three core components:1. A heat transfer model for simulating the temperature field at a given time.2. A surrogate model which takes scan pattern information and temperature data and predicts the solidification conditions of the bulk as well as the meltpool depths for a spot melt.3. A decision algorithm to decide which spot melt should be printed next based on the outputs of the surrogate model.Each of these components can be changed without changing the overall method. Within this work, this method is applied in the creation of an algorithm containing a semi-analytic heat transfer model to simulate the temperature field, a fully convolutional neural network (FCNN) as the surrogate model, and a greedy decision algorithm. The resulting algorithm produced complex scan patterns which gave strong results for simulated microstructure control.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Interplay of strain and phase evolution of laser powder bed fusion Ti–6Al–4V

While additive manufacturing (AM) provides a method of producing geometrically complex and highly detailed structures, the generation of residual strain in AM processes like laser powder bed fusion (L-PBF) can negatively impact performance-enabling properties. In applications such as orthopedic implants, specific performance windows require optimized microstructures in order to obtain desirable properties from multi-phase alloys like Ti–6Al–4V. This research aims to quantify the microscale origins of strain in L-PBF manufactured Ti–6Al–4V by understanding how strain is distributed at the grain and sub-grain scale, the interplay between phase evolution and strain, and examining post-processing strain relief strategies to control these features. Model spinal cage implants were manufactured from Ti–6Al–4V powder via L-PBF and then subjected to strain relieving heat treatment cycles above and below the Ti–6Al–4V β transus as a function of time and cooling rate. Residual strain was then studied via high resolution electron backscatter diffraction (HR-EBSD), and 2D strain maps with sub-micron resolution were generated for each post-processing state. It was found that macroscale thermal strains decreased with heat treatment time, but additional contributions from phase stabilizing residual strains retained primarily in the α' grains as lattice distortive strain remained. Additionally, the retention of β phase significantly changed the strain and dislocation distribution while reducing overall residual strain. In conclusion, these results were validated and reinforced with 3D mesoscopic micromechanical modeling of strain behavior across simulated microstructures, confirming that the local lattice dilation of α’ martensite is a primary contributor of microscale strain generation and retention in L-PBF Ti–6Al–4V.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling the thermodynamics of the FeTi hydrogenation under para-equilibrium: An ab-initio and experimental study

We report FeTi-based hydrides have recently re-attracted attention as stationary hydrogen storage materials due to favorable reversibility, good sorption kinetics and relatively low costs compared to alternative intermetallic hydrides. Employing the OpenCalphad software, the thermodynamics of the (FeTi) 1-x H x (0 ≤ x ≤ 1) system were assessed as a key basis for modeling hydrogenation of FeTi-based alloys. New thermodynamic data were acquired from our experimental pressure-composition-isotherm (PCI) curves, as well as first-principles calculations utilizing density functional theory (DFT). The thermodynamic phase models were carefully selected based on critical analysis of literature information and ab-initio investigations. Key thermodynamic properties such as dissociation pressure, formation enthalpies and phase diagrams were calculated in good agreement to our performed experiments and literature-reported data. This work provides an initial perspective, which can be extended to account for higher-order thermodynamic assessments and subsequently enables the design of novel FeTi-based hydrides. In addition, the assessed thermodynamic data can serve as key inputs for kinetic models and hydride microstructure simulations.

08 HYDROGEN↗
DOE OSTI2021

Metamaterial-like aerogels for broadband vibration mitigation

We present metamaterial-like aerogels that provide broadband low-frequency vibration mitigation through their self-assembled locally resonant microstructures. Simulations confirm the presence of a negative dynamic stiffness for the aerogels.

Malakooti, Sadeq↗
DOE OSTI2023

Fierro Version 2.x

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three-dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2023

Fierro

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗