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.

SourceResearch connectionFrom yearThrough year

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

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Dusty Gas Model for Solid Oxide Fuel Cell Fuel Electrode

This model applies the Dusty Gas Model simulate multi-component species transport in SOFC (solid oxide fuel cell) anodes which considers the pressure gradient across the fuel electrode. This studyhas been verified with the analytical solution for different fuel electrode thicknesses and with literature values. The model was developed using the VoronoiFVM platform in Julia which is a built in implicit and semi implicit solver that integrates electrochemical behavior, microstructural effects, and transient analysis for accurate prediction of species transport under varying conditions.

dusty gas model (DGM)↗

DOE OSTI2021

Columnar Grain-driven Plasticity and Cracking in Nanotwinned FCC Metals

The mechanisms of strengthening and plasticity in columnar-grained metals with preferentially oriented nano-sized twins have been examined traditionally by considering dislocation processes, but rarely from the perspective of grain boundary (GB) deformation. Here, the effects of GB strain accommodation on plastic deformation in four different columnar-grained nanocrystalline nanotwinned (nt) face-centered-cubic metals (Cu, Ag, Al, and Ni) were studied by large-scale molecular dynamics simulations. It is observed that in tensile deformation parallel to coherent twin boundaries (CTBs), the dislocation mechanisms in each metal are identical and associated with GB emissions of jog and threading dislocations at small and large CTB spacings, respectively. However, CTB strengthening effects are increasingly more pronounced in columnar-grained nt metals as their shear modulus increases, which is rationalized by the dependence of GB stress concentrations on twin size, metal type and strain rate. Also, while flow stresses in nt-Cu, nt-Ag, and nt-Al metals increase linearly with decreasing CTB spacing, a maximum strength limit is reached in nt-Ni below a critical CTB spacing of 6 nm. The strength limit in nt-Ni results from columnar GB cracking induced by prominent GB sliding. For columnar-grained microstructures, GB sliding is equivalent in nt-Ag and nt-Al and slightly lower in nt-Cu but markedly higher in nt-Ni. These findings underscore the importance of new GB deformation mechanisms on plasticity and fracture in columnar-grained nt metals and enrich our understanding of CTB strengthening in fcc metals synthesized in the literature.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Coupled THM modeling of bentonite heating and hydration in tank tests with a new temperature-dependent water retention model

This study presents a coupled thermo-hydro-mechanical (THM) model for simulating the heating and hydration behavior of bentonite, a buffer material in deep geological repositories (DGRs). The model incorporates a new temperature-dependent soil water retention curve which captures the thermal-induced shift in water retention behavior. It also distinguishes between liquid and gas permeability, modeling intrinsic gas permeability as a function of accessible porosity to improve vapor transport and desaturation predictions. The model was validated against two large-scale tank tests, demonstrating good agreement with measured temperature, relative humidity, and water inflow data. It revealed a complex porosity evolution driven by thermal expansion, vapor movement, vapor condensation, and hydration-induced swelling during heating and hydration processes. The simulation results also suggest that the permeability of the hydration layer plays a critical role in controlling water intake. Clogging of this layer can significantly reduce the volume of water inflow during the hydration phase. Furthermore, while the model effectively captures key THM behavior, further development of the mechanical constitutive law is required to account for possible thermo-elasto-plastic volume changes and microstructural effects. Overall, the model provides a robust tool for evaluating the evolution of bentonite-based barrier material in DGRs.

Guo, Guanlong [Lawrence Berkeley National Laborato↗

DOE OSTI2026

Origins of Enhanced Ion Transport in Nanostructured Anion-Conducting Polyelectrolytes

Ion-conducting polymer chemistry and microstructure profoundly impact membrane water uptake and ionic conductivity. Water uptake strongly impacts ionic conductivity; yet excess water uptake compromises ion-exchange membrane mechanical properties. Although nanophase separation has been proposed to overcome this trade-off, it is unclear how polymer backbone architecture governs ionic nanostructure and its subsequent impact on water uptake and conductivity. Here, we integrate experiments and molecular dynamics simulations to elucidate the role of backbone chemistry in governing ionic nanostructure, hydration behavior, and ion transport in anion-conducting polyelectrolytes (ACPs). We systematically investigate hydrocarbon polynorbornene (PNB)-based ACPs with three distinct backbone architectures: vinyl-addition polymerization (VAP), ring-opening metathesis polymerization (ROMP), and hydrogenated ROMP. While maintaining comparable ion exchange capacities (IECs) and identical side-chain chemistry, we isolate the effects of backbone structure. We show that nanophase-separated ionic nanostructures originate in the dry state and evolve upon hydration through heterogeneous water uptake, with water preferentially partitioning into ion-rich domains. This nanophase separation arises from a delicate interplay between ionic segregation propensity and the entropic barrier imposed by backbone stiffness. Specifically, flexible backbones intensify attractive ion–ion interactions by reducing the entropic penalty for backbone deformation, promoting nanophase separation, while rigid backbones suppress ionic nanostructure formation. Nanophase-separated ion domains locally concentrate water upon hydration, which in turn enables the connectivity required for fast transport at lower water concentration values. Furthermore, these findings demonstrate that backbone chemistry can be tuned as a design lever to promote nanophase separation and enhance ion transport without excessive water uptake.

Anions↗

DOE OSTI2021

Relativistically transparent magnetic filaments: scaling laws, initial results and prospects for strong-field QED studies

Abstract Relativistic transparency enables volumetric laser interaction with overdense plasmas and direct laser acceleration of electrons to relativistic velocities. The dense electron current generates a magnetic filament with field strength of the order of the laser amplitude (>10 5 T). The magnetic filament traps the electrons radially, enabling efficient acceleration and conversion of laser energy into MeV photons by electron oscillations in the filament. The use of microstructured targets stabilizes the hosing instabilities associated with relativistically transparent interactions, resulting in robust and repeatable production of this phenomenon. Analytical scaling laws are derived to describe the radiated photon spectrum and energy from the magnetic filament phenomenon in terms of the laser intensity, focal radius, pulse duration, and the plasma density. These scaling laws are compared to 3D particle-in-cell (PIC) simulations, demonstrating agreement over two regimes of focal radius. Preliminary experiments to study this phenomenon at moderate intensity ( a 0 ∼ 30) were performed on the Texas Petawatt Laser. Experimental signatures of the magnetic filament phenomenon are observed in the electron and photon spectra recorded in a subset of these experiments that is consistent with the experimental design, analytical scaling and 3D PIC simulations. Implications for future experimental campaigns are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Size-Dependent Nucleation in Crystal Phase Transition from Machine Learning Metadynamics

In this Letter, we present a framework that combines machine learning potential (MLP) and metadynamics to investigate solid-solid phase transition. Here, based on the spectral descriptors and neural networks regression, we develop a scalable MLP model to warrant an accurate interpolation of the energy surface where two phases coexist. Applying it to the simulation of B4–B1 phase transition of GaN under 50 GPa with different model sizes, we observe sequential change of the phase transition mechanism from collective modes to nucleation and growths. When the size is at or below 128 000 atoms, the nucleation and growth appear to follow a preferred direction. At larger sizes, the nuclei occur at multiple sites simultaneously and grow to microstructures by passing the critical size. The observed change of the atomistic mechanism manifests the importance of statistical sampling with large system size in phase transition modeling.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Reactive Carbide‐Based Synthesis and Microstructure of NASICON Sodium Metal All Solid‐State Electrolyte

Reactive carbide precursor-based synthesis of NASICON-type NZSP (Na 1+x Zr 2 Si x P 3-x O 12 ) solid-state electrolyte (SSE) is demonstrated, in contrast to the established oxide-based approach. Exothermic decomposition of ZrC and SiC in air homogenizes microstructure, yielding 98% compact density after conventional sintering at 1200 °C. Quantitative stereology demonstrates that significant microstructural differences are present. Compacts of carbide-derived Carb-NZSP are 98% dense with a secondary zirconium oxide (ZrO 2 ) volume fraction of 0.2% ± 0.3%, versus 93% dense and 3% ± 1% for oxide-derived baseline. For Carb-NZSP, the secondary glassy phosphate phase is agglomerated, while for baseline, it is dispersed and percolated. Electrochemical testing combined with post-mortem analysis demonstrates how microstructural control of secondary phases is critical for dendrite suppression: Carb-NZSP critical current density (CCD) is 3.1 ± 0.8 mA cm −2 at 0.1 mAh cm −2 , versus 1.0 ± 0.7 mA cm −2 at 0.1 mAh cm −2 . Cryogenic focused ion beam (cryo-FIB) analysis demonstrates that in both materials, the porous 2D sheet-like sodium metal dendrites propagate around and subsume NZSP grains, likely following a path enriched with glassy phase and with porosity. Dendrites also flow around isolated zirconia particles. Phase field simulation reveals deflection of dendrites by mechanically tough zirconia, while brittle glassy phase accelerates dendrite growth, especially when finely distributed.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Texture and microstructure evolution in thermomechanically processed Mg-Ca and Mg-Zn-Ca alloys

Mg-Zn-Ca alloys have the potential for producing Mg alloy sheets with weaker basal textures and therefore improved formability. Thermomechanical processing (TMP) also plays a substantial role in determining the final sheet texture. The interplay of these variables, TMP and alloying, complicates comparisons of alloys across the literature. This work systematically explores the texture evolution and recrystallization behavior in Mg-Ca and Mg-Zn-Ca alloys during plane strain compression (PSC) using a Gleeble thermomechanical simulator. It is demonstrated that the basal texture intensity and texture characteristics change significantly during post-deformation annealing, particularly in the ternary alloys. It is also shown that careful selection of the TMP processing variables used during PSC is essential to producing weak textures. In particular, it is important to limit the recrystallization that occurs between compressive hits. This is achieved by adjusting solute content, strain rate, and the duration of the soak between passes.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Intergranular Hotspots: A Molecular Dynamics Study on the Influence of Compressive and Shear Work

Numerous crystal- and microstructural-level mechanisms are at play in the formation of hotspots, which are known to govern high explosives initiation behavior. Most of these mechanisms, including pore collapse, interfacial friction, and shear banding, involve both compressive and shear work done within the material and have thus far remained difficult to separate. Here, we assess hotspots formed at shocked crystal–crystal interfaces using quasi-1D molecular dynamics simulations that isolate effects due to compression and shear. Two high explosive materials are considered (TATB and PETN) that exhibit distinctly different levels of molecular conformational flexibility and crystal packing anisotropy. Temperature and intramolecular strain energy localization in the hotspot are assessed through parametric variation of the crystal orientation and two velocity components that respectively modulate compression and shear work. The resulting hotspots are found to be highly localized to a region within 5–20 nm of the crystal–crystal interface. Compressive work plays a considerably larger role in localizing temperature and intramolecular strain energy for both materials and all crystal orientations considered. Shear induces a moderate increase in energy localization relative to unsheared cases only for relatively weak compressive shock pressures of approximately 10 GPa. These results help isolate and rank the relative importance of hotspot generation mechanisms and are anticipated to guide the treatment of crystal–crystal interfaces in coarse-grained models of polycrystalline high-explosive materials.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Design of high-pressure iron Rayleigh–Taylor strength experiments for the National Ignition Facility

Iron is an important metal, scientifically and technologically. It is a common metal on Earth, forming the main constituent of the planet's inner core, where it is believed to be in solid state at high pressure and high temperature. It is also the main component of many important structural materials used in quasistatic and dynamic conditions. Laser-driven Rayleigh–Taylor instability provides a means of probing material strength at high pressure and high temperature. The unavoidable phase transition in iron at relatively low pressure induces microstructural changes that ultimately affect its strength in this extreme regime. This inevitable progression can make it difficult to design experiments and understand their results. Here, we address this challenge with the introduction of a new approach: a direct-drive design for Rayleigh–Taylor strength experiments capable of reaching up to 400 GPa over a broad range of temperatures. We use 1D and 2D hydrodynamic simulations to optimize target components and laser pulse shape to induce the phase transition and compress the iron to high pressure and high temperature. At the simulated pressure–temperature state of 350 GPa and 4000 K, we predict a ripple growth factor of 3–10 depending on the strength with minimal sensitivity to the equation of state model used. The growth factor is the primary observable, and the measured value will be compared to simulations to enable the extraction of the strength under these conditions. These experiments conducted at high-energy laser facilities will provide a unique way to study an important metal.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Strengthening of tungsten by coherent rhenium precipitates formed during low fluence irradiation

Experimental data show that the accumulation of rhenium and osmium from transmutation reactions severely affect the microstructural evolution and property degradation of tungsten-based materials under neutron irradiation. Theory and modeling have confirmed that Re atom transport in W is by irradiation-produced migrating self-interstitial atoms. With this diffusion mode in operation, a specific microstructure evolution is realized when at relatively low neutron fluence the Re-rich precipitates are formed, while the void and interstitial loop population development is suppressed, affecting the mechanical properties. Here, this research shows the effect of small coherent Re-rich precipitates on the dislocation glide under stress, investigated using the molecular dynamics approach with empirical interatomic potentials. The results are compared with an earlier simulation of void hardening in W. It is demonstrated that small coherent Re-rich precipitates of less than 6 nm diameter represent relatively weak obstacles for moving edge dislocations. The implication of these results on the interpretation of experimental results is discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Measurement of Transport Properties of Woody Biomass Feedstock Particles Before and After Pyrolysis by Numerical Analysis of X-Ray Tomographic Reconstructions

Lignocellulosic biomass has a complex, species-specific microstructure that governs heat and mass transport during conversion processes. A quantitative understanding of the evolution of pore size and structure is critical to optimize conversion processes for biofuel and bio-based chemical production. Further, improving our understanding of the microstructure of biochar coproduct will accelerate development of its myriad applications. This work quantitatively compares the microstructural features and the anisotropic permeabilities of two woody feedstocks, red oak and Douglas fir, using X-ray computed tomography (XCT) before and after the feedstocks are subjected to pyrolysis. Quantitative analysis of the three-dimensional (3D) reconstructions allows for direct calculations of void fractions, pore size distributions and tortuosity factors. Next, 3D images are imported into an immersed boundary based finite volume solver to simulate gas flow through the porous structure and to directly calculate the principal permeabilities along longitudinal, radial, and tangential directions. The permeabilities of native biomass are seen to differ by three to four orders of magnitude in the different principal directions, but we find that this anisotropy is substantially reduced in the biochar formed during pyrolysis. The quantitative transport properties reported here enhance the ability of pyrolysis simulations to account for feedstock-specific effects and thereby provide a useful touchstone for the biorefining community.

09 BIOMASS FUELS↗

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.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

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

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↗

NASA NTRS2017

GEER Status Update

History: The Glenn Extreme Environments Rig (GEER) first became operational in the early part of 2015. Since that time GEER has completed a number of scientific tests and has undergone improvements in the chemical delivery system and analytics following a year of operations experience. Recent Updates: In June 2016, the GEER process system was rebuilt to provide a more robust system, higher accuracy and new capabilities. New insulation was installed on the exterior of the pressure vessel and gas lines. The newly revamped GEER plumbing system can provide extremely precise custom gas mixtures using any gas desired by the investigator in any combination. GEER can heat the resulting mixture up to 500 deg C and 1500 psia. The 304 stainless steel vessel walls were polished to reduce corrosion rate and reduce unwanted chemical reactions. The process lines were replaced with high purity Sulfinert coated tubing. The GEER team added the ability to individually boost specialty gases to GEER thus allowing operators to make very precise changes to the gas chemistry inside of GEER during a test while at high temperature and pressure. High accuracy mass flow meters were added to further improve gas mixing accuracy and precision. An in-line, integrated Inficon MicroGC Fusion was added for real time gas analysis along with a high purity gas sampling system, providing fully automated, real time analysis of the gas chemistry inside of GEER in minutes. This complements a co-located mass spectrometer and both are used for regular monitoring of the vessel chemistry. All internal vessel components were replaced with polished 304SS equivalents. Hot vent down capability was increased. Finally, an automated liquid injection system was added and is rated for max vessel operating conditions of (1500 psia, 500 C). Recent Results and Publications: In May 2016, GEER completed a test that exposed high temperature electronics to Venus surface conditions for 21.5 days. This demonstrated the potential for operating robotic spacecraft in the Venus environment without the need for thermal or environmental protection. Results from this test were published in December 2016 and received national media attention. In April 2017 GEER implemented an 80 day test at Venus surface conditions to simulate chemical weathering of expected Venus minerals. This test supported a ROSES award to a team led by Prof. Ralph Harvey of Case Western Reserve University. The test concluded in July 2017 and nearly doubled previous operation record of 42 days at Venus surface conditions. Preliminary results of these and previous experiments were presented at the recent Venus Modeling Workshop. In June 2017, NASA TM2017-219437 "Chemical and Microstructural Changes in Metallic and Ceramic Materials Exposed to Venusian Surface Conditions" was published. This report provides an extensive and valuable resource detailing the behavior of a variety of engineering materials at Venus surface conditions. Community Involvement: An external science advisory panel has been formed.

Kremic, Tibor↗

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↗