Search NASA⌕ Search

SEARCH · Search NASA

Results for “species 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 163 records · Page 9

Impact of phase state and non-ideal mixing on equilibration timescales of secondary organic aerosol partitioning

Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.

54 ENVIRONMENTAL SCIENCES↗

Mechanisms of Alkali Ionic Transport in Amorphous Oxyhalides Solid State Conductors

Amorphous oxyhalides have attracted significant attention due to their relatively high ionic conductivity (1 mS cm –1 ), excellent chemical stability, mechanical softness, and facile synthesis routes via standard solid‐state reactions. These materials exhibit an ionic conductivity that is almost independent of the underlying chemistry, in stark contrast to what occurs in crystalline conductors. In this work, we employ machine learning interatomic potentials to construct large‐scale molecular dynamics trajectories encompassing hundreds of nanoseconds to obtain statistically converged transport properties. We find that the amorphous state consists of chain fragments of metal‐anion tetrahedra of various lengths. By analyzing the residence time of alkali cations migrating around tetrahedrally‐coordinated metals, we find that oxygen anions limit alkali diffusion. By computing the full Einstein expression of the ionic conductivity, we demonstrate that the alkali transference number of these materials is strongly influenced by distinct‐particles correlations, while alkali transport is dictated by uncorrelated self‐diffusion. By extending this analysis to chemical compositions AMX 2.5 O 0.75 , spanning different alkaline (A = Li, Na, K), metallic (M = Al, Ga, In), and halogen (X = Cl, Br, I) species, we clarify why the diffusion properties of these materials remain largely insensitive to variations in atomic isovalent chemistry.

amorphous materials↗

Ni coarsening under humid atmosphere in the electrode of solid oxide cells: A combined study of density-functional theory and phase-field modeling

Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.

25 ENERGY STORAGE↗

Influence of the Soret effect on binary-species boundary layers at high pressure

Direct numerical simulations of binary-species temporal boundary layers at high pressure are performed. The main objective is to investigate the influence of the Soret effect on flow physics of binary-species boundary layers where the fluid has a uniform composition. The working fluid is a mixture of 25 % methane and 75 % nitrogen in mass fraction. Although the fluid composition is uniform at the initial condition, the mass fraction of methane increases near the wall when the wall temperature is hotter than the free stream temperature, whereas it decreases when the wall temperature is colder. The non-uniform mass fraction indicates that the uphill diffusion occurs near the wall. Investigation of fluctuations of the mass fraction reveals that the mass fraction fluctuates in the whole boundary layer, indicating that the uphill diffusion occurs even far from the wall. Examination of the species-mass diffusion balance for mean flow fields clarified that the Soret effect flux becomes large near the wall, and the large flux causes the non-uniform profile of the mass fraction near the wall.

Binary Diffusion Coefficient↗

Ab Initio Insights to Metal Passivation: Diffusion and Defect Formation in Amorphous Zirconia and Alumina

Metal passivation refers to the formation of protective oxide films on metals, which shield them from further corrosion and oxidation, playing a crucial role in maintaining their stability. The mesoscopic Point Defect Model has successfully predicted passivity as a steady state process where oxide growth from oxygen vacancies at the metal/film interface competes with oxide dissolution at the film/environment interface. In this work, informed by the Point Defect Model parameters, we use first-principles calculations to calculate defect formation and atomic diffusion in amorphous materials and correlate these descriptors with the behavior and growth of the oxide film. Focusing on amorphous zirconia and alumina, we demonstrate that defect formation energies exhibit significant variability in amorphous systems. In alumina, vacancies dominate, with cation and anion vacancies occurring at comparable concentrations. Diffusion calculations for stoichiometric amorphous alumina and zirconia, as well as oxygen-deficient zirconia, reveal faster diffusion in the oxygen-deficient case, highlighting the impact of defects on transport. Comparison of calculated self-diffusion coefficients for the dominant defect species with experimentally measured oxide thicknesses shows a clear correlation, suggesting that first-principles-derived diffusivity information can serve as a key descriptor for surface passivation film growth.

Farnell, Mackinzie S. [University of California, B↗

A Scaling Study for Incompressible Multispecies Solver in Vertex-CFD

Multispecies incompressible flows occur widely in engineering and environmental applications, such as chemical reactors, fuel cells, ocean mixing, and biomedical systems. However, accurately resolving the complex transport and mixing phenomena associated with multiple interacting species remains computationally challenging, especially for large-scale problems. In this study, we present a robust, high-performance computing--enabled multispecies incompressible Navier–Stokes solver integrated within the Vertex-CFD framework. Our solver employs a fully coupled, implicit, finite element--based formulation that accurately captures the advection, diffusion, and interaction of multiple species in incompressible flows by leveraging the Kokkos library for parallel computing to achieve high computational efficiency. For pressure coupling, the entropically damped artificial compressibility method is utilized. We validated the solver against canonical test cases, including multispecies advection, diffusion, and Bateman systems; the results demonstrate second- and third-order spatial accuracy and consistent convergence. Additionally, we demonstrated the strong and weak scaling study results obtained on the leadership-class high-performance computing system, Frontier at Oak Ridge National Laboratory.

Oz, Furkan [ORNL] (ORCID:0000000265831724)↗

Sodium diffusion in heterogeneous porous media: Connecting laboratory experiments and simulations

Sodium has been suggested as a tracer for brine in reservoir formations where a significant amount of sodium ion is found mainly in the aqueous phase. Fortunately, the most abundant sodium isotope, 23 Na, is Nuclear Magnetic Resonance (NMR) active, making it possible to study the structure and dynamical properties of both bulk and pore confined fluid systems. We investigated the diffusion of key dissolved sodium species in bulk solution and porous media as a function of composition, ionic strength, and permeability of the carbonate samples using NMR relaxometry and pulsed-field gradient (PFG) NMR techniques. We use Indiana limestone as an example of natural porous media and water as a freely diffusible tracer and carrier. We demonstrate apparent diffusion measurement of sodium based on changes in spin–spin relaxation time (T 2 ) signal in pore confinement. The diffusion rate of sodium decreases with increasing counter-ion size. This effect is greater at higher ionic strengths and lower chemical potentials in porous media. The reactive transport code, CrunchFlow, was used to complement the NMR experiments to simulate diffusion behavior in porous media. Combining these two methods provides a powerful approach to estimating effective diffusion coefficients in heterogeneous matrices. The modeling considers the influence of physical properties (porosity and tortuosity) and chemical properties (geochemical composition and chemical gradients). The results presented in this work highlight the advantage of measuring apparent diffusivity using NMR T 2 relaxometry in conjunction with numerical simulation to derive effective diffusivity and the corresponding matrix properties (i.e., tortuosity) of the system.

58 GEOSCIENCES↗

Small to Large Polaron Behavior Induced by Controlled Interactions in Perovskite Quantum Dot Solids

The polaron is an essential photoexcitation that governs the unique optoelectronic properties of organic–inorganic hybrid halide perovskites, and it has been subject to extensive spectroscopic and theoretical investigation over the past decade. A crucial but underexplored question is how the nature of the photogenerated polarons is impacted by the microscopic perovskite structure and what functional properties this affects. To tackle this question, we chemically tuned the interactions between perovskite quantum dots (QDs) to rationally manipulate the polaron properties. Through a suite of time-resolved spectroscopies, we find that inter-QD interactions open an excited-state channel to form large polaron species, which exhibit enhanced spatial diffusion, slower hot polaron cooling, and a longer intrinsic lifetime. At the same time, polaronic excitons are formed in competition via localized band-edge states, exhibiting strong photoluminescence but are limited by shorter intrinsic lifetimes. Furthermore, this control of polaron type and function through tunable inter-QD interactions not only provides design principles for QD-based materials but also experimentally disentangles polaronic species in hybrid perovskite materials.

14 SOLAR ENERGY↗

Novel Structural Motif To Promote Mg-Ion Mobility: Investigating ABO 4 Zircons as Magnesium Intercalation Cathodes

There is an increasing need for sustainable energy storage solutions as fossil fuels are replaced by renewable energy sources. Multivalent batteries, specifically Mg batteries, are one energy storage technology that researchers continue to develop with hopes to surpass the performance of Li-ion batteries. However, the limited energy density and transport properties of Mg cathodes remain critical challenges preventing the realization of high-performance multivalent batteries. In this work, ABO 4 zircon materials (A = Y, Eu and B = V, Cr) are computationally and experimentally evaluated as Mg intercalation cathodes. Remarkably good Mg-ion transport properties were predicted and Mg-ion intercalation was experimentally verified in sol-gel synthesized zircon YVO 4 , EuVO 4 , and EuCrO 4 . Among them, EuVO 4 exhibited the best electrochemical performance and demonstrated repeated reversible cycling. While we believe that the one-dimensional diffusion channels and redox-active species tetragonal coordination limit the value of many zircons as high-performance cathodes, their unique structural motif of overlapping polyhedra along the diffusion pathway appears instrumental for promoting good Mg-ion mobility. The motif results in a favorable "6-5-4" change in coordination that avoids unfavorable sites with lower coordination along the diffusion pathway and a structural design metric for future Mg cathode development.

25 ENERGY STORAGE↗

Larger than uncorrelated vacancy diffusion contributions in chemically disordered crystalline materials

Arising from a variational approach to compute diffusion coefficients, we introduce “superkinetic kinosons”, which are contributions to the diffusion flux of mobile defects and atomic species by single jump mechanisms that can exceed their uncorrelated values. In vacancy-mediated diffusion in crystalline random and high entropy materials, these contributions primarily arise from intermittent jumps of slow moving atoms that remove the vacancy out of correlation traps and enable a quantification of such effects on the diffusion coefficients. They can be significant contributors to diffusion, even when faster moving atoms have infinite percolation networks. Furthermore, their existence and importance provide fundamental insights about underlying aspects of vacancy diffusion such as formation of localized correlation traps and deviation of vacancy diffusion from percolation behavior.

36 MATERIALS SCIENCE↗

Flue Gas CO 2 Capture via Electrochemically Mediated Amine Regeneration: Desorption Unit Design and Analysis

We analyzed the performance of an electrochemical desorption unit to be implemented in an electrochemically mediated amine regeneration (EMAR) process, which separates CO 2 and regenerates the amine upon modulation of a targeted species that complexes with the amine molecules. Using a simplified one-dimensional boundary layer model, we combined the kinetic expressions of the faradaic reactions with mass transfer effects to analyze the polarization behavior of the flow cell. Mass transfer limitations become dominant at high current densities, leading to increases in energy consumption. An understanding of scaling laws with key parameters, such as the Sherwood number, is required to achieve high power operation with minimal mass transfer losses. Because the EMAR chemistry is sensitive to temperature, heat distribution in the unit is also examined. Our results indicate that the desorption performance increases with addition of a high-temperature heat source. Additionally, incorporation of membrane crossover suggests energy losses because of diffusion and migration of active species. Furthermore, these findings inform on scaling laws for at-scale desorber design and provide a model for scaling up the EMAR system for larger scale operations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigation of N 2 /O 2 plasma interaction with Pt-catalyst: effect of metastable adsorbates on product hysteresis

The coupling of catalysts and atmospheric-pressure plasma has the potential to improve the efficiency of certain catalytic reactions. Understanding the changes that the catalyst surface undergoes during exposure to plasma is key to improving plasma–catalytic performance. In this work, long term exposure of Pt–Al 2 O 3 powder catalyst to an Ar/N 2 /O 2 non-equilibrium atmospheric-pressure plasma-jet was investigated. Products produced by the interaction were analyzed downstream with Fourier-transform infrared spectroscopy while surface species were analyzed operandi with diffuse reflectance infrared Fourier transform spectroscopy. During exposure, the catalyst temperature was ramped cyclically between 100 °C and 350 °C to understand how substrate temperature affects the plasma–catalyst interaction. Long-lasting changes were revealed to take place on the catalyst surface during plasma exposure. At low temperatures, Pt–O and Pt–NO accumulate on the surface which react at elevated temperatures to form NO 2 . NO 2 initially appears to spill on to the Al 2 O 3 support as nitrites and nitrates instead of desorbing. Stable surface conditions are only achieved after prolonged plasma exposure, when nitrate sites on the Al 2 O 3 support are filled. By changing the catalyst temperature at various rates, the impact of total plasma species flux to the surface was analyzed. It was found that decreasing the heating rate increased the hysteresis in the pattern of NO 2 formation during thermal cycling. The variation with temperature demonstrates that plasma exposure results in a buildup of surface NO x and oxygen species which react or desorb at high temperatures. The observed changes are discussed from the generic viewpoint that a non-equilibrium plasma interacting with a catalyst at low temperature introduces metastable steady-state surface conditions. Upon heating above a threshold temperature, the introduced surface modifications can change either due to thermal effects, or, for a plasma environment, by additional interaction with the incident plasma species flux. The surface/material changes take place in a highly predictable fashion and after sufficient time above the threshold temperature reach a steady-state condition that is different from the transient behavior that is observed during initial heating. During cooling the plasma-surface interaction exhibits a different behavior than during heating, and this results in hysteresis of diverse observables. The metastability/hysteresis description appears quite generic and analogous to hysteresis behavior seen for different systems. Furthermore, it is expected to be useful for understanding the consequences of plasma–catalyst surface interactions for various systems.

36 MATERIALS SCIENCE↗

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{∞}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Unveiling the Origin of Morphological Instability in Topologically Complex Electrocatalytic Nanostructures

Coarsening and degradation phenomena in metals have largely focused on thermally driven processes, such as bulk and surface diffusion. However, dramatic coarsening has been reported in high-surface-area, nanometer-sized Pt-based catalysts during potential cycling in an electrolyte at room temperature─a temperature too low for the process to be explained purely by surface mobility values measured in both vacuum and electrolytes (∼10 –22 and ∼10 –18 cm 2 /s, respectively). This morphological evolution must be due to a different mechanism for mass transport that is sensitive to electrochemical conditions (e.g., electrolyte composition, potential limits, and scan rate). However, there have been no notable studies of electrochemically induced coarsening in nanometer-sized electrocatalysts. Here, we unveil the origins of coarsening in an electrolyte through coupled in situ experiments and atomistic kinetic Monte Carlo (kMC) simulations. Our work demonstrates electrochemical coarsening is driven by two concurrent mechanisms that can be explained at the atomistic level: (i) dissolution/redeposition during the reduction of an oxidized species and (ii) rapid surface diffusion of undercoordinated atoms.

gold↗

Study of Am electrodeposition from AmCl 3 -LiCl-KCl Molten Salt

In this work, the electrodeposition reaction of Am in AmCl 3 -LiCl-KCl molten salt was studied using voltammetry and chronoamperometry. Electrodeposition of Am was shown to proceed via a two – step reaction scheme, involving the one-electron transfer reduction of Am 3+ to Am 2+ , followed by the two-electron transfer reduction of intermediate Am 2+ to Am 0 . A low Coulombic efficiency of Am deposition was measured despite the Am deposition occurring within the thermodynamic stability window of the supporting LiCl-KCl electrolyte. We hypothesize that the low efficiency is due to out-diffusion of the intermediate Am 2+ species away from the electrode surface. Experimental observations provide evidence of a kinetically – limited electrodeposition reaction, which allows for the loss of intermediate Am 2+ via diffusion leading to Am deposition inefficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implementation and verification of a user-defined element (UEL) for coupled thermal-hydraulic-mechanical-chemical (THMC) processes in saturated geological media

Efficient and accurate modeling of the coupled thermal-hydraulic-mechanical-chemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. Here, this work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user-defined element (UEL) subroutine which is then combined with the general-purpose FE software Abaqus Standard to solve initial-boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic-mechanical (HM) and thermal-hydraulic-mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed-form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three-dimensional case study is performed with reference to the in-situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical-concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Method of Finite Averages: A rigorous upscaling methodology for heterogeneous porous media

Rigorous upscaling techniques offer accurate and computationally-efficient strategies for modeling the average behaviors of multi-physical, multiscale phenomena in geological porous media. However, such techniques often rely on a variety of methodological assumptions that prohibit their rigorous application to practical systems (e.g., systems involving heterogeneous porous media, system-scale boundary conditions, and fine-scale dynamics that are not diffusion-dominant). In this work, we aim to formulate an upscaling methodology with few methodological assumptions to provide high levels of model generality and foster the utilization of rigorously-derived upscaled models in practice. In particular, we introduce the Method of Finite Averages (MoFA), a novel upscaling methodology for rigorously modeling heterogeneous porous media and system-scale boundary conditions. We then detail MoFA’s implementation for the advective–diffusive transport of a single species and compare the methodology with classic numerical techniques, as well as other rigorous upscaling techniques, to highlight MoFA’s unique combination of rigor and generality. We then validate the derived model while demonstrating its benefits in three numerical experiments. The results suggest that (1.) the applicability and a priori error guarantees of MoFA models do not directly depend on system geometry, (2.) a model’s applicability and error guarantees can be can arbitrarily expanded and reduced, respectively, with further computational expense, and (3.) downscaling with MoFA provides an efficient strategy for generating accurate pore-scale solutions from upscaled results. Ultimately, the results evidence that upscaled models can be rigorously derived for heterogeneous porous media systems and resolved in a fraction of the time it takes to perform the equivalent pore-scale simulations.

58 GEOSCIENCES↗

Dynamic and static corrosion of chromium-alumina refractory in simulated nuclear waste glass

Monofrax® K-3, a chromium–alumina refractory, is widely used in nuclear waste glass melters due to its high durability. However, aggressive vitrification conditions still lead to corrosion that limits melter lifespan. Here, this study investigates two key corrosion modes, subsurface and melt-line corrosion, under static and dynamic conditions, using a custom setup enabling precise control of melt velocity and temperature. Experimental results, interpreted using diffusion-limited dissolution models, show that (i) subsurface corrosion increases by approximately 10% per mm/s increase in melt velocity, while melt-line corrosion remains virtually unaffected, (ii) corrosion rate dependence on temperature is inversely proportional to melt viscosity, and (iii) glass composition affects both corrosion modes similarly, based on Cr 2 O 3 , Al 2 O 3 , and alkali contents. The experimental data also allowed the determination of diffusion coefficients of major dissolving species, which were found to be consistent between the melt-line and subsurface corrosion models, indicating their suitability for future CFD studies.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗