Search NASA⌕ Search

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.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Se diffusion in CdTe thin films for photovoltaics

Manipulating CdSeTe bandgrading to enhance photocurrent and carrier lifetime is an essential part of high-performance CdTe photovoltaics (PVs). Here, we examine Se diffusion kinetics in single-crystal and polycrystalline CdTe during deposition, thermal annealing, and CdCl 2 treatments. Se distributions are determined by dynamic secondary-ion-mass spectroscopy and Auger electron spectroscopy depth profiling and coupled with electron backscatter diffraction images of the crystalline structure. Effective bulk and grain boundary diffusion coefficients are determined by analytical models and discussed in the context of processing and film morphology. Se is found to diffuse in CdTe at much higher rates during CdCl 2 treatments than with thermal processing alone. GB diffusion also occurs at a significantly faster rate than bulk diffusion. As a result of these two effects, the near interface bulk and GB Se diffusion during CdCl 2 treatments dominates the bandgrading profiles in CdTe PVs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Leveraging Sequential Doping of Semiconducting Polymers to Enable Functionally Graded Materials for Organic Thermoelectrics

With the ability to modulate electronic properties through molecular doping coupled with ease in processability, semiconducting polymers are at the forefront in enabling organic thermoelectric devices for thermal energy management. In contrast to uniform thermoelectric material properties, an alternative route focuses on functionally graded materials (FGMs) where one spatially controls and optimizes transport properties across the length of a thermoelectric material. While primarily studied in the context of inorganic materials, the concept of FGMs for organic thermoelectrics has not been explored. Herein, we introduce how molecular doping of semiconducting polymers enables spatial compositional control of thin-film FGMs. Specifically, we utilize sequential vapor doping of poly[2,5-bis(3-tetradecylthiophen-2-yl) thieno [3,2-b]thiophene] (PBTTT) with the small molecule acceptor 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ) to fabricate the simplest form of FGMs – double- segmented thin films. The two thin-film segments are of equal length (7.5 mm), but each set to different doping levels. Our study focuses on understanding the thermoelectric properties (Seebeck coefficient, α, and electronic conductivity, σ) and structural properties (through X-ray scattering, UV-vis-NIR spectroscopy, and Raman spectroscopy) within and across the two segments. We observe the presence of a small diffuse interfacial region of 0.5 – 1 mm between the two segments where the doping level and transport properties vary continuously. Despite the diffuse interface, the measured effective Seebeck coefficient (αeff) across the two segments is simply the average of α within each segment. Importantly, this experimental result is consistent with reported mathematical models describing αeff in graded thermoelectric materials. Finally, our results demonstrate the facile fabrication and characterization of functionally graded organic thermoelectric materials, providing guidelines for further development on more complex FGMs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamics of deuterium retention and desorption from plasma-facing materials in fusion reactor-relevant conditions

Hydrogen isotopes retention and desorption during and after discharges in fusion devices are still not well understood due to the complex device conditions and limitations of in-situ diagnostics and measurements. Here, we simulated well-diagnosed recent experiments at the DIII-D facility to benchmark our ITMC-DYN integrated package of modeling deuterium diffusion, retention, and desorption during and after D discharge irradiation. Modeling results were compared with detail experimental data of D desorption fluxes for various irradiation conditions. We predicted the temporal evolution of free and trapped D distribution in tungsten (W) plasma-facing material (PFM). Effects of key parameters namely diffusion coefficient, recombination rate, trapping energies against different defect types, were examined in these simulations. Existing experimental data of these parameters in literature varies significantly which makes it harder to identify key mechanisms and physics responsible for hydrogen isotope retention and desorption. The purpose of this work is to accurately simulate recent well-diagnosed reactor experiments given the uncertainties in such parameters and identify mechanisms responsible for the retention and desorption. We implemented the best identified diffusion, recombination, and trapping parameters in ITMC-DYN package that integrate both various collisional and thermal processes. We predicted, for example, that sample cooling between discharges in DIII-D operations can significantly affect the spatial distribution of trapped D in W under reactor irradiation conditions. Correct prediction of desorption spectra from samples irradiated during 10 DIII-D discharges showed that up to 35% of D can be retained in high binding energy defects such as vacancy clusters or voids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

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↗
DOE OSTI2021

A Pore‐Scale Investigation of Mineral Precipitation Driven Diffusivity Change at the Column‐Scale

Abstract Mineral precipitation affects the pore structure and thus transport properties of porous media. In this study, we investigated the pore‐scale dynamics of precipitation in diffusion controlled systems and the resulting impacts on the effective diffusivity, using a micro‐continuum reactive transport model. Forty two‐dimensional pore structures representing both idealized and realistic geometries were simulated with consideration of different precipitation scenarios and rates. A homogeneous nucleation scenario reproduced patterns observed in previous experimental study showing mixing‐induced precipitation, and a surface growth scenario captured the pattern for mineral precipitation on a substrate with the same or similar mineral structures. In all cases, local precipitation resulted in the reduction in the average porosity of the domain ( Φ ) until the cessation of diffusive transport and the termination of precipitation. The minimum porosity reached was referred to as the critical porosity ( Φ c ). The effective diffusivity ( D eff ) decreased with Φ and dropped sharply to effectively zero, that is, the critical effective diffusivity ( D eff c ), as Φ c was reached. These pore‐scale dynamics can be captured by a revised relationship that explicitly considers the critical porosity and the corresponding effective diffusivity, and the pre‐exponential coefficient and the exponent of the relationship varied with initial pore structure and the precipitation kinetics. Overall, the homogeneous nucleation scenario results in systematically larger Φ c and coefficients that give rise to a sharper decrease in diffusivity as Φ c is approached, compared to the surface growth scenario. The revised relationship was also implemented at continuum scale and used to examine column scale diffusivity change and reactions.

58 GEOSCIENCES↗
DOE OSTI2025

Choline Chloride-Based Water-in-Salt Electrolyte for Efficient Iron Electrodeposition

Electrochemical production of iron is a promising low-cost and modular approach to replace the traditional blast furnace. Aqueous electrolytes for iron electrolysis are advantageous as they can be operated at near-ambient temperatures, but they suffer from inefficiencies due to the parasitic hydrogen evolution reaction. In this work, we identify a new water-in-salt electrolyte (WiSE) based on choline chloride (ChCl) for high coulombic efficiency (>85%) iron deposition. Electrochemical analysis of the partial current densities of iron plating and hydrogen co-evolution revealed that, at optimal WiSE compositions, water reduction is kinetically suppressed resulting in an increase in the Fe plating efficiency. Decreased coordination of water and increased coordination of choline’s alcohol group with the Fe 2+ ion were observed through 1 H NMR providing evidence that water reduction is kinetically suppressed in WiSE. Additionally, Raman spectroscopy revealed that complexation effects (with Cl – ) reduce Fe 2+ diffusion coefficients and corresponding limiting currents as the ChCl concentration is increased. This results in an optimal WiSE composition (4 M ChCl + 1 M FeCl 2 ) that provides kinetic suppression of HER but also low transport resistance to Fe plating yielding 85% coulombic efficiency at high current densities.

Sinclair, Nicholas Scott [Case Western Reserve Uni↗
DOE OSTI2021

Taylor dispersion in osmotically driven laminar flows in phloem

Sucrose is among the main products of photosynthesis that are deemed necessary for plant growth and survival. It is produced in the mesophyll cells of leaves and translocated to different parts of the plant through the phloem. Progress in understanding this transport process remains fraught with experimental difficulties, thereby prompting interest in theoretical approaches and laboratory studies. The Münch pressure and mass flow model is one of the accepted hypotheses describing the physics of sucrose transport in the phloem. It is based on osmosis creating an energy potential difference between the source and the sink. The flow responding to this energy potential is assumed laminar and described by the Hagen–Poiseuille equation. This study revisits such osmotically driven flows in tubes with membrane walls by including the effects of Taylor dispersion on mass transport. This effect has been overlooked in phloem flow studies. Taylor dispersion can increase the effective transport of solutes by increasing the apparent diffusion coefficient. It is shown that, in addition to the conventional diffusive correction derived for impermeable tube walls, a new adjustment to the mean advective terms arises because of osmotic effects. Because the molecular Schmidt number is very large for sucrose in water, the sucrose front speed and travel times have a direct dependence on the Péclet number for different ranges of the Münch number. Fially, this study establishes upper limits on expected Taylor dispersion enhancement of sucrose transport.

42 ENGINEERING↗
DOE OSTI2021

Estimating the Diffusion Coefficient of Lithium in Graphite: Extremely Fast Charging and a Comparison of Data Analysis Techniques

Galvanostatic intermittent titration experiments were performed in three-electrode cells to characterize the effect of C/2, 2-C and 4-C charge rates on the observed lithium diffusion coefficient. As part of the data analysis process, we compared the classic Weppner-Huggins analysis of polarization data with a newer (Wang et al.) analysis method for depolarization data. At low values of x in Li x C 6 , both analysis methods showed the same general trend in the apparent lithium diffusion coefficient, 4-C > 2-C > C/2. The two techniques differed in the magnitude of the estimated diffusion coefficient by about a factor of 100. The observed increase in diffusion coefficient does not last over a large compositional range. Since the estimates from the method of Weppner and Huggins may contain artifacts due to the use of particulate electrodes and high charge rates, the method of Wang et al. may produce better values.

25 ENERGY STORAGE↗
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Connecting Hindered Transport in Porous Media across Length Scales: From Single-Pore to Macroscopic

Hindered mass transport is widely observed in various porous media, including filtration membranes and porous rocks. However, there is no universal model capable of predicting transport in porous media due to the heterogeneity of porous structures and the complexity of the underlying microscopic mechanisms, including electrostatic and hydrodynamic effects. Here, we used a highly ordered porous medium (i.e. a silica inverse opal film) as a model system to directly explore the effects of geometric parameters (i.e. pore size, pore throat size and tracer particle size) and microscopic interaction parameters (e.g. controlled by ionic strength) on nanoparticle transport in porous environments using three-dimensional single-particle tracking. We found a linear scaling relation between the macroscopic diffusion coefficient and microscopic diffusion behaviors and a combination of parameters associated with microscopic features and phenomena, including both geometric effects and particle-level interactions. The proportionality coefficient relating micro and macro behaviors was complex and related to the connectivity of the matrix and pore-size variation, which could lead to tortuous diffusion pathways, hindering macroscopic transport. Furthermore, these findings provide important insights into nanoparticles and macromolecules transport in various porous media.

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

Resolution of a few problems in the application of quasilinear theory to calculating diffusion coefficients in heliophysics [Slides]

Three problems with the historical literature were discussed. Kennel and Engelmann (1966): ‘notational error’ that implies that diffusion coefficients depend on the distribution of wave power with respect to ψ. This dependence is unphysical and is not present in Lyons+ 1971. Kennel and Engelmann (1966) assume that only a single frequency can be associated with each wave vector: not necessarily true for EMIC, e.g. Lyons (1974a, 1974b) incorrectly used the Jacobian to relate the power spectral density in (ω,χ) space to that into (κ ⊥ , κ ∥ space. The correct use of the Jacobian: (1) Eliminates undesirable properties associated with N(ω) : lack of integrability, lack of superposition; and (2) Can produce orders of magnitude changes to the diffusion coefficients, but depends on the wave distribution: bigger effect for more oblique waves, and depends on the plasma parameter: bigger effect for smaller $ω_{pe}/ω_{ce}$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Insights into Rotational and Translational Dynamics in Mixtures of Ethylene Glycol and Choline Chloride Using Nuclear Magnetic Resonance Techniques

This work examines molecular dynamics and interactions in ethylene glycol–choline chloride (EG–ChCl) mixtures across 0–33 mol % ChCl, spanning the true eutectic region near 17–20 mol % and the commonly used 1:2 formulation. We combine pulsed-field-gradient (PFG) diffusion, fast-field-cycling (FFC) relaxometry, temperature-dependent 13 C T 1 , and nuclear Overhauser effect spectroscopy (NOESY) to disentangle local from macroscopic dynamics. PFG and FFC show that both translational and average rotational motions largely track the strong increase in viscosity with ChCl content, with ethylene glycol consistently diffusing faster than the choline cation and no global dynamical anomaly at the eutectic composition. More subtle, site-specific composition effects nevertheless emerge. The ratio of the diffusion coefficient of the hydroxyl group of choline to the diffusion coefficient of the methyl group of choline displays a shallow minimum in the 17–25 mol % region, indicating a modest change in how the hydroxyl-bearing end of choline samples the underlying translational motion relative to the methyl groups. 13 C T 1 analysis shows that rotational correlation times at 25 °C generally increase with ChCl, reflecting viscosity-coupled slowing, while the CH 2 –N α site exhibits a small but reproducible deviation from this monotonic trend near the eutectic. NOESY spectra at similar compositions reveal enhanced cross-relaxation between EG and choline protons, consistent with increased headgroup–solvent contact density rather than a wholesale structural rearrangement. Overall, our multitechnique study demonstrates that EG–ChCl dynamics are predominantly viscosity-dominated, with the eutectic region acting as a subtle dynamical crossover where specific choline segments become maximally coupled to the hydrogen-bond network. These insights refine the structure–dynamics picture of choline-chloride DESs and provide practical guidance for tuning composition in electrochemical, separation, and catalytic applications.

diffusion↗
DOE OSTI2020

Data on the comprehensive first-principles diffusion study of the aluminum-magnesium system

First-principles calculation of diffusion coefficients between Mg and Al is investigated comprehensively using density functional theory (DFT). The effect of different uncertainty sources arising from first principles calculations has been investigated systematically. These sources include the diffusion model, energetic, entropic and attempt frequency calculations. Variation in self and impurity diffusion coefficients of Mg and Al in stable phases are quantified using different DFT settings and compared with the experiments. Using the optimal DFT settings, diffusion coefficients in metastable phases of Al and Mg are predicted. The dataset refers to “An integrated experimental and computational study of diffusion and atomic mobility of the aluminum-magnesium system” [1].

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗
DOE OSTI2024

Role of Short-Range Order on Diffusion Coefficients in the Li–Mg Alloy

Li–Mg alloys are important because of their beneficial role in fostering uniform plating and stripping of lithium in all-solid-state batteries. The alloy Li x Mg 1–x forms a solid solution on the BCC crystal structure when the lithium content is greater than x ≈ 0.3. The activation barriers of lithium and magnesium exchanges with a vacancy, crucial for substitutional diffusion, are predicted to be exceptionally low and almost identical, with negligible dependence on the alloy composition. The equilibrium vacancy concentration at room temperature is predicted to be very low, and it also remains almost constant with no dependence on Mg content in the alloy (for x Li ≥ 0.5). Nevertheless, both experiments and kinetic Monte Carlo simulations indicate that the tracer diffusion coefficients decrease by almost an order of magnitude with the addition of Mg to the alloy. In this contribution, the crucial role that chemical short-range order plays in affecting the diffusion coefficients is studied. Chemical short-range order is found to increase the effective activation barrier for lithium and magnesium diffusion by making successive atomic hops with vacancies correlated.

36 MATERIALS SCIENCE↗