Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

Analyzing the effect of pressure on the properties of point defects in $\gamma\mathrm{U–Mo}$ through atomistic simulations

Uranium–molybdenum (U–Mo) alloys in monolithic fuel foil are the primary candidate for the conversion of high-performance research reactors in the USA. Monolithic fuel is utilized in a plate-type design with a zirconium diffusion barrier and aluminum cladding. These fuel types are unique in that they contain no plenum for the release of fission gases, which, in conjunction with the aluminum cladding, can lead to large stress states within the fuel. The nature of how fundamental processes of radiation damage, including the evolution of point defects, under such stresses occur is unknown. In this work, we present molecular dynamics simulations of the formation energy of point defects under applied stress. We report this work will allow for the implementation of stress-dependent microstructural evolution models of nuclear fuels, including those for both fission gas bubble growth and creep, which are critical to ensure the stable and predictable behavior of research reactor fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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 OSTI2023

Nanoscale reshaping of resonant dielectric microstructures by light-driven explosions

Femtosecond-laser-assisted material restructuring employs extreme optical intensities to localize the ablation regions. To overcome the minimum feature size limit set by the wave nature of photons, there is a need for new approaches to tailored material processing at the nanoscale. Here, we report the formation of deeply-subwavelength features in silicon, enabled by localized laser-induced phase explosions in prefabricated silicon resonators. Using short trains of mid-infrared laser pulses, we demonstrate the controllable formation of high aspect ratio (>10:1) nanotrenches as narrow as ~ $λ / 80$. The trench geometry is shown to be scalable with wavelength, and controlled by multiple parameters of the laser pulse train, such as the intensity and polarization of each laser pulse and their total number. Particle-in-cell simulations reveal localized heating of silicon beyond its boiling point and suggest its subsequent phase explosion on the nanoscale commensurate with the experimental data. The observed femtosecond-laser assisted nanostructuring of engineered microstructures (FLANEM) expands the nanofabrication toolbox and opens exciting opportunities for high-throughput optical methods of nanoscale structuring of solid materials.

42 ENGINEERING↗
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 OSTI2020

Crystallographic orientation and spatially resolved damage in a dispersion-hardened Al alloy

The in-situ neutron diffraction technique, in combination with both the full-field crystal elasto-viscoplastic finite element model and microstructural characterization, was used to study the deformation-induced damage anisotropy in a commercial Al alloy, subjected to uniaxial tensile and cyclic loading. The simulations capture well the crystallographic-orientation-dependent lattice strain behavior. The hard grains, e.g. those orientated with the <111> and <422> orientations parallel with the uniaxial loading direction (LD), feature large Taylor factors and seem more prone to form damage-related band structures. Their effective elastic moduli decrease with the accumulation of damage, which are different from the soft grains orientated with the <200> orientation along the LD. Correlation between the distribution of voids and that of the residual lattice strain developed after failure may exist. The maximum tensile type residual lattice strain observed after failure may be resulted from the band structure formed in the hardest <111> grains. It was revealed that the band structure triggered by the hard particles could be one of sources of damage. In addition, while the specimen was obviously damaged, a fast stress relief was evidenced after unloading from the tension, especially at the beginning of unloading. Our present investigations provide a novel method for exploring the damage mechanisms of polycrystalline materials during plastic deformation.

36 MATERIALS SCIENCE↗
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 OSTI2023

Direct Measurements of Ionic Transport Behavior of Dual-Layer Porous Electrodes

To improve power and cycling performance of lithium-ion batteries, dual-layer or porosity-gradient electrodes have been proposed. By engineering a higher porosity close to the separator, the intention is to improve ion transport where it is most impactful. In this research, MacMullin numbers of two dual-layer anode samples are tested using an impedance measurement technique developed previously. To characterize the microstructure of each layer independently, we developed an improved transmission-line model that accounts for each layer's properties.Virtual experiments using COMSOL Multiphysics to simulate impedance measurements are used to examine and improve the accuracy of the numerical inversion procedure. The results for the two dual-layer anodes studied show that MacMullin numbers follow expected trends, though the anodes are quite different from each other.

Dual-Layer Electrodes↗
DOE OSTI2020

Imaging transient solidification behavior

Solidification processing offers the first opportunity to control microstructure, properties, and performance in metallic alloy components. Until recently, microstructural evaluations were limited to post-solidification characterization by destructive methods. We review the development of time-resolved, in situimaging techniques capable of capturing solid-liquid interfacial evolution in metallic alloys with high spatial and temporal resolution underdiverse solidification conditions relevant for applications ranging from conventional directional solidification, crystal growth, and casting,to welding and additive manufacturing. These experiments enable direct visualization of transient behaviors that would otherwise remain unknown, uniquely providinginsights into the physics that impact microstructure and defect development and strategies for microstructuralcontrol and defect mitigation. Understanding microstructural evolution and the characteristics that form under various solidification conditions is essential for the development of multiscale, experimentallyinformed predictive modeling. This is highlightedby solidification simulations that utilize in situmeasurements of solidification dynamics from state-of-the-art experimental techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Investigation on laser absorption and x-ray radiation in microstructured titanium targets heated by short-pulse relativistic laser pulses

The enhancement effect of a microstructured surface on laser absorption and characteristic Kα emission has been investigated by measuring K-shell emission from titanium (Ti) targets irradiated with high-intensity (~10 20 W cm -2 ), subpicosecond (500 fs) laser pulses. The experimental results indicate a modest enhancement (1.6×) of Kα emission from microstructured targets compared to flat foils, but with a similar intensity and profile of Heα and Li-like satellites. Particle-in-cell (PIC) simulations are implemented to further understand the underlying physical processes in the laser interaction with both targets, interpreting the mechanisms responsible for the Kα enhancement. The reasons for the lower-than-expected enhancement of Kα emission are discussed. The rapid heating of the bulk plasma might result in the premature shutdown of Kα emission before the thermalization of hot electrons or even the end of laser pulses, suggesting that the use of Kα emission as a diagnostic of the hot-electron yield or relaxation could lead to a misinterpretation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Characterization of Fe-6Si Soft Magnetic Alloy Produced by Laser-Directed Energy Deposition Additive Manufacturing

Abstract Commercial electrical steels, Fe-Si alloys with < 4 wt.% Si, are inexpensive and efficient materials for electrical power conversion. Further efficiency improvements require increasing the silicon concentration to 6 wt.%, at which point the material becomes brittle and difficult to form by conventional rolling and sheet fabrication methods. Additive manufacturing stands to overcome challenges with commercial manufacturing techniques by leveraging near-net-shape fabrication. The wide array of process conditions provides additive manufacturing with increased flexibility, enabling control over the microstructure and mechanical properties. This work explores the microstructures and magnetic properties of ring-shaped Fe-Si alloys produced using concentric and cross-hatch tool paths on a laser-directed energy deposition additive manufacturing system. Concentric-built samples exhibit elongated grain structures while cross-hatch-built samples comprise lower aspect ratio grain structures. Thermal finite element analysis simulations model the stress conditions produced by the different scan path geometries. Microhardness measurements probe the mechanical properties as a function of anneal temperature, providing a qualitative understanding of the intergranular defect density. Soft magnetic properties measured under quasistatic and AC conditions show frequency- and microstructure-dependent coercivity and permeability. Finally, analysis of the core loss quantifies how the build strategies and thermal treatments influence efficiency in electrical power conversion applications. Understanding the influences of scan path geometry and thermal treatment provides a pathway towards application of additively manufactured soft magnetic materials.

Adamczyk, Jesse M.↗
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 OSTI2020

Theoretical Correlation of Elemental Distribution of Nd and Pr in Ce-Fe-B Microstructure With Hard Magnetic Properties

Relatively resource-rich but property inferior-Ce-Fe-B magnet can be improved by partial replacement of Ce by Nd and/or Pr. In addition to the amount of Nd/Pr, their distribution profile in microstructure plays an important role. From our first principles density functional theory (DFT) calculation, the substitution energy of Ce by Pr/Nd is negative in Ce 2 Fe 14 B (2:14:1) while that for laves phase, CeFe 2 is positive, implying that Nd/Pr stabilize 2:14:1 and suppress the formation of the CeFe 2 phase. Further, micromagnetic simulation indicates that homogenized distribution of Nd/Pr improves squareness of demagnetization curve, while core (Ce-rich)-shell (Nd/Pr-rich) 2:14:1 grain structure enhances coercivity. Magnetic properties of Ce-Fe-B can be optimized by manipulating distribution profile of chemical element in microstructure based on their subtle difference in thermodynamic property, which is an effective pathway to design optimized chemical composition and processing route for high-performance magnet.

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 OSTI2022

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Active sampling for neural network potentials: Accelerated simulations of shear-induced deformation in Cu–Ni multilayers

Neural network potentials (NNPs) can greatly accelerate atomistic simulations relative to ab initio methods, allowing one to sample a broader range of structural outcomes and transformation pathways. In this work, we demonstrate an active sampling algorithm that trains an NNP that is able to produce microstructural evolutions with accuracy comparable to those obtained by density functional theory, exemplified during structure optimizations for a model Cu–Ni multilayer system. We then use the NNP, in conjunction with a perturbation scheme, to stochastically sample structural and energetic changes caused by shear-induced deformation, demonstrating the range of possible intermixing and vacancy migration pathways that can be obtained as a result of the speedups provided by the NNP. The code to implement our active learning strategy and NNP-driven stochastic shear simulations is openly available at https://github.com/pnnl/Active-Sampling-for-Atomistic-Potentials .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2024

Elucidating Abnormal Grain Growth in Thermomagnetic Processed Materials with Transfer Learning and Reinforcement Learning

The goal of this research program is to establish the mechanism governing local grain boundary motion, which is needed to design and process desirable microstructures for better performance, by identifying the relative contributions of grain boundary (GB) energy and mobility to grain growth. Classical models for grain growth assume that the primary mechanism for reducing the total interfacial energy is area reduction and that GB restructuring is not significant. This assumption implies that grain growth is locally driven by curvature. However, recent experimental observations using new non-destructive 3D x-ray diffraction microscopy techniques (3D-XRM) reveal that classic descriptors (i.e., curvature, number of neighbors, grain size) do not predict real grain growth. Instead, local GB motion appears to be governed by its energy relative to its neighbors such that low-energy boundaries replace those of higher energy. However, simulations that incorporate GB energy anisotropy still fail to reproduce these observations. These discrepancies suggest that the common assumption for grain growth theory must be re-examined to predict and, thus, control microstructure evolution in real polycrystals. A significant challenge to testing this assumption is due to anisotropic GB mobility. Mobility may cause abnormal grain growth or affect the final grain shapes or growth rate but its true contributions are unknown because it is difficult to measure. For example, observations in Fe have found that grains associated with high energy and high mobility boundaries tend to experience abnormal grain growth, whereas abnormal grain growth is associated with low energy and high mobility boundaries in alumina. As mobility and energy both control GB motion, it is challenging to isolate the local driving forces necessary to test the common assumption that the primary mechanism is area reduction. The novelty of this work is the use of machine learning tools to capture GB mobility and energy from 3D-XRM measurements in polycrystals to test the common assumption used in grain growth models. Machine learning can capture high-order correlations in dynamic systems like those found in the evolving GB topology. The PIs have developed a physics-regularized interpretable machine learning microstructure evolution (PRIMME) model that accurately replicates the grain growth behavior of its trained data set.

36 MATERIALS SCIENCE↗