Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion behavior”

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

Microstructure and mechanical behavior of laser remelted amorphous Al-Ni-La welds

To explore the feasibility of using laser powder bed fusion systems to print bulk amorphous Al alloys, a single autogenous weld produced by laser remelting on a cast Al-5La-9Ni (at.%) alloy was studied for its microstructure and mechanical behavior. The solidification rates experienced by material within the welds were high enough to produce entirely amorphous regions within the welds. Welds were characterized using SEM, STEM, and APT and were found to contain Ni-rich amorphous clusters. Micropillar compression tests were used to assess mechanical properties of the welds and found that all of the amorphous regions exhibited yield strengths around 1 GPa. In conclusion, these results are discussed in the context of diffusivity of Ni and rare earth elements in liquid, glass forming ability, thermodynamic driving force for phase formation during solidification, free volume concentrations in bulk metallic glasses, and cluster-related softening.

APT↗
DOE OSTI2025

Quantum Monte Carlo Benchmarking of Molecular Adsorption on Graphene-Supported Single Pt Atom

The precise understanding of adsorption energetics and molecular geometry at catalytic sites is fundamental for advancing catalysis, particularly under the constraints of resource efficiency and environmental sustainability. Here, this study benchmarks the performance of density functional theory (DFT) calculations against diffusion Monte Carlo (DMC) calculations for adsorption properties of small gas molecules relevant to CO oxidation—namely O 2 , CO, CO 2 , and atomic oxygen—on a single Pt atom supported by pristine graphene. Our findings reveal that DMC calculations provide a significantly different landscape of adsorption energetics compared to DFT results. Notably, DFT predicts different lowest-energy configurations and spin states, particularly for O 2 , which suggests potential discrepancies in predicting the catalytic behavior. Furthermore, this study identifies the critical issue of CO poisoning, highlighted by the large disparity between the DMC adsorption energies of O 2 (−1.23(2) eV) and CO (−3.37(1) eV), which can inhibit the catalytic process. These results emphasize the necessity for more sophisticated computational approaches in catalysis research, aiming to refine the prediction accuracy of reaction mechanisms and to enhance the design of more effective catalysts.

Ahn, Jeonghwan [University of Illinois at Urbana-C↗
DOE OSTI2024

Regioisomeric Spirobifluorene CANAL Ladder Polymers and Their Gas Separation Performance

For this work, we synthesized and characterized two isomeric microporous hydrocarbon ladder polymers from catalytic arene norbornene annulation (CANAL) of regioisomeric bis-norbornene fused spirobifluorenes, where the ladder chains are connected either through the same fluorene unit or across two different fluorene units in spirobifluorene. This pair of isomeric polymers was used to investigate the effect of ladder macromolecular structures on the microporosity and transport properties. Both polymers form mechanically intact films with thermal stability up to 480 °C and relatively high BET surface areas. The polymer formed from 2,7-dibromospirobifluorene showed higher BET surface area and higher gas permeability than the polymer from 2,2'-dibromospirobifluorene despite similar intersegmental spacing as indicated by X-ray scattering. The aging behavior for both polymers followed the same trend as the previously reported CANAL-fluorene polymers, with dramatically increased permselectivities over time, resulting in gas separation performance above the 2008 upper bounds for H 2 /CH 4 and O 2 /N 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Chromium segregation-induced oxide evolution in Ni-10Cr alloys during high-temperature oxidation

The oxidation behavior of a Ni-10(wt%)Cr alloy under high-temperature O 2 conditions is investigated using transmission electron microscopy and first-principles calculations. Results reveal that chromium segregation plays a central role in driving the evolution of complex oxide phase structures during oxidation. At low Cr concentrations, Cr preferentially segregates to NiO grain boundaries or internal pores, substituting for Ni atoms and forming Ni(Cr)O solid solutions. As Cr content increases, enhanced diffusion promotes Cr penetration into the NiO lattice, leading to the formation of multiphase oxide structures. First-principles modeling corroborates these findings: at low Cr concentrations, Cr atoms favor surface and grain-boundary segregation, while higher concentrations lead to Cr aggregation within the NiO bulk. Furthermore, the integrated experimental-theoretical approach provides atomistic insights into Cr-mediated mass transport mechanisms during alloy oxidation and offers valuable guidance for controlling oxide growth kinetics and phase stability in Ni-Cr alloys, with implications for improving oxidation resistance in high-temperature structural applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High-Mass Loading of Flower-like Ni-MoS 2 microspheres Toward Efficient Intercalation pseudocapacitive Electrode

This work reports the exploration of intercalation pseudocapacitance in a thicker electrode of flowerlike Ni-doped MoS 2 microspheres that features a mass loading of ~10 mg/cm 2 without sacrificing the gravimetric capacitance (~425 F/g at 5 mV/s). Integration of Ni atoms into MoS 2 microspheres not only stabilized the structural integrity but also ameliorated the rapid intercalation and deintercalation of electrolyte ions even at a commercial-level mass loading. The energy instability by Ni doping significantly changed the local bonding behavior and the overall electronic structure of MoS 2 , facilitating the breaking of the MoS 2 layer and generation of more active edge sites, which are responsible for faster reaction kinetics. The experiments attribute the overall capacitance enhancement in (Mo-Ni)S 2 to the increased rate of electrolyte ion insertion and extraction, which is confirmed by b-values close to 0.5, at different potentials, indicating that the current response predominantly depends on the diffusive mechanism for both MoS 2 and Ni-MoS 2 thicker electrodes. The symmetric device constructed with Ni-MoS 2 microspheres exhibited a capacitance value of 101 F/g in 1 mV/s, for which the energy density is 9 Wh/kg, as well as attained an outstanding cycling stability of 10000 cycles with 60% retention at 2 A/g. In addition to providing insights into the development of 2D TMDs, this work explores the design of robust and highly efficient intercalation electrode material for electrochemical energy storage devices.

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

Probing the electronic structure and spectroscopy of pyrrolyl and imidazolyl radicals using high-resolution photoelectron imaging of cryogenically cooled anions

High-resolution photoelectron imaging and photodetachment spectroscopy of cryogenically cooled pyrrolide and imidazolide anions are used to probe the electronic structure and spectroscopy of pyrrolyl and imidazolyl radicals. The high-resolution data allow the ground state vibronic structures of the two radicals to be completely resolved, yielding accurate electron affinities of 2.1433 ± 0.0008 eV and 2.6046 ± 0.0006 eV for pyrrolyl and imidazolyl radicals, respectively. Fundamental frequencies for eight vibrational modes of pyrrolyl and ten vibrational modes of imidazolyl are measured, including several nonsymmetric Franck–Condon-forbidden modes. Two electronic excited states are also observed for the two radicals, displaying diffuse spectral features in both systems. The observations of nonsymmetric vibrational modes in the ground states and the diffuse excited state features provide strong evidence for vibronic couplings between the ground state and the two close-by excited states. The 2-pyrrolide isomer is also observed as a minor species from the electrospray ionization source and the electron affinity of 2-pyrrolyl is measured to be 1.6690 ± 0.0030 eV along with five vibrational frequencies. Even though the HOMOs of both pyrrolide and imidazolide anions are p orbitals, photodetachment spectroscopy reveals completely different threshold behaviors for the two anions: a d-wave-dominated spectrum for pyrrolide and an s-wave-dominated spectrum for imidazolide. Furthermore, the current study provides a wealth of electronic and spectroscopic information, which is ideal to compare with more accurate vibronic coupling calculations for these two important radicals, as well as interesting information about the photodetachment dynamics of the two anions.

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

Martensitic transformation in superlattices of two non-transforming metals

A number of studies have shown that the coherent integration of a second non-transforming phase in martensitic materials results in unprecedented thermo-mechanical properties. We demonstrate the possibility of martensitic transformation by epitaxially combining two non-transforming materials, Mg (hexagonal closed packed) and MgLi (body centered cubic). Density functional theory calculations predict that MgLi/Mg superlattices undergo a transformation crystallographically equivalent to MgSc alloys and can be expected to result in finite temperature shape memory behavior. We predict the zero-temperature energy of the martensite phase to be 9.3 meV/atom lower than austenite for a 50–50 at. % MgLi/Mg nanolaminate, which is expected to result in a martensitic transformation temperature of approximately 150 K. This energy difference can be tuned between 4.9 and 18.1 meV/atom by varying the fraction of the pure Mg phase in the superlattice, potentially enabling the transformation from cryogenic to room and higher temperatures. While Li diffusion might limit the applications of the systems proposed, this work opens the door to a new approach to create martensitic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A Theory for the Balance between Warm Rain and Ice Crystal Processes of Precipitation in Mixed-Phase Clouds

Abstract Mixed-phase clouds contain both supercooled cloud liquid and ice crystals. In principle, precipitation may be initiated either by the liquid phase or by the ice phase. Ice crystals may grow by vapor diffusion to become snow (“ice crystal process”), forming “cold” precipitation. Equally, cloud droplets, when large enough, coalesce to form “warm” precipitation by the “warm rain process.” Warm rain could be supercooled and freeze as “warm” graupel. In the present paper, a new simplified theoretical analysis is provided to examine the microphysical system consisting of three species of hydrometeor, namely, cloud liquid, “cold ice” (crystals, snow), and “warm rain” (frozen or supercooled). This is obtained by nondimensionalizing and simplifying the evolution equations for the mass of each species. Analytical formulas are given for equilibria. Feedback analysis shows that the sign of the feedback is linked to the abundance of precipitation, with a neutral surface in the 3D phase space. The system’s precipitation amount explodes while in the initial unstable regime, crossing the neutral surface and approaching the equilibrium point that is a stable attractor. Positive and negative feedbacks are elucidated. In a standard case, the cold ice mass is about 1000 times larger than the warm rain mass. To illustrate the physical behavior of the theory, sensitivity tests are performed with respect to environmental conditions (e.g., aerosol, updraft speed) and microphysical parameters (e.g., riming and sedimentation rates for cold ice). Cold ice prevails, especially in fast ascent, due to its low bulk density, favoring slow sedimentation and a wide cross-sectional area for riming. Significance Statement The theory elucidates how the ice phase can prevail in the precipitation from any mixed-phase clouds with supercooled cloud liquid and crystals. The ice phase radically suppresses cloud liquid by riming when active and “wins” the competition against coalescence. This prevalence of ice is shown to arise from the low bulk density of snow. The cloud is viewed as a system of negative and positive feedbacks that prevail in realms of stability and instability in a 3D phase space.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

In Situ Oxidation Studies of High-Entropy Alloy Nanoparticles

Although high-entropy alloys (HEAs) have shown tremendous potential for elevated temperature, anticorrosion, and catalysis applications, little is known on how HEA materials behave under complex service environments. Herein, we studied the high-temperature oxidation behavior of Fe 0.28 Co 0.21 Ni 0.20 Cu 0.08 Pt 0.23 HEA nanoparticles (NPs) in an atmospheric pressure dry air environment by in situ gas-cell transmission electron microscopy. It is found that the oxidation of HEA NPs is governed by Kirkendall effects with logarithmic oxidation rates rather than parabolic as predicted by Wagner's theory. Further, the HEA NPs are found to oxidize at a significantly slower rate compared to monometallic NPs. The outward diffusion of transition metals and formation of disordered oxide layer are observed in real time and confirmed through analytical energy dispersive spectroscopy, and electron energy loss spectroscopy characterizations. Localized ordered lattices are identified in the oxide, suggesting the formation of Fe 2 O 3 , CoO, NiO, and CuO crystallites in an overall disordered matrix. Hybrid Monte Carlo and molecular dynamics simulations based on first-principles energies and forces support these findings and show that the oxidation drives surface segregation of Fe, Co, Ni, and Cu, while Pt stays in the core region. The present work offers key insights into how HEA NPs behave under high-temperature oxidizing environment and sheds light on future design of highly stable alloys under complex service conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

VARI3D & PERSENT: Perturbation and Sensitivity Analysis (Revision 5)

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry-based transport code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Probing phonon focusing, thermomechanical behavior, and moiré patterns in van der Waals architectures using surface acoustic waves

Surface acoustic waves (SAWs) propagate along solid-air, solid-liquid, and solid-solid interfaces. Their characteristics depend on the elastic properties of the solid. Combining transmission electron microscopy (TEM) experiments with molecular dynamics (MD) simulations, we probe atomic environments around intrinsic defects that generate SAWs in vertically stacked two-dimensional (2D) bilayers of MoS 2 . Our joint experimental-simulation study provides insights into SAW-induced structural and dynamical changes and thermomechanical responses of MoS 2 bilayers. Using MD simulations, we compute mechanical properties from the SAW velocity and thermal conductivity from thermal diffusion of SAWs. The results for Young’s modulus and thermal conductivity of an MoS 2 monolayer are in good agreement with experiments. The presence of defects, such as nanopores which generate SAWs, reduces the thermal conductivity of 2D-MoS 2 by an order of magnitude. We also observe dramatic changes in moiré patterns, phonon focusing, and cuspidal structures on 2D-MoS 2 layers.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Interfacial Hydrogen-Bond Dynamics in Transition Metal Compounds

Understanding how water behaves when confined within atomic layers of active transition-metal carbides, nitrides, and carbonitrides is essential for uncovering the fundamental principles needed to engineer solid–liquid interfaces at the atomic scale. Yet, how lattice element chemistry and surface termination groups collectively regulate the structure and mobility of such interlayer water remains poorly understood. Here, we present a composition-controlled investigation of interlayer water dynamics in layered transition-metal nitride, carbide, and carbonitride systems using a systematic integration of quasi-elastic neutron scattering (QENS), ab initio molecular dynamics (AIMD) simulations, and density functional theory (DFT) calculations. QENS measurements show that nitride-rich systems host mobile, translationally diffusing water with thermally activated self-diffusion coefficients on the order of 10 –10 m 2 s –1 , whereas mixed C/N lattices confine water to localized, nontranslational motion that is insensitive to temperature. AIMD and DFT reveal that lattice C/N chemistry and surface functional group composition reshape the first hydration layer by modulating the surface electronic structure and termination-dependent hydrogen-bond networks, leading to pronounced differences in water ordering and thermal resilience. On the other hand, fully carbide systems exhibit intermediate behavior, highlighting that water mobility is not primarily controlled by the hydration level alone but by the coupling between lattice composition and surface chemistry. Overall, this study establishes how surface chemistry and lattice composition jointly control interfacial hydrogen bond dynamics, offering a mechanistic framework for designing transition-metal layered materials with tailored interfacial transport properties.

Hydration↗
DOE OSTI2025

Dynamical Signatures of Thermotoga maritima Maltose-Binding Proteins Affected by Ligand Binding

Functional segregation among protein isoforms depends on the interplay of their overall structures and the molecular dynamics of these structures. Thermotoga maritima maltose-binding protein (tmMBP) isoforms show size-dependent differential binding of maltose and malto-oligomers while maintaining remarkable fold conservation. This differential behavior needs detailed characterization in native-like aqueous conditions to understand the effects of protein dynamics on ligand binding and recognition. Small-angle neutron scattering (SANS), neutron spin echo (NSE) spectroscopy, and dynamic light scattering (DLS) were used in conjunction with previously published computational molecular dynamics (MD) simulations to understand the dynamic behavior of tmMBPs experimentally. SANS provided information on the overall structure of the molecules, while NSE was used to determine the dynamics in the nanosecond time scale. Both tmMBP2 and tmMBP3 have a bidomain architecture linked with a flexible hinge, with the binding pocket sitting in the cleft between the two domains. tmMBP2 and tmMBP3 showed different solution dynamics, with the translational and rotational components dominating the dynamics of both systems, resulting in a clear differentiation of their diffusion pattern. A faster dynamics component was also observed and was attributed to segmental dynamics. Differences observed between the ligand-free (apo) and ligand-bound (holo) states of the two proteins are attributed to conformational entropy. Our results highlight the intricacies of how structure and dynamics can together shape binding to a repertoire of substrates in structurally similar proteins.

Diffusion↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
DOE OSTI2025

Porosity, swelling, and composition evolution in high-burnup monolithic U-Mo fuel

The microstructural progression of very high-burnup (>8 × 10 21 fissions/cm 3 ) monolithic uranium-molybdenum (U-10wt.%Mo) was analyzed, providing crucial insights into the behavior of post-recrystallized nuclear fuel, where scant data exists. Three focused ion beam cuboids sourced from a fuel plate with varying local burnups of 8.86 × 10 21 , 9.05 × 10 21 , and 9.36 × 10 21 fissions/cm 3 were characterized. The porosity and composition of the samples were evaluated to characterize the evolution of the microstructure as a function of fission density and different locations on the fuel plate, while simultaneously isolating plate-specific parameters such as Zr diffusion barrier thickness, hot-isostatic-pressing conditions, enrichment, and reactor conditions. The porosity was segmented, and the three-dimensional distribution of the porosity was extrapolated from the two-dimensional segmentation. The composition was assessed and quantified using energy-dispersive X-ray spectroscopy areal mapping. The porosity fraction increased as a function of the burnup from 27.77±0.51, 35.12±1.54, and 37.71±0.44 % for 8.86 × 10 21 , 9.05 × 10 21 , and 9.36 × 10 21 fissions/cm 3 , respectively. When compared to literature, the porosity volume fraction plateaus at burnups greater than 6 × 10 21 fissions/cm 3 , while the pore size grows linearly as a function of fission density. The number of large pores increased in number density as a function of burnup, while the smallest pores (<0.3 µm) increased up to 9.05 × 10 21 fissions/cm 3 , followed by a decrease at 9.36 × 10 21 fissions/cm 3 . The delamination and cracking in the fuel plate propagated through an interconnected porosity sublayer identified ∼5 µm from the diffusion barrier. The local swelling of the specimens was within or near the prediction bounds of the Robinson-Williams model for local swelling. The fission products, strontium, barium, cerium, and cesium, precipitated into the pores, while neodymium accumulated adjacent to the pores. Furthermore, these findings have direct implications for the development of fuel performance codes and the accurate documentation of the microstructure evolution in high burnup U-Mo, thus enhancing the safety and efficiency of nuclear fuel usage.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗