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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Hydrogen charging and desorption from microstructural viewpoint: A method for deconvoluting hydrogen desorption spectra and unveiling the hydrogen-microstructure interaction

Understanding the interaction of hydrogen with microstructural features in metallic materials is crucial for designing hydrogen-resistant alloys. Although thermal desorption spectroscopy (TDS) is widely used for investigating the hydrogen binding behavior of various microstructural features, its application to face-centered cubic (fcc) metals and alloys that exhibit low hydrogen diffusivity is limited due to the lumped TDS desorption signals. This paper shows that, by coupling a Sofronis–McMeeking type hydrogen transport model with a microstructure-informed finite-element model, TDS data can be deconvoluted to reveal the underlying adsorption–diffusion–desorption processes, hydrogen diffusivity, and trap-binding energies. In conclusion, the austenitic steel SS316L in solution-annealed condition is used as a demonstration material, and we focused on investigating the interaction of deuterium (hydrogen isotope) with grain boundaries, which is difficult to investigate from experiments alone but critical for design of alloys for hydrogen infrastructure.

Finite element simulation, Polycrystalline Microst↗

DOE OSTI2025

A computational investigation of high-flux, plate-and-frame membrane modules for industrial carbon capture

In this work, we study the application of membrane-based separation systems for carbon capture, considering plate-and-frame membrane modules. The successful deployment of membrane CO 2 capture system relies on high-performing membranes as well as effective membrane modules that can fully exploit the developed membranes. A plate-and-frame membrane module is especially attractive for CO 2 capture from industrial flue gas due to its lower pressure drop compared to its counterparts such as spiral wound modules and hollow fiber modules. To design better plate-and-frame modules, we investigate their basic unit - a single membrane stack through a combination of computational modeling and experimental investigations. The modeling approach is based on Computational Fluid Dynamics (CFD) to represent a multiphysics problem, including the fluid flow and diffusion processes within a membrane module. We use experimental data collected under different operating conditions to validate the CFD model. Numerical results suggest a good agreement between experiments and model outputs for the CO 2 recovery, CO 2 mole fraction in the retentate and permeate, and stage-cut. The CFD model is able to predict accurately the flow behavior, providing valuable insights on the effects of fluid dynamics on mass transfer of CO 2 . We also carry out a sensitivity analysis to identify the effect of key parameters on the CO 2 recovery and the CO 2 purity of the outlet streams.

CFD simulation↗

DOE OSTI2025

Modeling aerosol bolus inhalations in the human lung with the multiple path particle deposition model: Comparison with experimental data

Existing one-dimensional (1D) models of aerosol dosimetry often ignore mixing mechanisms of inhaled aerosols during their transport in the lung. This mixing or aerosol dispersion results from different physical mechanisms in different regions of the lung. It is a higher order effect, which cannot be directly captured in 1D modeling approaches, and thus is sometimes modeled as a diffusive process. Here, in this study, we improved our recently developed alveolar mixing module incorporated in the multiple path particle dosimetry model (MPPD) to account for flow irreversibility and particle trapping in the alveolar spaces, as well as mixing occurring in the tracheobronchial region. This new version of MPPD was coupled with CFPD-based predictions of aerosol bolus dispersion in the oral airway. The model was used to predict the deposition, dispersion, and mode shift of aerosol bolus inhaled at different penetration depths within the lung for breathing patterns and particle size matching those used in a previous experimental study (Darquenne et al., 2016). Even though a quite simplified approach was used, the computations appear to describe subject-specific and test-specific experimental data reasonably well. The proposed combined dispersion-deposition model can be a useful tool for targeted drug delivery and also for exposure health risk assessment.

MPPD↗

DOE OSTI2025

In Situ Observation of Ion Migration in a Ferroelectric Ionic Conductor Rb-KTP during Thermal Annealing

Ion exchange in Rb-doped KTiOPO 4 has facilitated significant advancements in ferroelectric domain engineering, yet understanding the underlying mechanisms remains in its infancy. We perform time-of-flight secondary ion mass spectrometry analysis on multiple periodically ion-exchanged and periodically poled Rb-doped KTiOPO 4 samples under different temperatures and annealing durations. The results are compared between annealing in air, which involved ex situ annealing before periodic poling, and vacuum annealing conducted in situ after periodic poling. The Rb + diffusion profile after periodic ion exchange forms a tooth-shaped pattern. We show that in situ annealing causes a surface pinning effect on the nonpolar face, limiting Rb + migration along the polar axis at the surface. Once the pinned layer is removed through milling, the underlying Rb + diffusion is distinctively different from the surface. Additionally, the rate of Rb + diffusion during in situ annealing is linear, while the periodic domain structures remain stable after annealing. These results contribute to understanding the ionic diffusion process in a ferroelectric ionic conductor and using ion exchange to tailor the linear and nonlinear properties of KTiOPO 4 .

36 MATERIALS SCIENCE↗

DOE OSTI2025

Chromium versus Aluminum: Impact of Nickel Alloy Composition and Interfacial Kinetics on High-Temperature Passivating Oxide Formation

High-temperature corrosion resistance depends critically on the formation of a passivating surface oxide, which is highly sensitive to alloy composition and structure. Such details often elude experimental investigation, and simplified analytical models fail to provide a truly chemical view of passivating oxide evolution. Here, we explicitly compare the fundamental chemistry of Cr and Al as prototypical passivating elements in Ni alloys by directly simulating competing reaction and diffusion processes within the oxide film using kinetic Monte Carlo and density functional theory. We find that the origin and expression of passivating behavior during early-stage thermal oxidation are qualitatively different between the two alloy systems. Ni–Cr alloys feature a sudden onset of passivation associated with a sharp phase transition upon Cr enrichment that directly couples oxidation kinetics to phase transformation behavior. In contrast, Ni–Al alloys display more continuous oxide phase variation with Al enrichment, ultimately resulting in a lower composition threshold for passivation and a thinner passivating layer. In addition, we elucidate the nonobvious role of metal exchange within the alloy near the oxide boundary, which fundamentally alters film composition and passivating behavior. Furthermore, our results have key implications for engineering improved corrosion-resistant alloys, both in terms of compositional variation and processing.

Alloys↗

DOE OSTI2025

Asymmetric rotations slow down diffusion under confinement

Translation and rotation are the two most fundamental forms of diffusion, yet their coupling mechanism is not clear, especially under confinement. Here, we provided evidence of the coupling between rotation and translation using a substituted benzene molecule as an example. A counterintuitive behavior was observed where the movement of the smaller molecule with an asymmetric shape was unexpectedly slower than the larger one with a symmetric shape in confined channels of zeolite. We showed that this diffusion behavior was caused by the presence of the specific and selective interaction of the asymmetric guest with the pores, which increased the local restricted residence time, thus inhibiting the translation under confinement, as further confirmed by dynamic breakthrough curves, uptake measurements, quasi-elastic neutron scattering, and 2 H solid-state NMR techniques. Our work correlated asymmetric rotation and diffusion under a confined environment, which enriched our understanding of the coupling between rotation and translation and could shed light on a fundamental understanding of the diffusion process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Neutron spin echo spectroscopy

Neutron spin echo (NSE) spectroscopy is a powerful technique used to probe the internal dynamics and diffusion processes within matter, allowing researchers to quantify how materials respond to changes in external conditions. To fully understand the bulk properties of a material, it is often crucial to investigate dynamic processes occurring on nanometre length scales and nanosecond timescales. NSE spectroscopy provides an exceptional capability to measure such dynamics across a broad spectrum of condensed matter systems, ranging from proteins and quantum magnets to catalysts. In this Primer, we present an overview of NSE spectroscopy, highlighting the instruments used, methods of data collection and representative applications in both soft-matter and hard-matter sciences.

Czakkel, Orsolya [Institut Max von Laue-Paul Lange↗

DOE OSTI2025

The origin of the Stokes–Einstein relation in simple dense liquids

Here, we investigate the origin of the universal relation between structural relaxation and diffusion in simple dense liquids, known as the Stokes–Einstein (SE) relation. The fact that this relation, originally derived from a hydrodynamic model of a macroscopic particle in a viscous medium, can describe the microscopic-scale liquid dynamics still eludes understanding. We introduce a new universal measure of structural relaxation in a system of N identical particles based on an explicit decomposition of the configuration space into N! congruent convex polyhedra. This measure makes it possible to quantify the correlation between two distinct particle configurations in terms of their minimal Euclidean distance, optimized with respect to particle permutations. Using this measure alongside a model of independent random walkers under the single-occupancy constraint, we derive a master equation that quantifies the SE relation. It allows us to demonstrate that the universal relation between structural relaxation and diffusion in simple dense liquids is caused by two conditions: (a) the confinement of the dominant density fluctuations to the first coordination shell, manifested by de Gennes narrowing, and (b) Gaussianity of the diffusion process; the former is shown to be violated in low-density fluids, and the latter is known to be violated in supercooled liquids.

Physics - Condensed matter physics↗

DOE OSTI2025

Combining Distributed Electrochemical Measurements with a Semi-Empirical Model to Identify Local Variations in Overpotentials for PEM Water Electrolyzers

Abstract A segmented cell and PCB approach was used to demonstrate an advanced technique for measuring and understanding the electrochemical behavior of PEM water electrolyzers. By employing a semi-empirical model, discrete modes of overpotential and their spatial distributions were identified in an operating cell. Distributed area-specific resistance (DASR measurements indicate that local dehydration of the polymer electrolyte membrane is the primary driver of increased polarization, behaving similarly to mass transport-limited currents at low flow rates due to ineffective gas-phase bubble removal. To improve quantification, the overall ohmic overpotential is divided into background (η Ohm,b ) and mass transport-coupled (η Ohm,mt ) components. Two types of anode diffusion media, a baseline PTL and a thin foil liquid/gas diffusion layer (LGDL), were examined. The results show significant variations in η Ohm,b and highlight that ineffective bubble removal leads to pronounced local dehydration. The effect of cell temperature on various overpotentials was also examined. Negligible differences in diffusion overpotential were observed when comparing PEWE performance using the PTL and LGDL, despite major differences in two-phase fluid transport within each diffusion media. This underscores the critical role of the catalyst layer's micro-porous structure in diffusion processes. The study highlights the complex nature and intricacies of mass transport limitations within PEWEs.

Roy, Anirban (ORCID:0000000163803386)↗

DOE OSTI2026

Conditional Latent Diffusion for High-Resolution Prediction of Electrochemical Surface Morphology

A conditionally guided generative latent diffusion process that is trained on a set of experimental processing parameters and their associated resulting electron microscope images of the electrodeposition process is able to interpolate between processing parameters in a physically consistent way. Electrodeposition of rhenium with pulse and pulse-reverse waveforms is used as a model system, and the process is adaptable to other electrodeposition, electropolishing, or corrosion processes. The method is able to extrapolate, predicting estimates of material morphologies for experimental setups unseen in the training data. The results are demonstrated with experimental data.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Non-Invasive Electrode Impedance Estimation for Optimized Charge Profile Parameterization of Lithium-Ion Batteries

This work presents a non-invasive method for parameterizing a physically motivated equivalent circuit model of lithium-ion batteries using operando electrochemical impedance spectroscopy and time-domain data. The proposed model consists exclusively of linear circuit elements, enabling computationally efficient simulation and real-time implementation on battery management system chips. By integrating frequency- and time-domain measurements, the model accurately estimates internal states such as the negative electrode potential, achieving a root mean square error of 12.3 mV during fast charging validation. Parameterization requires only rate tests with sinusoidal perturbations at three different ambient temperatures, making the approach experimentally accessible. The model reveals key insights into battery behavior, including rate-dependent overpotentials primarily governed by charge transfer kinetics at the positive electrode, and temperature-dependent impedance contributions from both charge transfer and solid-state diffusion processes. Validation using reference electrodes confirms the model’s ability to detect lithium plating onset and reproduce impedance behavior across a wide range of operating conditions. The approach enables in situ optimization of fast charging profiles and lays the foundation for future extensions incorporating aging effects and plating dynamics.

25 ENERGY STORAGE↗

DOE OSTI2025

The Galactic Annihilation Line Explorer (GALE): Balloon-borne Mission to study Galactic 511-keV annihilation emission

The Galactic Annihilation Line Explorer (GALE) mission will address the long-standing question of the sources of Galactic positrons: whether they are produced by unresolved astrophysical sources or created via diffuse processes, possibly due to dark matter decay and/or annihilation. The problem of Galactic positrons that produce 511-keV gamma-ray emission from the Galactic plane and Galactic Center has existed for years and is well-defined by the measurements of INTEGRAL/SPI. GALE is designed to achieve the needed angular resolution (fraction of a degree) and point source sensitivity (< 10−4 ph/cm2/s) required to understand the 511-keV emission structure in the Galactic Center. The GALE instrument combines a coded-aperture mask (CAM) with a cadmium-zinc-telluride imaging calorimeter (ImCal) to form a combined CAM/Compton telescope. This is mated with the Wallops Arc-Second Pointer (WASP) support system to achieve the precise telescope pointing control required for imaging the Galactic Center. In this paper, the science goals, GALE instrument concept, anticipated performance and mission results will be discussed.

Moiseev, Alex A [University of Maryland, College P↗

DOE OSTI2026

Conservation laws and effective hadronization models

Hadronization models based on local string-breaking dynamics are typically Markovian by construction, yet the physical ensemble of final states is shaped by global constraints that couple the entire fragmentation trajectory. Recasting hadronization as a conditioned stochastic diffusion process provides a precise mathematical resolution to this tension. In particular, this language reveals explicitly that constraints stemming from conservation laws induce non-Markovian correlations between otherwise independent fragmentation steps, and that these correlations can be absorbed exactly into a renormalization of the local dynamics through a Doob $h$-transform. We develop this formalism for a $q\bar{q}$ string in the chiral limit, where the longitudinal-transverse factorization of the Lund kernel becomes exact, enabling systematic power counting and clean ultraviolet (UV)/infrared (IR) separation. The dynamics organize naturally into a tower of effective theories distinguished by the remaining string mass, spanning a UV fixed point with scale-invariant transport coefficients, an intermediate regime where transverse phase space induces controlled running, and an IR boundary layer where non-local effects enter at leading order. The tower exhibits genuine Wilsonian structure, including $β$-functions, anomalous dimensions, and systematic matching conditions. The resulting framework achieves a clean factorization of universal microscopic fragmentation dynamics from infrared constraint effects, and opens new directions for both the theoretical analysis and practical simulation of hadronization.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗

DOE OSTI2025

A dynamic likelihood approach to filtering transport processes: advection-diffusion dynamics

A Bayesian data assimilation scheme is formulated for advection-dominated advective and diffusive evolutionary problems, based upon the Dynamic Likelihood (DLF) approach to filtering. The DLF was developed specifically for hyperbolic problems –waves–, and in this paper, it is extended via a split step formulation, to handle advection-diffusion problems. In the dynamic likelihood approach, observations and their statistics are used to propagate probabilities along characteristics, evolving the likelihood in time. The estimate posterior thus inherits phase information. For advection-diffusion the advective part of the time evolution is handled on the basis of observations alone, while the diffusive part is informed through the model as well as observations. We expect, and indeed show here, that in advection-dominated problems, the DLF approach produces better estimates than other assimilation approaches, particularly when the observations are sparse and have low uncertainty. The added computational expense of the method is cubic in the total number of observations over time, which is on the same order of magnitude as a standard Kalman filter and can be mitigated by bounding the number of forward propagated observations, discarding the least informative data.

97 MATHEMATICS AND COMPUTING↗

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↗