SEARCH · Search NASA

Results for “Microstructural dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Defect Production and Microstructural Feature Impact for Radiation Damage in Additively Manufactured 316 Stainless Steel

This milestone presents multi-scale modeling research results for additively manufactured 316 stainless steel. A combination of phase field, cluster dynamics, molecular dynamics, and density functional theory with machine learning is used, allowing for predictions of radiation-driven microstructural evolution in additively manufactured 316 stainless steel over a range of temperatures, damage rates, neutron spectra, and microstructures, and supporting the development of combined ion and neutron irradiation for material qualification. Informed by ion irradiation and neutron irradiation results across the Advanced Materials and Manufacturing Technologies program, we investigate the unique aspects of radiation-driven microstructure evolution in additively manufactured 316 stainless steel. In particular, we focus on understanding the impact of carbon concentration (varying, for example, between the 316L and 316H standards) on void formation; radiation-induced segregation at dislocation cells and grain boundaries; and the evolution of dislocation loops and network dislocation populations. We find that the unique characteristics of the additively manufactured microstructures must be accounted for in understanding the evolution of dislocation populations under thermal and irradiation conditions, such as the variation in sink strengths arising due to the variation in dislocation density. We also find that the radiation-induced segregation of Cr and Ni to grain boundaries and cell walls differs due to the differences in their defect sink biases. We also find that increasing the Ni content can slow vacancy diffusion, which may provide a mechanism for the observed reduction in transient swelling rate for austenitic Fe-Cr-Ni alloys with increasing Ni content. In addition, ion irradiations have shown that increasing carbon content in 316 SS results in a larger population of smaller voids, suggesting reduced vacancy diffusion. Our results show that the carbon content of additively manufactured 316 SS has a significant impact on the migration rate of defect clusters. The presence of carbon atoms results in carbon-vacancy trapping, significantly reducing the diffusion rate of vacancies. Carbon atoms may also be trapped near the surface of a void, which may reduce void growth by trapping vacancies that diffuse toward the void.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Shock propagation in aerogel and TPP foams for inertial fusion energy target design

Achieving practical inertial fusion energy (IFE) requires the development of target designs with well-characterized microstructure and compression response. We measured shock dynamics in low-density (17.5–500 mg/cm 3 ) aerogel and two-photon polymerization (TPP) foams using x-ray phase contrast imaging (XPCI) methods and the Velocity Interferometer System for Any Reflector. By analyzing shock front evolution, we examined how target type and density influence shock propagation and energy dissipation. Talbot-XPCI shows that aerogels support a smooth, bowed shock front due to their homogeneous nanometer-scale pore network. In contrast, TPP foams exhibit irregular, stepwise propagation driven by interactions with their periodic micrometer-scale lattice. Shock velocity follows a power-law relation: aerogels deviate from classical ρ −1/2 scaling due to pore-collapse dissipation, while TPP foams follow the trend with larger uncertainties from density variations. Comparisons with xRAGE simulations reveal systematic underestimation of shock speeds. These results provide the first experimental constraints on shock propagation in TPP foams over a wide density range and highlight the influence of internal structure on anisotropic shock behavior. Our findings support improved benchmarking of EOS and hydrodynamic models and inform the design of foam architectures that promote implosion symmetry in IFE capsules.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Development of a MOOSE-based crystal plasticity model with irradiation defect evolution for irradiation creep in 316

Irradiation creep and irradiation swelling are two of the lifetime limiting factors for structural materials in nuclear reactors. These mechanical effects are driven by irradiation defect evolution and the interaction of those defects with dislocations in the microstructure. We present here a coupled cluster dynamics – crystal plasticity approach to model irradiation swelling and creep behavior in 316 SS. The time-dependent evolution of irradiation defects is calculated with a cluster dynamics approach and passed to the crystal plasticity model to compute the dislocation evolution. We show the impact of the irradiation defect evolution on the stress state in the material, which drives inelastic deformation through dislocation motion. The inelastic deformation in the 316 SS is dependent on the dose rate, where the inelastic deformation driven by the early-stage irradiation defect evolution determines the mechanical behavior of the 316 SS.

316 Stainless Steels↗

DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Direct Ab Initio Simulation of the Synthesis of BaZrO 3 and the Microstructure Impacts on Proton Transport

Controlling and predicting the processing-structure-performance relationship in functional materials is a grand challenge in materials science, with important implications for a wide range of emerging applications; a high fidelity understanding of the performance impact of microstructures formed under synthesis conditions is required to develop advanced materials, such as solid-state fuel cells and electrolyzers. Using the ceramic BaZrO 3 as a case study, we directly simulate the synthesis and investigate how proton transport is dictated by microstructures. We develop a framework that couples density functional theory (DFT), machine-learning interatomic potential (MLIP) driven molecular dynamics, and grand canonical Monte Carlo to perform large-scale, microstructure-resolved, atomistic simulations of proton transport in experimentally representative polycrystalline structures. Our fully ab initio approach, using a MLIP as a proxy for DFT, allows us to quantify the competition between two distinct diffusion mechanisms: one associated with grain-boundary regions and another within grains. When the impacts of grain boundaries are taken into account, proton transport exhibits substantial deviation from the bulk oxide limit. This addresses long-standing discrepancies between theory and experiments. Our integrated approach provides atomistic insight into microstructure-dependent proton pathways in BaZrO 3 and establishes a general protocol for predicting processing-structure-performance relationships.

organic↗

DOE OSTI2024

Microstructure Clones

Microstructure drives component behavior. Contemporary crystal plasticity studies compare strain measurements of polycrystal specimens to models. Because each specimen is unique, it is impossible to know which differences are significant. In this project, we invented microstructure clones and explored their use in understanding crystal plasticity. Microstructure clones are specimens with nearly identical microstructures, which allows for multiple destructive tests of a microstructure, insight into how a specimen will deform, variability quantification, and the ability to measure the effects of microstructural changes. Several sets of microstructure clones, pure nickel tensile bars, were tested. The techniques of digital image correlation, crystal plasticity finite element analysis, high resolution electron backscatter diffraction, transmission electron microscopy, and dislocation dynamics were used to understand the structural behavior of these microstructures. This work reshapes the fields of crystal plasticity and structure-property relationships by providing a technique to control for specific variables, quantify microstructural stochasticity, and replicate experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Correlated Ion Transport Governed by Dynamic Local Structure in High Concentration and Localized High Concentration Electrolytes

Understanding the dynamics of cluster formation and network percolation provides the mechanistic link between microscopic solvation structure and transport in concentrated electrolytes, including localized high-concentration electrolytes (LHCEs). Although recent studies have shown LHCEs to form micelle-like aggregates at specific compositions, a quantitative understanding of how solvation structures and transport properties depend on salt–solvent–diluent ratios remains limited. Here, we integrate molecular dynamics with Onsager transport analyses to chart the evolution of solvation microstructure and associated ionic transport in LiFSI/DMC/TTE, achieving good agreement with experimental conductivities across composition. We show that micelle formation arises from a dynamic instability of the cation–anion network, and that network percolation is the primary determinant of conductivity. Two compositional thresholds emerge: A critical network concentration (CNC) and a critical micelle concentration (CMC) that delineate transitions from extended percolating networks to micelle-like clusters and then to fragments. These structural transitions rationalize the nonintuitive conductivity decrease at intermediate dilution despite monotonically decreasing viscosity and provide composition-level design rules for LHCEs.

Mohanakrishnan, Rohith Srinivaas [University of Ca↗

DOE OSTI2025

Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations

Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale molecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.

X-ray imaging↗

DOE OSTI2025

Bridging the Gap in Carbon Free Iron Making: How Hydrogen Affects the Reduction of Iron Ore between 900 and 1590 °C

Hydrogen-based reduction of iron ore for iron and steel production has emerged as a promising alternative to coal and natural gas. Unlike other hydrogen-based iron ore reduction studies, this research focuses on a wide temperature range across 900–1590 °C, encompassing reduction in solid, mixed, and liquid (slag) phases. For a 20 min exposure to hydrogen, the reduction degree increased monotonically from ∼35% at 900 °C to >90% at 1550 °C, except between 1100 °C and 1400 °C, where it stagnated ∼60%. This experimental work challenges the widely accepted notion that higher temperatures enhance the reduction process. Instead, it reveals an overlooked kinetic bottleneck, suggesting complex thermodynamic and mass transfer limitations influenced by phase transformations, diffusion barriers, and microstructural changes. Density functional theory-based molecular dynamics simulations indicate that oxygen diffusivity in BCC iron is 3.88 × 10 –5 cm 2 /s which is ∼5–10 times higher than that in FCC iron. This study reports that this stagnant reduction degree in the mixed solid–liquid phase is due to competition of multiple mechanisms, such as surface- and bulk-diffusion, pore collapse mechanisms, and crystallographic transitions.

08 HYDROGEN↗

DOE OSTI2025

Evolution of dislocations during the rapid solidification in additive manufacturing

Materials processed by fusion-based additive manufacturing (AM) typically exhibit relatively high dislocation densities, along with cellular structures and elemental segregation. This representative structural feature significantly influences material performance; however, post-mortem microstructure characterizations of AM materials cannot capture the dynamic evolution of dislocations during the manufacturing process, thereby offering limited mechanism-based guidance for further advancing AM techniques and facilitating the qualification and certification of AM products. In this study, we conduct operando high-energy synchrotron X-ray diffraction experiments on wire-laser directed energy deposition of 316 L stainless steel. Through a unique configuration, our operando synchrotron experiments semi-quantitatively probe the dislocation density in solid phases and their dynamic changes during solidification and subsequent cooling. By integrating this advanced synchrotron technique with multi-physics simulation, in-situ neutron diffraction, and multi-scale electron microscopy characterization, our mechanistic study aims to elucidate the effects of rapid cooling and subsequent thermal cycling on the dislocation generation and evolution.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Effect of Thermal Oxidation on the Structure, Surface Texturing, and Microstructure Evolution in Nanocrystalline Ga-O-N Films

An extensive examination of the nanoscale, crystallographic growth dynamics of the system, which is impacted by the thermal energy given to the GaN, is carried out to derive a deeper understanding of the growth kinetics, morphology and microstructure evolution, chemical bonding, and optical properties of Ga-O-N films. Thermal annealing of GaN films is performed in the temperature range of 900–1200 °C. Crystal structure, phase formation, chemical composition, surface morphology, and microstructure evolution of Ga-O-N films are investigated as a function of temperature. Increasing temperature induces surface oxidation, which results in the formation of stable β-Ga2O3 phase in the GaN matrix, where the overall film composition evolves from nitride (GaN) to oxynitride (Ga-O-N). While GaN surfaces are smooth, planar, and featureless, oxidation induced granular-to-rod shaped morphology evolution is seen with increasing temperature to 1200 °C. The considerable texturing and stability of the nanocrystalline Ga-O-N on Si substrates can be attributed to the surface and interface driven modification because of thermal treatment. Corroborating with structure and chemical changes, Raman spectroscopic analyses also indicate that the chemical bonding evolution progresses from fully Ga-N bonds to Ga-O-N. While the GaN oxidation process starts with the formation of β-Ga 2 O 3 at an annealing temperature of 1000 °C, higher annealing temperatures induce structural distortion with the potential formation of Ga-O-N bonds. The structure-phase-chemical composition correlation, which will be useful for nanocrystalline materials for selective optoelectronic applications, is established in Ga-O-N films made by thermal treatment of GaN.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Rethinking materials simulations: Blending direct numerical simulations with neural operators

Abstract Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Low Field NMR Relaxometry Characterization of Water Adsorption in Corn Stover Anatomical Fractions

Low magnetic field 2 MHz NMR relaxometry is applied to measure water adsorption processes of intact corn stover anatomical fractions. Comparison with high magnetic field 250 MHz relaxometry of milled corn stover fractions as a function of water activity, i.e., relative humidity, provides insight into the microstructural changes of the biomass and varying water molecular dynamics during adsorption. As a result, the data presented establish a basis for low field portable NMR of biomass in situ in field and processing environments.

09 BIOMASS FUELS↗

DOE OSTI2025

Assessment of Microstructure Prediction Capabilities for Powder Bed Fusion Stainless Steel 316

The Advanced Materials and Manufacturing Technologies program aims to accelerate the development, qualification, demonstration, and deployment of advanced materials and manufacturing technologies to enable reliable and economical nuclear energy. However, the characteristic process-structure-property relationships of additive manufacturing (AM) materials pose challenges for the qualification and certification of AM nuclear components. In particular, component-scale variations in microstructure and properties can be driven by localized changes in melt pool dynamics due to how process parameters interact with different part geometries. Computational modeling tools can play a crucial role in predicting and controlling this variability. This report presents final results on process modeling tools designed to predict microstructure variability in additively manufactured stainless steel 316 parts. It details the software packages and physical modeling approaches employed to simulate an AM component within an automated process modeling workflow. Results are demonstrated through comparisons between predicted microstructures and experimental measurements across various representative processing conditions. The report concludes by discussing identified challenges and future opportunities for connecting the developed simulation workflow with mechanics simulations for prediction of part performance.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Experimental and Theoretical Thermokinetic Analysis of the Na 3 VO 4 –CO 2 Reaction System under Different Physicochemical Conditions

Sodium vanadate (Na 3 VO 4 ) was synthesized, structural and microstructurally characterized as well as tested as possible carbon dioxide (CO 2 ) captor through thermodynamic calculations in addition to thermogravimetric dynamic and isothermal analyses. Ceramic structural characterization evidenced the formation of the Na 3 VO 4 crystal phase, while microstructural features evidenced dense agglomerates with a poor specific surface area, although synthesis temperature (600 °C) was not as high as that for other sodium ceramics. Na 3 VO 4 was investigated for the CO 2 capture process through dynamic and isothermal experiments, in the presence or absence of oxygen. All these experiments showed that Na 3 VO 4 possesses interesting CO 2 capture properties in a specific temperature range (520 and 600 °C). Isothermal products characterization allowed to elucidate the reaction pathway, implying the Na 4 V 2 O 7 formation, as possible reactive intermediate. In fact, reaction path evolution was supported with theoretical thermodynamic data. Moreover, the kinetic analysis, performed using the Avrami-Erofeev mathematical model and subsequent thermodynamic data obtention, probed the positive influence of the temperature, in both CO 2 superficial and bulk chemical sorption processes, while oxygen addition did not enhance the superficial process. Additionally, as it would be expected, kinetics were enhanced as a function of the CO 2 concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Multiscale Modeling Framework for Lithium Nucleation in 3D Porous Carbon Anodes

Porous carbon scaffolds offer a promising route for mitigating non-uniform lithium (Li) plating to enhance the safety and longevity of Li metal batteries. However, the influence of microstructural morphology on Li nucleation is not well understood. Here, we present a multiscale modeling framework to investigate how the porous microstructure of carbon materials affects Li nucleation behavior. Ab initio molecular dynamics simulations quantify the nucleation energy barriers of Li on graphene as a function of Li content, surface curvature, and applied potential, providing key parameters for a classical nucleation theory (CNT) model. From macroscale half-cell simulations, we obtained Li concentration and electrical potential profiles to define boundary conditions for mesoscopic simulations. At the mesoscale, three distinct synthetic 3D microstructures with different porosities and characteristic feature sizes are generated to resolve local distributions of Li flux, current density, and mechanical stress. These outputs are integrated into the CNT model to map spatial variation in nucleation rates. Our findings reveal trade-offs between suppressing nucleation rates and achieving spatial uniformity, offering design guidelines for optimizing porous carbon anodes to balance nucleation control and mechanical integrity.

Materials science↗

DOE OSTI2024

High plasticity in refractory composite fabrication by ultrasonic additive manufacturing

Refractory metal composites are desirable for use in extreme environments that require materials with high specific strengths and resilience to external environments such as that found in nuclear and aerospace. However, due to the high melting temperature of refractories, liquid state joining processes such as welding remains difficult. Ultrasonic additive manufacturing (UAM) provides a potential route for processing refractory composites because it is a solid-state (i.e., no-melting) process and allows for intermittent machining operations to be performed between welds. Further, to demonstrate refractory composite fabrication, this study utilized UAM to machine a cavity to locate and sequester a Mo foil in a Zircaloy-4 (Zry-4) baseplate and additively build over the top with Zry-4 foils, thereby embedding the Mo in Zry-4 matrix. Significant deformation of the Zry-4 microstructure was observed: this deformation caused adiabatic heating and subsequent dynamic recrystallization through the transformation from α→β and then back to α as the material cooled. Flexural testing of the Zr–Mo composite revealed no delamination or failure, but the strength was not as expected, falling lower than a cold-worked Zry-4 sample. Finite element analysis supported that some bonding must have existed between the Zry-4 and Mo. There was indeed an interdiffusion zone at the Zry-4 foil–Mo foil interface, observing a metastable body-centered cubic β-Zr lathe. It was determined that sufficient strain energy was present to encourage the nucleation of the β-Zr grain along α-Zr grains. Future work is warranted to investigate UAM for refractory composite fabrication.

36 MATERIALS SCIENCE↗