Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The Riddle of Dark LLZO: Cobalt Diffusion in Garnet Separators of Solid‐State Lithium Batteries

Abstract Solid‐state batteries (SSBs) with a Li 7 La 3 Zr 2 O 12 (LLZO) garnet electrolyte are attracting much attention as robust and safe alternative to conventional lithium‐ion batteries. Technical challenges in the practical implementation of garnet SSBs are related to the need for high‐temperature sintering, which often leads to undesirable chemical reactions with the cathode material. While these reactions are well understood for composite cathodes, very little is known about similar processes between cathode and separator during battery fabrication. This work focuses on understanding the processes between the composite LiCoO 2 ‐LLZO cathode and the LLZO separator and how they affect the battery performance. The extensive diffusion of Co‐ions within LLZO, which leads to the often‐observed LLZO darkening, is shown to have a significant impact on ionic conductivity, electronic conductivity, and dendrite stability of the separator. Experimental data coupled with large‐scale molecular dynamics simulations uncover the diffusion mechanism for Co‐ions and identify secondary phases that form during these interactions. In addition to extensive Co‐ion diffusion within the grains, a non‐uniform segregation of Co‐ions at grain boundaries is found leading to the formation of three distinct Co‐containing phases. This work offers a general approach to studying the fundamental ion diffusion processes that occur during the fabrication of oxide SSBs.

electrochemical properties↗
DOE OSTI2024

Interpenetrated Structures for Enhancing Ion Diffusion Kinetics in Electrochemical Energy Storage Devices

The architectural design of electrodes offers new opportunities for next-generation electrochemical energy storage devices (EESDs) by increasing surface area, thickness, and active materials mass loading while maintaining good ion diffusion through optimized electrode tortuosity. However, conventional thick electrodes increase ion diffusion length and cause larger ion concentration gradients, limiting reaction kinetics. We demonstrate a strategy for building interpenetrated structures that shortens ion diffusion length and reduces ion concentration inhomogeneity. This free-standing device structure also avoids short-circuiting without needing a separator. The feature size and number of interpenetrated units can be adjusted during printing to balance surface area and ion diffusion. Starting with a 3D-printed interpenetrated polymer substrate, we metallize it to make it conductive. This substrate has two individually addressable electrodes, allowing selective electrodeposition of energy storage materials. Using a Zn//MnO 2 battery as a model system, the interpenetrated device outperforms conventional separate electrode configurations, improving volumetric energy density by 221% and exhibiting a higher capacity retention rate of 49% compared to 35% at temperatures from 20 to 0 °C. Our study introduces a new EESD architecture applicable to Li-ion, Na-ion batteries, supercapacitors, etc.

25 ENERGY STORAGE↗
DOE OSTI2022

Anisotropic thermal diffusivity and conductivity in SiC/SiC tubes studied by infrared imaging and X-ray computed tomography

Traditional flash diffusivity evaluation of thermal diffusivity/conductivity of composite tubes require machining of specimens. For a thin-wall tube, this method can only be used to obtain through-thickness transport property. A novel method to evaluate anisotropic thermal diffusivity in a composite tube has been developed. Braided SiC/SiC composite tubes were subjected to a xenon flash heating pulse. A high-speed, high-sensitivity infrared camera was used to measure surface temperature changes as a function of time and nondestructively detect subsurface defects/damages, such as macroscopic pores. Standard reference material (Pyroceram 9606) and curved SiC/SiC composite tube specimens were used to validate thermal diffusivity obtained from infrared imaging. Unlike the traditional method, there is no need prepare special specimens, and thermal diffusivity values in three orientations are obtained after a single flash. A finite element analysis model based on x-ray computed tomography scans was developed to simulate the heat transfer. In conclusion, this technique is significant in assessing thermal conductivity and inspecting the health of ceramic tubes during and after service.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Migration energy barriers and diffusion anisotropy of point defects on tungsten surfaces

We report tungsten is one of the most promising candidates for plasma-facing materials in future fusion devices, owing to its high performance under extreme irradiation conditions. However, irradiation-induced surface morphology varies significantly, depending on the irradiation type, fluence, and flux. Therefore, it is critical to examine the dynamics of point defects on tungsten (W) surfaces to understand how irradiation affects surface morphology. In this study, we employ the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method to search for potential migration paths of point defects on W (1 0 0), (1 1 0) and (1 1 1) surfaces. The first-principles calculations are then used to accurately determine the migration energy barriers. The obtained paths and barriers are incorporated into a kinetic Monte Carlo (KMC) model to determine trajectories and diffusivities, which are described by diffusion tensors to demonstrate their anisotropic features. Multiple diffusion mechanisms have been identified on different surfaces, with various anisotropy factors. Particularly, point defects on the W (1 1 0) surface have the highest diffusivities, with anisotropy factors independent of temperatures. In comparison, the anisotropy of the W (1 0 0) surface decreases as temperature increases, while the W (1 1 1) surface is isotropic for point defects. This study provides insights into defect transport properties on different surfaces, which are essential for understanding the early stages of irradiation-induced microstructural evolutions and surface morphology of tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Extended defects-enhanced oxygen diffusion in ThO 2

Oxygen self-diffusion is key to understanding stoichiometry and defect structures in oxide nuclear fuels. Experimentally, low activation-barrier oxygen migration was found in ThO 2 , a candidate nuclear fuel, possibly due to short-circuit diffusion mechanisms. Here, we perform extensive molecular dynamics simulations to show that various types of extended defects can enhance oxygen self-diffusion with a much-reduced activation barrier in ThO 2 . In this work, we consider extended defects including 1D (dislocation), 2D (grain boundary), and 3D (void) defects. Due to the distinct characteristics of each type of extended defect, the modulation of oxygen diffusion varies. Further, these results provide a quantitative description of oxygen transport, which is significantly enhanced within a close distance (nanometer scale) from the extended defects. Among all these considered defects, grain boundary, particularly the low-energy $Σ$3 twin boundary, exhibits the strongest effect on increasing oxygen transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Diffusion behavior of lanthanide fission products in bcc Fe cladding: A first-principles study

Fuel-cladding chemical interaction poses significant challenges in nuclear reactors, where fission products generated from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of lanthanide fission products, Lanthanum (La), Cerium (Ce), Praseodymium (Pr), and Neodymium (Nd), within body-centered cubic (bcc) Fe cladding using the density functional theory, nudged elastic band method, and self-consistent mean field theory. Our results reveal significant vacancy binding energies, particularly with the 1st and 2nd nearest neighbors, which diminish beyond the 5th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The nudged elastic band calculations indicate significant high barriers for the dissociation of 1st nearest neighbor vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products were derived in an Arrhenius form, with a magnetic correction that accounts for the high-temperature paramagnetic state. The significant trapping effect of vacancies caused by a very dilute concentration of fission products reduces vacancy mobility, potentially leading to modifications in point defect supersaturation, void nucleation, and swelling under irradiation. These represent critical challenges for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. Furthermore, this study provides essential insights for understanding fission product transport in cladding materials and informs future design strategies to mitigate fuel-cladding chemical interaction, ultimately enhancing nuclear reactor safety and performance.

Diffusion↗
DOE OSTI2025

Radiation-reduced diffusion of Nd in bcc Fe

Here, this study combines ab initio calculations and kinetic Monte Carlo simulations to investigate the bulk diffusion of Nd, an exemplar lanthanide element, in bcc Fe. Contrary to the conventional belief that irradiation accelerates diffusion processes, our study reveals that irradiation can unexpectedly reduce the diffusion of Nd due to the strong binding between substitutional Nd atoms and vacancies in bcc Fe. Combined with the high concentration of nonequilibrium vacancies, this promotes the abundant formation of Nd+3Va defect complexes containing a single Nd substitutional atom and three vacancies, accompanied by the depletion of Nd-vacancy pairs that mediate the diffusion of Nd.

Density functional theory (DFT)↗
DOE OSTI2023

Accelerating Markov Chain Monte Carlo sampling with diffusion models

Global fits of physics models require efficient methods for exploring high-dimensional and/or multimodal posterior functions. We introduce a novel method for accelerating Markov Chain Monte Carlo (MCMC) sampling by pairing a Metropolis-Hastings algorithm with a diffusion model that can draw global samples with the aim of approximating the posterior. We briefly review diffusion models in the context of image synthesis before providing a streamlined diffusion model tailored towards low-dimensional data arrays. We then present our adapted Metropolis-Hastings algorithm which combines local proposals with global proposals taken from a diffusion model that is regularly trained on the samples produced during the MCMC run. Our approach leads to a significant reduction in the number of likelihood evaluations required to obtain an accurate representation of the Bayesian posterior across several analytic functions, as well as for a physical example based on a global fit of parton distribution functions. Our method is extensible to other MCMC techniques, and we briefly compare our method to similar approaches based on normalising flows. A code implementation can be found at https://github.com/NickHunt-Smith/MCMC-diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2023

High-temperature chromium diffusion in austenitic stainless steel: Ab initio molecular dynamics simulations

Chromium self-diffusion through stainless steel (SS) matrix and along grain boundaries is an important mechanism controlling SS structural materials corrosion. For this work, Cr diffusion in austenitic SS was simulated using canonical ab initio molecular dynamics with realistic models of type-316 SS bulk, with and without Cr vacancies, and a low-energy Σ3 twin boundary typically observed at active corrosion sites. Cr self-diffusion coefficients at 750 and 850 °C calculated using Einstein’s diffusion equation are 4.2 × 10 –6 and 8.1 × 10 –6 Å 2 ps –1 in pristine bulk, 3.8 × 10 –3 and 5.5 × 10 –3 Å 2 ps –1 in bulk including Cr vacancies, and 9.5 × 10 –2 and 1.0 × 10 –1 Å 2 ps –1 at a Σ3[1 1 1]60° twin boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulating policy interventions in the interfirm diffusion of low-carbon technologies: An agent-based evolutionary game model

Policy interventions are frequently used by authorities around the world to mitigate carbon emissions. However, economic systems are essentially complex adaptive systems, which often exhibit unexpected responses to exogenous interventions and leave intervenors in a dilemma, even if the interventions are elaborately designed. To explore the behaviors of the diffusion system regarding low-carbon technologies, this study builds an agent-based model (ABM) to simulate enterprises' reactions to multiple policy interventions aimed at spurring low-carbon technology diffusion. The simulated enterprises are in a complex network where they play evolutionary games with their neighbors, which enables the model to possess two critical features of economic systems: adaptiveness and equilibrium. The model reveals a dilemma of policy interventions: intuitively, carbon taxes, asymmetric penalties (only imposed on un-low-carbon enterprises) and subsidies can improve the diffusion, which is consistent with previous research; counter-intuitively, all these policies turn out to be inefficient or even harmful to low-carbon enterprises because of the diffusion system's high adaptiveness. Specifically, when carbon taxes and penalties increase, both the low-carbon and un-low-carbon enterprises end up with the equilibrium of equivalent but lower profits. In contrast, all the enterprises earn equivalent but more profits even though subsidies are only given to the low-carbon enterprises, which implies that the un-low-carbon enterprises (indirectly) grab a portion of the subsidies and accordingly weaken the incentive effect of subsidies. These system behaviors are summarized as "equalizing effect" because the system tends to equalize the impact of both positive and negative interventions among all enterprises even though the policy interventions are asymmetrically imposed on one type of the enterprises. The findings also indicate that policies implemented to enlarge green market sizes can help policymakers bypass the dilemma.

Agent-based modeling↗
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

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗
DOE OSTI2023

A diffusion–deformation model with damage for polymer undergoing rapid decompression failure

Green hydrogen is emerging as one of the most promising alternatives to replace fossil fuels. While hydrogen gas has a good energy density by weight, its poor energy density by volume requires it to be stored under high pressure for commercial use. The hydrogen infrastructure developed to handle this high pressure hydrogen contains multiple components made of polymers as sealing agents. Although polymers do not react chemically with hydrogen gas, they undergo a mechanical failure when high pressure hydrogen gas is suddenly depressurized. This phenomenon known as 'rapid decompression failure' occurs due to the diffusion of hydrogen through polymer and getting trapped inside pre-existing cavities or voids. Here, in this paper, a continuum mechanics-based coupled diffusion-deformation-damage model was developed to predict the hydrogen distribution, stress distribution, and damage propagation inside the polymer while it undergoes rapid decompression failure. The polymer was modeled as a hyperelastic material because it represents the nonlinear material response observed in uniaxial tensile tests perfectly. The effects of hydrogen diffusivity, pre-existing cavity size, cavity location, applied hydrogen pressure, and depressurization rate on damage initiation were studied. It was found that the coefficient of diffusion plays an important role in damage initiation and damage was mostly concentrated in the inside areas rather than near the surface. Experiments were conducted with EPDM polymer which agreed well with the predicted trends using the given model. The effect of adding carbon black and silica filler particles and plasticizer to the pure EPDM polymer was also studied. It was found that damage during RDF decreases with the addition of fillers, but increases with the addition of the plasticizer. Finally, the damage evolution in the presence of two cavities was also studied, and was found that the interaction of stress fields around the cavities alters the damage occurring during RDF.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Micro-tensile testing of neutron-irradiated Al/Zr and Zr/U-Mo diffusion bonds

This study focuses on micro-tensile testing of neutron-irradiated Al/Zr and Zr/U-Mo diffusion bonds, which are integral to the structure and performance of nuclear fuel plates used in the fuel being developed for U.S. high-power research reactors. During fabrication, two distinct diffusion bonds occur—one at the Zr/U-10Mo interface via roll bonding, and the Zr/Al interface via hot isostatic pressing at 833 K. This work utilizes micro-tensile testing specimens approximately 7 × 7 × 18 μm 3 , which are sufficiently large to encompass the identified diffusion zones to evaluate the failure strength and identify the failure location. For the Zr/Al interface, all failures occurred in the bulk aluminum. The average yield strength, ultimate strength, and strain at failure for the two parent samples were: 125 ± 10 MPa, 139 ± 24 MPa, 156 ± 21 MPa, 185 ± 14 MPa, 37.9 ± 5.6 %, and 29.7 ± 3.4 %, respectively. For specimens across the Zr/U-Mo interface, failures occurred either in the bulk Zr or U-Mo, with no failures observed in the diffusion region. For specimens that failed in the bulk Zr the failure strength and strain for the two parent samples were 570 MPa, 710 ± 131 MPa, and 43.3 %, 27.3 ± 6.6 % respectively. The specimens that failed in the bulk U-Mo had failure strengths and strains of 596 ± 32 MPa, 548 ± 124 MPa, and 2.2 ± 1.3 %, 1.5 ± 0.6 % respectively. In conclusion, these findings support the development of accurate thermo-mechanical models for irradiated fuel performance by identifying the mechanical limits and failure locations within the bonded regions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A Robust Numerical Treatment of Solid-Phase Diffusion in Pseudo Two-Dimensional Lithium-Ion Battery Models

Solid-phase diffusion in active materials of lithium-ion batteries significantly affects charging and safety-related behavior of lithium-ion batteries. Therefore, it is essential to develop an efficient and robust numerical algorithm for solving solid-phase diffusion equations in physics-based battery models. In this work, we discuss the origins of numerical instabilities that can occur when solving the solid-phase diffusion equations using iterative methods. Then, in order to resolve such issues, we propose a simple numerical treatment to the surface flux term of discretized solid-phase diffusion equations. To demonstrate its numerical robustness, the proposed method is implemented into a pseudo two-dimensional (P2D) physics-based battery model and simulations are conducted at wide ranges of operating conditions. Even with extremely poor initial guesses for the Li+ concentrations of the active materials, computations using the proposed method do not diverge and the their computational speeds are comparable to those with conventional initial guesses. Comprehensive tests of the proposed method are also performed with a dynamic current profile based on US06 driving profile and a multi-stage charging profile with very high initial C-rate (12C).

battery modeling↗
DOE OSTI2025

Effect of high scandium doping in barium zirconate on nickel diffusion and performance of proton-conducting solid oxide electrolyzer cells

Proton-conducting solid oxide electrolyzer cells (p-SOECs) are emerging but promising technologies for hydrogen production. However, due to the lack of a robust electrolyte, p-SOECs struggle simultaneously to display high performance, Faradaic efficiency, and durability. Motivated by its high proton concentrations and stability as a barium zirconate, we have investigated BaZr 0.6 Sc 0.4 O 3-δ (BZSc40) as a potential next-generation p-SOEC electrolyte. Here, we found elevated levels of NiO diffusion through BZSc40 electrolytes during high-temperature sintering, attributed to the large oxygen vacancy concentrations present in BZSc40, as revealed by first-principle computational results. Controlling NiO diffusion is critical, as it can facilitate densification and grain size growth, but it may also detrimentally impact performance by causing electronic leakage. By optimizing sintering temperature when fabricating BZSc40 cells, we successfully controlled NiO diffusion, achieving sufficient electrolyte densification along with high performance and Faradaic efficiency. BZSc40 cells reached −0.99 A/cm 2 at 1.3 V and 600 °C and exhibited enhanced durability with a 3.37 mV/kh degradation rate at −0.8 A/cm 2 over a 200-h testing period. BZSc40 electrolytes demonstrated superior performance over BaZr 0.8 Y 0.2 O 3-δ (BZY20). In addition to elevated current densities and grain sizes, BZSc40 cells achieved Faradaic efficiencies of 76 % compared to 54 % for BZY20 at −0.2 A/cm 2 and 600 °C. This work lays the foundation for BZSc40 as a potential electrolyte due to its advantages over BZY20 while demonstrating the significance of controlling NiO diffusion when fabricating p-SOECs.

Electrolyzer↗
DOE OSTI2024

Electric field enhanced diffusion welding of alloy 617: Microstructural characteristics and mechanical properties

This study investigated the microstructural characteristics and mechanical behavior of diffusion welded nickel-based Alloy 617 obtained by electric field-assisted sintering (EFAS) using various parameters. The interfacial microstructure exhibited different characteristics including good grain boundary (GB) migration across the interface in the samples diffusion-welded at 1100 °C and a flat interface in the samples joined at 1000 °C and 1050 °C. The interface consisted of fine Al 2 O 3 oxides, while precipitation of interfacial M 23 C 6 carbides was not observed. Grain boundaries migrated across the Al 2 O 3 oxides, leaving these oxides within the grains. Graded grain size was observed, with grain coarsening being more significant near the sample surface due to the temperature gradient induced by EFAS. Tensile testing revealed that the specimens fractured in the matrix away from the interface, indicting strong diffusion-welded joints. Further, the peak tensile strength of 807 MPa was obtained in the samples welded at 1000 °C due to minimal grain growth. The materials obtained at 1100 °C exhibited reduced tensile strength but improved ductility. Strain maps revealed by digital image correlation showed alternating high and low strain segments in the samples produced at 1000 °C and 1050 °C, indicating that the flat interfaces with no GB migration were less ductile compared to the matrix. A greater strain uniformity was observed along the bond interfaces with improved GB migration. The hardness reduced near the sample surfaces due to enlarged grains induced by temperature gradient. This study demonstrates that GB migration and enhanced mechanical strength can be achieved in diffusion-welded Alloy 617.

36 MATERIALS SCIENCE↗