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 469 records · Page 26

A Deep Learning Pipeline for Optimizing Large-scale Phase Field Simulations

Phase field (PF) simulations are computationally expensive but remain a key analysis tool to understand the complex mechanisms of additive manufacturing (AM) processes. Each PF simulation-aided analysis requires thousands of node hours on leadership-class supercomputers. One of the main goals of these analyses is the study of microstructure evolution during the build process which begins with the onset of nucleation. Nucleation occurs under certain thermomechanical conditions which are not known a priori and many PF simulations are required to identify ranges of input thermo-mechanical parameters that can result in the onset of nucleation. Since many of the simulations do not result in nucleation, an analysis campaign often ends up wasting tremendous amounts of precious computing resources executing nucleation-absent simulations. The goal of this work is to design and train deep learning models to inform a PF simulation about the likelihood of the occurrence of nucleation in a future simulation time-step based on the state summary over a finite number of past time-steps of a running simulation. If the prediction determines that the running simulation is unlikely to reach nucleation in the allotted time, then its execution is stopped immediately ultimately resulting in vast reduction in wasted computations when accrued over all the PF simulations typically performed in a single or multiple analysis campaign(s). The paper presents the performance of a machine learning pipeline that uses a convolutional neural network (CNN) model to learn an embedding which is then used with a self-attention network to build a multi-task deep learning model to predict the likelihood of nucleation. The model also predicts the input parameters used in a simulation. Performance is compared with a baseline pipeline that uses an off-the-shelf LeNet-5 model to learn the initial embedding. Despite their smaller size, performance results indicate significant improvement in accuracy of the proposed models compared to the larger baseline models.

Kannan, Ramakrishnan {ramki}↗

Microstructure-Based Estimation of Strength and Ductility Distributions for \(\alpha +\beta \) Titanium Alloys

Titanium alloys are processed to develop a wide range of microstructure configurations and therefore material properties. While these properties are typically measured experimentally, a framework for property prediction could greatly enhance alloy design and manufacturing. In this work, a microstructure-sensitive framework is presented for the prediction of strength and ductility as well as estimates of the bounds in variability for these properties. The framework explicitly considers distributions of microstructure via new approaches for instantiation of structure in synthetic samples. The parametric evaluation strategy, including the finite element simulation package FEpX, is used to create and test virtual polycrystalline samples to evaluate the variability bounds of mechanical properties in Ti-6Al-4V. Critical parameters for the property evaluation framework are provided by measurements of single crystal properties and advanced characterization of microstructure and slip system strengths in 2D and 3D. Property distributions for yield strength and ductility are presented, along with the validation and verification steps undertaken. Comparisons between strain localization and slip activity in virtual samples and in experimental grain-scale strain measurements are also discussed.

36 MATERIALS SCIENCE↗

Phase-field simulations of fission gas bubbles in high burnup UO2 during steady-state and LOCA transient conditions

To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized this is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory's Marmot application for phase-field simulation of nuclear materials. Simulations of bubble growth in the high burnup region during steady-state conditions showed that initially overpressurized bubbles decreased in pressure during growth, but still remained above equilibrium pressure. During a LOCA transient, simulations of bubbles in the high burnup region showed that bubble size did not change significantly, and the pressure increase due to the transient was calculated and passed to a phase-field model of fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mirostructure Characterization of Friction Consolidated Copper-Nickel using a Machine Learning Approach: Developing Process to Microstructure Associations

Friction consolidation (FC) is a solid phase processing approach where discrete material forms such as powders, chips, nuggets, etc. are densified via shear deformation. The precursors are placed in a billet container and brought in contact with a rotating tool that applying the desirable amount of normal force. Under the combined action of the rotation and normal pressure, the discrete precursor is consolidated through porosity reduction and shear deformation. FC is increasingly being studied as an attractive approach to manufacturing fully dense parts from powder forms owing to its ability to mix, alloy and consolidate difficult-to-process precursors in minimal number of process steps. Material consolidation and deformation in shear consolidation processes have been studied extensively previously for different material combinations previously. However, despite the extensive research in this area, understanding of the mechanistic processes in pore consolidation, deformation-induced mixing and material solubility during FC is still evolving. Material development using solid phase processing approaches such as FC is often performed based on research experience/education, which can be biased. Conventional analysis and simulation tools in this area tend to be successful only when material thermodynamic pathways and microstructural evolution sequences resulting from processing are clearly defined or known. They are not as effective for emerging advanced manufacturing technologies where material evolution pathways are not well established. The ability to predict optimal process parameters based on material chemistry and bulk properties is essential to accelerate materials design and processing, as are an understanding of the relevant structure-processing-property relationships. These structure-processing-property-performance relationships are at the core of materials science research. Microstructure characterization provides the link to these four core areas, often through visualizing material microstructure using imaging techniques. However, linking microstructure image data (i.e., micrographs) to variables of interest (e.g., processing parameters, material chemistry) in a reproducible, generalizable, and quantitative manner is a significant challenge. Typically, quantitatively linking image data to processing history relies on significant domain knowledge and manual or subject matter expert (SME)-heuristic based image analysis. Such an approach to image analysis has the potential to be biased, inefficient, and difficult to replicate.

36 MATERIALS SCIENCE↗

Stochastic microstructural analysis of ceramic matrix composites using a high-fidelity multiscale framework

This paper presents a high-fidelity multiscale framework to simulate the mechanical behavior of ceramic matrix composites (CMCs), accounting for the complex micromorphology captured using detailed material characterization. First, high-resolution micrographs are obtained for the specific carbon fiber silicon-carbide-nitride matrix (C/SiNC) CMC to characterize the variability of the architectural features and manufacturing-induced defects at the microscale. An image processing algorithm is then used to precisely estimate the size and distribution of all subscale features and defects from the micrographs. The information is then used to generate a three-dimensional stochastic representative volume element (SRVE) to reconstruct microscale constituents accounting for the variability. As a result, the generated SRVEs are simulated using the high-fidelity generalized method of cells (HFGMC) micromechanics theory to investigate the effects of defects on the elastic properties of C/SiNC CMCs.

Ceramic Matrix Composites↗

Stochastic parametric skeletal dosimetry model for humans: General approach and application to active marrow exposure from bone-seeking beta-particle emitters

The objective of this study is to develop a skeleton model for assessing active marrow dose from bone-seeking beta-emitting radionuclides. This article explains the modeling methodology which accounts for individual variability of the macro- and microstructure of bone tissue. Bone sites with active hematopoiesis are assessed by dividing them into small segments described by simple geometric shapes. Spongiosa, which fills the segments, is modeled as an isotropic three-dimensional grid (framework) of rod-like trabeculae that “run through” the bone marrow. Randomized multiple framework deformations are simulated by changing the positions of the grid nodes and the thickness of the rods. Model grid parameters are selected in accordance with the parameters of spongiosa microstructures taken from the published papers. Stochastic modeling of radiation transport in heterogeneous media simulating the distribution of bone tissue and marrow in each of the segments is performed by Monte Carlo methods. Model output for the human femur at different ages is provided as an example. The uncertainty of dosimetric characteristics associated with individual variability of bone structure was evaluated. An advantage of this methodology for the calculation of doses absorbed in the marrow from bone-seeking radionuclides is that it does not require additional studies of autopsy material. The biokinetic model results will be used in the future to calculate individual doses to members of a cohort exposed to 89,90 Sr from liquid radioactive waste discharged to the Techa River by the Mayak Production Association in 1949–1956. Further study of these unique cohorts provides an opportunity to gain more in-depth knowledge about the effects of chronic radiation on the hematopoietic system. In addition, the proposed model can be used to assess the doses to active marrow under any other scenarios of 90 Sr and 89 Sr intake to humans.

59 BASIC BIOLOGICAL SCIENCES↗

A Discrete Dendrite Dynamics Model for Epitaxial Columnar Grain Growth in Metal Additive Manufacturing with Application to Inconel

Epitaxial columnar grain growth is a prevalent microstructural feature in the additive manufacturing (AM) of metal components such as Inconel, with cubic unit cell crystal lattice structure (face centered cubic (FCC) or body centered cubic (BCC)). These columnar grains evolve from the partly molten grains in the substrate or the solidified metal. This work proposes an efficient model to simulate the competitive growth of epitaxial columnar dendritic grains. The proposed model tracks the dynamic changes in the dendrites emanating from discrete points along the solid/liquid interface of a quasi-steady melt pool (MP). These dynamic changes include convergence and divergence of growing dendrites. The model is extended to predict the microstructure of large 3D parts and experimentally validated by comparing the simulation results for laser powder bed fusion (L-PBF) and wire-arc additive manufacturing (WAAM) processes. The microstructure and pole figures are predicted for Inconel 718 samples produced by L-PBF and Inconel 740H samples produced by WAAM processes. Finally, the model predictions compare well with the observed microstructure and pole figures results for both the L-PBF and WAAM processes.

20 FOSSIL-FUELED POWER PLANTS↗

Enhancing microsegregation during rapid directional solidification through ternary microalloying

The nano-cellular dendritic microstructure formed during rapid directional solidification in powder bed fusion additive manufacturing creates unique properties such as simultaneous improvement in strength and ductility. However, process control of microsegregation features remains challenging due to low sensitivity of critical solidification mechanisms to process parameters. This study leverages microalloying to achieve large changes in dendrite composition, microstructure, and interdendritic zone width during laser powder bed fusion without modifying process parameters. CALPHAD simulations predict that the addition of Zr significantly steepens the solidus line of the dilute Cu-Cr alloy system, leading to enhanced Cr rejection into the melt and greater than 95% reduction in solubility of Cr in the solidified Cu matrix. Experimental validation using time-of-flight secondary ion mass spectrometry and Kelvin probe force microscopy reveals that the ternary alloy containing 0.01 wt% Zr exhibited wider interdendritic regions compared to the binary, a significantly higher number of Cr-rich particles within interdendritic regions, near-complete ejection of oxygen impurities from the matrix, and greater nanoscale work function contrast. These features indicate more aggressive Cr segregation in the presence of Zr and a purer Cu matrix and provide a potentially robust method for engineering the nano-cellular dendritic solidification microstructure.

CALPHAD↗

Simulating Catalysis with Realistic Pellet Geometries Using Mesoflow: A Case Study of Catalytic Propane Dehydrogenation

We present a case study of catalytic propane dehydrogenation with our open-source multiphysics solver, Mesoflow. The solver was developed to simulate reactive flow coupled to heterogeneous catalytic reactions and deactivation in the context of complex, mesoscale geometry. The method leverages cartesian block-structured adaptive mesh refinement to capture realistic catalyst microstructural features acquired directly from X-ray computed tomography data. A kinetic model for propane dehydrogenation and catalyst deactivation was developed based on temporal analysis of products (TAP) reactor experiments. The TAP reactor experiments allow for precise characterization of intrinsic kinetic reaction steps which are implemented into Mesoflow simulations to model the spatial and temporal evolution of reactants, products, and catalyst active sites. The short-term and long-term deactivation behavior is studied by using XCT data collected from fresh and aged catalyst pellets, which exhibit different microstructural features. This study employs time-splitting algorithms to connect disparate reaction and flow timescales, enabling the simulations to achieve realistic deactivation timescales on the order of minutes while the flow time-scales for small particles (100 microns) are several milliseconds. We also introduce a flexible automated python script that writes the necessary files to construct a Mesoflow simulation from user-created chemical mechanisms. We will also introduce a few new features that are added to Mesoflow such as higher order schemes, implicit chemistry integrators and the ability to run on AMD and NVIDIA graphics-processing-units.

AMReX↗

The kinetics of static recovery by dislocation climb

Abstract The initial microstructure of a wide range of structural materials is conditioned by thermo-mechanical treatments such as hot-working, tempering, or solution annealing. At the elevated temperatures associated with these treatments the dislocation microstructure evolves, usually decreasing in density through a process known as static recovery. Despite its technological relevance, static recovery is not fully characterized from a theoretical standpoint, with even the controlling mechanisms subject to debate. In this study, a climb-enabled discrete dislocation dynamics (DDD) capability is leveraged to explore the kinetics of static recovery in pure Fe when controlled by dislocation climb. Quantitative data from these simulations is used to develop a revised static recovery law, and provides the parameters appropriate for predictive microstructure models in Fe. This law differs from previous analytical derivations invoking climb of dislocations, following the logarithmic trends typical of experimental observations where prior work did not. Direct comparison between the recovery law derived from DDD to experimental recovery data in alpha Fe shows strong agreement across a range of temperatures, and suggests that climb is the controlling mechanism for static recovery in pure metals.

36 MATERIALS SCIENCE↗

Systematic discovery of new nano-scale metastable intermetallic eutectic phases in laser rapid solidified Aluminum-Germanium alloy

Laser surface remelting of as-cast Al-Ge eutectic alloy is shown to produce ultrafine lamellar eutectic morphology with interlamellar spacing refined up to ∼60 nm and composed of FCC Al solid solution and unusual Al x Ge y intermetallic phases that do not form during near-equilibrium solidification. The microstructures are characterized and analyzed using a combination of selected area electron diffraction, high-resolution scanning transmission electron microscopy, energy dispersive X-ray spectroscopy to obtain high-resolution elemental maps, and atomistic modeling using density functional theory followed by atomic-scale image simulation. Depending on the local solidification conditions, the crystallography of the Al x Ge y intermetallic phases in the eutectic microstructure is either monoclinic (C 2/c) or monoclinic (P 2 1 ), with high densities of defects in both cases. This is in sharp contrast to the as-cast alloys that showed nominally pure Al and Ge phases with significant solute partitioning and equilibrium FCC and diamond cubic crystal structures, respectively. Corresponding kinetic phase diagrams are proposed to interpret the evolution of nano-lamellar eutectic morphologies with equilibrium Al and metastable Al x Ge y phases, and to explain increased solid solubility in the Al phases manifested by precipitation of ultrafine clusters of Ge. Furthermore, the reasons for the formation of these metastable eutectics under laser rapid solidification are discussed from the perspective of the competitive growth criterion.

Al-Ge eutectic↗

Elucidating the Onset of Plasticity in Sliding Contacts Using Differential Computational Orientation Tomography

Abstract Depending on the mechanical and thermal energy introduced to a dry sliding interface, the near-surface regions of the mated bodies may undergo plastic deformation. In this work, we use large-scale molecular dynamics simulations to generate “differential computational orientation tomographs” (dCOT) and thus highlight changes to the microstructure near tribological FCC alloy surfaces, allowing us to detect subtle differences in lattice orientation and small distances in grain boundary migration. The analysis approach compares computationally generated orientation tomographs with their undeformed counterparts via a simple image analysis filter. We use our visualization method to discuss the acting microstructural mechanisms in a load- and time-resolved fashion, focusing on sliding conditions that lead to twinning, partial lattice rotation, and grain boundary-dominated processes. Extracting and laterally averaging the color saturation value of the generated tomographs allows us to produce quantitative time- and depth-resolved maps that give a good overview of the progress and severity of near-surface deformation. Corresponding maps of the lateral standard deviation in the color saturation show evidence of homogenization processes occurring in the tribologically loaded microstructure, frequently leading to the formation of a well-defined separation between deformed and undeformed regions. When integrated into a computational materials engineering framework, our approach could help optimize material design for tribological and other deformation problems. Graphic Abstract .

36 MATERIALS SCIENCE↗

Phase field modeling of discontinuous dynamic recrystallization in hot deformation of magnesium alloys

We have developed a quantitative phase field model of discontinuous dynamic recrystallization (PF-DDRX) for investigating grain structure evolution and its impact on mechanical response in magnesium alloy during hot deformation. A set of isothermal compression tests were conducted by Gleeble thermo-mechanical simulator and used to determine PF-DDRX model parameters such as critical stresses of recrystallization, activation energies of dislocation, grain boundaries and recrystallization grain nucleation. In the simulations, the experimentally measured average grain size and distribution were used to generate initial microstructures. The results show that the typical ‘necklace’ microstructures observed in DDRX forms at the earlier stage of recrystallization. The inhomogeneity of microstructures first increases rapidly and then decreases gradually and reaches a constant with the increase of dynamic recrystallization volume fraction. Both predicted grain structure and mechanical response are in good agreement with experimental results. This demonstrates that the model parameters determined by experimental data are reliable for the developed PF-DDRX model to predict the grain structure evolution and flow behavior of AZ80 magnesium alloy under hot deformation. This work presents a path to combine experiments and modeling to give quantitative prediction and minimize experimental cost.

magnesium alloy, Phase field method, dynamic recry↗

Anisotropic colloidal interactions & assembly in AC electric fields

In this work, we match experimental and simulated configurations of anisotropic epoxy colloidal particles in high frequency AC electric fields by identifying analytical potentials for dipole-field and dipole-dipole interactions. We report an inverse Monte Carlo simulation algorithm to determine optimal fits of analytical potentials by matching simulated and experimental distribution functions for non-uniform liquid, liquid crystal, and crystal microstructures in varying amplitude electric fields. Two potentials that include accurate particle volume and dimensions along with a concentration dependent prefactor quantitatively capture experimental observations. At low concentrations, an effective ellipsoidal point dipole potential works well, whereas a novel stretched point dipole potential is found to be suitable at all concentrations, field amplitudes, and degrees of ordering. The simplicity, accuracy, and adjustability of the stretched point dipole potential suggest it can be applied to model field mediated microstructures and assembly of systematically varying anisotropic particle shapes.

42 ENGINEERING↗

Low-voltage magnetoelectric coupling in membrane heterostructures

Strain-mediated magnetoelectric (ME) coupling in ferroelectric (FE)/ferromagnetic (FM) heterostructures offers a unique opportunity for both fundamental scientific research and low-power multifunctional devices. Relaxor-FEs, such as (1-x)Pb(Mg 1/3 Nb 2/3 )O 3 -(x)PbTiO 3 (PMN-xPT), are ideal FE layer candidates because of their giant piezoelectricity. However, thin films of PMN-PT suffer from substrate clamping, which substantially reduces piezoelectric in-plane strains. Here, we demonstrate low-voltage ME coupling in an all-thin-film heterostructure that uses the anisotropic strains induced by the (011) orientation of PMN-PT. We completely remove PMN-PT films from their substrate and couple with FM Ni overlayers to create membrane PMN-PT/Ni heterostructures showing 90° Ni magnetization rotation with 3 V PMN-PT bias, much less than the bulk PMN-PT ~100-V requirement. Scanning transmission electron microscopy and phase-field simulations clarify the membrane response. These results provide a crucial step toward understanding the microstructural behavior of PMN-PT thin films for use in piezo-driven ME heterostructures.

36 MATERIALS SCIENCE↗

ExaCA grain structure predictions for laser powder bed fusion processing

This dataset contains simulated cross-sections of laser powder bed fusion grain structures produced using the microstructure model ExaCA, in turn using time-temperature history data produced using the heat transport model AdditiveFOAM. These predictions show the grain structure using various permutations of hatch spacing and nucleation density. Also contained in the dataset are the input files necessary to reproduce the results. AdditiveFOAM (https://github.com/ORNL/AdditiveFOAM) and ExaCA (https://github.com/LLNL/ExaCA), both open source software, are required to reproduce the results in this dataset using the given input files. This DOI was updated 12/09/2025.

36 MATERIALS SCIENCE↗

Understanding and Strategies for Controlled Interfacial Phenomena in Lithium-Ion Batteries and Beyond

Electrolyte chemistry and properties and electrode structure and chemical properties particularly at the interfaces are crucial for the development of advanced battery components. This project focused on elucidating the roles of the electrolyte and electrode on the formation and evolution of the SEI layer and cell electrochemical performance on silicon and on Li metal anodes. For Si nano/microstructures we evaluated lithiation, volume changes, reactivity, and chemo-mechanical transformations as functions of nanoparticle size, shape, presence of coatings, electrolyte composition, electron leakage to the electrolyte, and cycling. For Li metal anodes we provided detailed analyses of the electrode and electrolyte effects on interfacial reactivity and as a function of microstructural evolution and underlying stochasticity as dendrites nucleate and grow. We thoroughly investigated and simulated how the chemistry of the various components of the electrolyte and the electrode architecture may affect the electrochemical reactions as well as cell degradation. A multiscale modeling approach was utilized, where atomistic simulations informed about the microscopic behavior of the system, and intermediate time and length scales were investigated with mesoscopic models. Results were also tested against selected experiments from our collaborators. The combined theoretical-experimental strategy saves costs by utilizing predictions from computational analyses to guide experimentation.

25 ENERGY STORAGE↗

NSUF Report EBSD analyses on metallic fuel samples: Tips and improvements in sample preparation

EBSD (Electron Backscatter Diffraction) is an analytic technique used to evaluate the microstructure of a specimen. It can be used to evaluate interesting crystallographic data ranging from grains size to grain orientation, from grain boundaries to stresses. Moreover, in the field of nuclear materials EBSD can aid in estimating properties important for reactor performance, such as dislocation and residual strain [1]. EBSD patterns are generated by focusing a stationary electron gun on the specimen and detecting the backscatter electrons that satisfy Bragg conditions, and thus are diffracted by the crystal planes forming the so called Kikuchi bands. These bands correspond to each of the lattice diffracting crystal planes and are detected on a phosphor screen via a low-light CCD camera [2]. The patterns are generated within a small interaction volume located at a depth of less than 50-100 nm. Because of this, the data quality is extremely sensitive to the integrity of the crystallographic lattice order at the surface of the sample and thus to sample preparation and to oxidation [3]. While conventional sample preparation methods are widely used in materials science, they often do not lead to quality electron backscatter patterns needed for the characterization of nuclear materials [4]. This has been related by many authors to the challenges associated with handling radiological specimens, their rapid oxidation, and artifacts in such samples (e.g. high number of defects, precipitates and/or porosities) [4-8]. EBSD, however, can be an invaluable technique to evaluate microstructural evolution during fabrication and after irradiation at the required mesoscale for model and simulation. It also provides important information on grain structure necessary to understand irradiation effect, such as grain subdivision. Thus, it is of paramount importance to further the application of this technique to nuclear materials and especially irradiated fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗