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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effects of Internal and External Heat Sources on Cladding Microstructure and Rupture Performance

In the event of a Loss of Coolant Accident (LOCA), the primary supply of cooling water for a nuclear reactor is lost, leading to a significant pressure differential across the cladding wall. Without adequate cooling, the fuel rods continue to heat as a result of fission reactions. Research at Oak Ridge National Laboratory’s (ORNL’s) Severe Accident Test Station (SATS) is currently focused on evaluating fuel cladding performance using an external infrared lamp as a heat source, whereas legacy testing primarily utilized an internal heating approach. While external heating may better simulate the effect of neighboring fuel rods heating a central rod, it may not accurately represent the internal heat absorption from the fuel during accident transients. The heating dynamics depend greatly on the fuel rod's position within the bundle and the reactor. To fully understand the implications of a LOCA event and assess the influence of internal heating on cladding performance, combined internal heating and pressurization capability was developed at ORNL. Tests were conducted to compare cladding segments heated internally (representing heat from the fuel within the rod) and externally (representing heat from adjacent fuel rods). The findings indicate that both internal and external heating result in comparable rupture temperatures during 5°C/s laboratory LOCA tests and also agree with legacy test data. Axial temperature gradients and internal heat source dispersal were found to significantly impact cladding deformation and rupture geometry. Clear modifications were outlined to further improve the capability with heating rates above 5°C/s.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Dislocation loop evolution in Kr-irradiated ThO 2

Here, the early stage of microstructural evolution of ThO 2 , under krypton irradiation at 600, 800, and 1000°C, was investigated using in situ transmission electron microscopy (TEM). Dislocation loops grew faster, whereas their number density decreased with increasing irradiation temperature. Loop density was found to decrease with ion dose. Interstitial dislocation loops, including Frank loops with Burgers vector of $\textit{a}/3\langle111\rangle$ and perfect loops with Burgers vector of $\textit{a}/2\langle110\rangle$, were determined by traditional TEM and atomic resolution–scanning TEM techniques. Atomistic and mesoscale level modeling are performed to interpret experimental observations. The migration energy barriers of defects in ThO 2 were calculated using density-functional theory. The energetics of different dislocation loop types were studied using molecular dynamics simulations. Loop density and diameter were analyzed using a kinetic rate theory model that considers stoichiometric loop evolution. This analysis reveals that loop growth is governed by the mobility of cation interstitials, whereas loop nucleation is determined by the mobility of anion defects. Lastly, a rate theory model was used to extract the diffusion coefficients of thorium interstitials, oxygen interstitials, and vacancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

The evolution of grain boundary energy in textured and untextured Ca–doped alumina during grain growth

The role of anisotropic grain boundary energy in grain growth is investigated using textured microstructures that contain a high proportion of special grain boundaries. Textured and untextured Ca–doped alumina was prepared by slip casting inside and outside a high magnetic field, respectively. At 1600°C, the textured microstructure exhibits faster growth than the untextured microstructure and its population of low–angle boundaries increases. Atomic force microscopy (AFM) is employed to measure the geometry of thermal grooves to assess the relative grain boundary energy of these systems before and after growth. In the textured microstructure, the grain boundary energy distribution narrows and shifts to a lower average energy. Conversely, the energy distribution broadens for the untextured microstructure as it grows and exhibits abnormal grain growth. Further analysis of the boundary networks neighboring abnormal grains reveals an energy incentive that facilitates their growth. These results suggest that coarsening is not the only dominant grain growth mechanism and that the system can lower its energy effectively by replacing high energy boundaries with those of low energy. The faster growth of lower energy boundaries suggests that isotropic simulations do not adequately account for anisotropic grain growth mechanisms or anisotropic mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Pore Resolved Simulations of Joule Heating in Fibrous Media using an Embedded Boundary Method

Joule heating has been regarded as an energy-efficient and sustainable method for heating materials and gases at large scales. The modeling of local temperature effects at pore-resolved scales for such systems, however, has been difficult to achieve due to challenges in coupling thermo-chemical processes in complex porous media and in large representative volume elements (RVEs). To this end, we developed an electro-thermal model at the pore scale to study Joule heating effects in large heterogeneous systems with different microstructures. This was achieved using the level set method to implicitly delineate distinct regions within the domain, and an embedded boundary method to facilitate heat exchange across the fluid-solid interface. Moreover, we applied this method to investigate unsteady non-linear electro-thermal effects in non-woven fibrous graphite conductors for RVEs with characteristic lengths of 2 mm, with different fiber orientations, porosity (80% – 90%) and fiber diameters (10 – 20µm). The coupled equations were solved numerically and they produced peak temperatures greater than 2000 K resulting in heating rates as high as 80,000 K/s. Moreover, the results depended strongly on the microstructure of the fiber skeleton and current density. Geometries with large fibers (∼ 20µm) had the highest average and peak temperatures with the mean temperature increasing by 3.9 % while the peak temperature increased by 9.9 %. Anisotropic domains on the other hand had the lowest mean and peak temperatures with peak and mean temperatures of 2293 K and 1437.7K respectively representing a corresponding 12.1% and 5.1% drop in the temperatures. An increase in porosity from 80% to 90%, however, led to an increase in the peak temperature by 5.1%.

Joule heating↗
DOE OSTI2026

Data‐Efficient Generation of Synthetic Microstructures of Polymer‐Bonded Energetic Material With Fine‐Tuned Stable Diffusion

Among current deep learning approaches for synthetic image generation, diffusion-based models stand out in terms of algorithmic stability and ability to retain high-fidelity image features with detailed resolution. Here, in this work, we employ Dreambooth, a method for fine-tuning Stable Diffusion, on X-ray CT images of microstructure of the polymer-bonded form (PBX) of a commonly used high explosive, Pentaerythritol tetranitrate (PETN), which yields generative models for creating synthetic PBX images. The models developed here represent five classes (or ‘lots’) of microstructures and demonstrate successful generation of images of each class with high fidelity, as verified by computed classification accuracy of ∼ 94% or higher. Data augmentation afforded by such image synthesis can be used to more reliably decipher underlying statistics, build processing-structure correlations, recognize off-normal structural anomalies, and identify age-related changes. Ideas related to converting image data into appropriate density mapping and performing mesoscale simulation or surrogate modeling of detonation are also discussed.

Dreambooth↗
DOE OSTI2022

Nanotwin assisted reversible formation of low angle grain boundary upon reciprocating shear load

Severe plastic deformation of metals is known to lead to superior properties that cannot be achieved by any traditional metallurgic process. Origin of the superior properties is perceived to be closely associated with grain refinement, a fundamental process during the severe plastic deformation, which is essentially the formation of new grain boundaries. However, the atomistic mechanism of grain boundary formation remains largely obscure. Here, by using in-situ transmission electron microscopy and molecular dynamic simulation, we reveal, for the first time at atomic level, the shear-induced low-angle grain boundary (LAGB) formation processes in Au nanocrystal. In this work, we discover the LAGB formation is accomplished through inward propagation of nanotwins accompanied by dislocations gliding on the twin boundaries, a nanotwin-mediated dislocation slip mechanism, which shows reversible characteristic under reciprocating shear load and is affected by the nanocrystal microstructure and orientation. Our result unveils unprecedented atomistic insights on shear driven grain refinement towards nanostructure of superior properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Crystal plasticity modeling of strain-induced martensitic transformations to predict strain rate and temperature sensitive behavior of 304 L steels: Applications to tension, compression, torsion, and impact

This paper advances crystallographically-based Olson-Cohen (direct γ → α’) and deformation mechanism (indirect γ→ε→α’) phase transformation models for predicting strain-induced austenite to martensite transformation. Here, the advanced transformation models enable predictions of not only strain-path sensitive, but also of strain-rate and temperature sensitive deformation of polycrystalline stainless steels (SSs). The deformation of constituent grains in SSs is modeled as a combination of anisotropic elasticity, crystallographic slip, and phase transformation, while the hardening is based on the evolution of dislocation density and explicit shifts in phase fractions. Such grain-scale deformation is implemented within the meso-scale elasto-plastic self-consistent (EPSC) homogenization model, which is coupled with the implicit finite element (FE) method to provide a constitutive response at each FE integration point for solving boundary value problems at the macro-scale. Parameters pertaining to the hardening and transformation models within FEEPSC are calibrated and validated on a suite of data including flow curves and phase fractions for monotonic compression, tension, and torsion as a function of strain-rate and temperature for wrought and additively manufactured (AM) SS304L. To illustrate the potential and accuracy of the integrated multi-level FE-EPSC simulation framework, geometry, mechanical response, phase fractions, and texture evolution are simulated during gas-gun impact deformation of a cylinder and quasi-static tension of a notched specimen made of AM SS304L. Details of the simulation framework, comparison between experimental and simulation results, and insights from the results are presented and discussed.

304L steels↗
DOE OSTI2025

Bubble Transport through a Porous Lattice with an Applied Inlet Flow

Within gas-evolving electrochemical systems, bubbles negatively impact performance by covering electrode active sites for reactions, blocking electric field lines, and obstructing liquid electrolyte flow causing pressure buildup. Recent additive manufacturing advances have enabled tuned porous electrode microstructures to be created, but producing systems that maximize electrochemical throughput and minimize bubble impact remains challenging. Thus, improved physical understanding of and modeling capabilities for bubble behavior are critical to improve electrolyzer design. To address this need, this study examines rising stage bubbles within a lattice with an applied liquid flow—an underexplored regime that strongly influences an electrochemical bubble’s fate. Notably, theoretical predictions and resolved bubble simulations are complemented by experiments from a 3D-printed visualization cell that matches the simulation geometry. The minimum threshold flow rate to achieve bubble breakthrough is found to be larger for higher porosities and for smaller bubbles. Different-sized bubbles decrease expected electrochemical performance in different ways; smaller bubbles tend to stay stuck but cover less solid surface, while larger bubbles more readily break through but cover more surface while in the lattice. The bubble trajectory, deformation, and contact area provide insight into these different behaviors. These findings provide design guidelines toward creating more effective electrolyzers.

Guo, Jack [Lawrence Livermore National Laboratory ↗
DOE OSTI2026

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗
DOE OSTI2026

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2022

In situ study on heavy ion irradiation induced microstructure evolution in single crystal Cu with nanovoids at elevated temperature

Understanding the response of nanovoids to irradiation damage advances our capability to design radiation tolerant materials. Recent in situ studies reveal the shrinkage of hexagon-shaped nanovoids in (110) textured Cu. However, the evolution of voids with irregular geometries under irradiation is less well understood. Here, in situ Kr ion irradiation was performed at 100 °C on single-crystal Cu (112) that possesses nanovoids with high aspect ratio. In situ studies show these elongated voids fragment into smaller voids and shrink gradually with increasing dose. Phase field simulations confirm that fragmentation of the void is governed by the competing kinetics between atoms diffusing toward and away from the elongated nanovoid. In conclusion, cost irradiation analysis reveals the formation of high-density defect clusters.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Mechanisms of Three-Dimensional Solid-Phase Epitaxial Crystallization of Strontium Titanate

Strontium titanate (SrTiO 3 , STO) is a complex metal oxide with a cubic perovskite crystal structure. Due to its easily described and understood crystal structure in the cubic phase, STO is an ideal model system for exploring the mechanistic details of solid-phase epitaxy (SPE) in complex oxides. SPE is a crystallization approach that aims to guide crystal growth at low homologous temperatures to achieve targeted microstructures. Beyond planar thin films, SPE can also exploit the addition of a chemically inert, noncrystallizing, amorphous obstacle in the path of crystallization to generate complex three-dimensional structures. The introduction of this mask fundamentally alters the SPE process, inducing a transition from two- to three-dimensional geometries and from vertical to lateral crystal growth under the influence of the crystal/mask/amorphous boundary. Using a combination of molecular dynamics simulations and experiments, we identify several unique phenomena in the nanoscale growth behaviors in both conventional (unmasked) and masked SPE. Examining conventional SPE of STO, we find that crystallization at the interface is strongly correlated to, and potentially driven by, density fluctuations in the region of the amorphous STO near the crystalline/amorphous interface with a strong facet dependence. In the masked case, we find that the crystalline growth front becomes nonplanar near contact with the mask. We also observe a minimum vertical growth requirement prior to lateral crystallization. Both phenomena depend on the relative bulk and interfacial free energies of the three-phase (crystal/mask/amorphous) system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Experimental Insights into the Coupling of Methane Combustion and Steam Reforming in a Catalytic Plate Reactor in Transient Mode

The microstructured reactor concept is very promising technology to develop a compact reformer for distributed hydrogen generation. In this work, a catalytic plate reactor (CPR) is developed and investigated for the coupling of methane combustion (MC) and methane steam reforming (MSR) over Pt/Al 2 O 3 -coated microchannels in cocurrent and counter-current modes in transient experiments during start-up. A three-dimensional (3D) computational fluid dynamics (CFD) simulation shows uniform velocity and pressure distribution profiles in microchannels. For a channel velocity from 5.1 to 57.3 m/s in the combustor, the oxidation of methane is complete and self-sustainable without explosion, blow-off, or extinction; nevertheless, flashbacks are observed in counter-current mode. In the reformer, the maximum methane conversion is 84.9% in cocurrent mode, slightly higher than that of 80.2% in counter-current mode at a residence time of 33 ms, but at the cost of three times higher energy input in the combustor operating at ~1000 °C. Nitric oxide (NO) is not identified in combustion products, but nitrous oxide (N 2 O) is a function of coupling mode and forms significantly in cocurrent mode. Furthermore, this research would be helpful to establish the start-up strategy and environmental impact of compact reformers on a small scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Residual Stress Measurements of Cr-coated Cladding

Accident-tolerant fuel concepts for light-water reactor applications have been developed and tested in diverse research programs around the world. Industry teams have developed coated cladding concepts using different coating application methods, such as physical vapor deposition and coldspraying. Recently, variations in thickness and application method were investigated, and their effects on the microstructure and mechanical properties were reported. The results suggest an increase in performance in an unirradiated state, however, it remains unclear what phenomenon drove this performance. One explanation to this performance is associated with the coating process imparting a residual stress. Residual stress analysis was performed on coatings produced with various process parameters to understand how a thin layer of Cr impacts the mechanical properties after coating application and after simulated operating conditions. The gathered information is intended to address safety concerns or support margin identification related to accident-tolerant fuels and to support subsequent modeling and simulation efforts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Neutron thermalization in nuclear graphite: A modern story of a classic moderator

Graphite has been of great interest since the dawn of the nuclear era. In this work, we will focus on the thermal neutron scattering of two types of nuclear graphite that have different microstructures, namely, the historic PGA which has a very anisotropic medium-grained needle coke, and the new G347A which is a near-isotropic fine-grained graphite. We report the results of room temperature inelastic neutron scattering measurements for the scattering function, double differential scattering cross-section, and phonon density of states. The obtained results are discussed in conjunction with results extracted from first-principles lattice dynamics simulations of crystalline graphite, using the incoherent approximation and by also including the coherent-one phonon contribution. Excellent agreements have been found between the measured and calculated scattering functions and phonon densities of states. The inclusion of the coherent-one phonon scattering shows a significant improvement of the calculated scattering functions, double differentials, and inelastic scattering cross-sections. We also highlight the physical concerns about the current thermal neutron scattering ENDF/B-VIII of nuclear graphite library and its corresponding phonon density of states. This work aims to provide a new presentation of neutron thermalization by graphite and offers a new set of high quality experimental and computational data for consideration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

In situ deformation of antigorite-olivine two-phase mixtures: Implications for dynamics and seismic anisotropy in the mantle wedge

Water released from hydrous minerals in subducting slabs reacts with the overlying plate, resulting in widespread serpentinization in the mantle wedge. Deformation of serpentinized peridotites has been invoked to explain forearc seismic anisotropy, yet studies of the mechanical properties and deformation behaviors of serpentine-bearing multiphase aggregates remain limited. Here we deformed olivine-antigorite mixtures containing 70, 50, and 20 vol.% of antigorite at 2.5 – 7.6 GPa, 673 K and strain rates of ∼10 –5 –10 –4 s –1 . Elasto-viscoplastic self-consistent simulations, constrained by synchrotron X-ray diffraction (XRD) data, were used to estimate lattice strain, stress–strain partitioning, crystallographic preferred orientations (CPO) development, and aggregate strength. Selected run products were also analyzed by electron backscatter diffraction for comparison with the CPO results obtained from XRD experiments. We found olivine transitions from A- or B-type to C-type when antigorite fraction drops to 20 vol.%, coinciding with a microstructural change from interconnected weak layers to a load-bearing framework (LBF). An additional run on a sample Atg50/Ol50 with preexisting microstructures suggested the formation of LBF was promoted by these microstructures, although the preexisting antigorite CPO has been overprinted at 20.8 % strain and could be erased completely by subsequent deformation in nature. Estimated viscosity of the two-phase mixtures suggests that low-degree serpentinization (≤20 %) in the mantle wedge may increase the strength of olivine-rich peridotite and hinder slab-mantle decoupling, whereas high-degree serpentinization (≥50–70 %) weakens the peridotite and favors decoupling if sufficient viscosity contrast (>10) develops. Seismic anisotropy shows a nonlinear dependence on antigorite fraction: antigorite CPO governs the anisotropy of the mixtures with ≥50 vol.% antigorite, whereas olivine CPO dominates at low fractions (∼20 vol.%). The presence of pre-existing microstructures reduces seismic anisotropy of the deformed mixtures, however the persistence of pre-existing CPO in actively subducting slabs remains uncertain, making their significance over geological timescales questionable.

Crystallographic preferred orientation↗
DOE OSTI2025

Molecular-resolution imaging of ice crystallized from liquid water by cryogenic liquid-cell TEM

Despite the ubiquity of ice, a molecular-resolution image of nanoscopic defects or microstructures in ice crystallized from liquid water has never been obtained. This is mainly due to the difficulties in preparing and preserving crystalline ice samples that can survive under high-resolution imaging conditions. Here, we report the stabilization and Å-resolution electron imaging of ice Ih crystallized from liquid water by developing cryogenic liquid-cell transmission electron microscopy (CRYOLIC-TEM). We combine lattice mapping with molecular dynamics simulations to reveal that ice formation is highly tolerant to nanoscale defects such as misoriented subdomains and trapped gas bubbles, which are stabilized by molecular-scale structural motifs. Importantly, bubble surfaces adopt low-energy nanofacets and create negligible strain fields in the surrounding crystal. These bubbles can dynamically nucleate, grow, migrate, dissolve, and coalesce under electron irradiation and be monitored in situ near a steady state. This work improves our understanding of water crystallization behaviors at a molecular spatial resolution.

Materials science↗
DOE OSTI2025

High Energy Density Physics of Inertial Confinement Fusion Ablator Materials (Final Technical Report)

The goal of this project was to conduct dynamic compression experiments and predictive simulations to reveal the fundamental high-energy-density (HED) physics of amorphous carbon. These results are essential for assessing amorphous carbon as a potential ablator material for next generation of inertial confinement fusion (ICF) capsules. We made significant progress in exploring the HED properties of amorphous carbon through experiments at Omega EP Laser and the European XFEL, in addition to billion-atom, quantum-accurate molecular dynamics (MD) simulations. Through our joint experimental and simulation program, we mapped the phase diagram of amorphous carbon, uncovering its range of metastability and identifying phase transitions to diamond and liquid carbon along the Hugoniot and at higher pressures using double shock compression pathways. Our findings indicate that amorphous carbon melts at significantly lower shock pressures than high-density carbon (diamond). However, nanocrystalline diamond nucleates across a broad range of pressures and temperatures. This emergence of the nanocrystalline microstructure during compression can negatively impact the planarity of the shock front and potentially trigger ablator/fuel mixing during Inertial Fusion Energy (IFE) applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗