Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion”

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 613 records · Page 34

Double Diffusion Maps and their Latent Harmonics for scientific computations in latent space

In this work, we introduce a data-driven approach to building reduced dynamical models through manifold learning; the reduced latent space is discovered using Diffusion Maps (a manifold learning technique) on time series data. A second round of Diffusion Maps on those latent coordinates allows the approximation of the reduced dynamical models. This second round enables mapping the latent space coordinates back to the full ambient space (what is called lifting); it also enables the approximation of full state functions of interest in terms of the reduced coordinates. In our work, we develop and test three different reduced numerical simulation methodologies, either through pre-tabulation in the latent space and integration on the fly or by going back and forth between the ambient space and the latent space. The data-driven latent space simulation results, based on the three different approaches, are validated through (a) the latent space observation of the full simulation through the Nyström Extension formula, or through (b) lifting the reduced trajectory back to the full ambient space, via Latent Harmonics. Latent space modeling often involves additional regularization to favor certain properties of the space over others, and the mapping back to the ambient space is then constructed mostly independently from these properties; here, we use the same data-driven approach to construct the latent space and then map back to the ambient space.

97 MATHEMATICS AND COMPUTING↗

Exponential Runge-Kutta Parareal for non-diffusive equations

Parareal is a well-known parallel-in-time algorithm that combines a coarse and fine propagator within a parallel iteration. It allows for large-scale parallelism that leads to significantly reduced computational time compared to serial time-stepping methods. However, like many parallel-in-time methods it can fail to converge when applied to non-diffusive equations such as hyperbolic systems or dispersive nonlinear wave equations. Here, this paper explores the use of exponential integrators within the Parareal iteration. Exponential integrators are particularly interesting candidates for Parareal because of their ability to resolve fast-moving waves, even at the large stepsizes used by coarse propagators. This work begins with an introduction to exponential Parareal integrators followed by several motivating numerical experiments involving the nonlinear Schrödinger equation. These experiments are then analyzed using linear analysis that approximates the stability and convergence properties of the exponential Parareal iteration on nonlinear problems. The paper concludes with two additional numerical experiments involving the dispersive Kadomtsev-Petviashvili equation and the hyperbolic Vlasov-Poisson equation. These experiments demonstrate that exponential Parareal methods offer improved time-to-solution compared to serial exponential integrators when solving certain non-diffusive equations.

97 MATHEMATICS AND COMPUTING↗

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗

A score-based diffusion model approach for adaptive learning of stochastic partial differential equation solutions

In this paper, we propose a novel framework for adaptively learning the time-evolving solutions of stochastic partial differential equations (SPDEs) using score-based diffusion models within a recursive Bayesian inference setting. SPDEs play a central role in modeling complex physical systems under uncertainty, but their numerical solutions often suffer from model errors and reduced accuracy due to incomplete physical knowledge and environmental variability. To address these challenges, we encode the governing physics into the score function of a diffusion model using simulation data and incorporate observational information via a likelihood-based correction in a reverse-time stochastic differential equation. This enables adaptive learning through iterative refinement of the solution as new data becomes available. To improve computational efficiency in high-dimensional settings, we introduce the ensemble score filter, a training-free approximation of the score function designed for real-time inference. Numerical experiments on benchmark SPDEs demonstrate the accuracy and robustness of the proposed method under sparse and noisy observations.

97 MATHEMATICS AND COMPUTING↗

Diffusion coefficient and nucleation density studies on electrochemical deposition of aluminum from chloroaluminate ionic liquid electrolytes

In this study, aluminum was electrodeposited from ionic liquids comprising a melt of alkyl imidazolium chlorides and aluminum chloride (AlCl 3 ). The ionic liquids utilized in this study were 1-ethyl-3-methyl-imidazolium chloride (EMIC)-AlCl 3 , 1-butyl-3-methylimidazolium chloride (BMIC)-AlCl 3 , and 1-hexyl-3-methyl-imidazolium chloride (HMIC)-AlCl 3 at the AlCl 3 mole fraction of 0.667 (molar ratio of 1:2). The electrochemical behavior of chloroaluminate species in the three ionic liquids was investigated by cyclic voltammetry (CV) and chronoamperometry (CA) techniques at different temperatures. The cyclic voltammograms indicated that the reduction of Al 2 Cl 7 – species to metallic aluminum followed a diffusion-controlled phenomenon. Moreover, even at a less negative applied potential, the higher current density was obtained for EMIC-AlCl 3 ionic liquid. That indicates that EMIC-AlCl 3 favors less energy consumption during the electrodeposition. The chronoamperometric analysis revealed that the onset of aluminum deposition from such ionic liquids proceeds via a three-dimensional instantaneous nucleation process. The concentrations of Al 2 Cl 7 – ions in the ionic liquids were calculated based on the thermodynamic data, which were 2520, 2271, and 2035 mol m –3 for EMIC-AlCl 3 , BMIC-AlCl 3 , and HMIC-AlCl 3 , respectively. The diffusion coefficient (D) values of Al 2 Cl 7 – species in such ionic liquids were also calculated at various temperatures. The D values determined from the CA technique at 363 K are 2.16 × 10 –11 , 1.03 × 10 –11 , and 0.87 × 10 –11 m 2 s –1 for EMIC-AlCl 3 , BMICAlCl 3 , and HMIC-AlCl 3 , respectively. In addition, the D value increased as temperature increased and decreased as the hydrocarbon group in the ionic liquid increased. As a result, the calculated grain sizes of nucleation range from 2 to 4 μm for these ionic liquids and are in good agreement with experimental data obtained from SEM micrographs.

42 ENGINEERING↗

Evolving radon diffusion through earthen barriers at uranium waste disposal sites

Field measurements of Rn-222 fluxes from the tops and bottoms of compacted clay radon barriers were used to calculate effective Rn diffusion coefficients (D Rn ) at four uranium waste disposal sites in the western United States to assess cover performance after more than 20 years of service. Values of D Rn ranged from 7.4 × 10 -7 to 6.0 × 10 -9 m 2 /s, averaging 1.42 × 10 -7 . Water saturation (S W ) from soil cores indicated that there was relatively little control of D Rn by S W , especially at higher moisture levels, in contrast to estimates from most steady-state diffusion models. Further, this is attributed to preferential pathways intrinsic to construction of the barriers or to natural process that have developed over time including desiccation cracks, root channels, and insect burrows in the engineered earthen barriers. A modification to some models in which fast and slow pathway D Rn values are partitioned appears to give a good representation of the data; 4% of the fast pathway was needed to fit the data regression. For locations with high S w and highest D Rn (and fluxes) at each site, the proportion of fast pathway ranged from 1.7% to 34%, but for many locations with lower fluxes, little if any fast pathway was needed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Non-equilibrium molecular dynamics studies of thermal diffusion of hydrogen isotopes in low concentration zirconium hydrides

Tritium permeability in zirconium-based tritium getter critically impacts tritium storage and environmental safety during operation of tritium-producing burnable absorber rods (TPBARs). Previous experiments indicated that during irradiation operation, the hydrogen equilibrium pressured is increased. Further experimental and modeling studies suggested that the enhanced tritium release observed for reactor scale assemblies might be related to a thermal diffusion known as the Soret effect. A direct measurement of the Soret factor, however, has not been performed. To improve TPBAR and other nuclear applications, here we have applied two non-equilibrium molecular dynamics methods to study thermal diffusion of hydrogen isotopes in low-concentration zirconium hydrides. One of the methods produces sufficiently converged results to distinguish crystal orientation, isotope type, and concentration effects. In conclusion, with this method, crystal orientation, isotope type, and concentration effects are discussed.

36 MATERIALS SCIENCE↗

Microstructural features and deuterium diffusion in lithium penta-aluminate pellets under He + and D + ion irradiation

Lithium (Li) penta-aluminate (LiAl 5 O 8 ) is investigated as a potential tritium (T) breeding material, with a focus on microstructural response to ion irradiation and deuterium (D) diffusion behavior. Under high-fluence ion irradiation (2 x 10 17 (He + +D + )/cm 2 ) at 773 K, LiAl 5 O 8 exhibits significant disorder on the Li sublattice, as revealed by atomic-resolution scanning transmission electron microscopy, while the Al and O sublattices remain stable, demonstrating strong resistance to structural amorphization. Irradiation induces the formation of platelet-shaped antiphase boundaries (APBs), which may serve as effective D trapping sites. Atom probe tomography suggests the presence of 6 LiD clusters in the mass spectra, though definite conclusions regarding APB composition are hindered by signal overlap and limited data statistics. Time-of-flight secondary ion mass spectrometry reveals that D retention approaches to saturation at 3 x 10 17 (He + +D + )/cm 2 . Isothermal and isochronal annealing studies determine an average diffusivity of 1.6 x 10 -13 at 773 K and an effective activation energy of 0.8 ± 0.1 eV for D migration. Compared to γ-LiAlO 2 , LiAl 5 O 8 demonstrates superior irradiation resistance, minimal Li loss, and enhanced D retention, underscoring its potential as a durable breeder material for T production. In conclusion, these findings provide key insights into the microstructural evolution, defect dynamics, and D retention mechanisms in LiAl 5 O 8 under reactor-relevant conditions.

42 ENGINEERING↗

Redox-stable symmetrical solid oxide fuel cells with exceptionally high performance enabled by electrode/electrolyte diffuse interface

Here, in this study, we report a high performance and redox-stable symmetrical solid oxide fuel cell (SOFC) based on (Ba 0.5 Sr 0.5 ) (Mo 0.1 Fe 0.9 )O 3-δ (BSMF) electrode and La 0.8 Sr 0.2 Ga 0.8 Mg 0.2 O 3-δ (LSGM) electrolyte. BSMF is able to operate both as anode and cathode. Excellent electrocatalytic activity has been achieved on BSMF towards hydrogen oxidation and oxygen reduction. Due to its closely matched lattice parameter to LSGM electrolyte, a unique diffuse interface is formed between BSMF and LSGM. Compared to a clean interface, e.g. BSMF/gadolinium doped ceria interface, this diffuse interface promotes the performance of BSMF electrode 1–1.8 times in 600–800 °C. Polarization resistance of the BSMF/LSGM specimen is as low as 0.047 and 0.007 Ωcm 2 in humidified H 2 and in air at 800 °C, respectively. On the BSMF/LSGM/BSMF symmetrical cell, a maximum power density of 2.28 W/cm 2 is achieved at 800 °C, the highest among with redox-stable ceramic electrodes to the best of our knowledge. Redox stability of this cell is confirmed. The role of anode and cathode is reversed back and forth in different operation modes. No apparent degradation is observed through 4 cycles within a 110 h operation period. These findings demonstrate that (Ba 0.5 Sr 0.5 ) (M o0.1 Fe 0.9 )O 3-δ coupled with LSGM electrolyte is an excellent choice to build a high performance, redox-stable SOFC.

25 ENERGY STORAGE↗

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗

Solid-State Calcium-Ion Diffusion in Ca 1.5 Ba 0.5 Si 5 O 3 N 6

Rechargeable batteries based on multivalent working ions are promising candidates for next-generation high-energy-density batteries. Development of these technologies, however, is largely limited by the low diffusion rate of multivalent ions in solid-state materials, thereby necessitating a better understanding of the design principles that control multivalent-ion mobility. We report Ca 1.5 Ba 0.5 Si 5 O 3 N 6 as a potential calcium solid-state conductor and investigate its Ca migration mechanism by means of ab initio computations and neutron diffraction. This compound contains partially occupied Ca sites in close proximity to each other, providing a unique mechanism for Ca migration. Nuclear density maps obtained with the maximum entropy method from neutron powder diffraction data provide strong evidence for low-energy percolating one-dimensional pathways for Ca-ion migration. Ab initio molecular dynamics simulations further support a low Ca-ion migration barrier of ~400 meV when Ca vacancies are present and reveal a unique "vacancy-adjacent"concerted ion migration mechanism. This work provides a new understanding of solid-state Ca-ion diffusion and insights into the future design of novel cation configurations that utilize the interactions between mobile ions to enable fast multivalent-ion conduction in solid-state materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thin-Film Paradigm to Probe Interfacial Diffusion during Solid-State Metathesis Reactions

We highlight a paradigm for studying complex reaction mechanisms that guide the synthesis of materials. Thin films of FeCl 2 and Na 2 S 2 were deposited to study the metathesis reaction to form FeS 2 and NaCl. In situ X-ray reflectivity was used to monitor the interface between materials, which revealed a slow, impeded reaction at high temperatures as compared to previous studies using powder samples. AC impedance and X-ray photoelectron spectroscopy provided insight into distribution of elements and conductivity of the phases present during the reaction, and phase-field modeling was used to elucidate the diffusion of ions throughout the thin-film bilayers. Here, the use of thin-film bilayers provides a simplified system to study solid-state metathesis reactions and highlights the complexity of diffusion at solid-state interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Defect Diffusion Graph Neural Networks for Materials Discovery in High-Temperature Energy Applications

Here, the migration of crystallographic defects dictates material properties and performance for a plethora of technological applications. Density functional theory (DFT)-based nudged elastic band (NEB) calculations are a powerful computational technique for predicting defect migration activation energy barriers, yet they become prohibitively expensive for high-throughput screening of defect diffusivities. Without introducing hand-crafted (i.e., chemistry- or structure-specific) descriptors, we propose a generalized deep learning approach to train surrogate models for NEB energies of vacancy migration by hybridizing graph neural networks with transformer encoders and simply using pristine host structures as input. With sufficient training data, computationally efficient and simultaneous inference of vacancy defect thermodynamics and migration activation energies can be obtained to compute temperature-dependent vacancy diffusivities and to down-select candidates for more thorough DFT analysis or experiments. Thus, as we specifically demonstrate for potential water-splitting materials, candidates with desired defect thermodynamics, kinetics, and host stability properties can be more rapidly targeted from open-source databases of experimentally validated or hypothetical materials.

14 SOLAR ENERGY↗

Reaction–Diffusion Coupling Facilitates the Sequential Precipitation of Metal Ions from Battery Feedstock Solutions

Here, the development of new technologies for chemical separations is urgently needed to meet the surging demand for critical materials that has strained resources and caused environmental challenges. Inspired by the classic Liesegang experiment, we demonstrated the separation of critical metal ions based on the coupling of ion diffusion and precipitation kinetics. For this purpose, a model feedstock solution simulating dissolved battery electrodes was placed on top of a hydrogel loaded with a precipitating agent, namely sodium hydroxide. As the lithium, manganese, cobalt, and nickel ions diffused into the gel, a gradient of precipitates formed along the length of the reactor. Elemental analysis of the spatially distributed precipitates showed the enrichment of nickel near the gel-solution interface, followed by the formation of an almost pure (>96%) manganese product further along the reactor. Optimization experiments revealed that a sodium hydroxide concentration of 10 mM and a gel/solution volume ratio of 2:1 favored efficient separations. The robustness of the method was demonstrated in four out of five feedstock compositions of typically used battery cathodes. Our proof-of-concept experiments present a paradigm for critical materials separations that does not require specialty chemicals, binding agents, membranes, or toxic solvents.

25 ENERGY STORAGE↗

Substrate-Dependent Exciton Diffusion and Annihilation in Chemically Treated MoS 2 and WS 2

Atomically thin semiconductors such as monolayer MoS 2 and WS 2 exhibit nonlinear exciton–exciton annihilation at notably low excitation densities (below ~10 excitons/μm 2 in exfoliated MoS 2 ). In this work, we show that the density threshold at which annihilation occurs can be tuned by changing the underlying substrate. When the supporting substrate is changed from SiO 2 to Al 2 O 3 or SrTiO 3 , the rate constant for second-order exciton–exciton annihilation, k XX [cm 2 /s], is reduced by 1 or 2 orders of magnitude, respectively. Using transient photoluminescence microscopy, we measure the effective room-temperature exciton diffusion coefficient in bis(trifluoromethane)sulfonimide-treated MoS 2 to be in the range D = 0.03–0.06 cm 2 /s, corresponding to a diffusion length of L D = 350 nm for an exciton lifetime of τ = 18 ns, which does not depend strongly on the substrate. We discuss possible mechanisms for the observed behavior, including substrate permittivity, long-range exciton–exciton or exciton–charge interactions, defect-mediated Auger recombination, and spatially inhomogeneous exciton populations arising from substrate-induced disorder. Exciton annihilation limits the overall efficiency of 2D semiconductor devices operating at high exciton densities; the ability to tune these interactions via the underlying substrate is an important step toward more efficient optoelectronic technologies featuring atomically thin materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Confinement on the Density and Diffusivity of Organic Molecules in Single-Digit Nanopores Relative to Bulk Fluid Behavior

Prediction of fluid properties in microporous and mesoporous confinement is non-trivial due to distinct influences from pore surface functionality and shape and fluid molecular characteristics. The vast majority of literature focuses nearly exclusively on a small number of components that are generally low molecular weight gases and vapors. Missing are studies that utilize a relatively large number of compounds that are liquid at ambient conditions to probe the impact of molecular shape and size on confined fluid properties. Such materials are of great interest in the extraction of oil and gas from shale reservoirs or for the design of selective adsorbents for chemical separation. In this work, we reveal non-intuitive relationships that emerge when comparing bulk fluid properties to those when confined in amorphous silica single-digit nanopores (approx. 5 nm SBA-15), specifically the confined density and self-diffusivity. A counterintuitive result was obtained where larger molecules had relatively faster self-diffusivities; however, this behavior was revealed to be caused by these larger compounds exhibiting a relatively low density when confined. As the prediction of confined fluid properties impacts a diversity of disciplines, the results highlight the need for additional investigations where consistent and large-scale experimental procedures are employed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Gold Catalyst Surface Morphology on Wetting Behavior and Electrochemical CO 2 Reduction Performance in a Large-Area Zero-Gap Gas Diffusion Electrolyzer

We report catalyst surface area and wetting behavior are key factors in determining the performance of gas diffusion electrode (GDE) electrolyzers for electrochemical CO 2 reduction. In this work, we report the integration of sub-1 μm thick nanoporous gold (npAu) catalyst coatings into a large-area (25 cm 2 ) zero-gap electrolyzer. The npAu coatings were prepared by magnetron sputtering (MS) of thin AgAu alloy films on the microporous carbon layer of a gas diffusion layer (GDL) followed by Ag leaching. Compared to MS Au films of the same thickness, npAu catalyst coatings enable higher Faradaic efficiencies and improved catalyst stability for CO 2 -to-CO reduction with Faradaic efficiencies of up to 88% at 100 mA/cm 2 . For a 800 nm npAu coating, the device level energy efficiency for CO 2 to CO conversion reaches 45% (52% for CO + H 2 ) at 100 mA/cm 2 with a single pass CO 2 conversion efficiency of ~12%. Contact angle measurements reveal that npAu coatings provide a more hydrophobic electrode interface compared to MS Au coatings, suggesting that the more hydrophobic interfacial environment of npAu coatings helps mitigating electrode flooding which is associated with performance deterioration over time.

30 DIRECT ENERGY CONVERSION↗

Diffusion and Adsorption of 2-Methylpentane and 3-Methylpentane in MFI-Type Zeolite Crystals

Due to the intimate contact of guest molecules with the inner surface of microporous adsorbents, slight changes in the structure of the guest molecules may lead to remarkable changes in their microdynamics. As an impressive example, in this work, diffusion measurements via IR microimaging (IRM) show that the transport diffusivity of 2-methylpentane (2MP) exceeds the value attained for 3-methylpentane (3MP) by a factor of 2−4.5 in silicalite- 1, an MFI-type zeolite crystal and a technically relevant porous host system. This experimental finding is rationalized by comparison with the outcome of dynamically corrected transition state theory (dcTST) simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗