Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective diffusivity”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Puzzling 𝐵⁡(𝐸⁢2;0 + → 2 + ) strength in the proton dripline nucleus 36 Ca

Recent measurements of the 𝐸⁢2 transition rate from the ground state to the first 2 + excited state of the proton dripline nucleus 36 Ca show an unusual pattern when compared to its isotopic neighbor 38 Ca: despite having a higher 𝐸 𝑥 ⁡(2$^{+}_{1}$) excitation energy, the 𝐵⁡(𝐸⁢2;0$^{+}_{1}$ → 2$^{+}_{1}$) rate in 36 Ca is larger. The question that naturally arises is to what extent this observation can be attributed to the unbound character of the 2$^{+}_{1}$ state. To understand the influence of the continuum space on the low-energy properties of 36 Ca, we carried out Gamow shell model calculations that can account for the continuum coupling effects associated with the occupation of unbound 𝑓⁡𝑝 shells. We found that in the threshold 2 + state, 36 Ca is spatially diffused, which impacts the observed 𝐵⁡(𝐸⁢2) trend.

20 ≤ A ≤ 38↗
DOE OSTI2022

3D Printed Carbon Aerogels for Polymer-Electrolyte Fuel Cells

Carbon aerogels (CA) are known for ultra-low density, very high surface area, fine porous structure, and excellent electrical conductivity. Recent advances in high resolution 3D printing at LLNL enable the manufacturing of complex structured CA electrodes which have been demonstrated in electrolysis, flow battery, and supercapacitor applications. Polymer-electrolyte fuel cells (PEFC) may be an excellent candidate for these materials/methods as designers seek to optimize membrane-electrode assemblies to achieve higher efficiencies. A major challenge to PEFC design is two-phase flow due to water production which causes pore flooding and leads to cell instability and degradation. To mitigate these effects, groups have investigated the creation of liquid specific pathways in the porous layers to alleviate liquid buildup and enhance gas diffusion. While some progress has been made, traditionally manufactured GDLs, which consist of randomly distributed carbon fibers, are limiting. 3D printed CAs are a promising solution as complex macroporous electrode structures can be manufactured out of micro and nano porous filaments with a high level of control. Additionally, the process may allow for the consolidation of multiple layers (MPL, GDL, & flow channels) into a single printed architecture. In this work various CA lattice structures were printed using a direct ink write (DIW) method and liquid breakthrough pressures were measured to demonstrate the ability to tune the water pathways to achieve specific pressure thresholds. Furthermore, a discussion of the prospects for printed CA electrodes in PEFC will focus on material properties, synthesis techniques, and current challenges.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tusqh

SAND2025-00675O Tusqh is a software tool that generates cubical meshes in 2D and 3D and computes the homology of these meshes using persistent homology. It includes a grid cell in the output if its volume-fraction is above a selectable threshold, estimated by sampling points within the cell. Tusqh incorporates anti-aliasing algorithms to mitigate grid orientation and scale effects. It is designed for creating finite element meshes for simulations and can be used in various applications such as heat diffusion, mechanical simulations, and computer graphics rendering. The software outputs meshes in an open format compatible with downstream software. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗
DOE OSTI2026

Formation of Non-Doped Cubic Lithium Lanthanum Zirconium Oxide Nanofibers: Insights from In Situ Synchrotron X-Ray Scattering

This study investigates the formation mechanism of non-doped cubic lithium lanthanum zirconium oxide (c-LLZO) nanofibers using in situ synchrotron X-ray scattering techniques. Electrospun polymer precursor nanofibers were annealed at temperatures up to 800 °C, enabling real-time tracking of phase transitions via simultaneous small-angle X-ray scattering (SAXS), wide-angle X-ray scattering (WAXS), and evolved CO 2 gas analysis. The results reveal a three-step transformation pathway: polymer decomposition, formation of La 2 Zr 2 O 7 (LZO), and direct conversion of LZO to c-LLZO without intermediate tetragonal phases detected within the sensitivity of our in situ WAXS measurement. Cryo-electron energy loss spectroscopy (EELS) further elucidates the role of lithium diffusion, showing Li enrichment at fiber surfaces and Li deficiency in the interior, which stabilizes the cubic phase. This Li segregation effect in nanostructured LLZO materials extends beyond the previously reported size effect. This work advances the understanding of c-LLZO formation mechanisms and provides practical insights for optimizing synthesis routes to achieve phase-pure c-LLZO for solid-state battery applications.

LLZO phase stability↗
DOE OSTI2023

Molecular Simulations of the Chain Length Dependent Adsorption of C7‐C14 n ‐Alkanes in ZIF‐8

Recent experiments show that the diffusivities of C8 and C10 n-alkanes in ZIF-8 are higher than those of C7 and C9, respectively. We investigated this unusual ‘odd-even' effect by simulating the adsorption of C7-C14 n-alkanes in ZIF-8 using hybrid Monte Carlo molecular simulations. The resultant adsorption isotherms, guest-host energies, isosteric heats, and chain length distributions are analyzed for trends among the n-alkanes. ZIF-8 cages filled with n-alkanes are characterized using a combination of image processing and data science to quantify the differences in packing that occur with chain length. In conclusion, results indicate that packing changes drastically from C7 to C8 and from C9 to C10, which is consistent with the diffusion trends previously reported in experiments.

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

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Energy-momentum-conserving stochastic differential equations and algorithms for the nonlinear Landau-Fokker-Planck equation

Coulomb collision is a fundamental diffusion process in plasmas that can be described by the Landau-Fokker-Planck (LFP) equation or the stochastic differential equation (SDE). While energy and momentum are conserved exactly in the LFP equation, they are conserved only on average by the conventional corresponding SDEs, suggesting that the underlying stochastic process may not be well defined by such SDEs. Here, in this study, we derive new SDEs with exact energy-momentum conservation for the Coulomb collision by factorizing the collective effect of field particles into individual particles and enforcing Newton's third law. These SDEs, when interpreted in the Stratonovich sense, have a particularly simple form that represents pure diffusion between particles without drag. To demonstrate that the new SDEs correspond to the LFP equation, we develop numerical algorithms that converge to the SDEs and preserve discrete conservation laws. Simulation results are presented in a benchmark of various relaxation processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Precise Linker Length and Dynamic Bond Exchange Control Penetrant Diffusion in Dense Vitrimers

Polymer networks with dynamic covalent bonds have been investigated for their self-healing ability, recyclability, and potential as more sustainable materials. Recent results have indicated that in some cases, bond exchange can enhance the transport of penetrants in dense networks, pointing to their potential for separations of membranes. Here, imine dynamic bonds in ethylene oxide (EO) networks with precise linker lengths were synthesized to investigate the transport of N,N′-bis(2,5-di-tert-butylphenyl)-3,4,9,10-perylenedicarboximide (BTBP), a large, anisotropic dye molecule. Networks with mesh sizes smaller than, comparable to, and greater than the size of the penetrant axes were investigated to probe the effects of bond exchange and network confinement on transport. Mesh sizes, which ranged from 0.5 to 1.62 nm, were determined from shear rheology, glass transitions by calorimetry, and probe diffusion coefficients by fluorescence recovery after photobleaching. Permanent networks with identical EO chain lengths were prepared as control samples, and up to a 3 orders of magnitude increase in diffusion coefficient is observed in the dynamic systems for short linkers containing 13 backbone atoms. The longest linkers with 71 backbone atoms show no difference between the permanent and dynamic networks. Linkers shorter than 11 backbone atoms, corresponding to a mesh size smaller than the penetrant small axis, diffusion is no longer observable on the experimental time scale, indicating a sharp cutoff attributed to the precise linkers and narrow mesh size distribution. The dynamic imine exchange time scales were compared to the diffusive hopping times of penetrants and indicate that exchange can occur during a diffusive displacement. Furthermore, these findings provide insights into the factors affecting penetrant transport in dense polymers and inspire the development of next-generation selective polymer membranes.

Diffusion↗
DOE OSTI2023

Mechanical behavior of polyamide-6 after combined photo-oxidative and hygrothermal aging

The effect of moisture on the photo-oxidative degradation of polyamide-6 (PA-6) was studied by analyzing the mechanical response after two different accelerated aging procedures. In the first aging procedure, the PA-6 was only exposed to ultra-violet (UV) radiation at 60 °C. In the second procedure, the same duration of UV radiation was periodically interrupted while the relative humidity was raised to 100%. Diffusion-limited and nominally homogeneous degradation conditions were investigated using bulk and film specimens, respectively. Accelerated UV aging reduced the ductility of PA-6, but the additional hygrothermal exposure had no effect on the ductility or strength, indicating that humidity did not influence the photo-oxidation of PA-6. In conclusion, this finding contrasts with previous studies that found thermo-oxidation of PA-6 was accelerated by moisture.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Jet wake from linearized hydrodynamics

We explore how to improve the hybrid model description of the particles originating from the wake that a jet produced in a heavy ion collision leaves in the droplet of quark-gluon plasma (QGP) through which it propagates, using linearized hydrodynamics on a background Bjorken flow. Jet energy and momentum loss described by the hybrid model become currents sourcing linearized hydrodynamics. By solving the linearized hydrodynamic equations numerically, we investigate the development of the wake in the dynamically evolving droplet of QGP, study the effect of viscosity, scrutinize energy-momentum conservation, and check the validity of the linear approximation. We find that linearized hydrodynamics works better in the viscous case because diffusive modes damp the energy-momentum perturbation produced by the jet. We calculate the distribution of particles produced from the jet wake by using the Cooper-Frye prescription and find that both the transverse momentum spectrum and the distribution of particles in azimuthal angle are similar in shape in linearized hydrodynamics and in the hybrid model. Their normalizations are different because the momentum-rapidity distribution in the linearized hydrodynamics analysis is more spread out, due to sound modes. Since the Bjorken flow has no transverse expansion, we explore the effect of transverse flow by using local boosts to add it into the Cooper-Frye formula. After including the effects of transverse flow in this way, the transverse momentum spectrum becomes harder: more particles with transverse momenta bigger than 2 GeV are produced than in the hybrid model. Although we defer implementing this analysis in a jet Monte Carlo, as would be needed to make quantitative comparisons to data, we gain a qualitative sense of how the jet wake may modify jet observables by computing proxies for two example observables: the lost energy recovered in a cone of varying open angle, and the fragmentation function. We find that linearized hydrodynamics with transverse flow effects added improves the description of the jet wake in the hybrid model in just the way that comparison to data indicates is needed. Our study illuminates a path to improving the description of the wake in the hybrid model, highlighting the need to take into account the effects of both transverse flow and the broadening of the energy-momentum perturbation in spacetime rapidity on particle production.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Manipulating topological Hall-like signatures by interface engineering in epitaxial ruthenate/manganite heterostructures

Topologically protected non-trivial spin textures (e.g. skyrmions) give rise to a novel phenomenon called the topological Hall effect (THE) and have promising implications in future energy-efficient nanoelectronic and spintronic devices. Here, we have studied the Hall effect in SrRuO 3 /La 0.42 Ca 0.58 MnO 3 (SRO/LCMO) bilayers. Our investigation suggests that pure SRO has hard and soft magnetic characteristics but the anomalous Hall effect (AHE) in SRO is governed by the high coercivity phase. Here, we have shown that the proximity effect of a soft magnetic LCMO on SRO plays a critical role in interfacial magnetic coupling and transport properties in SRO. Upon reducing the SRO thickness in the bilayer, the proximity effect becomes the dominant feature, enhancing the magnitude and temperature range of THE-like signatures. The THE-like features in bilayers can be explained by a diffusive Berry phase transition model in the presence of an emergent magnetic state due to interface coupling. This work provides an alternative understanding of THE-like signatures and their manipulation in SRO-based heterostructures, bilayers and superlattices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Effect of thermal oxidation on helium implanted pure iron

The impact of thermal oxidation on Helium (He) implanted pure iron (Fe) was investigated to experimentally evaluate how thermal oxidation influences the diffusion and distribution of He within the material. In case of the sample with the lowest dose (2 × 10 17 ions/cm 2 ), the thinnest oxide layer was observed compared to the non-implanted pure Fe. It is due to He nano bubbles implanted in the material, which hinder the diffusion of Fe ions necessary for forming the oxide layer. For the medium dose sample (5 × 10 17 ions/cm 2 ), the oxide layer was slightly thicker than that of the lowest dose, however it had more porous structure due to the numerous nano He bubbles. In addition, a large bubble region around 150–200 nm depth was generated. Here, the highest dose sample (1 × 10 18 ions/cm 2 ) was observed that the escape of He was accelerated by the large amount of nano He bubbles from within the metal, forming a porous Fe matrix on the top surface and an oxide layer that is very porous but thicker than those of lower doses.

He implantation↗
DOE OSTI2025

An unexpectedly shrunken bandgap in V2O5 nanoparticles

Synchrotron x-ray spectroscopy was employed to determine the effects of nanostructuring on electronic band structure in V2O5, a promising cathode material and widely used catalyst. V2O5nanoparticle and bulk powders were characterized via P-XRD, electron microscopy, and diffuse reflectance ultraviolet/visible/near-infrared spectroscopy to confirm the optical bandgap. X-ray emission spectroscopy revealed the nanoparticle valence band O 2pstates to be upshifted relative to the bulk, while x-ray absorption spectroscopy and resonant inelastic x-ray scattering showed the lowest V 3dconduction band states to be static. Together, these changes (in conjunction with an increased density of unoccupied lower conduction band states) produce a shrunken bandgap in the V2O5nanoparticles that defies the Burstein-Moss effect. Changes in nanoparticle band structure are generally attributed to oxygen vacancy defects. While nanostructure bandgap reduction is in line with much previous computational work, it is unexpected from most previous experimental results. To our knowledge, this is the first synchrotron x-ray spectroscopy study of a shrunken bandgap achieved in pure V2O5nanoparticles.

Anquillare, E L↗
DOE OSTI2025

100 kHz OH PLIF of High-Pressure H2/NH3 and H2/CH4 Flames

The reaction zones in a multi-element, multi-stage, micromix injector are investigated using planar laser-induced fluorescence (PLIF) imaging of the hydroxyl (OH) radical at a 100 kHz repetition rate. The combustor is operated at a pressure of ∼10 bar. Multiple fuel blends of hydrogen, methane, and ammonia are considered to study the effect of changing fuel composition on the flame structure. Compact regions of high OH signal are observed in each flame for fuel mixtures with high hydrogen content, indicating rapid reactant consumption in the shear layers. Instantaneous images exhibit wrinkled flame fronts and the formation of unreacted pockets due to highly turbulent burning. The addition of methane is observed to reduce the reactivity and flame speed, resulting in more diffuse regions of OH between each flame. The shape of the reactant core is altered, and the length scale of the flame surface wrinkling is increased. Ammonia addition has a more severe effect on the distribution of OH and the local flame structure than methane addition. A time series of images shows the presence of OH vortical structures and significant turbulent fluctuations in the flames near the wall, suggesting the presence of a recirculation zone that could drive the observed disparity in flame lengths.

Hodge, Alexander J. [Purdue University]↗
DOE OSTI2022

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

The effect of stress on the migration of He gas bubbles under a thermal gradient in Fe by phase-field modeling

Here a phase-field model is parameterized to study the effect of elastic stresses on the migration of He gas bubbles in Fe under a temperature gradient. Stresses caused by the gas bubble pressure and residual stress in the Fe matrix are considered. The dependence of He bubble migration velocity on the magnitude of the residual stress, average temperature, temperature gradient, and bubble size is measured. In agreement with a theoretical model based on surface diffusion, simulation results demonstrate that He bubbles move towards the high temperature region with velocities in Fe that are orders of magnitude faster than previously reported in UO 2 . It is found that local stresses in the matrix caused by the He bubble have negligible effect on the bubble migration process; however, residual stresses in the Fe matrix, potentially caused by processing or irradiation, can modestly modify bubble kinetics through pressure dependence of the He diffusion coefficients. Compressive residual stress decreases diffusion coefficients for bulk and surface diffusion mechanisms, thus reducing the migration velocity of the gas bubble. In contrast, tensile residual stress increases the diffusion coefficients, resulting in an increase in the gas bubble migration velocity. This pressure dependence is also consistent with a theoretical model. This phase field model lays the foundation for analysis of bubble coalescence-induced fracture in He bubble-containing steels.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evaluation of Sb-Nd and Te-Nd phases within the U-Zr fuel matrix and their interactions with HT9 alloy

Antimony (Sb) and tellurium (Te) were investigated as potential additives for U-10Zr (wt.%) metallic fuel to limit the fuel-cladding chemical interaction (FCCI) with HT-9 alloy. Neodymium (Nd) was utilized to simulate the formation of lanthanide-based solid fission products which are known to play a detrimental role in FCCI. Fuel alloys of U-Zr-Sb-Nd and U-Zr-Te-Nd were evaluated in their annealed condition and compared against their as-cast conditions. Isothermal diffusion couple experiments were performed between U-Zr-Nd, U-Zr-Sb-Nd, and U-Zr-Te-Nd against the cladding alloy HT9 to evaluate the effectiveness of the additives to stabilize Nd within the fuel alloys, as well as investigate the interaction regions that form between the different fuel alloys and HT9. Further, SbNd and Sb 3 Nd 4 , and TeNd are found to be the primary neodymium-based phases formed in the U-Zr-Sb-Nd and U-Zr-Te-Nd alloys, respectively. The zirconium-based phase, Zr 2 Sb, is also found to form within the former alloy. All phases were found to remain stable through the diffusion experiments and exhibited no interaction with HT9 constituent elements. Preferential interaction between Nd with additivities Te and Sb compared to constituting elements in HT9 was further verified based on density functional theory (DFT) calculated enthalpy of mixing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗