Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multicomponent 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Cross-field electron diffusion due to the coupling of drift-driven microinstabilities

In this paper, the nonlinear interaction between kinetic instabilities driven by multiple ion beams and magnetized electrons is investigated. Electron diffusion across magnetic field lines is enhanced by the coupling of plasma instabilities. Here, a two-dimensional collisionless particle-in-cell simulation is performed accounting for singly and doubly charged ions in a cross-field configuration. Consistent with prior linear kinetic theory analysis and observations from coherent Thomson scattering experiments, the present simulations identify an ion-ion two-stream instability due to multiply charged ions (flowing in the direction parallel to the applied electric field) which coexists with the electron cyclotron drift instability (propagating perpendicular to the applied electric field and parallel to the ExB drift). Small-scale fluctuations due to the coupling of these naturally driven kinetic modes are found to be a mechanism that can enhance cross-field electron transport and contribute to the broadening of the ion velocity distribution functions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

The CALPHAD approach for HEAs: Challenges and opportunities

Phases are key microstructural features that determine material properties. The most intuitive way of representing phase stability in a material is by its phase diagram. The vastness of the composition space is the most compelling reason for investigating high-entropy alloys (HEAs), but it also creates major challenges in controlling phases and microstructures of HEAs due to lack of experimental phase diagrams in highly concentrated multicomponent space. The CALPHAD (Calculation of Phase Diagram) approach, although it needs improvements, is the only viable way to calculate multicomponent phase diagrams. The two major property databases in CALPHAD, Gibbs energy and mobility, when coupled, can provide not only the equilibrium phase information, but also information about diffusion and kinetics. In this article, we will review the use of the CALPHAD approach in understanding phase stability and diffusion in high-entropy alloys. Finally, we will discuss the notable trend of coupling the high-throughput CALPHAD approach with machine learning in HEA design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2023

Fundamental study of gas species transport in the oxygen electrode of solid oxide fuel and electrolysis cells

A fundamental analysis of multicomponent gas transport models was performed in application to the oxygen electrodes of solid oxide cells. It is common practice to neglect the effect of pressure gradients within oxygen electrodes, even though a net molar flux at the electrolyte surface implies that a pressure gradient must exist. The influence of both Darcy velocity and Knudsen flux are considered in the context of ordinary (Fickian) diffusion, the dusty gas model, and the binary friction model. Comparisons between the models and different sets of assumptions are made via parametric studies on operating load, oxygen partial pressure, microstructural properties, and electrode thickness. Results show that the pressure gradient will have a significant impact on the oxygen concentration distribution and therefore the concentration overpotential. In electrolysis mode, pressure increases up to 1 atm are predicted, indicating that pressure at the electrode/electrolyte interface could contribute to electrode delamination. Additionally, it is found that Darcy's law is insufficient for calculating the pressure distribution without accounting for the flux due to Knudsen diffusion. Additionally, it is found that for the range of properties typical of oxygen electrodes, there is a negligibly small difference between the dusty gas model and binary friction model from a practical standpoint.

08 HYDROGEN↗
DOE OSTI2025

Mesoporous Thin Film Architectures: Addressing Material Demands through Molecular Self-Assembly

Mesoporous thin films spark interest across a wide range of disciplines due to their tunable nanostructures, large internal surface areas, and strong compatibility with planar optical, electronic, and microfluidic devices. While attention in the porous materials community has shifted toward macroporous or disordered nanoporous systems, a resurgence in mesoporous thin film research is underway, driven by new molecular self-assembly methods, advanced materials chemistry, and improved characterization techniques. The integration of high-χN block copolymer design, kinetically persistent micelle templating, and postdeposition processing protocols now allows control over structural parameters such as pore size, wall thickness, porosity, and connectivity. These advances have overcome many of the thermodynamic and processing constraints that previously limited widespread adoption. Rather than serving only as high-surface-area supports, mesoporous thin films are engineered as active interfaces where responsive chemistries and nanoscale confinement act in tandem. Embedding switchable ligands, thermoresponsive polymers, redox mediators, or ion-selective groups directly within the pore walls enables real-time control over transport, optical, and electrochemical properties. These capabilities open up new directions in adaptive coatings, gated membranes, and fast-response biosensors. To further expand their functional scope, mesoporous films are integrated into hierarchical and multicomponent architectures. Techniques such as triblock terpolymer templating, crack-directed assembly, and nanoimprint lithography allow for control over spatial organization on the micron and submicron scale and pore system orientation. This enables programmable anisotropy, enhanced molecular diffusion, and wavelength-selective photonic behavior, essential for next-generation sensing, catalysis, and energy applications. Such structural and functional complexity requires equally sophisticated characterization. Multimodal and in situ techniques can track material dynamics under operational conditions. Recent progress includes extended-range ellipsometric porosimetry (EP) for hierarchical architectures, vacuum EP for interface energetics, time-resolved EP for diffusion kinetics, and correlative AFM-SAXS mapping. The introduction of advanced neutron-based spectroscopies, particularly quasielastic neutron scattering (QENS), promises to provide real-time access to ion transport dynamics and segmental motion under nanoscale confinement, offering a path toward deeper mechanistic understanding of structure-performance correlations in mesoporous systems. This Account reflects the technical advances made and the interdisciplinary collaborations that have shaped our collective vision. The particular dimensions of mesopores enable us to subtly tune interactions at the molecular, interfacial, and mesoscopic levels that permit us to harness nanoconfinement. What emerges is a versatile, modular platform capable of chemical gating, energy transduction, and sensing with a level of tunability unmatched by other porous materials. We highlight critical challenges including the need for more robust large-area processing, a deeper understanding of dynamic behavior under cycling, and better integration with device-level architectures. Our strategies support the transition of mesoporous thin films into active high-performance components in next-generation energy, environmental, and biomedical systems.

oxides↗
DOE OSTI2020

Solar-metallicity gas in the extended halo of a galaxy at z ∼ 0.12

ABSTRACT We present the detection and analysis of a weak low-ionization absorber at z = 0.121 22 along the sightline of the blazar PG 1424+240, using spectroscopic data from both HST/COS and STIS. The absorber is a weak Mg ii analogue, with an incidence of weak C ii and Si ii, along with multicomponent C iv and O vi. The low ions are tracing a dense (nH ∼ 10−3 cm−3) parsec-scale cloud of solar or higher metallicity. The kinematically coincident higher ions are either from a more diffuse (nH ∼ 10−5–10−4 cm−3) photoionized phase of kiloparsec-scale dimensions or are tracing a warm (T ∼ 2 × 105 K) collisionally ionized transition temperature plasma layer. The absorber resides in a galaxy overdense region, with 18 luminous (>L*) galaxies within a projected radius of 5 Mpc and velocity of 750 km s−1. The multiphase properties, high metallicity, and proximity to a 1.4L* galaxy, at ρ ∼ 200 kpc and separation |Δv| = 11 km s−1, favour the possibility of the absorption tracing circumgalactic gas. The absorber serves as an example of weak Mg ii–O vi systems as a means to study multiphase high-velocity clouds in external galaxies.

Pradeep, Jayadev↗
DOE OSTI2022

Geochemical Modeling of Celestite (SrSO 4 ) Precipitation and Reactive Transport in Shales

Celestite (SrSO 4 ) precipitation is a prevalent example of secondary sulfate mineral scaling issues in hydraulic fracturing systems, particularly in basins where large concentrations of naturally occurring strontium are present. Herein, we present a validated and flexible geochemical model capable of predicting celestite formation under such unconventional environments. Simulations were built using CrunchFlow and guided by experimental data derived from batch reactors. These data allowed the constraint of key kinetic and thermodynamic parameters for celestite precipitation under relevant synthetic hydraulic fracturing fluid conditions. Effects of ionic strength, saturation index, and the presence of additives were considered in the combined experimental and modeling construction. This geochemical model was then expanded into a more complex system where interactions between hydraulic fracturing fluids and shale rocks were allowed to occur subject to diffusive transport. We find that the carbonate content of a given shale and the presence of persulfate breaker in the system strongly impact the location and extent of celestite formation. The results of this study provide a novel multicomponent reactive transport model that may be used to guide future experimental design in the pursuit of celestite and other sulfate mineral scale mitigation under extreme conditions typical of hydraulic fracturing in shale formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Cations impact radical reaction dynamics in concentrated multicomponent aqueous solutions

Ultraviolet (UV) photolysis of nitrite ions (NO 2 – ) in aqueous solutions produces a suite of radicals, viz., NO · , O – , · OH, and · NO 2 . The O – and NO · radicals are initially formed from the dissociation of photoexcited NO 2 – . The O – radical undergoes reversible proton transfer with water to generate · OH. Both · OH and O – oxidize the NO 2 – to · NO 2 radicals. The reactions of · OH occur at solution diffusion limits, which are influenced by the nature of the dissolved cations and anions. Here, we systematically varied the alkali metal cation, spanning the range from strongly to weakly hydrating ions, and measured the production of NO · , · OH, and · NO 2 radicals during UV photolysis of alkaline nitrite solutions using electron paramagnetic resonance spectroscopy with nitromethane spin trapping. Comparing the data for the different alkali cations revealed that the nature of the cation had a significant effect on production of all three radical species. Radical production was inhibited in solutions with high charge density cations, e.g., lithium, and promoted in solutions containing low charge density cations, e.g., cesium. Through complementary investigations with multinuclear single pulse direct excitation nuclear magnetic resonance (NMR) spectroscopy and pulsed field gradient NMR diffusometry, cation-controlled solution structures and extent of NO 2 – solvation were determined to alter the initial yields of · NO and · OH radicals as well as alter the reactivity of NO 2 – toward · OH, impacting the production of · NO 2 . The implications of these results for the retrieval and processing of low-water, highly alkaline solutions that comprise legacy radioactive waste are discussed.

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

Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study

The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Transport phenomena in electrolyte solutions: Nonequilibrium thermodynamics and statistical mechanics

The theory of transport phenomena in multicomponent electrolyte solutions is presented here through the integration of continuum mechanics, electromagnetism, and nonequilibrium thermodynamics. The governing equations of irreversible thermodynamics, including balance laws, Maxwell's equations, internal entropy production, and linear laws relating the thermodynamic forces and fluxes, are derived. Green–Kubo relations for the transport coefficients connecting electrochemical potential gradients and diffusive fluxes are obtained in terms of the flux–flux time correlations. Furthermore, the relationship between the derived transport coefficients and those of the Stefan–Maxwell and infinitely dilute frameworks are presented, and the connection between the transport matrix and experimentally measurable quantities is described. To exemplify the application of the derived Green–Kubo relations in molecular simulations, the matrix of transport coefficients for lithium and chloride ions in dimethyl sulfoxide is computed using classical molecular dynamics and compared with experimental measurements.

42 ENGINEERING↗
DOE OSTI2025

Understanding the impact of SPMAK and PEGPEA in crosslinked PEGDA membranes: Methanol-carboxylate co-transport behavior and correlating structure-physicochemical-transport properties

Investigating multicomponent transport in dense, hydrated polymer membranes is necessary in applications such as fuel cells, electrolyzers, and desalination systems. Of particular interest are photoelectrochemical CO 2 reduction cells (PEC-CRC) which produce liquid products such as alcohols (methanol) and carboxylates (formate, acetate). Flux coupling and competitive sorption behavior of these products (solutes) have been found to affect permselectivity, thereby motivating us to investigate fundamental membrane structure-physicochemical-transport relationships. Past research has shown that systematic tuning of crosslinked cation exchange membranes (CEMs) with charge-neutral monoacrylate monomers containing alkyl and phenyl groups (i.e., poly(ethylene glycol) methacrylate (PEGMA), and poly(ethylene glycol) phenyl ether acrylate (PEGPEA)) can suppress acetate transport (permeability) in co-permeation with methanol. To further investigate this transport behavior and to enhance membrane ionic conductivities, 3-sulfopropyl methacrylate potassium (SPMAK), a sulfonated monoacrylate monomer is incorporated here. SPMAK content is varied with neutral PEGPEA and diacrylate crosslinker, poly (ethylene glycol) diacrylate (PEGDA) of two different chain lengths (n = 10 and 13) to prepare membranes of various compositions. Electrochemical and physicochemical properties including ionic conductivity, ion exchange capacity, and water uptake increase with increasing charged SPMAK content. Different states of water within hydrated membranes are probed using differential scanning calorimetry (DSC), where increasing intermediate (loosely bound) water is observed with increasing hydrophilic SPMAK content. The transport behavior of methanol and carboxylates (formate, acetate, propionate) are investigated, where permeabilities vary as methanol > formate > acetate ≈ propionate. Interestingly, permeabilities decrease with increasing PEGPEA content and are more dependent on solute diffusion than sorption. Permeabilities and diffusivities decrease while permselectivities increase with decreasing PEGDA chain length.

25 ENERGY STORAGE↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2020

Hydration and reaction mechanisms on sodium silicate glass surfaces from molecular dynamics simulations with reactive force fields

Recent development of reactive force fields have enabled molecular dynamics simulations of interactions between silicate glasses and water at the atomistic scale. While multicomponent silicate glasses encompass a wide variety of compositions and properties, one common structural feature in these glasses is the combination of the network structure that is made up of silica tetrahedra linked through corner sharing interspersed with network modifiers like alkali and alkaline-earth ions that break up the Si–O–Si linkages by forming nonbridging oxygen. In reactions with water, ion exchange between alkali ions in the glass and proton or hydronium in the solution, as well as hydrolysis reaction of the Si–O–Si linkages and subsequent silanol formation, is observed and well documented. Furthermore, we have used a set of recently developed reactive force field to investigate the reactions between water and the surfaces of silica and sodium silicate glasses of different compositions for reactions up to 8 nanoseconds. Our results indicate sodium leaching into water and diffusion of water molecules up to 25 Å into the glass surface. We examined the structural and compositional changes inside the glass and around the diffused ions and use these to explain the rates of silanol formation at the surface. We also observed proton transport in the glass which has an indirect influence on the silanol formation rates. While the surface of the glass was rough to start with, it undergoes further modification into a hydrated gel-like structure in the glass for up to 5 Å in the higher alkali containing glasses. It was found that the leached sodium ions remain close to the interface and that fragments of silicate network from the surface is capable of dislodging from the bulk glass and enter the aqueous solution. These simulations thus provide insights into the formation and structure of an alteration layers commonly observed in multicomponent silicate glasses corroded in aqueous solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Extended kinetic lattice grand canonical Monte Carlo simulation method for transport of multicomponent ion mixtures through a model nanopore system

Abstract An extended version of the original kinetic lattice grand canonical Monte Carlo simulation method combined with mean field theory (KLGCMC/MF) ( J. Chem. Phys. 2007 , 127, 024706) is presented for the study of transport of multicomponent ion mixtures through a model nanopore. Comparison of the extended KLGCMC/MF (eKLGCMC/MF) simulation results with Poisson–Nernst–Planck (PNP) calculations is also made to confirm the validity of the extended simulation approach. Unlike the original version of KLGCMC/MF simulation method that treats only a binary ionic solution with one cation and one anion species, this extended version can deal with a system that includes ternary ion mixtures. A diffusion probability algorithm is also added to the extended version of the simulation method to describe the inhomogeneous diffusivity of ions that is often observed in the ion permeation through nanopores. Both Legendre and Chebyshev polynomials of the second kind were tested as a basis set for the basis set expansion (BSE) method with which to calculate the reaction field energy in the eKLGCMC/MF simulation. It turned out that the Legendre polynomials perform better than the Chebyshev polynomials, and as a result, the Legendre polynomials were implemented in the current version of eKLGCMC/MF simulation algorithm. The presented eKLGCMC/MF simulation method with new features finds its potential applications in nanopore systems where the correlation between ion species with the same sign of charges plays a key role such as oscillating ion currents or anomalous mole fraction effects.

Choi, Inhyeok↗
DOE OSTI2025

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Coupling between bulk thermal defects and surface segregation dynamics

Surface segregation is a phenomenon that depends on the delicate interplay between thermodynamic driving forces and kinetic obstacles, for which elevated temperature is often needed to enhance the atom mobility and reach equilibrium. Using the classic system of Cu 3 Au ( 100 ) under nonisothermal conditions, in this study, we show an adatom process underlying transient surface segregation dynamics through the temperature-change-driven creation and annihilation of thermal vacancies in the bulk and the resulting bulk-surface mass exchanges. This is demonstrated by monitoring the surface composition evolution of Cu 3 Au ( 100 ) with temperature changes between 250 and 500 ° C , showing that the increase in temperature decreases monotonically the surface Au concentration as a result of the transfer of more Cu than Au from the bulk to the surface to form Cu-rich clusters of adatoms. Such a bulk thermal defect effect is expected to be universal in inducing the disparity in the bulk-surface mass exchanges of dissimilar atoms in multicomponent materials because of the inherent differences in the vacancy formation energies of the constituent atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ion Atmosphere of Wormlike Micelles Profiled by Contrast Variation Small-Angle Neutron Scattering

Structural studies of wormlike micelles have so far mostly focused on the conformational properties of surfactant aggregates. The diffuse ionic atmosphere, which has a profound influence on various micellization phenomena such as thermodynamic stability and structural polymorphism, remains largely unexplored experimentally. In this report a strategy of contrast variation small-angle neutron scattering for this crucial structural study is outlined. Underlined by a general criterion established for unbiasedly identifying the length scale relevant to charge association from the spectral evolution, our analytical framework can provide a quantitative description of counterion distribution in a mathematically tractable manner. Overall, our method can be conveniently extended to facilitate structural studies of complex multicomponent systems using contrast variation neutron scattering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗