Search NASA⌕ Search

SEARCH · Search NASA

Results for “microstructure prediction”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Substitution-Mediated Calcination of Nickel-Based Cathodes: Decoupling Lithiation and Crystallization

Nickel-based layered cathodes such as LiNiO 2 offer high energy density for lithium-ion batteries, yet improvements in cycling performance and safety are required for practical use–often achieved through manganese and cobalt substitution as in LiNi 0.80 Mn 0.10 Co 0.10 O 2 (NMC811). However, how such substitution impacts calcination, the key process that governs lithiation, structural ordering, crystallization, and ultimately the resulting material properties, remains unclear. Here, in this study, we investigate substitution-mediated calcination dynamics in NMC811 compared to LiNiO 2 using multiscale-correlated in situ spectroscopy and atomistic-to-mesoscale modeling. While both systems progress through the same sequence of intermediates toward the thermodynamically favored layered phase, NMC811 exhibits an earlier onset of layering, concurrent with hydroxide decomposition followed by sluggish crystallization. Modeling reveals that Mn and Co lower the energy barrier for lithium incorporation and ordering but increase the penalty for interlayer gliding, thereby slowing crystal growth at elevated temperatures. This substitution-mediated decoupling of lithiation and crystallization explains the fine-grained microstructure observed in NMC811 versus coarsened particles in LiNiO 2 and establishes a mechanistic framework for predictive microstructure engineering of Ni-based cathodes.

Calcination↗
DOE OSTI2025

Mesoscale Modeling for Restructuring and Fragmentation in High Burnup UO 2

This report summarizes the mesoscale modeling conducted in fiscal year 2025 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, focusing on the microstructural evolution and restructuring in high burnup UO 2 nuclear fuel and its impact on fuel fragmentation. We developed a pioneering phase-field model to simulate restructuring behavior across different regions of high burnup fuel, including the dark zone and rim region. A grand-potential-based phase-field model is employed to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of temperature and burnup rate were studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to initiate around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. The restructuring occurs faster with higher burnup rate. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. The predicted microstructures are consistent with experimental observations of the restructured regions. Based on the mesoscale simulations, a mechanistic model for restructuring and grain size evolution was implemented in BISON. Thus, this work provides a first-of-its-kind restructuring model for different regions of high-burnup fuel to BISON. It was shown that the model predicts appropriate grain size evolution along fuel radius as those observed in experiments. This model enables BISON to predict the effect of restructuring on the fission gas release. Finally, the phase-field fracture simulation with dark zone specific microstructures were presented to provide the fragmentation criteria for the dark zone. This work introduces a first-of-its-kind restructuring model for high burnup fuel, significantly enhancing BISON's predictive capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Phase-field modeling for restructuring in the dark zone of high burnup UO 2

This report summarizes the mesoscale modeling work performed in fiscal year 2024 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstructural evolution and restructuring observed in the dark regions of high burnup UO 2 nuclear fuel. This is the first attempt to realistically simulate the restructuring behavior observed in different region of a high burnup fuel. We employ a grand-potential based phase-field model to concurrently evaluate the formation of subgrains and growth of fission bubbles within the fuel. A energy-based subgrain formation criteria is introduced to simulate the restructuring process. Effect of different initial conditions and different modeling parameters are studies systematically to capture how each of these parameters influence the characteristics of the restructured fuel. It is observed that the subgrain formation begins around existing fission gas bubbles and then proceeds towards triple junctions, grain boundaries and grain interiors. It is demonstrated that restructuring is influenced by a combination of initial dislocation densities, subgrain formation rate, and temperature. Rate of restructuring increases with increase in fuel temperature. A restructuring bias is observed within the microstructure due to variation in defect accumulation among different grains. Furthermore, bubble sizes and distribution does not have a significant effect on rate of restructuring. The predicted microstructures resembles the characteristics of the restructured regions as observed in experiments. Finally, a correlation is presented that demonstrates the evolution of the restructuring volume fraction as a function of local effective burnup. This work provides a first of its kind restructuring model for darkzone that can be used by BISON for performance prediction of high burnup UO 2 fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Symposium MT02: Statistical Mechanics-Based Computational Tools for the Study of Phase Transformation in Complex Materials (Final Report)

Symposium MT02 brought together a diverse and interdisciplinary community of scientists specializing in Statistical Mechanics-based computational modeling to investigate phase transformations in materials exhibiting complex disordered structures. As the demand for materials with extreme performance metrics grows—from aerospace components to next-generation optical fibers—the ability to predict microstructural evolution under non-equilibrium conditions has become paramount. The primary goal of this symposium was to identify, evaluate, and discuss advanced computational tools capable of designing precise manufacturing conditions to tailor material properties efficiently. By fostering a dialogue between computational theorists and experimentalists, the symposium sought to establish new protocols for predicting how processing history—such as cooling rates or strain paths—dictates the final microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Comparison of hardening and microstructures of ferritic/martensitic steels irradiated with fast neutrons and dual ions

Ferritic/martensitic steels T91 and HT9 were irradiated with neutrons (BOR-60 reactor) and dual ions (9 MeV Fe 3+ and 3.42 MeV energy degraded He 2+ ) from 369 to 520 °C and damage levels of 16.6 to 72 dpa to quantify the possibility of using ion irradiation to simulate neutron irradiation in terms of microstructures and mechanical properties. Nanoindentation testing was performed to obtain the bulk equivalent hardness of the dual-ion irradiated samples. For the neutron irradiated samples, both nanoindentation and Vickers hardness testing were conducted. Transmission Electron Microscopy (TEM) characterizations of the cavities, dislocation loops and precipitates were conducted to account for the strengthening contribution of each microstructure element. The good agreement between the microstructure-predicted (dispersed barrier hardening) and measured strength of the irradiated specimens demonstrated the accuracy of the strengthening model and the nanoindentation tests. Furthermore, the comparison of mechanical property and microstructure changes in ion and neutron irradiated structural materials indicated that ion irradiation replicated many neutron irradiation features. However, a single 70 °C temperature shift is insufficient to match all complex microstructures of neutron vs. ion irradiation over the irradiation temperature range of 369–520 °C.

36 MATERIALS SCIENCE↗
NASA NTRS2022

Analysis of Base Metal Microstructures and Mechanical Properties of a Single-Piece, Spin-Formed Forward Pressure Vessel Bulkhead

In 2012, a pathfinder forward pressure vessel bulkhead (FPVBH) for the Orion Multi-Purpose Crew Vehicle (MPCV) was fabricated using aluminum (Al) alloy 2219. The demonstration article was offered as a replacement for the baseline configuration, which at that time was a multi-piece welded construction using aluminum-lithium (Al-Li) alloy 2195. In 2014, the NASA Engineering and Safety Center (NESC) funded a Phase I effort to explore spin-forming as a manufacturing method to produce a complex-shaped, single-piece FPVBH for the crew module (CM) [1]. The primary focus of the NESC effort was to expand spin-forming technology to the fabrication of an Al-Li 2195 FPVBH. As a result, limited testing and analysis was performed on the original Al 2219 article. In this study, a sufficiently large plate of Al 2219 was not available for the fabrication of a FPVBH. Therefore, the forming blank comprised two plates (from different lots) butt welded together using a single friction stir weld. The welded forming blank was then spin formed into the FPVBH configuration and heat treated to the T62 temper. Mechanical property testing of the spin-formed 2219-T62 FPVBH revealed that the tensile properties were comparable to the Metallic Material Properties Development and Standardization (MMPDS) A-basis allowable for Al 2219-T6 wrought products. However, post-spin forming application of the standard heat treatment produced large, recrystallized grains in the material. The response varied within the two plates; one plate lot exhibited nearly complete recrystallization, while the other lot contained bands of recrystallization. Although no reduction in tensile properties was observed due to these microstructural variations, a more detailed investigation was recommended to better understand the effect of spin forming on the recrystallization mechanism and the impact on other material properties. The current study examined the microstructure and mechanical properties of the fully processed Al 2219-T62 FPVBH with emphasis on comparing differences between the two plates. Testing was limited to the cone region of the FPVBH. The two specific tasks in this study were: a) Characterization of the microstructure and mechanical properties in the two plates used in the FPVBH. b) Isolation of the impact of each step in the thermo-mechanical processing (TMP) sequence of the FPVBH on recrystallization and grain growth. Mechanical property testing sampled regions containing both large, recrystallized and small, unrecrystallized grains. Tensile strengths and elongations were equivalent in both regions and exceeded MMPDS A-basis allowables. Fracture toughness was equivalent in both regions and fracture morphology was typical of aluminum alloys. The fatigue precrack region of the fracture toughness specimens appeared more faceted in the recrystallized plate. This suggests that property testing that is more sensitive to grain size, such as fatigue crack growth rate and stress corrosion cracking, may be warranted. Alcoa reviewed the production records for the two plate lots and confirmed that composition of both plates was within the allowable range for Al 2219 and that plate processing steps were all performed within production limits. However, the Fe content was notably higher in the fully recrystallized plate, resulting in a greater volume fraction of large Fe-bearing inclusions and Al-Cu constituent particles. Alcoa did not provide details regarding plate processing variables for the two plate lots. It was surmised that deformation by hot rolling was the most accurate TMP simulation of commercial spin forming within a laboratory setting. Both hot rolling and spin forming comprise combinations of shear and compressive deformation. The through-thickness strain caused by hot rolling is symmetrical about the mid-plane, as the rollers contact both sides of the material. Shear strains are maximum at the inner and outer surfaces and trend toward zero at the mid-thickness. In contrast, spin forming utilizes a single roller on the outer surface, producing a gradient in shear stress through the thickness. Shear strains in spin-formed material are highest at the outer surface due to direct contact with the roller, and trend toward zero at the inner surface. Consequently, there will undoubtedly be some discrepancies in the mechanical response of hot-rolled and spin-formed products in the through-thickness direction. Interrupted TMP simulations were performed on remnant plate from the FPVBH forming blank in order to isolate the effects of a post-weld anneal, spin-forming deformation/thermal cycling, and solution heat treatment (SHT) on development of the recrystallized microstructures. Analysis confirmed that large, recrystallized grains formed exclusively during SHT and only after a certain deformation level was exceeded. Remnant plate from both lots that received thermal processing only did not exhibit recrystallization but exhibited fully recrystallized microstructures when processed by hot rolling. This indicates that (1) the level of deformation exceeded the critical level for recrystallization in both plates and (2) that the level of deformation during hot rolling was greater than that which accumulated during spin forming of the FPVBH. Comparison of the hot-rolled plate with the spin-formed material showed that the recrystallized grain size was larger in the FPVBH, providing further evidence that deformation levels were lower in the FPVBH. The most plausible explanation for the differing degrees of recrystallization in the two plates is particle stimulated nucleation (PSN) at the Fe-bearing inclusions and large Al-Cu constituent particles. In the plate with higher solute content, the greater volume fraction of these particles reduces the deformation level required to promote recrystallization. During spin forming the critical deformation level was likely exceeded for the higher Fe content plate resulting in a fully recrystallized microstructure. The lower Fe content plate developed a microstructure that exhibited bands of recrystallization due to the combination of a non-uniform distribution of particles and deformation. The critical deformation required for recrystallization was likely only exceeded during spin forming in regions of higher particle content. While there were differences in Fe content between these lots of Al 2219 plate, the composition of both was within defined alloy limits and Alcoa reported no anomalies during plate production. The observed differences in recrystallization after spin forming may reflect the inherent variability between these two lots of plate. The spin-forming vendor, Spincraft, reported that similar variations in recrystallized microstructures have been observed during examination of multiple spin-forming trials of a different 2xxx series aluminum alloy. Hence, tighter chemical specification limits on impurity levels in Al 2219 plate may help ensure uniform, predictable microstructures. Finally, the tensile and fracture toughness properties were not affected by the microstructural variations noted.

single-piece↗
DOE OSTI2025

Mesoscale modeling of restructuring in high burnup UO 2 fuel

Here, this work aims to simulate the restructuring behavior observed in different regions of high burnup fuel, providing a first-of-its-kind restructuring model for the dark zone and rim region of high-burnup UO 2 fuel. We employed a grand-potential-based phase-field model to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of different initial conditions and different modeling parameters were systematically studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to begin around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Restructuring was demonstrated to be influenced by a combination of initial dislocation densities, burnup rate, subgrain formation rate, and temperature. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. Furthermore, bubble size and distribution do not significantly affect the rate of restructuring. The predicted microstructures are consistent with experimental observations of the restructured regions. Finally, we present a correlation demonstrating the evolution of the restructuring volume fraction as a function of local burnup.

UO2↗
DOE OSTI2025

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

ExaCA v2.0: A versatile, scalable, and performance portable cellular automata application for additive manufacturing solidification

The previously established ExaCA software for performance portable alloy grain structure simulation has been updated to better represent the solidification behavior during complex alloy processing conditions, such as those encountered during metal additive manufacturing (AM), and for improved performance and scalability. Here, an extension to the time–temperature history input data format and the core ExaCA algorithm to include an arbitrary number of melting and solidification events yielded improved prediction of texture for various melt pool geometries, expanding the range of AM-relevant conditions that can be accurately simulated. Improved heat transport process simulation coupling, including the creation of large raster datasets from single track time–temperature history data and in-memory coupling with the new, performance portable finite difference code Finch, were also demonstrated in example studies on the effect of multilayer AM microstructure predictions on hatch spacing and cell size, respectively. Additional new features are detailed and demonstrated, including the ability to perform simulations using various interfacial response function forms, execute simulations on state-of-the-art hardware, improved usability through post-processing versatility, and improved strong and weak scaling performance. The performance, physics, and versatility improvements demonstrated here will further enable large-scale studies on AM process–microstructure relationships that were not previously possible. Furthermore, the usability improvements and ability to run coupled AM process–microstructure simulations using the Finch-ExaCA workflow will facilitate broader use of this open-source software by the computational materials community.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An atomistic study connecting underlying dislocation behavior with superior mechanical properties of NiCoCr medium entropy alloy

NiCoCr-based medium-entropy alloy (MEA) with a simple face-centered cubic crystal phase exhibits excellent mechanical properties, often attributed to the synergy of multiple deformation mechanisms. However, the atomistic origin of their outstanding mechanical response, including microstructural evolution and dislocation behavior under varying strain-rates and orientation, remains unclear. In this work, we employ large-scale molecular dynamics (MD) simulations to investigate the changes in deformation mechanisms along three distinct orientations ([110], [111], [100]) under varying strain rates (1 ×10 8 /sec, 1 ×10 10 /sec, 1 ×10 12 /sec) in the NiCoCr MEA. The presence of the stair-rod and the Shockley partial dislocations under uniaxial tensile strain are found to play a key role in the formation of deformation twinning and ε-martensite, which positively correlates with strain-rate dependent dislocation analysis. These findings further establish the role of the dislocations in controlling the superior mechanical response and excellent fracture toughness of the NiCoCr MEA. Systematic transmission-electron microscopy tests performed on the [111]-oriented crystals, deformed at different strain levels, at room temperature provide clear evidence of both the extended stacking-fault and the stair rods, confirming the predicted microstructural features. Finally, this study offers key insights into the complex nucleation mechanisms of deformation twinning and ε-martensite, such as twinning – and transformation–induced plasticity (TWIP-TRIP), providing valuable guidelines for studying similar material classes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Numerical framework for integrated additive manufacturing-compression molding (AM-CM) of thermoplastic composites

Additive manufacturing-compression molding (AM-CM) has emerged as a transformative technology in advanced composite manufacturing. Additive manufacturing (AM) offers high design flexibility and the ability to produce complex geometries with precisely aligned fibers in the preferred orientation. Compression molding (CM) enhances composite materials by providing excellent dimensional stability, reduced porosity, high production rates, and a smooth surface finish. Despite these advantages, extensive integrated analysis is required to optimize processing conditions for improved fiber orientation distribution (FOD) and porosity control. Here, this study develops a comprehensive numerical model to simulate the AM-CM manufacturing process. The model isolates the effects of both the AM and CM phases while also capturing their integration. Additionally, it accounts for heat transfer, temperature-dependent viscosity, and fiber orientation in the extruded fiber-filled polymer, accurately representing material behavior during processing. This approach enables the analysis of interactions between deposited beads of complex strand shapes and their interface regions after full compression. Moreover, the model predicts key parameters such as polymer flowability, fiber orientation, and temperature evolution in AM-CM parts. By optimizing processing conditions, it facilitates a controlled and predictable microstructure.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Dislocation mechanism based model for stage II fatigue crack propagation rate

Repeated plastic deformation, which of course depends on dislocation mechanism, at or near the crack tip leads to the fatigue crack propagation. By involving the theory of thermally activated flow and the cumulative plastic strain criterion, an effort is made here to model the stage II fatigue crack propagation rate in terms of the dislocation mechanism. The model, therefore, provides capability to ascertain: (1) the dislocation mechanism (and hence the near crack tip microstructures) assisting the crack growth, (2) the relative resistance of dislocation mechanisms to the crack growth, and (3) the fracture surface characteristics and its interpretation in terms of the dislocation mechanism. The local microstructure predicted for the room temperature crack growth in copper by this model is in good agreement with the experimental results taken from the literature. With regard to the relative stability of such dislocation mechanisms as the cross-slip and the dislocation intersection, the model suggests an enhancement of crack growth rate with an ease of cross-slip which in general promotes dislocation cell formation and is common in material which has high stacking fault energy (produces wavy slips). Cross-slip apparently enhances crack growth rate by promoting slip irreversibility and fracture surface brittleness to a greater degree.

Mazumdar, P. K.↗
NASA NTRSDate not supplied

Convolution-Based Numerical Solutions of Transient Temperature Fields during Powder Bed Fusion Additive Manufacturing: Theory, Accuracy, and Computational Cost

Powder bed fusion (PBF) additive manufacturing has found numerous applications in the aerospace domain. However, components fabricated via PBF have a complex time-temperature history that significantly impacts subsequent mechanical performance. This study examines convolution-based numerical solutions of transient temperature fields that support simulations involving arbitrary beam shapes and paths during PBF. The convolutional approach is verified through comparisons with analytical solutions of the temperature field. The computational speed and accuracy of the method are assessed through comparisons with other explicit and implicit numerical techniques. In addition, the straightforward translation of the approach from a CPU to a GPU implementation and the resultant performance improvement are presented. The role of the technique in predicting microstructure evolution during PBF (for a greater process-structure-property-performance framework) is also demonstrated. This work supports the development of computational materials methods for understanding and controlling the time-temperature history during PBF.

powder bed fusion↗
DOE OSTI2025

Microstructurally validated stable and predictable swelling in low-enriched uranium monolithic U-10Mo fuel mini-plates

Qualification of the low-enriched uranium (LEU) monolithic U-10 wt%Mo (U-10Mo) plate-type fuel system requires a demonstration of a stable and predictable fuel swelling behavior over the anticipated operating conditions of the United States high-performance research reactors (USHPRRs) selected for conversion to LEU operation. This will allow each reactor to develop appropriate safety margins that will retain fuel element lifetime coolability. Additionally, the fuel system must maintain performance attributes when fabricated at a commercial scale. The Mini-plate 1 experiment represents the first irradiation test of commercially fabricated miniaturized monolithic LEU U-10Mo fuel plates. Here, the swelling behavior within this experiment was compared against that of historical fuel developmental tests to reveal that the commercially fabricated fuel performed within the current recommended U-10Mo swelling model's predictions. Additionally, the fuel microstructural evolution was evaluated to link initial conditions to subtle variations detected in the swelling response, providing validation and confidence that the fuel system is robust.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

FFTF HT9 Cladding Microstructure Characterization

The sodium-cooled fast reactor (SFR) is a promising candidate for next generation nuclear reactors, operating at extreme conditions which include high temperatures (>500?C core outlet temperature) and significant neutron damage. High-Cr martensitic HT9 steel is an excellent candidate for SFR cladding and duct material due to its compatibility with liquid sodium, good thermal conductivity, resistance to void swelling, and strong creep rupture strength [1-4].However, the harsh in-core environment of SFRs can cause complex microstructural changes and mechanical property degradation in HT-9. Ensuring the safe use of HT9 cladding for metallic fuel requires both a thorough understanding of its mechanical response to microstructure evolution as well as reliable microstructure-sensitive modeling predictions. Microstructure-sensitive modeling of high temperature creep behavior in HT9 cladding for SFR applications currently lack experimental data to model the phenomena accurately. To fill this need, methods to perform microstructural characterization have been developed and performed on HT9.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Coupling Microstructural Evolution Simulations to Material Property Degradation Predictions for Plasma-Facing Materials

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Predicting Mechanical Properties from Microstructure Images in Fiber-Reinforced Polymers Using Convolutional Neural Networks

Evaluating the mechanical response of fiber-reinforced composites can be extremely time-consuming and expensive. Machine learning (ML) techniques offer a means for faster predictions via models trained on existing input–output pairs and have exhibited success in composite research. This paper explores a fully convolutional neural network modified from StressNet, which was originally used for linear elastic materials, and extended here for a non-linear finite element (FE) simulation to predict the stress field in 2D slices of segmented tomography images of a fiber-reinforced polymer specimen. The network was trained and evaluated on data generated from the FE simulations of the exact microstructure. The testing results show that the trained network accurately captures the characteristics of the stress distribution, especially on fibers, solely from the segmented microstructure images. The trained model can make predictions within seconds in a single forward pass on an ordinary laptop, given the input microstructure, compared to 92.5 h to run the full FE simulation on a high-performance computing cluster. These results show promise in using ML techniques to conduct fast structural analysis for fiber-reinforced composites and suggest a corollary that the trained model can be used to identify the location of potential damage sites in fiber-reinforced polymers.

Sun, Yixuan (ORCID:0000000311093380)↗
DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗