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 433 records · Page 24

Scan strategies in EBM-printed IN718 and the physics of bulk 3D microstructure development

Three-dimensional (3D) characterization provides opportunities for understanding processing-structure relationships in additively manufactured (AM) materials. Bulk samples of Inconel 718 were fabricated via electron beam melting (EBM) in order to study microstructural development as a function of energy input and beam scan strategy. TriBeam tomography of bulk Inconel 718 microstructures built under steady-state growth conditions reveals the sensitivity of microstructure formation and evolution to machine process parameters. Here, samples manufactured using a narrow range of energy input per unit build area result in varied grain morphologies and crystallographic textures. Using TRUCHAS, a thermal simulation software, the thermal history of bulk scan strategies was predicted, and combined with a calibrated microstructure-processing map to accurately predict bulk grain morphologies. The solidification parameters and the 3D measured nucleation density are used to predict the transition between columnar and equiaxed grain morphologies, providing a process map to guide AM parameter choices to locally control as-printed microstructure. A two-dimensional metric for characterizing bulk grain morphology was also found to agree well with predictions from the process map calibrated by 3D data. Combined with 3D tomography and thermal modelling, the physics of structure development were understood at a new level of detail with respect to the competing processes of grain nucleation and epitaxial growth.

36 MATERIALS SCIENCE↗

New Mechanism for Yield Point Phenomena

Yield point phenomena (YPP) are widely attributed to discrete dislocation locking by solute atmospheres. An alternate YPP mechanism was recently suggested by simulations of Ta single crystals without any influence of solutes or discrete dislocations. The general meso-scale (GM) simulations consist of crystal plasticity (CP) plus accounting for internal stresses of geometrically necessary dislocation content. GM predicted the YPP while CP did not, suggesting a novel internal stress mechanism. The predicted YPP varied with crystal orientation and boundary conditions, contrary to expectations for a solute mechanism. The internal stress mechanism was probed by experimentally deforming oligocrystal Ta samples and comparing the results with independent GM simulations. Strain distributions of the experiments were observed with high-resolution digital image correlation. A YPP stress–strain response occurred in the 0–2% strain range in agreement with GM predictions. Shear bands appeared concurrent with the YPP stress–strain perturbation in agreement with GM predictions. At higher strains, the shear bands grew at progressively slower rates in agreement with GM predictions. It was concluded that the internal stress mechanism can account for the existence of YPP in a wide variety of materials including ones where interstitial-dislocation interactions and dislocation transient avalanches are improbable. Furthermore, the internal stress mechanism is a CP analog of various micro-scale mechanisms of discrete dislocations such as pile-up or bow-out. It may operate concurrently with strain aging, or either mechanism may operate alone. A suggestion was made for a future experiment to answer this question.

42 ENGINEERING↗

Highly Ordered Eutectic Mesostructures via Template‐Directed Solidification within Thermally Engineered Templates

Template-directed self-assembly of solidifying eutectics results in emergence of unique microstructures due to diffusion constraints and thermal gradients imposed by the template. Here, the importance of selecting the template material based on its conductivity to control heat transfer between the template and the solidifying eutectic, and thus the thermal gradients near the solidification front, is demonstrated. Simulations elucidate the relationship between the thermal properties of the eutectic and template and the resultant microstructure. The overarching finding is that templates with low thermal conductivities are generally advantageous for forming highly organized microstructures. When electrochemically porosified silicon pillars (thermal conductivity < 0.3 Wm −1 K −1 ) are used as the template into which an AgCl-KCl eutectic is solidified, 99% of the unit cells in the solidified structure exhibit the same pattern. In contrast, when higher thermal conductivity crystalline silicon pillars (≈100 Wm −1 K −1 ) are utilized, the expected pattern is only present in 50% of the unit cells. The thermally engineered template results in mesostructures with tunable optical properties and reflectances nearly identical to the simulated reflectances of perfect structures, indicating highly ordered patterns are formed over large areas. This work highlights the importance of controlling heat flows in template-directed self-assembly of eutectics.

36 MATERIALS SCIENCE↗

Quantification of grain rotation in bulk nanostructured copper via in situ heating white beam synchrotron X-ray diffraction

Grain rotation during microstructural relaxation under heating is conventionally studied extensively through transmission electron microscopy and simulations. However, there is a shortage in examining grain rotations at larger length and volume scales in bulk materials. It is critical to understand the thermal stability of bulk nanostructured metals, since those enhanced mechanical properties have been well recognized. Here, this study aimed to employ a white beam microdiffraction X-ray technique under in situ heating from 300 K to 1073 K at 12 K/min on a nanostructured copper processed by high-pressure torsion, yielding an initial grain size of ∼260 nm before the heating. Evaluation across a wide range of temperatures reveals transition temperatures associated with microstructural relaxation processes. By tracking separate Laue diffraction peaks stemming from individual grains, changes in their orientations can be estimated and quantified. This approach becomes particularly effective when the number of grains within the probed volume is reasonably small. At temperatures above 940 K, about 7.5 % of Cu grains are rotating at speeds of ∼2 × 10 –3 °/s. The radial direction of the disk specimen is found to be the preferred direction of grain rotation, with the rotation axis along the shear direction. Further analysis found no correlation between peak intensity and peak size, related to grain sizes and deviatoric strain, respectively, with rotation conditions or speed. These findings demonstrate that the diffraction technique utilizing white beam X-rays is a complementary and novel characterization tool for understanding microstructural evolution, especially in grain rotations, of bulk crystalline materials.

Grain orientation↗

Machine learning insights into microstructural origins of transport and mechanical properties in porous microstructures

Multifunctional porous materials are increasingly needed across various fields, but their complex microstructures create significant challenges due to the intricate microstructure-property relationships. This complexity, combined with limitations of traditional analysis methods, hinders efforts to understand and optimize microstructure–property relationships. Here, to address this, we integrate physics-based mesoscale modeling with interpretable machine learning (ML) to uncover how microstructural features govern effective diffusivity and elastic modulus. At constant porosity, we show diffusivity varies by over 150 × and modulus by ∼50 ×, highlighting the power of microstructure engineering. Statistical analysis reveals bimodal behavior in diffusivity and unimodal in modulus. ML identifies connectivity as the dominant factor, while modulus is also sensitive to domain size and feature interactions. Controlled simulations further highlight domain shape as a critical feature for modulus. This framework enables efficient exploration of microstructure-property correlations, offering new insights to guide the design of advanced porous materials.

Bicontinuous microstructure↗

Effect of decades of corrosion on the microstructure of altered glasses and their radiation stability

Understanding the microstructural evolution of glasses during their interaction with water and radiation is of fundamental importance in addressing the corrosion of nuclear waste forms under geological disposal conditions. Here we report the results of more than 21 years of corrosion of two borosilicate glasses showing the formation of mesoporous C–S–H gels in Ca-bearing glasses and a mainly microporous microstructure in Al-bearing glasses. These porous corroded glasses were then irradiated with heavy ions to simulate the effects of recoil nucleus damage and monitored in real time using transmission electron microscopy with in situ ion irradiation. The ballistic collisions remarkably healed the porous corroded glasses to a pore-free homogeneous microstructure. Besides providing new insights and predictions about how doped glasses and actual waste forms may evolve under corrosion and irradiation, the results highlight the non-universal nature of the existing corrosion models and the important role that the glass composition and radiation damage play in the evolution of the microstructure during corrosion.

36 MATERIALS SCIENCE↗

The effect of temperature and burnup on U-10Zr metallic fuel chemical interaction with HT-9: A SEM-EDS study

The fuel cladding chemical interaction (FCCI) between Uranium-Zirconium-based metallic fuel and cladding materials during in-pile service is one of the most constraining phenomena affecting the performance of this fuel system. In this study, we investigated the effect of temperature and burnup on the FCCI development in two U-10 wt.% Zr (U-10Zr) fuel samples with HT-9 cladding irradiated as part of the MFF-3 irradiation test in the Fast Flux Test Facility (FFTF). One sample achieved a burnup of 13.1 at.% and operated with an average inner cladding temperature of 530°C, while the other achieved a burnup of 8.5 at.% and was subjected to an average inner cladding temperature of 615°C. Automated scanning electron microscopy (SEM) back-scattered electron (BSE) imaging of entire fuel cross-sections and SEM energy dispersive x-ray spectroscopy (EDS) analysis on specific fuel-cladding interface regions successfully provided a comprehensive characterization of the depth and type of interaction happening under different irradiation conditions. Further, our analysis shows that FCCI development on both fuel and cladding side is strongly influenced by the inner cladding temperature and, to some extent, the formation and integrity of Zr-rich layers between the fuel and cladding, while the impact of burnup and power is negligible. Measured FCCI thicknesses were compared to BISON simulations using both an empirical model based upon legacy data from the Experimental Breeder Reactor II (EBR II) irradiations and a mechanistic model currently under development for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, showing satisfactory agreement. Nonetheless, this comparison supports the need for additional microstructural characterization in intermediate ranges of temperature, power, and burnups in prototypic-length pins.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dendrite Growth Morphology Modeling in Liquid and Solid Electrolytes

The main goal of this project is to develop a multi-scale modeling approach that connects micron-scale phase-field models and atomic-scale density functional theory (DFT)-based simulations via parameter- and relationship-passing in order to predict Li-metal dendrite morphology evolution, in both liquid and solid electrolytes. The key hypothesis of the DFT-informed phase-field multiscale modeling approach is that it can capture the electrochemical-mechanical driving forces and incorporate the roles of nano-meter-thin solid electrolyte interphase (SEI) in liquid electrolytes as well as of the microstructures of micro-meter-thick solid electrolytes (SEs) for all-solid-state batteries. In this project, we have formulated and implemented phase-field models to incorporate the electrochemical driving forces in liquid electrolytes and then incorporate mechanical driving forces to simulate dendrite growth in solid electrolytes with resolved microstructures. We have implemented two treatments for the SEI: an explicit model to include the microstructure of the SE or SEI in the phase field model and an implicit model to simulate the impact of nano-meter thick SEI in liquid electrolytes by varying the electrode/electrolyte interfacial properties. The key interfacial properties, including the electronic and ionic transport properties, the charge transfer reaction kinetics, and mechanical properties, were computed by DFT-based calculations. At the DFT-based model, one key advancement is to directly predict the charge transfer reaction kinetics at a complex Li/SEI/electrolyte interface by linking DFT with density functional tight binding (DFTB) calculations. As the main accomplishments, we have demonstrated two successful predictions in both solid electrolyte and liquid electrolyte based on this multiscale approach. The predicted intergranular Li dendrite growth in LLZO revealed the importance of trapped electrons at internal interfaces in the microstructure of LLZO. The predicted electroplating morphology of mossy Li and faceted Mg agreed well with experiments. The insights provided by the multiscale model and the model enabled electrolyte and SEI design will accelerate the development of Li-metal electrode for high energy density batteries, that meet DOE’s target on cell density (>350 Wh/kg) and cost below $100/kWhuse for EV applications.

25 ENERGY STORAGE↗

Micropolar deep material network

This study extends the Deep Material Network (DMN), a physics-informed machine learning framework, to predict the homogenized mechanical response of composite materials with micropolar (Cosserat-type) constitutive behavior. This extension incorporates microstructure-dependent size effects, enabling accurate, efficient, and size-aware predictions for composites with complex internal architectures. While traditional, direct numerical simulation micropolar models effectively capture size effects by introducing extra local degrees of freedom, they bring significant computational challenges, particularly for multiscale analyses relevant to engineering applications. The micropolar DMN developed in this paper achieves high accuracy while significantly reducing computation time compared to micropolar direct numerical simulations. This advancement enables multiscale analyses and parameter studies that were previously impractical, such as high-cycle fatigue simulations and comprehensive investigations of internal length scale effects notably in size-dependent plastic response and the optimization of lattice structures. By uniting microstructure-sensitive modeling, physics-driven learning, and scalable surrogate modeling, the micropolar DMN paves the way for accelerated material design, large-scale parametric studies, and the reliable incorporation of size-dependent effects across a wide range of engineering applications, including optimization and next-generation composite design.

36 MATERIALS SCIENCE↗

Aging heat treatment design for Haynes 282 made by wire-feed additive manufacturing using high-throughput experiments and interpretable machine learning

Wire-feed additive manufacturing (WFAM) produces superalloys with complex thermal cycles and unique microstructures, often requiring optimized heat treatments. To address this challenge, we present a hybrid approach that combines high-throughput experiments, precipitation simulation, and machine learning to design effective aging conditions for the WFAM Haynes 282 superalloy. Our results demonstrate that the γ’ radius is the critical microstructural feature for strengthening Haynes 282 during post-heat treatment compared with the matrix composition and γ’ volume fraction. New aging conditions at 770°C for 50 hours and 730°C for 200 hours were discovered based on the machine learning model and were applied to enhance yield strength, bringing it on par with the wrought counterpart. This approach has significant implications for future AM alloy production, enabling more efficient and effective heat treatment design to achieve desired properties.

CALPHAD↗

Microstructure-based modeling of inner oxygen pressure in solid oxide electrolysis cells: Analysis of electrode delamination and mitigation

One major degradation mechanism during long-term operation of solid oxide electrolysis cells (SOECs) is delamination of oxygen electrodes (OEs). The driving force for the electrode delamination could be the generated high inner oxygen pressure near the electrode-electrolyte interface during operation. However, the effects of transport properties and electrode thickness on the inner oxygen partial pressure are not well understood. Here a microstructure-based electrochemical model, which includes the conduction of electrons and oxygen ions coupled with Butler-Volmer-type chemical reactions at triple-phase-boundaries (TPBs), is employed to investigate the oxygen pressure in lanthanum strontium manganate (LSM)-based SOECs. The model is applied to both two-dimensional (2D) prototype microstructures and three-dimensional (3D) realistic microstructures, and the oxygen pressure is analyzed as a function of transport properties and electrode thickness under both potentiostatic and galvanostatic operations. The simulation results suggest strategies to suppress electrode delamination. The simulation results are compared to an analytical solution, and the discrepancies are attributed to the Butler-Volmer-type kinetics included in the microstructure-based model.

25 ENERGY STORAGE↗

Microstructure-Based Modeling of Inner Oxygen Pressure in Solid Oxide Electrolysis Cells

One major degradation mechanism during long-term operation of solid oxide electrolysis cells (SOECs) is delamination of oxygen electrodes (OEs). The driving force for the electrode delamination could be the generated high inner oxygen pressure near the electrode-electrolyte interface during operation. However, the effects of transport properties and electrode thickness on the inner oxygen partial pressure are not well understood. Here, a microstructure-based electrochemical model which includes the conduction of electrons and oxygen ions coupled with Butler-Volmer-type chemical reactions at triple-phase-boundaries (TPBs), is employed to investigate the oxygen pressure in lanthanum strontium manganate (LSM)-based SOECs. The model is applied to both two-dimensional (2D) prototype microstructures and three-dimensional (3D) realistic microstructures, and the oxygen pressure is analyzed as a function of transport properties and electrode thickness under both potentiostatic and galvanostatic operations. The simulation results suggest strategies to suppress electrode delamination. The simulation results are compared to an analytical solution, and the discrepancies are attributed to the Butler-Volmer-type kinetics included in the microstructure-based model.

Xue, Fei↗

Phase transformation kinetics model for metals

We develop a new model for phase transformation kinetics in metals by generalizing the Levitas–Preston (LP) phase field model of martensite phase transformations (see Levitas and Preston (2002a,b) and Levitas et al. (2003)) to arbitrary pressure. Furthermore, we account for and track: the interface speed of the pressure-driven phase transformation, properties of critical nuclei, as well as nucleation at grain sites and on dislocations and homogeneous nucleation. The volume fraction evolution of each phase is described by employing KJMA (Kolmogorov, 1937; Johnson and Mehl, 1939; Avrami, 1939, 1940, 1941) kinetic theory. We then test our new model for iron under ramp loading conditions and compare our predictions for the α → ϵ iron phase transition to experimental data of Smith et al. (2013). In conclusion, more than one combination of material and model parameters (such as dislocation density and interface speed) led to good agreement of our simulations to the experimental data, thus highlighting the importance of having accurate microstructure data for the sample under consideration.

36 MATERIALS SCIENCE↗

Phase-field simulations to inform nuclear fuel performance modeling

Software tools to simulate nuclear fuel performance at the engineering scale, such as Idaho National Laboratory (INL)’s BISON code, are increasingly relied upon in regulatory and economic decision-making. However, accurate results from these tools depends on the availability of materials parameters that are used as input. In recent years, atomistic and mesoscale simulation methods have emerged as a cost-effective, expedient means to obtain such input parameters. Phase-field simulations using INL’s Marmot application have been used to obtain microstructure-level parameters and to improve material models for fuel performance modeling using BISON. In this talk, recent examples of this process are given, including applications in UO2, U3Si2, and UZr fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multimodal Nanoscale Tomographic Imaging for Battery Electrodes

Accurate representations of the 3D structure within a lithium-ion battery are key to understanding performance limitations. However, obtaining exact reconstructions of electrodes, where the active particles, the carbon black and polymeric binder domain, and the pore space are visualized is challenging. In this work, it is shown that multimodal imaging can be used to overcome this challenge. High-resolution ptychographic X-ray computed tomography are combined with lower resolution but higher contrast transmission X-ray tomographic microscopy to obtain 3D reconstructions of pristine and cycled graphite-silicon composite electrodes. This cross-correlation enables quantitative analysis of the surface of active particles, including the heterogeneity of carbon-black and binder domain and solid-electrolyte interphase coverage. Capturing the active particles as well as the carbon black-binder domain allows using these segmented structures for electrochemical simulations to highlight the influence of the particle embedding on local state of charge heterogeneities.

25 ENERGY STORAGE↗

Modeling Microstructural Effects on Heterogeneous Temperature Fields within Polycrystalline Explosives

The paper addresses the role of crystal anisotropy on the evolution of heterogeneous temperature fields in plastic-bonded explosives (PBXs) under conditions of weak shock. The modeling approach is based on simplified idealizations of PBX microstructure including RDX grains and estane binder regions subjected to velocity boundary conditions representative of impact conditions in situ. The constitutive description of the microstructure constituents includes a dislocation-based, anisotropic, single crystal plasticity model for the explosive grains and a linear viscoelastic model to represent the estane binder. Large suites of simulations were used to systematically study the correlation in heating of the grains with local wave dynamics, crystal orientation, and the microstructural neighborhood. These correlations were identified through the selection of characteristics of crystal anisotropy including oriented wave speeds and Taylor factor. A key observation is that the wave dispersion within a certain distance from the impact interface plays a dominant role in the temperature field. Beyond this distance, individual crystal orientations play a more dominant role, but cannot entirely account for the observed heterogeneity without consideration of the local microstructure neighborhood.

42 ENGINEERING↗

Chemical vapor infiltration of additively manufactured preforms: Pore-resolved simulations and experimental validation

The densification of additively manufactured porous preforms by chemical vapor infiltration (CVI) is studied using pore-resolved simulations and experiments. Experimentally, 3D printed silicon carbide (SiC) preforms are subject to CVI synthesis using methyltrichlorosilane (MTS) precursor to obtain high purity SiC/SiC composites. Optical images of the cross sections of the processed preforms are analyzed to obtain the spatial porosity distribution. The numerical method is based on a level set formulation to capture the spatial distribution and time evolution of the pore scale microstructural characteristics. The coupled transport and kinetic effects are represented using a dimensionless Thiele modulus. Simulations are initialized using representative synthetic preform geometries comprising of packed particles based on the size distribution of the powder used for 3D printing. The simulation results are validated against the experimental observations in terms of total density and the distribution of residual porosity. The densification characteristics, porosity classification, concentration profiles, and structure functions are analyzed as functions of processing temperature and Thiele modulus.

36 MATERIALS SCIENCE↗

Identifying the microstructural features associated with void nucleation during elevated–temperature deformation of copper

The microstructural-scale mechanisms that produce cracks in metals during deformation at elevated temperatures are relevant to applications that involve thermal exposure. Prior studies of cavitation during high-temperature deformation, for example, creep, suffered from an inability to directly observe the microstructural evolution that occurs during deformation and leads to void nucleation. Here the current study takes advantage of modern high-speed electron backscatter diffraction (EBSD) detectors to observe cavitation in oxygen-free, high-conductivity copper in situ during deformation at 300°C. Most voids formed at the triple junction between a twin boundary and a high-angle grain boundary (HAGB). This finding does not contradict previous studies that suggested that twins are resistant to cracking—it reveals that cracks in HAGBs originate at twin/HAGB triple junctions and that cracks preferentially grow along HAGBs rather than the accompanying twins. Atomistic simulations explored the origins of this observation and suggest that twin/HAGB triple junctions are microstructural weak points.

36 MATERIALS SCIENCE↗