SEARCH · Search NASA

Results for “Effective diffusion coefficient”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Investigation of Rare Earth Elements Mobility Using Molecular Dynamics

Rare Earth Elements (REE) are critical many applications such as electronics, permanent magnets, agricultural etc. Effective separation of REEs is a challenge due to their similarities in chemical and physical properties. Current separation processes are costly and environmentally detrimental. Therefore, research on REEs mobility in water-based solvents is necessary to design effective separation processes. In this research, the key kinetic and thermodynamic parameters, including diffusion coefficient and hydration shell are determined under the influence of external electric field. The diffusion coefficient of metal ions increases due to electric field in pure water.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2022

Investigation of Rare Earth Elements Mobility Using Molecular Dynamics

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2020

Machine Learning-Based Upscaling of Finite-Size Molecular Dynamics Diffusion Simulations for Binary Fluids

Molecular diffusion coefficients calculated using molecular dynamics (MD) simulations suffer from finite-size (i.e., finite box size and finite particle number) effects. Results from finite-sized MD simulations can be upscaled to infinite simulation size by applying a correction factor. For self-diffusion of single-component fluids, this correction has been well-studied by many researchers including Yeh and Hummer (YH); for binary fluid mixtures, a modified YH correction was recently proposed for correcting MD-predicted Maxwell–Stephan (MS) diffusion rates. In this study we use both empirical and machine learning methods to identify improvements to the finite-size correction factors for both self-diffusion and MS diffusion of binary Lennard-Jones (LJ) fluid mixtures. Using artificial neural networks (ANNs), the error in the corrected LJ fluid diffusion is reduced by an order of magnitude versus existing YH corrections, and the ANN models perform well for mixtures with large dissimilarities in size and interaction energies where the YH correction proves insufficient.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Spinel high-entropy oxides (FeNiCrMnZnX) 3 O 4 (X = Al, mg) as anode materials for high-performance lithium-ion batteries

To address the high cost, cobalt dependency, and resource constraints typical of conventional high-entropy oxide (HEO) anodes, this study reports the successful synthesis of two Co-free, six-component spinel-type HEOs(FeNiCrMnZnAl) 3 O 4 (HEO-Al) and (FeNiCrMnZnMg) 3 O 4 (HEO-Mg), via a sol-gel method. The distinct effects of Al 3+ and Mg 2+ incorporation on the electrochemical performance and lithium storage kinetics were systematically investigated. XRD, Raman, and TEM characterizations confirm that both materials possess a pure spinel phase, uniform particle size, and homogeneous elemental distribution. Notably, electrochemical evaluations reveal that HEO-Al delivers a superior reversible capacity of 480.7 mAh g −1 after 100 cycles at 0.1 A g −1 , and maintains 354.5 mAh g −1 after 1000 long-term cycles at 1 A g −1 , significantly outperforming HEO-Mg. Kinetic analysis indicates that HEO-Al exhibits lower charge transfer resistance, a higher Li + diffusion coefficient, and a pseudocapacitive contribution of up to 82%. Furthermore, these findings demonstrate that Al substitution effectively optimizes the structural stability and lithium storage kinetics of Co-free HEOs, providing a viable strategy for designing low-cost, highly stable HEO anode systems.

Anode materials↗

DOE OSTI2023

Defect-mediated diffusion of implanted Mg in GaN: Suppressing dopant redistribution by sequential thermal and microwave annealing

The diffusion behavior of Mg in Mg/N co-implanted GaN is investigated in response to a set of annealing conditions and methodologies, namely, 1000 °C/30 min thermal anneal, by high-temperature pulsed gyrotron microwave annealing at 1420 or 1500 °C, or by thermal and microwave annealing, sequentially. After 1000 °C annealing, the diffusion of Mg in GaN is found to be negligible, as measured by secondary ion mass spectrometry. Annealing by gyrotron microwave annealing alone induces the diffusion of Mg at a rate on the order of 10 −12 cm 2 /s. However, the use of a thermal anneal before microwave gyrotron annealing reduces this rate by an order of magnitude to 10 −13 cm 2 /s. We find that a model that considers Mg diffusion from an inhomogeneous medium that contains a defect-rich implanted region near-surface to a relatively pristine region below the implant range better explains the observed diffusion behavior than a conventional model that assumes a homogeneous medium. By analyzing the diffusion behavior using the Boltzmann–Matano method, we present a discussion of reduction in [V Ga ] by thermal annealing at 1000 °C, leading to a suppressed diffusion coefficient during subsequent high-temperature annealing relative to diffusion after 1420/1500 °C annealing alone. This effect holds potential for improvement in the precision of selectively doped regions for future applications based on the (Al)GaN material system. An improved field profile control in real devices can increase the breakdown and current-handling capabilities in power electronic applications.

Meyers, V.↗

DOE OSTI2020

Bridging Scales to Model Reactive Diffusive Transport in Porous Media

Published by ECS. Two novel scale-bridging algorithms to model reaction-diffusion transport in porous media are presented. The algorithms are based on direct numerical simulations and couple the information of a micro-scale model, which accounts for the large field of view provided by micro X-ray computed tomography (X-ray CT), and a nano-scale model, which locally resolves transport in the fine structure extracted from nano X-ray CT. The micro-scale model is discretized in the through-plane direction into a 1D grid, where effective properties and internal boundaries are determined based on the results from the nano-scale model. The validated algorithms are used to examine transport of oxygen in precious group metal-free electrodes considering both zero- and first-order kinetics. Unlike conventional methods, the results show that the effective diffusivity is not a passive property but increases in regions where the reaction-rate coefficient is large. The proposed algorithms account for the multiscale coupling of reaction-diffusion transport and material microstructure, thus improving the predictions compared to conventional methods.

42 ENGINEERING↗

DOE OSTI2022

Runaway electron deconfinement in SPARC and DIII-D by a passive 3D coil

The operation of a 3D coil—passively driven by the current quench (CQ) loop voltage—for the deconfinement of runaway electrons (REs) is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear magnetohydrodynamic (MHD) modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code to calculate advection and diffusion coefficients for REs based on the NIMROD-calculated fields, and the DREAM code to compute the runaway evolution in the presence of these transport coefficients. Compared with similar modeling in Tinguely et al (2021 Nucl. Fusion 61 124003), considerably more conservative assumptions are made with the ASCOT5 results, zeroing low levels of transport, particularly in regions in which closed flux surfaces have reformed. Of three coil geometries considered in SPARC, only the n = 1 coil is found to have sufficient resonant components to suppress the runaway current growth. Without the new conservative transport assumptions, full suppression of the RE current is maintained when the thermal quench MHD is included in the simulation or when the RE current is limited to 250kA, but when transport in closed flux regions is fully suppressed, these scenarios allow RE beams on the order of 1–2 MA to appear. Additional modeling is performed to consider the effects of the close ideal wall. In DIII-D, the CQ is modeled for both limited and diverted equilibrium shapes. In the limited shape, the onset of stochasticity is found to be insensitive to the coil current amplitude and governed largely by the evolution of the safety-factor profile. In both devices, prediction of the q-profile evolution is seen to be critical to predicting the later time effects of the coil.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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 OSTI2020

The impact of nanoparticle softness on its tracer diffusion coefficient in all polymer nanocomposites

The diffusion of nanoparticles in a polymer matrix is an area of current interest. However, a complete understanding is still limited as it is often difficult to quantify the much slower motion of nanoparticles in a polymer matrix. To combat this problem, we have developed a protocol to measure the diffusion coefficient of soft nanoparticles in a linear polymer matrix. Recently developed synthetic control over soft nanoparticle structures combined with this protocol provides a pathway to separately elucidate the effects of the molecular weight and nanoparticle softness on its diffusive behavior. These results indicate that the nanoparticle softness and deformability dictate its motion. Increasing the cross-linking density of the nanoparticle for all molecular weights increases its hardness and suppresses its motion in the linear matrix. Additionally, the nanoparticle molecular weight dependence deviates from the exponential dependence for star polymers suggesting that these nanoparticles benefit from the cooperative motion of the matrix to open pathways for the nanoparticle. Finally, comparison of these experimentally determined values to the Stokes–Einstein theory demonstrates that the nanoparticles diffuse much slower than a hard sphere. This is interpreted to indicate that there exist additional interactions between the nanoparticle and polymer matrix that are not captured by Stokes–Einstein, including threading or entanglement of the linear chain with the nanoparticle.

Rostom, Sahar↗

DOE OSTI2023

Study of Tritium Diffusivity in Pure and Sn-Defective Zr: A First-Principles Density Functional Theory Approach

Zirconium alloys (e.g., zircaloy-4) are used as tritium ( 3 H) getter materials in tritium-producing burnable absorber rods (TPBARs) owing to their ability to capture 3 H and chemically convert them into metal hydrides. Understanding of 3 H diffusion mechanisms in zircaloy is crucial for the optimal design of material performance in nuclear technology. Here, for this work, we perform first-principles density functional theory calculations to study the 3 H diffusion mechanism in pure and impure Zr with a low concentration of tin (Sn) atoms to determine the impact of the presence of Sn on the movement of 3 H atoms through the material. First, we calculated the diffusion barriers for 3 H in pure Zr by taking different migration pathways. We then introduced a low concentration of Sn impurity and systematically explored the impurity effect on the diffusion barriers for 3 H. Using our calculated diffusion energy barriers, we further obtained the diffusion coefficients and analyzed the results by comparing them with the experimental and previously calculated values. A diffusion coefficient of the order of 10 –8 m 2 /s is predicted. We also found that the presence of a Sn impurity could reduce the diffusivity up to 4 orders of magnitude. In conclusion, our results could serve as guidelines for further experimental investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Mesoscopic analyses of the impact of morphology and operating conditions on the transport resistances in a proton-exchange-membrane fuel-cell catalyst layer

Exploring the origins of local transport resistance and characterizing the oxygen transport resistances in the catalyst layer (RCL) are critical for cost reduction. In this paper, a comprehensive mesoscopic model for simulating coupled transport processes of oxygen and water vapor for different structural parameters under different operating conditions in reconstructed microstructures is proposed. The local transport resistance is calculated after achieving the limiting current density and the effective diffusivity of oxygen. The results demonstrate that RCL increases greatly with decreasing platinum loading (LPt) and the transport resistances in other components of the cell dominate for high-loadings. Both the reduced oxygen permeation coefficient in the ionomer thin-film and the adsorption resistance account for the origins of local transport resistance. The local transport resistance increases with the bare carbon ratio for a constant LPt and Pt/C ratio due to the decreased effective ionomer surface, and increases with the I/C ratio due to the increased ionomer thickness and decreased Knudsen diffusivity. Due to the presence of liquid water, a slight decrease followed by an increase of the local transport resistance versus relative humidity is obtained. The contribution of ionomer thin-films to RCL is more sensitive to liquid saturation compared with that of pores. This journal is

Mu, YT↗

DOE OSTI2024

Permeability and Induced Polarization of Mudstones

Electrical measurements can be used to estimate hydraulic properties such as permeability ( k ) in sedimentary rocks. Previous work has focused on sandstones, siltstones, and carbonates, while investigations on mudstones have rarely been reported. In this study, we report on electrical geophysical measurements for 23 mudstone samples using an experimental approach designed to reliably saturate these low permeability mudstones. The modified Hagen-Poiseuille model linking permeability to the formation factor ( F ) and an effective pore radius ( r ) provides an excellent fit to the data set with a near-constant pore radius, indicating that the effective porosity (1/ F ) is the controlling factor on k . In these samples, the surface area normalized to pore volume ( S por ), frequently used in permeability estimation models, varies by 1–2 orders of magnitude and is thus not a reliable proxy of the inverse effective hydraulic radius. The formation factor also exerts the primary control on induced polarization (IP) parameters, whereas Spor shows no relation to the IP parameters. A strong linear relationship is found between IP parameters (imaginary conductivity and normalized chargeability) and surface conductivity, although the proportionality factor is significantly lower than those observed in more permeable rocks and sediments. Apparent relationships between the polarization strength-derived and time constant-derived geophysical length scales and the effective hydraulic radius appear to be driven by variations in the electrochemical parameters (i.e., specific polarizability and diffusion coefficient). Overall, these findings emphasize that predicting hydraulic properties from electrical measurements in fine-grained rocks remains challenging and requires further investigation into the electrochemical properties involved.

58 GEOSCIENCES↗

DOE OSTI2026

Effects of confinement and pressure on the structure and dynamics of carbon dioxide in silica slit pores

An understanding of carbon dioxide fluid properties within geological mesopores is important in applications ranging from carbon sequestration to shale oil recovery. Here, a molecular dynamics study is presented that aims to shed light on these systems by simulation of CO 2 fluid confined in β-cristobalite silica slit pores with different pore widths and pressure conditions. The weakly associating nature of carbon dioxide leads to little difference in structural and interfacial dynamical properties for different pore sizes and pressures. Rather, the behavior is found to be dominated by the entropic effects, namely, how the CO 2 organizes next to the silica surface. The CO 2 self-diffusion coefficient shows the strongest pore size dependence. It is strongly diminished in small pores and does not reach the bulk fluid value even in 6 nm pores. It also decreases with pressure, yielding an activation volume that increases with pore size. These results provide new insight into the behavior of a compressible fluid in nanoscale confinement.

Godahewa, Sahan M. [Univ. of Kansas, Lawrence, KS ↗

DOE OSTI2023

Disentangling the effects of non-adiabatic interactions upon ion self-diffusion within warm dense hydrogen

Warm dense matter is a material state in the region of parameter space connecting condensed matter to classical plasma physics. In this intermediate regime, we investigate the significance of non-adiabatic electron-ion interactions upon ion dynamics. To disentangle non-adiabatic from adiabatic electron-ion interactions, we compare the ion self-diffusion coefficient from the non-adiabatic electron force field computational model with an adiabatic, classical molecular dynamics simulation. A classical pair potential developed through a force-matching algorithm ensures the only difference between the models is due to the electronic inertia. Here, we implement this new method to characterize non-adiabatic effects on the self-diffusion of warm dense hydrogen over a wide range of temperatures and densities. Ultimately we show that the impact of non-adiabatic effects is negligible for equilibrium ion dynamics in warm dense hydrogen.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Polysulfide entrapment and retardation in gel electrolyte Li–S batteries: experiments and modeling

Gelling the electrolyte in lithium–sulfur (Li–S) batteries provides advantages in terms of safety and mitigating polysulfide shuttling leading to reduced side reactions and dendrite growth on the anode side. In the current study, organic and ceramic crosslinkers have been used to improve the entrapment of polysulfides and prevent them from shuttling to the anode side. However, the performance of gel electrolytes is still far from being on par with liquid electrolytes in terms of raw capacity. The discharge curves exhibit an extended first plateau for gel electrolyte batteries compared to liquid electrolytes, showing that polysulfides are better trapped in the gel system. Yet, in the conversion of soluble polysulfides into insoluble polysulfides, the second plateau was delayed. Qualitative and quantitative analysis of the polysulfide diffusion in the liquid electrolyte and gel electrolyte systems were also conducted to predict the relative diffusivity ratio between the liquid electrolyte and gel electrolyte for simulation. To probe what happens at the end of the discharge, where the mass transport in the gel electrolyte appears to be slow, a modified discharge procedure was used where the cell is first paused after discharge and then discharged a second time before recharging normally. Coin cell tests for two gel electrolytes in addition to liquid electrolyte show up to 7% capacity recovery after the pause, which indicates that the diffusion of species during the pause is responsible for the failure recovery. Additionally, the recovery is higher for gel electrolytes compared to liquid, for higher C-rates, and for longer pause times. To better understand this behavior, a Li–S numerical model with mass transport limited reactions was used to examine the polysulfide diffusion coefficients expected in different electrolyte systems. The model is able to reproduce the trends seen in experiments and yields a higher recovery for smaller Li 2 S 4 diffusion coefficients, suggesting that insufficient Li 2 S 4 mass transport is responsible for failure at the end of discharge. By identifying the effect that gel electrolytes have on the state of the battery during discharge using both experimental and simulation methods, this work is a step towards leveraging the advantages of gel electrolytes and mitigating their weaknesses.

25 ENERGY STORAGE↗

DOE OSTI2024

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

In situ investigation of high-pressure hydrogen-induced swelling in elastomers and its correlation with material properties

The resistance of elastomeric materials to high-pressure hydrogen-induced damage is essential for ensuring the reliability of hydrogen infrastructure. Here, in this study, we systematically investigated the swelling behavior and hydrogen transport properties of four elastomer types – EPDM, NBR, FKM, and HNBR – using a custom in-situ view cell system capable of real-time monitoring during decompression from pressures up to 96.5 MPa. Each elastomer was formulated with and without fillers and plasticizers to assess the effects of formulation on swelling response. Thermal desorption analysis (TDA) was employed to determine equilibrium hydrogen content and diffusion coefficients, providing insight into gas uptake and mobility within each material. Correlation analyses using Pearson and Spearman coefficients revealed that the diffusion coefficient showed a stronger relationship with swelling behavior than hydrogen content, highlighting the dominant role of hydrogen mobility. Filled elastomers, particularly those with carbon black, consistently showed reduced swelling due to enhanced stiffness and reduced diffusivity. These results deepen our understanding of diffuso-mechanical interactions in elastomers and support the rational design of sealing materials for high-pressure hydrogen systems.

EPDM↗