Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microstructural evolution”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

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↗
DOE OSTI2026

Scaling kinetic Monte-Carlo simulations of grain growth with combined convolutional and graph neural networks

Graph neural networks (GNN) have emerged as a promising machine learning method for microstructure simulations such as grain growth. However, accurate modeling of realistic grain boundary networks requires large simulation cells, which GNN has difficulty scaling up to. To alleviate the computational costs and memory footprint of GNN, we suggest a hybrid architecture combining a convolutional neural network (CNN) based bijective autoencoder to compress the spatial dimensions, and a GNN that evolves the microstructure in the latent space of reduced spatial sizes. Our results demonstrate that the new design significantly reduces computational costs with using fewer message passing layer (from 12 down to 3) compared with GNN alone. The reduction in computational cost becomes more pronounced as the spatial size increases, indicating strong computational scalability. For the largest mesh evaluated (160 3 ), our method reduces memory usage and runtime in inference by 117× and 115×, respectively, compared with GNN-only baseline. More importantly, it shows higher accuracy and stronger spatiotemporal capability than the GNN-only baseline, especially in long-term testing. Such combination of scalability and accuracy is essential for simulating realistic material microstructures over extended time scales. The improvements can be attributed to the bijective autoencoder’s ability to compress information losslessly from spatial domain into a high dimensional feature space, thereby producing more expressive latent features for the GNN to learn from, while also contributing its own spatiotemporal modeling capability. Training data are generated from stochastic grain growth simulations, providing realistic variability for learning robust microstructure evolution. Comprehensive system validation confirms that the model is accurate, robust, and scalable.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Nonequilibrium defect-phase nanostructures stabilized by irradiation in undersaturated Ni-Si nanocrystalline alloy

Nanocrystalline thin films of the undersaturated alloy Ni-8.5 at% Si were subjected to 2 MeV Ti irradiation at temperatures ranging from 450˚C to 550˚C. Correlative microscopy combining transmission electron microscopy (TEM), scanning-TEM and atom probe tomography (APT revealed that large dose irradiation at 450˚C of samples with initial grain sizes below 100 nm stabilized a novel nanostructure which surprisingly contained three co-existing phases, the γ face-centered-cubic (FCC) matrix, γ' L12 ordered precipitates on intragranular dislocation loops and Ni 31 Si 12 precipitates at triple junctions (TJs). In contrast, irradiation at 550˚C and irradiation of larger grain-size samples at 450˚C only produced a γ-γ' two-phase coexistence. Analysis of the three-phase nanostructure and phase field simulations indicates that radiation-induced segregation is most pronounced at TJs, thus triggering the formation of Ni 31 Si 12 precipitates. These incoherent precipitates, in turn, are expected to stabilize the grain size under irradiation. The results are generalized using the concept of driven defect-phases. It is suggested that the stabilization of driven defect-phases may impart radiation resilience by providing localized relaxation modes to the microstructure evolution during and after temporary perturbations in irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Feasibility of in situ controlled heat treatment (ISHT) of Inconel 718 during electron beam melting additive manufacturing

A novel technique was developed to control the microstructure evolution in Alloy 718 processed using Electron Beam Melting (EBM). In situ solution treatment and aging of Alloy 718 was performed by heating the top surface of the build after build completion scanning an electron beam to act as a planar heat source during the cool down process. Results demonstrate that the measured hardness (478 ± 7 HV) of the material processed using in situ heat treatment similar to that of peak-aged Inconel 718. Large solidification grains and cracks formed, which are identified as the likely mechanism leading to failure of tensile tests of the in situ heat treatment material under loading. Despite poor tensile performance, the technique proposed was shown to successively age Alloy 718 (increase precipitate size and hardness) without removing the sample from the process chamber, which can reduce the number of process steps in producing a part. Lastly, tighter controls on processing temperature during layer melting to lower process temperature and selective heating during in situ heat treatment to reduce over-sintering are proposed as methods for improving the process.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Characterization of the Fe-Co-1.5V soft ferromagnetic alloy processed by Laser Engineered Net Shaping (LENS)

Processing of the low workability Fe-Co-1.5V (Hiperco® equivalent) alloy is demonstrated using the Laser Engineered Net Shaping (LENS) metals additive manufacturing technique. As an innovative and highly localized solidification process, LENS is shown to overcome workability issues that arise during conventional thermomechanical processing, enabling the production of bulk, near net-shape forms of the Fe-Co alloy. Bulk LENS structures appeared to be ductile with no significant macroscopic defects. Atomic ordering was evaluated and significantly reduced in as-built LENS specimens relative to an annealed condition, tailorable through selection of processing parameters. Fine equiaxed grain structures were observed in as-built specimens following solidification, which then evolved toward a highly heterogeneous bimodal grain structure after annealing. The microstructure evolution in Fe-Co is discussed in the context of classical solidification theory and selective grain boundary pinning processes. In conclusion, magnetic properties were also assessed and shown to fall within the extremes of conventionally processed Hiperco® alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Review of Computer-Aided Manufacturing (CAM) strategies for hybrid directed energy deposition

Hybrid additive manufacturing intertwines both additive and subtractive manufacturing layer by layer to digitally fabricate parts with complex geometries, improved surface finish and tight dimensional accuracies, the sum of which is difficult to obtain with any single process. Computer-Aided Manufacturing (CAM) software is required to orchestrate the machine toolpathing for both the deposition as well as the machining processes and is crucial for the successful fabrication of high quality structures. Additionally, CAM requires substantial operator input to account for challenging aspects of each fabricated structure. For example, deciding at which layer of deposition will the machining process continue to maintain access to complex cavities for finishing - internal features that would otherwise be unfinishable due to reach limitations or obstructions. Moreover, of the many commercially-available hybrid systems, each has a unique kinematic environment which can benefit from specific optimization of toolpath planning and substantial research has directly correlated toolpathing with the microstructure evolution, mechanical properties, porosity and residual stress state of the final fabricated part. This review explores the available strategies for CAM in the context of hybrid direct energy deposition, discusses the advantages and disadvantages of each and considers future CAM trends for this transformational digital manufacturing technology.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Improving fatigue performance of additively manufactured Ti-6Al-4V using sulfur-based self-terminating etching processes

The adoption of Additively Manufactured (AM) Ti-6Al-4V (Ti64) is continually increasing in the industry due to its benefits over the conventional manufacturing techniques. However, AM Ti64 parts are prone to poor surface finish, which significantly affects their fatigue performance. Therefore, AM Ti64 parts need post-processing to enhance surface finish and fatigue performance. Prior work with a sulfur-based self-terminating etching process for AM Ti64 demonstrated support removal and surface roughness reduction. To further this work, this manuscript details the impact of this sulfur-based post-processing technique on the mechanical properties and the fatigue life of AM Ti64 by characterizing the surface roughness, microstructural evolution, and fracture surfaces. Finally, this post-processing technique decreased the surface roughness parameter - S v by 55 % from 84 ± 11 µm to 38 ± 19 µm and increased the average fatigue life by 340 % from about 7000 cycles to 30,000 cycles.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Deformation behavior and phase transformation of nanotwinned Al/Ti multilayers

We report nanotwinned Al/Ti multilayers have exhibited size-dependent microstructure evolution and high strength. However, their deformation mechanisms are less well understood. In this work, we investigated the deformation mechanisms of nanotwinned Al/Ti multilayers with FCC/HCP layer interfaces by using in situ micropillar compression tests. Nanotwinned Al/Ti multilayers exhibit compressive strength up to 2.4 GPa and good work hardening capability. Post-compression TEM analyses reveal high-density stacking faults and the HCP-to-FCC phase transformations in Ti. Molecular dynamics simulations elucidate the mechanisms of deformation induced phase transformation in Ti and the influence of collective movement of partial dislocations on the deformability of Al/Ti multilayers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Properties of Nb3Sn films fabricated by magnetron sputtering from a single target

Superconducting Nb3Sn films were fabricated on sapphire and fine grain Nb substrates by magnetron sputtering from a single stoichiometric Nb3Sn target. The structural, morphological and superconducting properties of the films annealed for 24 h at temperatures of 800-1000 {\deg}C were investigated. The effect of the annealing time at 1000 {\deg}C was examined for 1, 12, and 24 h. The film properties were characterized by X-ray diffraction, scanning electron microscopy, atomic force microscopy, energy dispersive X-ray spectroscopy, and Raman spectroscopy. The DC superconducting properties of the films were characterized by a four-point probe measurement down to cryogenic temperatures. The RF surface resistance of films was measured over a temperature range of 6-23 K using a 7.4 GHz sapphire-loaded Nb cavity. As-deposited Nb3Sn films on sapphire had a superconducting critical temperature of 17.21 K, which improved to 17.83 K when the film was annealed at 800 {\deg}C for 24 h. For the films annealed at 1000 {\deg}C, the surface Sn content was reduced to ~11.3 % for an annealing time of 12 h and to ~4.1 % for an annealing time of 24 h. The Raman spectra of the films confirmed the microstructural evolution after annealing. The RF superconducting critical temperature of the as-deposited Nb3Sn films on Nb was 16.02 K, which increased to 17.44 K when the film was annealed at 800 {\deg}C for 24 h.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Integral boundary conditions in phase field models

Modeling the chemical, electric and thermal transport as well as phase transitions and the accompanying mesoscale microstructure evolution within a material in an electronic device setting involves the solution of partial differential equations often with integral boundary conditions. Employing the familiar Poisson equation describing the electric potential evolution in a material exhibiting insulator to metal transitions, we exploit a special property of such an integral boundary condition, and we properly formulate the variational problem and establish its well-posedness. Next, we compare our method with the commonly-used Lagrange multiplier method that can also handle such boundary conditions. Numerical experiments demonstrate that our new method achieves optimal convergence rate in contrast to the conventional Lagrange multiplier method. Furthermore, the linear system derived from our method is symmetric positive definite, and can be efficiently solved by Conjugate Gradient method with algebraic multigrid preconditioning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Order parameter re-mapping algorithm for 3D phase field model of grain growth using FEM

Phase field modeling (PFM) is a well-known technique for simulating microstructural evolution. To model grain growth using PFM, typically each grain is assigned a unique non-conserved order parameter and each order parameter field is evolved in time. Traditional approaches using a one-to-one mapping of grains to order parameters present a challenge when modeling large numbers of grains due to the computational expense of using many order parameters. This problem is exacerbated when using an implicit finite element method (FEM), as the global matrix size is proportional to the number of order parameters. While previous work has developed methods to reduce the number of required variables and thus computational complexity and run time, none of the existing approaches can be applied for an implicit FEM implementation of PFM. Here, we present a modular, dynamic, scalable reassignment algorithm suitable for use in such a system. Polycrystal modeling with grain growth and stress require careful tracking of each grain’s position and orientation which is lost when using a reduced order parameter set. In conclusion, the method presented in this paper maintains a unique ID for each grain even after reassignment, to allow the PFM to be tightly coupled to calculations of the stress throughout the polycrystal. Implementation details and comparative results of our approach are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Phase field benchmark problems for dendritic growth and linear elasticity

We present the second set of benchmark problems for phase field models that are being jointly developed by the Center for Hierarchical Materials Design (CHiMaD) and the National Institute of Standards and Technology (NIST) along with input from other members in the phase field community. As the integrated computational materials engineering (ICME) approach to materials design has gained traction, there is an increasing need for quantitative phase field results. New algorithms and numerical implementations increase computational capabilities, necessitating standard problems to evaluate their impact on simulated microstructure evolution as well as their computational performance. We propose one benchmark problem for solidifiication and dendritic growth in a single-component system, and one problem for linear elasticity via the shape evolution of an elastically constrained precipitate. We demonstrate the utility and sensitivity of the benchmark problems by comparing the results of 1) dendritic growth simulations performed with different time integrators and 2) elastically constrained precipitate simulations with different precipitate sizes, initial conditions, and elastic moduli. As a result, these numerical benchmark problems will provide a consistent basis for evaluating different algorithms, both existing and those to be developed in the future, for accuracy and computational efficiency when applied to simulate physics often incorporated in phase field models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Migration energy barriers and diffusion anisotropy of point defects on tungsten surfaces

We report tungsten is one of the most promising candidates for plasma-facing materials in future fusion devices, owing to its high performance under extreme irradiation conditions. However, irradiation-induced surface morphology varies significantly, depending on the irradiation type, fluence, and flux. Therefore, it is critical to examine the dynamics of point defects on tungsten (W) surfaces to understand how irradiation affects surface morphology. In this study, we employ the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method to search for potential migration paths of point defects on W (1 0 0), (1 1 0) and (1 1 1) surfaces. The first-principles calculations are then used to accurately determine the migration energy barriers. The obtained paths and barriers are incorporated into a kinetic Monte Carlo (KMC) model to determine trajectories and diffusivities, which are described by diffusion tensors to demonstrate their anisotropic features. Multiple diffusion mechanisms have been identified on different surfaces, with various anisotropy factors. Particularly, point defects on the W (1 1 0) surface have the highest diffusivities, with anisotropy factors independent of temperatures. In comparison, the anisotropy of the W (1 0 0) surface decreases as temperature increases, while the W (1 1 1) surface is isotropic for point defects. This study provides insights into defect transport properties on different surfaces, which are essential for understanding the early stages of irradiation-induced microstructural evolutions and surface morphology of tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The impact of anisotropic thermal expansion on the isothermal annealing of polycrystalline α -uranium

Although grain growth impacts microstructural evolution in a wide variety of materials systems, the effect of anisotropic thermal expansion on grain boundary mobility and texture evolution has not been widely studied. Anisotropic thermal expansion occurs in multiple non-cubic metals, and the thermomechanical processing behavior of these materials can be better understood with further study into the impact of thermal expansion on grain boundary mobility and texture evolution. In this work, we develop a mesoscale phase field model of grain growth that includes the effect of anisotropic thermal expansion, which is applied to study polycrystalline α-uranium, a highly anisotropic metal. Three-dimensional simulations on polycrystalline a-uranium with and without thermal expansion eigenstrains are performed to study the grain boundary mobility and texture evolution as a function of temperature. A strain-free temperature of 933 K is selected, and the system is studied at within the range of 873 – 933 K at intervals of ten degrees, resulting in increasing thermal eigenstrain with decreasing temperature. We also estimate a grain boundary mobility prefactor and activation energy based on existing experimental data of isothermal annealing of α-uranium. The grain boundary mobility is found to display significant deviation from Arrhenius behavior with the inclusion of thermal expansion eigenstrain as the amount of thermal eigenstrain (and thus elastic strain energy within the system) increases. This result explains an experimentally observed grain boundary mobility deviation from Arrhenius behavior. Furthermore, the texture evolution is affected, such that the grain orientations become less random with increasing thermal eigenstrain, which could explain experimentally observed texture behavior. These results indicate that the effect of thermal expansion should be considered when predicting the thermomechanical processing behavior of α-uranium and other materials with anisotropic thermal expansion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Massively parallel phase-field simulations targeting exascale

The interface thickness in the phase-field (PF) method limits its simulation scales. Consequently, large-scale PF simulations become prohibitively expensive for resolving the extremely fine microstructures that typically form during rapid solidification processing. This challenge is significant in predicting microstructure evolution in metal additive manufacturing and has been identified by the United States Department of Energy’s Exascale Computing Project. Here, to address this, we develop a multi-GPU and MPI-based massively parallel simulation code, utilizing state-of-the-art algorithms, software, and libraries, for large-scale three-dimensional (3D) PF simulations. We report the first GPU-parallel PF simulations on Frontier (currently the second TOP500 exascale cluster) and Summit machines, taking dendritic growth as an example problem. We evaluate the parallel performance of our implementation using scaling studies with more than 24 000 GPUs (among the largest known computations to date) and the acceleration performance using large-scale simulations of dendritic growth in 3D. Finally, massively parallel GPUs in these supercomputers enabled the first coupled multiscale simulations of laser melting and subsequent dendritic solidification on the scale of a full melt-pool, demonstrating the feasibility of performing PF simulations with a point total over 2 billion grid points within an acceptable time.

Exascale↗
DOE OSTI2025

Development of Predictive Model for Accurate Rupture Time from Multi-Axial Creep in Alloy 709 with Physics-Based Simulations

A physics-based model is developed to predict multiaxial creep behavior in Alloy 709 (A709), an advanced austenitic stainless steel intended for high-temperature applications such as Sodium Fast Reactors (SFRs). Compared to conventional stainless steels like 316H, A709 offers superior high-temperature performance; however, comprehensive data on its multiaxial creep response remain limited. To address this gap, a crystal plasticity finite element (CPFE) framework is used to simulate the deformation and failure mechanisms of A709 under multiaxial loading conditions. The model incorporates an extended Hu-Cocks dislocation creep formulation that accounts for precipitation effects, along with the Sham–Needleman model to capture grain boundary cavitation-driven failure. These advanced constitutive models enable a detailed understanding of the interplay between microstructural evolution and macroscopic creep response. Furthermore, the study evaluates the predictive accuracy of various effective stress measures in estimating creep rupture life, leveraging simulated multiaxial creep data. The findings provide critical insights into the applicability of different stress measures for engineering design and life prediction of A709 components operating under complex loading conditions. This work contributes to improving the reliability of high-temperature structural components by advancing predictive modeling capabilities for advanced austenitic steels.

Alloy 709↗