Data assimilation for microstructure evolution: Grain growth in kinetic Monte Carlo.
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
Explore the source record for details and available documents.
High fidelity simulations of grain formation in alloys are an indispensable tool for process-to-mechanical-properties characterization. Such simulations, however, can be computationally expensive as they require fine spatial and temporal discretizations. Their cost becomes an obstacle to parametric studies and ensemble runs and ultimately makes downstream tasks like optimal control and uncertainty quantification challenging. To enable such downstream tasks, we introduce GrainNN, an efficient and accurate reduced-order model for epitaxial grain growth in additive manufacturing conditions. GrainNN is a sequence-to-sequence long-short-term-memory (LSTM) deep neural network that evolves the dynamics of manually crafted features. Its innovations are (1) an attention mechanism with grain-microstructure-specific transformer architecture; and (2) an overlapping combination of several clones of the network to generalize to grain configurations that are different from those used for training. This design enables GrainNN to predict grain formation for unseen physical parameters, grain number, domain size and geometry. Furthermore, GrainNN not only reconstructs the quantities of interest but also can be pointwise accurate. In our numerical experiments, we use a polycrystalline phase field method to both generate the training data and assess GrainNN. For multiparametric, ensemble simulations with many grains, GrainNN can be orders of magnitude faster than phase field simulations, while delivering 5%–15% pointwise error. Additionally, this speedup includes the cost of the phase field simulations for generating training data.
Abstract SPPARKS is an open-source parallel simulation code for developing and running various kinds of on-lattice Monte Carlo models at the atomic or meso scales. It can be used to study the properties of solid-state materials as well as model their dynamic evolution during processing. The modular nature of the code allows new models and diagnostic computations to be added without modification to its core functionality, including its parallel algorithms. A variety of models for microstructural evolution (grain growth), solid-state diffusion, thin film deposition, and additive manufacturing (AM) processes are included in the code. SPPARKS can also be used to implement grid-based algorithms such as phase field or cellular automata models, to run either in tandem with a Monte Carlo method or independently. For very large systems such as AM applications, the Stitch I/O library is included, which enables only a small portion of a huge system to be resident in memory. In this paper we describe SPPARKS and its parallel algorithms and performance, explain how new Monte Carlo models can be added, and highlight a variety of applications which have been developed within the code.
This project addresses two critical and intertwined challenges in fusion energy, namely the shortage of a broadly trained scientific workforce and the lack of scalable manufacturing solutions for plasma-facing components (PFCs). Through a collaboration among Florida International University (FIU), Miami Dade College (MDC), and Purdue University, the project established structured, reproducible educational and research pathways that recruit and advance students from institutions historically outside the fusion energy enterprise, building the human capital that this field urgently needs. The project integrates the complementary research strengths of FIU and Purdue to investigate flash sintering as a transformative processing route for tungsten-based PFCs. Unlike conventional sintering approaches, flash sintering offers rapid densification at significantly reduced thermal budgets, making it a compelling candidate for fabricating complex tungsten geometries that must withstand extreme plasma-facing environments. Systematic experimental and modeling efforts will elucidate the fundamental mechanisms governing microstructure evolution, grain boundary chemistry, and thermomechanical response during flash sintering — knowledge that is presently lacking but essential for translating this technology into reliable manufacturing practice. The convergence of workforce development and cutting-edge manufacturing research positions this project to deliver measurable, durable impact: a pipeline of fusion-ready researchers cultivated through expanded institutional partnerships, and a validated materials processing framework that accelerates domestic readiness for next-generation fusion reactor construction.
In this study, prismatic geometries of GRCop-84 [Cu-8Cr-4Nb (at. pct)] were built with laser-powder bed fusion (L-PBF) process. The samples were sectioned parallel or perpendicular to the build direction and characterized in the as-built and after post-processing with a hot-isostatically pressing (HIP) treatment. The microstructure and phase evolutions were evaluated with optical microscopy, scanning electron microscopy (SEM), electron backscattered diffraction (EBSD), and high-temperature X-ray diffraction (HTXRD) up to 1223 K. The samples in the as-built conditions exhibited predominantly columnar epitaxial and misoriented Cu-FCC grains. The microstructure evolutions are discussed based on locations within the overall build geometry, the dynamics of small melt pool shape and sectioning effects. The above grain structure did not change significantly during post-process HIP treatment. The stability of this FCC grain structure is attributed to the formation of primary stable Cr 2 Nb (Laves phase) during L-PBF, even before the emergence of FCC grains from liquid. The stability of Cr 2 Nb in both as-built and HIPed samples were evaluated using high-temperature X-ray diffraction measurements and compared with that of gas-atomized powder. The significance of these results is discussed with reference to aerospace applications.
This project is aimed to develop a chromium (Cr) tolerant, highly active, and stable coating layer on the internal surfaces of the porous composite cathode from commercially available SOFCs. Such coating layer was developed using the additive manufacturing process of Atomic Layer Deposition (ALD) and has been applied on the cathode consisting of either an electronic conductor of LaxSr 1-x MnyO 3-δ (LSM) or mixed ionic and electronic conducting La x Sr 1-x Co y Fe 1-y O 3-δ (LSCF). PI's work has demonstrated that the internal surface of cathode from the commercial cells, can be further tailored using ALD coating to dramatically enhance the cell performance. For instance, ALD layer consisting heterostructured nano composite of nano-Pt and nano-(Mn 0.8 Co 0.2 ) 3 O 4 oxide on the internal surface of porous LSM/YSZ cathode from SOFCs, has resulted in the large reduction of the cell polarizations resistance by up to 53%, and enormous increase of power density over 370%. For the cells with LSCF/Sm 2 O 3 doped CeO 2 (SDC) cathode, the conformal layer of nano-composite consisting of superjacent CoOx and subjacent minimum amount of Pt nano-grains has resulted in the power density enhancement by 126% for the large scale industry tubular cells at 750°C, and both the performance enhancement and nanostructure of the ALD layer are stable over ~ 2000 h continuous operation performed at industry test station. In the meanwhile, those ALD coating layer developed by PI's work is also inherently Cr-tolerant, and could act as physical barrier for preventing Cr diffusion into the cathode backbone, so as to mitigate the Cr poisoning effect on the cathode. In this project, the impact of Cr on the performance of those ALD coated commercial cells has been evaluated. Based on evolution of the cell performance, the ALD coating layer chemistry and ALD coating layer thickness has been optimized to maximize the overall Cr tolerance, cell power density and cell longevity. Different ALD coating has been applied onto the internal surface of LSM/YSZ and LSCF/SDC backbone respectively. The architecture/scaffold structures on the internal surface of different cathode, designed by this project, was catalogued and analyzed using High Resolution Transmission Electron Microscopy (HRTEM), and cell power/durability performance are assured via comprehensive electrochemical performance testing in the industry operation relevant conditions. The impact of the electrochemical operation current density, the water humidity, the cell operation temperature, and cell operation duration on the Cr tolerance of ALD coated cells has been systematically investigated. There is completely different nanostructure degradation mechanisms between LSM and LSCF cells induced by Cr contamination. For the LSCF/SDC baseline cell, With the Cr source, there is no apparent Sr surface segregation phase even for the baseline cell operated for 3000 h at 750 °C. With the Cr source, there is significant amorphous (SrCr)Ox phase accumulated in the original pore region. For the commercial baseline cells, Cr contaminants on the LSM electrode severely impacted the entire cell's electrochemical performance and nanostructure degradation. Those degradations include (1). Peak power density loss of 64 % after 109 h of operation. The dramatic increase in Rp (2). They are cracking at LSM/SSZ interface, LSM grains. SSZ remains intact but with (CrMn)Ox. By contrast, ALD coating (MnCo)Ox/Pt dramatically improves the Cr resistance, as follows (1). ALD-coated cell with a power density is 280-380 % of the baseline cell, depending on the ALD layer thickness. (2). For a cell with a 20 nm thick ALD layer, there is a large performance enhancement (> 200 % power density) induced by ALD coating of Cr-tolerant Mn 0.8 Co 0.2 Ox. (3). For a cell with a 20 nm thick ALD layer, after 168 h at 750 °C power density of the ALD-coated cell is ~ 600% of that baseline cell upon operation with Cr contamination for 109 h. The ALD coating on the internal surface of cathode developed by this project integrated multi-functions. Those multi-functions include (1). Dramatically improving the cell power density for the commercial cells; (2). Dramatically improving contamination resistance of the cathode, for being an excellent protection coating layer sealing off Cr contamination. (3). Dramatically increasing the cell longevity by potentially preventing the microstructure evolution and grain coarsening of the cathode. Overall, this project will provide a simple solution to simultaneously enhance power density and increase the reliability, robustness, and endurance of commercial SOFCs, over the entire operating temperature range of 650-800 °C. For the inherently functional SOFC, the ALD coating of LSM based cathode mitigate the Cr-contamination. Power density of ALD-coated cell is ~ 600% of that baseline cell upon operation with Cr contamination. In addition to SOFCs, the novel on-demand design approach and creation of multifunctional heterogeneous architecture on the electrode surface presented in this work opens further research for their application in other types of fuel cells, batteries, and sensors for which electrochemical reactions on the surface are similarly critical.
We propose GrainGNN, a surrogate model for the evolution of polycrystalline grain structure under rapid solidification conditions in metal additive manufacturing. High fidelity simulations of solidification microstructures are typically performed using multicomponent partial differential equations (PDEs) with moving interfaces. The inherent randomness of the PDE initial conditions (grain seeds) necessitates ensemble simulations to predict microstructure statistics, e.g., grain size, aspect ratio, and crystallographic orientation. Here, currently such ensemble simulations are prohibitively expensive and surrogates are necessary.In GrainGNN, we use a dynamic graph to represent interface motion and topological changes due to grain coarsening. We use a reduced representation of the microstructure using hand-crafted features; we combine pattern finding and altering graph algorithms with two neural networks, a classifier (for topological changes) and a regressor (for interface motion). Both networks have an encoder-decoder architecture; the encoder has a multi-layer transformer long-short-term-memory architecture; the decoder is a single layer perceptron.We evaluate GrainGNN by comparing it to high-fidelity phase field simulations for in-distribution and out-of-distribution grain configurations for solidification under laser power bed fusion conditions. GrainGNN results in 80%–90% pointwise accuracy; and nearly identical distributions of scalar quantities of interest (QoI) between phase field and GrainGNN simulations compared using Kolmogorov-Smirnov test. GrainGNN's inference speedup (PyTorch on single x86 CPU) over a high-fidelity phase field simulation (CUDA on a single NVIDIA A100 GPU) is 150×–2000× for 100-initial grain problem. Further, using GrainGNN, we model the formation of 11,600 grains in 220 seconds on a single CPU core.
Friction stir processing (FSP) has been increasingly adopted for joining and processing materials in automotive, aerospace and industrial construction. During FSP a dynamic competition between high-speed shear deformation and deformation-induced heating brings about a complex competition between multiple dynamic microstructural evolution mechanisms making it difficult to predict the microstructural evolution pathway. Hence, improved understanding of microstructural evolution mechanisms during FSP can be beneficial for continued growth in the adoption of FSP for demanding applications of future. Towards this goal, this study uses a model binary Cu – 50 at.% Ni alloy to clarify the effect of single and double pass FSP on the microstructural evolution of a coarse grained and compositionally heterogeneous cast microstructure. High energy synchrotron X-ray diffraction, electron backscatter diffraction, and nanoindentation are used to clarify the microstructural evolution due to FSP. The process of compositional homogenization of as-cast segregations is studied by energy dispersive spectroscopy and atom probe tomography. Our results show that a single fast FSP pass at 30 mm.s -1 produces a 100 µm deep layer of submicrometric and hall-petch hardened CuNi grains. The initial cast compositional heterogeneities in a micrometric scale is rapidly transformed to nano-sized domains, mainly confined at grain boundaries. Double pass FSP increases the penetration depth of the processed layer and leads to a 2.9 times grain growth relative to single pass FSP. Grain fragmentation, discontinuous dynamic recrystallization, grain growth, and twinning mechanisms are discussed. In conclusion, these results highlight the value of FSP for ultrafast grain refinement and compositional homogenization of cast alloys.
Metal is most commonly used structural material in a wide range of spacecraft, and welding is the principal method for joining metal components into functional systems. However, conducting welding experiments under extreme environments—such as microgravity or vacuum conditions in space—is prohibitively expensive and experimentally challenging. To overcome these limitations, multi-physics computational welding models provide a cost-effective and versatile alternative. In this work, the authors have developed a coupled thermal (fluid) microstructure simulation framework to model metal welding under varying gravity conditions. The framework integrates a mixed-mode heat transfer formulation (conduction, convection, and radiation) with molten pool fluid dynamics, enabling accurate prediction of temperature fields and weld-pool geometry. A grain growth model is further incorporated to capture the spatial and temporal evolution of microstructure, including grain size distribution and morphological transitions during solidification. This approach provides detailed insight into molten pool evolution and grain-level microstructure development throughout the welding process. By explicitly parameterizing environmental conditions, the model supports extrapolation to off-Earth manufacturing scenarios such as welding on the lunar surface. Tantalum—chosen in this study due to its high melting point, oxidation resistance, and mechanical stability at elevated temperatures—serves as the material system for model demonstration. Beyond Tantalum, the integrated multi-physics framework offers broad applicability for predictive welding simulations of various structural and refractory metals or alloys used in extreme terrestrial or extraterrestrial environments.
The alpha phase of uranium exhibits the orthorhombic crystal structure, resulting in significant physical property anisotropy. These complex properties make the material challenging to work with in engineering settings but provide a rich arena to investigate the fundamental, multi-scale effects of anisotropy on physical behaviors such as grain growth and microstructure evolution under irradiation. We apply molecular dynamics and phase field modeling to study mass transport and grain growth in alpha-uranium. We characterize the mobilities of crystalline defects and the interfacial energy, and We find that the grain boundary energy is highly variable depending on the interface structure and that the relative rates of defect transport change with temperature. We also find that anisotropic thermal expansion has a significant impact on grain growth kinetics and texture development. Our results are used to explain experimental observations and to provide a basis for further hypotheses into the complex physical behavior of alpha-uranium.
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.
Laser powder bed fusion (L-PBF) has enabled the fabrication of geometrically complex metallic structures and components that are challenging to producing using conventional manufacturing approaches. The site-specific and far from equilibrium thermal conditions of L-PBF offer the potential to facilitate multi-length scale design of structure and properties across the atomic-through macro-levels. However, L-PBF systems face scalability challenges due to throughput constraints. Laser rotary powder bed fusion (L-RPBF) systems are being investigated as a solution to enhance the deposition rates compared to conventional L-PBF. Rotary systems also offer additional flexibility for controlling the time structure of melting through laser interleaving on alternating layers. Here, in this study, IN718 test samples were printed using single-laser or interleaved dual-laser configuration in a L-RPBF system to investigates the effect of process settings and melt-interleaving on as-fabricated microstructure. The microstructural evolution, such as grain size and crystallographic texture, was assessed by determining variations in the melt-pool shapes. Laser interleaving leads to a reduction in average grain size compared to single laser by ∼ 40 % at high power (400 W) and by ∼36 % at medium power (370 W). Results presented here identify key challenge for obtaining uniform microstructures and barriers for the broader adoption of high-deposition rate L-RPBF.
An experimental investigation of failure mechanisms in relation to the initial microstructure in high-purity titanium samples is presented. The initial microstructure of the specimens was measured in 3D using diffraction-contrast tomography (DCT). Ex-situ X-ray computed tomography (XCT) was used to visualize voids in specimens at various levels of deformation attained by interrupted tensile testing. Electron backscatter diffraction (EBSD) mapping was used to characterize the deformed microstructures per specimen. The investigation showed that the location of failures was unrelated to the formation of the first voids, but rather the nucleation and growth of cracks were from heavily deformed surface regions. The areas surrounding the failure-causing cracks had the highest propensity for void development of any area in the specimen. This facilitated crack propagation through the specimen via “linking up”. The failure of the specimens by microvoid coalescence was bolstered by “dimpled” fracture surfaces. Paint composed of gallium embedded in rubber was used to track the displacement of points along the loading direction of the specimens, to relate the location of failure to the material's initial microstructure. However, the correlation proved difficult, due to complex microstructural evolution involving voids, grain fragmentations, and profuse twinning. Nevertheless, it was shown that necking occurred near the largest grain in the gauge section of the specimens. This is likely caused by the decrease in strength of larger grains, causing strain to concentrate in their vicinity. Here, a comprehensive picture of the failure processes in high-purity titanium under axial loading is presented and discussed in the paper.
The welding process is a potential way of repairing a damaged metal component, especially cavity damage caused by a harsh environment like the lunar environment, which is characterized by large temperature differences and reduced gravity. The adjustment of welding parameter (e.g., hatch spacing) can improve production efficiency in the repair process. Seen from the microstructural level, hatch spacing sensitivity affects the metallic grain evolution and morphology in the welding process, which can further influence a repaired part’s mechanical properties; however, the study of hatch spacing’s effect on microstructure is challenging. Traditional experimental procedures are costly and time-consuming, and any change in hatch spacing value needs roll-back of experimental procedure. A modeling study can address the above challenges in experimental observation. In this research, a modeling approach based on the Kinetic Monte Carlo (KMC) Potts theory was used to simulate grain evolution and morphology with three hatch spacings. Through quantifying and analyzing the predicted grain morphologies, the effect of hatch spacing on microstructure in a welding-fabricated part was investigated. The predicted grain morphologies were validated with an EBSD image of welding microstructure, which has been published before. The primary grain morphologies were columnar grains with a small amount of fine equiaxed grains formed in the scanning path centerline. When increasing the hatch spacing, the columnar grains become larger and more lengthy, while the effect of hatch spacing on the equiaxed grains is not obvious.
The microstructural evolution at and near preexisting grain boundaries (GBs) and dislocations in materials under high radiation doses is still poorly understood. In this work, we use the creation relaxation algorithm (CRA) developed for atomistic modeling of high-dose irradiation in bulk materials to probe the athermal limit of saturation of GB and dislocation core regions under irradiation in fcc Ni. We find that, upon continuously subjecting a single dislocation or GB to Frenkel pair creation in the athermal limit, a local steady-state disordered defect structure is reached with excess properties that fluctuate around constant values. Case studies are given for a straight screw dislocation which elongates into a helix under irradiation and several types of low- and high-angle GBs, which exhibit coupled responses such as absorption of extrinsic dislocations, roughening, and migration. A positive correlation is found between the initial GB energy and the local steady-state GB energy under irradiation across a wide variety of GB types. Metastable GB structures with similar density in the defect core region but different initial configurations are found to converge to the same limiting structure under CRA. The mechanical responses of pristine and irradiated dislocations and GB structures are compared under an applied shear stress. Irradiated screw and edge dislocations are found to exhibit a hardening response, migrating at larger flow stresses than their pristine counterparts. Mobile GBs are found to exhibit softening or hardening responses depending on GB character. Although some GBs recover their initial pristine structures upon migration outside of the radiation zone, many GBs sustain different flow stresses corresponding to altered mobile core structures.
Not Available
Recent experimental studies suggest the use of spatially extended laser beam profiles as a strategy to control the melt pool during laser powder bed fusion (LPBF) additive manufacturing. However, linkages connecting laser beam profiles to thermal fields and resultant microstructures have not been established. Herein, we employ a coupled thermal transport-Monte Carlo model to predict the evolution of temperature fields and grain microstructures during LPBF using Gaussian, ring, and Bessel beam profiles. In addition, simulation results reveal that the ring-shaped beam yields lower temperatures compared with the Gaussian beam. Owing to the small melt pool size when using the Bessel beam, the grains are smaller in size and more equiaxed compared to those using the Gaussian and ring beams. Our approach provides future avenues to predict the impact of laser beam shaping on microstructure development during LPBF.