Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion process”

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

Remarks on the thesis of Aitor Bengoechea Fernández: Generalization of Doppler effect processing using S(α,β) thermal diffusion laws

The thesis titled “Generalization of Doppler effect processing using S(α,β) thermal diffusion laws” addresses a critical gap in the modeling of neutron scattering phenomena in nuclear materials, particularly where crystal lattice effects intersect with low-energy resonances. Conducted within the Doctoral School ED 352 at Aix-Marseille Université and hosted at CEA Cadarache, the research presents a theoretical and numerical framework for extending the applicability of Thermal Scattering Laws (TSLs) using higher-order correlation functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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

Parallel diffusion operator for magnetized plasmas with improved spectral fidelity

Diffusive transport processes in magnetized plasmas are highly anisotropic, with fast parallel transport along the magnetic field lines sometimes faster than perpendicular transport by orders of magnitude. This constitutes a major challenge for describing non-grid-aligned magnetic structures in Eulerian (grid-based) simulations. Here, the present paper describes and validates a new method for parallel diffusion in magnetized plasmas based on the anti-symmetry representation [Halpern and Waltz, Phys. Plasmas 25, 060703 (2018)]. In the anti-symmetry formalism, diffusion manifests as a flow operator involving the logarithmic derivative of the transported quantity. Qualitative plane wave analysis shows that the new operator naturally yields better discrete spectral resolution compared to its conventional counterpart. Numerical simulations comparing the new method against existing finite difference methods are carried out, showing significant improvement. In particular, we find that combining anti-symmetry with finite differences in diagonally staggered grids essentially eliminates the so-called “artificial numerical diffusion” that affects conventional finite difference and finite volume methods.

Anisotropic diffusion↗
DOE OSTI2025

Bayesian learning with Gaussian processes for low-dimensional representations of time-dependent nonlinear systems

This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. Furthermore, we demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.

Data-driven model reduction↗
DOE OSTI2025

Structural evolution and magnetic hardness of (Sm,Zr)(Fe,Co,Ti) 12 alloy particles via reduction-diffusion

The quest for achieving high coercivity in Sm(Fe,Co,Ti) 12 alloys, despite their inherent strong magnetocrystalline anisotropy, has posed significant challenges. Recently, (Sm,Zr)(Fe,Co,Ti) 12 monocrystalline particles have exhibited coercivity μ 0 H c > 1.2 T, showcasing promising prospects and significant potential for both manufacturing and research endeavors. This study delves into the structural evolution of (Sm,Zr)(Fe,Co,Ti) 12 (1:12) alloy particles made via the calciothermic reduction-diffusion synthesis process as influenced by the molar ratios of Ca atoms to O 2- ions (Ca/O), annealing time and annealing temperature. Critical insight that informs conditions to optimize the magnetic response is gained via systematic experimentation and advanced electron microscopy. Complex structural features, including core-shell morphologies and intricate multiphase compositions within individual particles, are unveiled. An optimal Ca/O ratio of 1.30 produces particles with a coercivity up to μ 0 H c = 1.63 T, while higher Ca/O ratios induce the formation of a Sm-rich TbCu 7 -type (1: 7 ) phase, which only partially transforms into the desired 1:12 phase during annealing. Persistent remnants of the 1:7 phase locally impact atomic structure, particle morphology, and coercivity. Furthermore, these findings underscore the complex interplay between synthesis parameters, resulting structures, and magnetic properties, informing the design and optimization of high-performance permanent magnets comprised of the (1:12) compound.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Radiolytic degradation of 240 Plutonium and 242 Plutonium oxalates

Raman, FTIR, and diffuse reflectance spectroscopy were used to study the auto-radiolytic degradation of 240 Pu and 242 Pu oxalates. The significant differences in the lifetimes of 240 Pu and 242 Pu enabled the differentiation between environmental and radiolytic mechanisms. 240 Pu oxalates were observed to decompose to PuOCO 3 at intermediate times (~ 20 weeks) followed by partial conversion to PuO 2 at times greater than one year. Atmospheric oxidation was shown to be the primary decomposition mechanism for 242 Pu(IV) oxalate, and the alpha radiolysis of aquo and oxalate ligands serves as a secondary decomposition mechanism. In conclusion, this study offers a fresh perspective on radiolytic aging, which is crucial for long-term storage applications.

Analytical Techniques in Art Conservation↗
DOE OSTI2025

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Carbon mineralization pathways in interfacial adsorbed water nanofilms

Carbon mineralization in humidified carbon dioxide offers a promising route to mitigate anthropogenic emissions in a world stressed by water security. Despite its technological importance, our understanding of carbonation in water-poor environments lags, as traditional dissolution-precipitation pathways struggle to explain the adsorbed water nanofilm-mediated reactivity. Here, we utilize in operando X-ray diffraction (XRD) and advanced molecular simulations to investigate nanoconfined reactions driving forsterite carbonation, the magnesium-rich olivine. By examining magnesium ion dissolution and transport in atomistic simulations of the forsterite-water-carbon dioxide interface and comparing these with the in operando XRD activation energies, we identify both processes as rate-limiting at saturation. Our simulations reveal a mechanistic view of interfacial carbonation, where dissolution and precipitation are mediated by anomalous quasi two-dimensional diffusion. The transport process involves intermittent diffusive hopping in the desorbed state, separated by crawling events that are spatially short but temporally long. This understanding transcends carbon mineralization, with implications for understanding the transport of contaminants in geosystems, the design of multifunctional materials, water desalination, and molecular recognition systems.

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

Ultrafast (1‐5 sec) Lamination of Perovskite Solar Cells With Self‐Encapsulation Using Rapid Joule Heating

Perovskite solar cells (PSCs) are traditionally fabricated using sequential layer‐by‐layer deposition, in which each layer of the device is processed on top of the preceding layer. This constrains the processing techniques and selection of transport layer materials that can be used in the solar cell. To overcome these challenges, two half‐cells can be processed independently and then diffusion‐bonded through a lamination process. However, current lamination processes for perovskite solar cells suffer from relatively long process times, which can limit throughput when moving toward high‐volume manufacturing. In this study, a custom platform was designed for rapid‐joule heating of perovskite materials and devices. This enabled more than a 99% reduction in lamination time from 26 min to 1 s. Perovskite samples that were laminated in 1 s exhibited comparable values of percent bonded area, interfacial toughness, grain domain size, and X‐ray diffraction spectra to those laminated in greater than 10 min. As a proof‐of‐concept, 18.3% efficient devices were successfully laminated in 5 s. A transient heat transfer model was developed to describe the relationship between the perovskite temperature and the electrical power supplied to the heaters, establishing a baseline for predicting processing conditions in large‐scale manufacturing systems. Ultra‐fast lamination provides a pathway toward scalable roll‐to‐roll or sheet‐to‐sheet manufacturing of PSCs.

heat transfer↗
DOE OSTI2025

A GPU Accelerated Mixed‐Precision Finite Difference Informed Random Walker (FDiRW) Solver for Strongly Inhomogeneous Diffusion Problems

In nature, many complex multi‐physics coupling problems exhibit significant diffusivity inhomogeneity, where one process occurs several orders of magnitude faster than others temporally. Simulating rapid diffusion alongside slower processes demands intensive computational resources due to the necessity for small time steps. To address these computational challenges, we have developed an efficient numerical solver named Finite Difference informed Random Walker (FDiRW). In this study, we propose a GPU‐accelerated, mixed‐precision configuration for the FDiRW solver to maximize efficiency through GPU multi‐threaded parallel computation and lower precision computation. Numerical evaluation results reveal that the proposed GPU‐accelerated mixed‐precision FDiRW solver can achieve a 117× speedup over the CPU baseline, while an additional 1.75× speedup is achieved by employing lower precision GPU computation. Notably, for large model sizes, the GPU‐accelerated mixed‐precision FDiRW solver demonstrates strong scaling with the number of nodes used in simulation. When simulating radionuclide absorption processes by porous wasteform particles with a medium‐sized model of 192 × 192 × 192, this approach reduces the total computational time to 10 min, enabling the simulation of larger systems with strongly inhomogeneous diffusivity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Irreversible Multielement Diffusion and the Resulting Compositional and Processing Flexibility in the Synthesis and Densification of Lithium Aluminum Lanthanum Zirconium Oxide

Despite the broad interest in Al-doped lithium lanthanum zirconium oxide (LLZO), the wide range of reported ionic conductivities suggests that the effect of processing parameters on the resulting phase purity and bulk conductivity is not well-understood. Here, in this work, we synthesized six separate series of LLZO with variations in Al concentration, Li excess, Li addition order, Li source, densification method, and mother powder to determine the effect on the composition, phase purity, density, and bulk conductivity. We found that a wider range of compositions (6.08–7.61 mol of Li and 0.06–0.23 mol of Al) than previously reported can result in cubic phase stability and that nearly all elements (Li, Al, Zr, and La) are lost to the MgO crucible during calcination and sintering. The manner in which different elements are lost is affected by the processing parameters. We observed that Al-doped LLZO shows great compositional flexibility in stabilizing the cubic phase, offering an explanation for the range of electrochemical performance metrics reported in the literature.

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

Matrix Diffusion Controls Mountain Hillslope Groundwater Ages and Inferred Storage Dynamics

Groundwater age distributions provide fundamental insights on coupled water and biogeochemical processes in mountain watersheds. Field-based studies have found mixtures of young and old-aged groundwater in mountain catchments underlain by bedrock; yet, the processes that dictate these groundwater age distributions are poorly understood. In this work, we use the coupled ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models to simulate groundwater age distributions on a lower montane hillslope in the East River Watershed, Colorado (USA). We develop a convolution-based approach to propagate fracture-matrix diffusion processes to the EcoSLIM advection-dominated age distributions. We compare observed 3 H and 4 He concentrations from two groundwater wells against model predictions that have varying advective transport times and matrix diffusion magnitudes. Based on a Monte Carlo analysis that considers uncertain matrix and fracture parameters, we find that matrix diffusion is needed to jointly predict 3 H and 4 He observations at both wells. The advection-dominated age distributions lack adequate mixing of young and old-aged water to capture the observed co-occurrence of 3 H and 4 He. The model scenario that best matches the 3 H, 4 He, and water level observations when considering both advective flowpath and matrix diffusion mixing processes has a dynamic bedrock groundwater reservoir that is susceptible to considerable storage losses during low-snow periods. This dynamic groundwater system amplifies the need to assimilate deeper bedrock groundwater into watershed hydro-biogeochemical predictions. This work further highlights the importance of considering matrix diffusion when interpreting environmental tracers in bedrock groundwater systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗
DOE OSTI2026

Diffusion and Solvation Dynamics of Ions in Water: Beyond the Brownian Approximation

The coupled dynamics of ions and water molecules in their first hydration shell impact a variety of processes including ion diffusion, selective ion transport in water-filled nanopores, and the kinetics of ion-pairing, ion adsorption, and metal-ligand binding reactions. In this work, we study these coupled dynamics for alkali metals (Li, Na, K, Rb, Cs), alkaline Earth metals (Mg, Ca, Sr, Ba), and chloride through the lens of their dependence on ion isotopic mass. Results are validated against previous measurements of the isotopic mass-dependence of ion diffusion coefficients in water and previous ab initio calculations of ion high-frequency dynamics in water. We find that the vibrational power spectra of ions in water consistently exhibit either two or three peaks, i.e., ions have several rattling frequencies within their solvations shells as previously reported for a subset of the species examined here. These frequencies have different sensititivies to isotopic mass that may serve as signatures of ion solvation processes (such as the tendency of ions to orient their first-shell water molecules) and that also may relate to Hofmeister-like effects including the relative affinity of different metals for ribonucleic acid (RNA).

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

Phase-field modeling of diffusion bonding in 316H stainless steel: Impact of processing conditions on grain morphology and bonding quality

A novel multi-phase, multi-component phase‐field model is presented to study the diffusion bonding of 316H stainless steel. Combined with targeted experimental investigations, this model simulates the bond-growth process and predicts the bonding quality. Unlike previous models, our approach captures the simultaneous evolution of voids and grain structures, while quantifying bonding quality using defined bonding ratio. A comprehensive analysis of bond process control is performed by changing temperature, pressure and surface roughness observing the resulting bond structure, which is consistent with experimental observations and analytical predictions. Temperature is determined to be the dominant factor, with the transition from a flat to a robust bond occurring between 1000 °C and 1050 °C. At the ideal bonding temperature of 1050 °C, a surface roughness exceeding 0.6 μm or an applied stress below 4 MPa results in poor bonding quality. Beyond this, higher pressures and smoother surfaces reduce void size, accelerate void shrinkage, and lead to improved bond integrity. This diffuse-interface model can be extended to other material systems if supplied with appropriate thermodynamic and kinetic data. In conclusion, this makes it an effective modeling platform for optimizing high-temperature diffusion bonding and developing reliable bonded components such as compact heat exchangers.

Diffusion bonding↗
DOE OSTI2025

Unveiling the Role of Critical Impurities in Spent LiFePO 4 Cathodes for Scalable Direct Regeneration

Direct regeneration offers a promising alternative to recycling End-of-Life (EoL) batteries by restoring metal elements and preserving the material structure, yet scaling these technologies to handle practical cathode black mass (CBM) with impurities remains challenging. Here, this study investigates the evolution of impurities, including aluminum (Al), polyvinylidene difluoride (PVDF) binder, and residual carbon (C), during direct recycling of spent LiFePO 4 (LFP) cathodes and their impact on electrochemical performance. Using various ex situ and in situ analyses, it is shown that the formation of lithium fluoride (LiF) during the traditional direct recycling process hinders lithium diffusion and deteriorates the reversible capacity. To address this major challenge, the combination of pH-controlled hydrothermal purification and the two-step sintering process is proposed effectively to regenerate spent LFP cathodes, eliminating the negative effect of Al and fluorine (F) impurities while mitigating any potential impacts of carbon residuals. The regenerated LFP from spent CBM achieves superior performance, retaining 152.5 mAh g −1 at 0.1 C and 133 mAh g −1 at 1 C with 98.7% capacity retention after 200 cycles. This approach is further validated using three distinct waste feedstocks from battery modules, enhancing impurity management and scalability in direct recycling. These findings present a sustainable and economically viable solution for large-scale LFP regeneration.

25 ENERGY STORAGE↗
DOE OSTI2026

Hydrogen uptake in graphite matrix at high temperature

Tritium management is a critical challenge for the next generation of nuclear reactors, such as Fluoride Salt Cooled High Temperature Reactors (FHRs) and High Temperature Gas-cooled Reactors (HTGRs), due to the higher production rate (up to 10,000 times) than conventional Light Water Reactors (LWRs). Graphitic materials employed as moderator, reflector, and fuel pebbles offer a potential pathway for tritium recovery by serving as a sink for tritium. Prediction of uptake capacity under reactor relevant conditions remains a challenge due to a lack of low partial pressure data and significant inter-grade variability of graphite. This study addresses these gaps by providing a comprehensive characterization of hydrogen (as a tritium surrogate) uptake and release behavior in the A3-3 graphite matrix (GM) used in fuel pebbles. Uptake measurements are performed at reactor relevant temperatures of 600- 800 °C, 1-200 Torr hydrogen pressure, and 15-120 min equilibration time, followed by thermal desorption spectroscopy up to 1100 °C. Uptake experiments at different equilibration times demonstrate the role of kinetics in hydrogen uptake, which can be modeled as a diffusion-with-trapping process. In the thermodynamics limit, the Sips adsorption model is shown to capture the uptake in A3-3 GM well. Our campaign provides a set of new results for hydrogen uptake in A3-3, including limiting uptake capacity at 600 °C, apparent diffusion coefficient at 600 °C, and the first estimates of the FHR/HTGR relevant (600 °C, 20 Pa partial pressure) equilibrium uptake capacity and time to saturation. Desorption data highlights a new site for hydrogen uptake, not observed in nuclear graphite, which we attribute to the non-graphitized binder. Using the Kissinger method, we estimate activation energy for release from the desorption peaks, confirming the activation energy for release from the basal planes and providing the first estimate for the activation energy of release from the binder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Computing Nonequilibrium Responses with Score-Shifted Stochastic Differential Equations

Using equilibrium fluctuations to understand the response of a physical system to an externally imposed perturbation is the basis for linear response theory, which is widely used to interpret experiments and shed light on microscopic dynamics. For nonequilibrium systems, perturbations cannot be interpreted simply by monitoring fluctuations in a conjugate observable and general response results rely on path ensemble averaging. Furthermore, these techniques do not apply to perturbations that affect the diffusion tensor in a stochastic system. Here, we introduce an “effective” physical process that represents the diffusion perturbed dynamics and enables accurate calculations of responses to a change in the diffusion. Interestingly, the effective dynamics contain an additional drift involving the instantaneous “score” of the system, and we leverage score matching algorithms to carry out nonequilibrium response calculations on systems for which the exact stationary distribution is unknown.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗