Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microstructure simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A sensing material-free and simple readout MEMS sensor for detecting Helium

Abstract In this work, we report a method that enables a standard electrostatic MEMS device to perform complex sensing functionalities, such as detecting the presence of helium without a sensing material or a conditioning circuit. Helium is a noble, odorless, non-reactive gas that is very challenging to detect. It is used in critical applications such as storing nuclear fuel waste inside a dry cask. In these applications, its leakage from the dry cask may indicate the cask's safe operation's degradation. A departure from the common practice of exciting the MEMS around its mechanical resonance, the method is based on exciting the MEMS around its electrical resonance circuit. This method shows that the tiny difference between the air dielectric constant (1.00059) and helium (1.000067) corresponding to only a few Femtofarad level capacitances produces a 25 mV difference without a conditioning circuit. Simulation results confirmed those findings and explored the sensor response at different operation conditions. This method eliminates the need for a heated microstructure and the need for absorption material. This method is not limited to gas sensing. It can be applied to other sensing mechanisms, such as acceleration and pressure measurements, and eliminate the complex circuit to read small capacitance in these applications.

Mohaidat, Sulaiman↗
DOE OSTI2025

Metallurgical Analysis and Forward Modeling of Weld Distortion in SMR Containment Vessels

This work aimed to apply Sandia’s expertise in metallurgy and modeling to enable the use of hybrid laser arc welding for building nuclear reactor containment structures, via a collaboration with Holtec International. Experimental observations were coupled with finite element analysis to resolve microstructure development, mechanical properties, distortion, and residual stress in welds relevant to the production of the Holtec SMR-160. High residual stresses were observed in welds that were not subjected to preheat. Meanwhile, the microstructure of the welds generally exhibited a narrow heat affected zone relative to conventional arc welds. FEA appeared to be effective in simulating the thermal/mechanical conditions that occur during hybrid laser arc welding of simplified and instrumented test welds. Subsequently, FEA was used to perform sensitivity analyses for various weld geometries that would be prohibitively costly to assess with physical experiments. Insights from the study were used to inform Holtec’s welding process, and successful production welds were performed in 2025.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Process-driven roadmap for depositing super duplex stainless steel via wire Arc additive manufacturing

Here, this study systematically investigates the effects of shielding gas, bead spacing, weld mode, and travel speed on the phase balance, porosity, and hardness of wire arc additively manufactured (WAAM) super duplex stainless steel (ER2594). Robotic WAAM was employed to fabricate multilayer walls under systematically varied process conditions, followed by phase transformation simulations, X-ray computed tomography (XCT), electron backscatter diffraction (EBSD), and microhardness evaluation. Thermodynamic simulations predicted rapid cooling of AM process can suppress the potential formation of deleterious precipitates which was later validated via cross-sectional microstructure analyses of printed samples. XCT revealed porosity levels below 0.2% for all deposits, with reduced travel speed significantly lowering defect volume. Microstructural analyses revealed the evolution of various austenite precipitates, including grain boundary austenite (GBA), Widmanstätten austenite (WA), and intergranular austenite (IGA), sequentially upon cooling of the ferrite phase. Among all process parameters, weld transfer mode exhibited the strongest influence on phase balance; pulsed mode promoted higher ferrite retention (~ 36%) compared to RapidX mode. No consistent relationship between stabilized phase fraction and captured microhardness was observed. This work provides critical insights for optimizing WAAM parameters to control phase balance and mechanical performance, which is essential for producing high-integrity super duplex stainless-steel components for nuclear and marine applications.

Grain orientation↗
DOE OSTI2025

Effect of pore fluid chemistry on the mechanical behavior of a divalent compacted bentonite, an experimental and constitutive study

Ongoing research in isolating high-level nuclear waste and spent fuel has highlighted compacted bentonite as a suitable material for engineered barrier systems in deep geological repositories due to its extraordinary swelling and retention properties. This research focuses on the chemo-mechanical behavior of compacted bentonite exposed to different pore fluids with different concentrations and loading conditions. The study involves swelling pressure and compressibility experiments along with mineralogy analysis employing X-ray diffraction (XRD) and Cation exchange. The tests were conducted on BCV (a Mg/Ca- bentonite) compacted at a dry density of 1.48 ± .02 Mg/m 3 . An advanced chemical-mechanical constitutive model for unsaturated highly expansive clays was adopted to simulate the material response and better understand its behavior. The model is able to account for the main phenomena at both macro and microstructural levels and the interactions between them. The model successfully replicated experimental observations. The XRD analyses support the macroscopic observation, indicating that salinity impacts crystalline swelling as demonstrated by the reduction of basal spacing from 19.27 Å to 15.68 Å when the osmotic suction increases from 0 MPa to 33 MPa. The results suggested that the osmotic pressure generated by the concentration in the pore fluids promotes a reduction in swelling pressures, swelling strains, and crystalline swelling of clay minerals. Also, it affects the pre-consolidation stress and the compressibility of the compacted samples. In conclusion, it was also observed that both solution type and solution concentration impact the clay swelling pressure.

Chemo-mechanical constitutive model↗
DOE OSTI2024

Efficient mapping between void shapes and stress fields using Deep Convolutional Neural Networks with sparse data

Establishing fast and accurate structure-to-property relationships is an important component in the design and discovery of advanced materials. Physics-based simulation models like the finite element method (FEM) are often used to predict deformation, stress, and strain fields as a function of material microstructure in material and structural systems. Such models may be computationally expensive and time intensive if the underlying physics of the system is complex. This limits their application to solve inverse design problems and identify structures that maximize performance. In such scenarios, surrogate models are employed to make the forward mapping computationally efficient to evaluate. However, the high dimensionality of the input microstructure and the output field of interest often renders such surrogate models inefficient, especially when dealing with sparse data. Deep convolutional neural network (CNN) based surrogate models have shown great promise in handling such high-dimensional problems. In this paper, a single ellipsoidal void structure under a uniaxial tensile load represented by a linear elastic, high-dimensional and expensive-to-query, FEM model. We consider two deep CNN architectures, a modified convolutional autoencoder framework with a fully connected bottleneck and a UNet CNN, and compare their accuracy in predicting the von Mises stress field for any given input void shape in the FEM model. Additionally, a sensitivity analysis study is performed using the two approaches, where the variation in the prediction accuracy on unseen test data is studied through numerical experiments by varying the number of training samples from 20 to 100.

surrogate modeling; convolutional neural networks;↗
DOE OSTI2024

Phase-Field Modeling of Mechanical Damages in Ceramic Matrix Composites

Developed a phase-field model for mechanical damages in CMCs, which incorporates the CMC microstructures, matrix cracking, fiber breakage, and interfacial sliding. Two types of mechanisms, fiber bridging and fiber pull-out, are considered. The obtained simulation results agree with experimental observations and an analytical solution. Simulation results suggest that the performance of CMCs would be enhanced with thicker fibers, longer fibers, and higher fiber density. Opposite trends of interfacial sliding resistances are suggested for the two types of situations. In reality, a mixture of the two situations may exist, and then an intermediate interfacial sliding resistance may be optimal.

Xue, Fei↗
DOE OSTI2025

Meso-scale modeling of UO2 nuclear fuel to high burnup

To improve the economics of light water reactors for commercial nuclear energy generation, utility operators are seeking to obtain regulatory approval to run UO2 fuel to higher levels of burnup. One potential impediment to obtaining this approval is the phenomenon of fuel fragmentation, relocation, and dispersal (FFRD). FFRD can result when fuel experiences a rapid temperature transient, such as that occurring during a Loss Of Coolant Accident (LOCA). FFRD has historically been most associated with the rim region in UO2 fuel pellets, where the phenomenon of fragmentation is also referred to as pulverization due to the small size of the fragments. More recent evidence suggests that the so-called “dark zone” (due to its appearance in micrographs) that can be observed in the mid-radial regions of high burnup fuel is also susceptible to FFRD. Although empirical fuel performance models have been developed that can adequately predict pulverization in the rim region under typical LWR conditions, a scientific understanding of what underlies fuel restructuring and subsequent FFRD is lacking even in the rim region, and no models are currently available for the behavior the dark zone. To address these challenges, the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has employed a multi-scale modeling approach to improve scientific understanding and develop new fuel performance models. In this talk, I will focus on meso-scale efforts, which form a crucial link between atomic-scale and engineering-scale models. Phase-field modeling combined with cluster dynamics is used to predict the restructuring process in the rim region. Phase-field fracture modeling, informed by atomistic simulations, is used to predict the onset of pulverization in the rim region. Combining these techniques together allows the extent of rim pulverization to be predicted. The formation and evolution of the dark zone has also been simulated with the phase-field method, using an improved approach to vacancy source term parameterization. The work shows the important impact of microstructure on fuel performance.

fracture↗
DOE OSTI2024

Virtual texture analysis to understand microstructure effects on deformation twinning and detwinning behavior in BCC metals

Understanding and predicting deformation twinning contributions to plastic deformation in BCC metals has been a long-standing challenge due to the interplay with dislocation slip and non-Schmid effects that render an asymmetry under tension and compression. This paper uses molecular dynamics simulations to understand the effect of unloading and grain orientation on deformation twinning in a nanocrystalline Fe (nc-Fe) system as a model BCC metal. The nc-Fe system is loaded under uniaxial stress tension and compression to understand the effect of grain orientation (Schmid effects) on deformation twinning behavior and the tension/compression asymmetry (non-Schmid effects). A new virtual texture analysis “VirTex” tool is used to understand the role of grain orientation on the nucleation of twins and their contributions to the observed stress–strain response. For certain grain orientations, the twinnability is observed to be different in tension and compression. In addition, the flow stress accommodation from twins in certain grains is observed to be different in tension and compression and different from that for the grains. Subsequent unloading leads to detwinning in the deformed microstructures, where the extent of detwinning depends on the strain from which the system is unloaded and on the morphology of the twin. Lastly, the simulations are carried out to analyze the role of the Schmid factor on the twinnability and asymmetry in tension and compression.

Kannan, Aadhithyan [University of Connecticut, Sto↗
DOE OSTI2025

Laser‐Powder Melt Pool Solidification Dynamics and Microstructural Engineering of Ti‐5553 Microlattices

The fine geometric and topological control afforded by additive manufacturing technologies has enabled the manufacture of architected materials across length scales, and enabling tunable mechanical performance as a function of local and global design. Progress has been made to tune the mechanical response of architected materials through geometry, but understanding how the geometry and processing conditions will inform the microstructure remains a challenge due to the rapid solidification in laser powder bed fusion. This study uses in situ X-ray imaging and electron backscatter diffraction microscopy to demonstrate that the melt pool size, microstructure morphology, and elastic strain distribution is influenced by a combination of lattice geometry and laser processing conditions. These results indicate that within larger melt pools the local thermal gradients are sufficient to enable a columnar-to-equiaxed transition across the melt pool. Furthermore, the solidification mechanisms producing these microstructures are examined across the first 5 ms of melting and solidification, described via in situ high-speed X-ray imaging and mirrored via multiphysics simulation.

additive manufacturing↗
DOE OSTI2025

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗
DOE OSTI2025

SOC Microstructural Property Estimator

This pre-trained ML model is a tool that uses basic compositional parameters for porous solid oxide cell (SOC) electrodes - the phase fractions and mean particle/pore diameters – as inputs and uses them to estimate additional electrochemical performance parameters: active (i.e., connected) TPB density, all tortuosity factors, and phase pair specific interfacial areas. The electrode is assumed to be composed of two solid phases and a pore phase. The property calculations are performed using neural network regression models trained on a large bank of synthetic electrode microstructural data that NETL has generated using the program DREAM3D (that bank is also hosted on EDX: https://edx.netl.doe.gov/dataset/soc-synthetic-microstructure-bank). This means the generated parameters are based on training from actual measured properties from 3D microstructures, not estimated from geometric simplifications. This tool was developed and is intended to replace percolation theory calculations in models that use hypothetical electrode properties. An example use case would be running SOC performance simulations across a parametric sweep of electrode designs (e.g., varying phase fractions and particle sizes) and assessing how it impacts the electrochemical performance of the SOC. Within the parameter space of the training data (statistics of that parameter space is provided in the readme file), this model achieves sub-5% mean absolute percent errors, an order of magnitude less error than percolation theory across the same parameter space. However, be aware that this tool was developed with parametric simulations in mind, and users are encouraged to assess accuracy for their own specific use case rather than taking accuracy metrics at face value. More info, including a usage guide, is in the included readme file. This tool should be cited with the DOI number provided.

Electrode Microstructure↗
DOE OSTI2025

In situ irradiation of spent nuclear fuels

To improve the economics of commercial nuclear reactors, nuclear vendors and utilities are seeking approval for increased burnup usage of the existing nuclear fleet. This is especially critical for meeting the clean energy initiative by the United States Government, Department of Energy (DOE). However, one of the key challenges the nuclear industry faces in this regard is that the regions exposed to high burnup and low temperatures exhibit a fine-grained microstructure with large bubbles known as high-burnup structure (HBS) [1]. The formation of HBS has been correlated to the diminished performance of the reactor, as well as fuel fragmentation and pulverization during transient and accidental conditions [2]. Therefore, it is paramount to understand the mechanisms for HBS formation along with its impact on the properties and performance of nuclear fuels. While existing programs, such as Nuclear Energy Advanced Modeling and Simulation (NEAMS) and Advanced Fuel Campaign (AFC) are focusing on evaluating the performance impact of HBS, the physical mechanisms contributing to HBS formation are still not fully understood. In addition, having predictive capabilities and sound understanding of the microstructural evolution of nuclear fuel is essential for accelerated development, qualification, and deployment of new nuclear materials and novel reactor designs for advanced nuclear reactors. There is a lack of consensus among researchers regarding the mechanisms leading to such restructuring observed in HBS. Grain subdivision due to polygonization versus recrystallization, continuous versus discrete recrystallization occurring in tandem or conjunction, etc., have been proposed and debated. In general, it is hypothesized that defect accumulation and dislocation interaction within the grains cause the realignment of dislocations into grain boundaries, leading to the new subgrain formation, which over time transforms into new grains. However, due to the lack of transient data, the importance of fission rate, irradiation, thermal, and stress history of the fuel on the restructuring could not be assessed. In situ microstructural evolution under various irradiation conditions is desired to bridge this gap. Alternatively, phase-field-based models have been developed to capture HBS formation via discrete recrystallization utilizing the classical nucleation approach [3–5]. However, in these models, the grain nucleation criteria are often defined based on empirical relations for burnup and fission gas density leading to dislocation density change. A mechanistic approach to capture the dislocation interaction with the microstructural features leading to grain subdivision is lacking.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Scale Invariance of Hot Spot Formation in TATB High Explosives

Shock-induced detonation of insensitive high explosives based on 1,3,5-triamino-2,4,6-trinitrobenzene starts with formation of hot spots at microstructural defects but has eluded atomistic modeling treatment at micron length scales. To this end, we performed multimicron scale all-atom molecular dynamics (MD) simulations of hot spots that form during the collapse of cylindrical pores with diameters between 10 and 300 nm. Our MD simulations show that hot spots formed at pores larger than 20 nm exhibit temperature fields with scale-invariant features for sizes up to at least 300 nm. Through a continuum-based grain-scale modeling framework, we span and extend beyond the size scales currently accessible to MD and find that hot spot scale invariance is a general feature that arises when the mechanical strength is insensitive to strain rate. Finally, our results demonstrate the applicability of all-atom MD to simulate the complicated dynamical evolution of micron-sized systems and bolster confidence in insights from MD simulations of materials that exhibit strength with negligible rate dependence over the relevant intervals.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of mechanistic, microstructure-informed BISON models for fission product-induced failure mechanisms in advanced nuclear fuels

The US continues to prioritize commercial deployment of advanced reactors—particularly those utilizing U-Zr metallic and TRISO particle fuels. Fission product-induced failure mechanisms for these systems include fuel–cladding chemical interaction in metallic fuel rods, which threatens cladding integrity, and Pd penetration in TRISO particles, which degrades SiC layer properties. Empirical models can be applied within the bounds of existing irradiation databases, but their utility is limited when considering new designs or investigating fundamental material behaviors. The Nuclear Energy Advanced Modeling and Simulation Program is therefore developing mechanistic models for these behaviors, which are expected to aid in fuel design, assist in development of failure mitigation strategies, and provide support for qualification and licensing. These efforts leverage multiscale capabilities to capture the microstructural features and processes that govern these behaviors. This talk summarizes past and ongoing efforts to develop and validate these models for the BISON fuel performance code.

FCCI↗
DOE OSTI2025

Microstructural evolution and phase stability in Nb-containing interstitial Fe-Mn-Co-Cr-C high-entropy alloys: An in-situ synchrotron X-ray diffraction study during laser melting

The influence of Nb on phase stability and microstructural evolution in an interstitial Fe-Mn-Co-Cr-C high-entropy alloy was investigated using in-situ synchrotron X-ray diffraction (SXRD) during laser melting. Scheil-Gulliver simulations predict the formation of σ and γ-f.c.c. phases in all three alloys, along with NbC in Nb-containing compositions. SXRD confirmed the presence of most predicted phases, but the σ phase was absent. Nb promotes crystallite refinement and increases dislocation density, though excessive additions reduce refinement efficiency due to solubility limits and secondary phase formation. Furthermore, Nb addition also enhances ε-h.c.p. phase formation by reducing stacking fault energy through NbC-induced carbon depletion. Analysis of intensity peak evolution reveals that Nb alters preferred grain orientations, reducing {111} γ intensity while enhancing {220} γ , leading to a more isotropic grain distribution. Overall, Nb plays a key role in phase selection, microstructure refinement, and preferred orientation evolution, allowing the tailored microstructure of high-entropy alloys via rapid solidification.

Interstitial high entropy alloys↗
DOE OSTI2026

Formation of Bimetallic Nanoparticles via Exsolution Using a Reducible Metal Oxide Capping Layer

Bimetallic nanoparticles are promising catalysts that can improve performance in heterogeneous catalysis and solid-state electrochemistry. Exsolution is a useful method for forming such nanoparticles; however, it is limited by the elements present within the host oxide lattice. Here, in this work, we develop and demonstrate a strategy to form bimetallic particles from La 0.5 Sr 0.5 Ti 0.94 Ni 0.06 O 3 (LSTN) exsolution and using a reducible SnO 2 capping layer, expanding the range of elements available for bimetallic nanoparticle formation. Using this capping layer strategy, we formed nickel–tin (Ni 0 –Sn 0 ) bimetallic nanoparticles via exsolution. We used in situ near-ambient pressure X-ray photoelectron spectroscopy to monitor surface chemical changes during exsolution, showing that first, SnO 2 volatilized. This SnO 2 loss exposed the perovskite surface of LSTN to reducing conditions, which induced Ni exsolution, and compounded with SnO 2 reduction led to the formation of bimetallic Ni 0 –Sn 0 particles. To evaluate the associated microstructural evolution, we measured grazing incidence small-angle X-ray scattering (GISAXS), which confirmed the loss of the SnO 2 capping layer, and scattering simulations suggested the formation of bimetallic particles. We confirmed the bimetallic nanoparticle composition and morphology by Auger spectroscopy and scanning transmission electron microscopy. The resulting bimetallic nanoparticles were smaller and more thermally stable than the monometallic Ni counterparts on LSTN. This capping layer and exsolution approach allow synthesizing multimetallic nanoparticles and can be applied to other reducible metal oxides and perovskite hosts, broadening the compositional space for advanced catalytic materials.

36 MATERIALS SCIENCE↗