Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solute 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Quantifying Species Populations in Multivalent Borohydride Electrolytes

Multivalent batteries represent an important beyond Li-ion energy storage concept. The prospect of calcium batteries, in particular, has emerged recently due to novel electrolyte demonstrations, especially that of a ground-breaking combination of the borohydride salt Ca(BH 4 ) 2 dissolved in tetrahydrofuran. Recent analysis of magnesium and calcium versions of this electrolyte led to the identification of divergent speciation pathways for Mg 2+ and Ca 2+ despite identical anions and solvents, owing to differences in cation size and attendant flexibility of coordination. To test these proposed speciation equilibria and develop a more quantitative understanding thereof, we have applied pulsed-field-gradient nuclear magnetic resonance and dielectric relaxation spectroscopy to study these electrolytes. Additionally, concentration-dependent variation in anion diffusivities and solution dipole relaxations, interpreted with the aid of molecular dynamics simulations, confirms these divergent Mg 2+ and Ca 2+ speciation pathways. These results provide a more quantitative description of the electroactive species populations. We find that these species are present in relatively small quantities, even in the highly active Ca(BH 4 ) 2 /tetrahydrofuran electrolyte. This finding helps interpret previous characterizations of metal deposition efficiency and morphology control and thus provides important fundamental insight into the dynamic properties of multivalent electrolytes for next-generation batteries.

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

Modeling the electrochemical behavior and interfacial junction profiles of bipolar membranes at solar flux relevant operating current densities

A 1-dimensional, multi-physics model that accounts for the migration and diffusion of solution species, electrostatics, and chemical reactions, in particular water dissociation (WD), at bipolar membrane (BPM) interfaces was developed to study the electrochemical behavior of bipolar membranes (BPMs) at solar flux relevant operating current densities (tens of mA cm -2 ). Significant partial current densities for WD were observed at BPM voltages much less than the equilibrium voltage, e.g., 59 mV × ΔpH from both experiments and modeling. The co-ion leakage across the BPM at pH differentials accounted for the early presence of the partial current density for WD. Two distinctive electric field dependent WD pathways, the un-catalyzed pathway and the catalyzed pathway, were quantitatively and parametrically studied to improve the turn-on potential of the BPM. The catalyzed pathway accounted for the majority of the partial current density for WD at low voltages, while the un-catalyzed pathway dominated the WD at relatively high voltages. Significant WD was observed only within the interfacial CL (<5 nm), in which a large electric field was present. To improve the electrochemical behavior and the turn-on potential of BPMs, the impacts of the pKa of the immobilized WD catalysts, the electric-field dependent rate constant, the thickness of the catalyst layer and fixed charge density in BPMs on the partial current densities for WD were studied systematically. In addition, the electrochemical behavior and concentration profiles of BPMs in a buffered electrolyte were studied and contrasted with those in an un-buffered electrolyte from both modeling and experiments.

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

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Comparing Constant and Transient Membrane Transport Parameters for Use in Wave Desalination Models

Directly pressurizing seawater for desalination with reverse osmosis membranes via wave motion is a promising and sustainable method for producing freshwater in coastal regions. However, such a system could result in significant pressure fluctuations and a departure from conventional steady-state desalination operations. This study sought to assess if membrane transport parameters (apparent water and salt permeability) should be modeled as transient or constant in solution–diffusion-based modeling efforts of dynamically operated desalination systems, such as those coupled to wave power. Two approaches were used to model membrane transport parameters: one considered each parameter to be a function of the net driving pressure of the system, and the other assumed they were constant across all conditions. A pilot-scale system was used to conduct steady-state and controlled ramping experiments. Data from steady-state experiments were used to calculate transient and constant transport parameters. Parameter combinations were used in a simulation model to predict water flux and effective permeate salinity, and simulation outcomes were compared against experimental ramping results. The transient relationships for both water and salt permeability produced the most accurate results for water flux and comparable results for effective permeate salinity. Development of such relationships would be unique to a specific system but could be valuable in modeling wave-driven desalination systems across the wide range of operating conditions they experience.

13 HYDRO ENERGY↗
DOE OSTI2020

OECD/NEA COUPLED NEUTRONIC/THERMAL-FLUIDS BENCHMARK OF THE MHTGR-350 MW CORE DESIGN: RESULTS FOR PHASE-I EXERCISE 1

This report presents the results from code-to-code comparisons for a prismatic modular high-temperature gas-cooled reactor 350-MW design for steady-state neutron transport calculations. Integral parameters as well as power and flux distributions are provided. Differences between solutions from diffusion and transport operators are analysed. Effects of control rod homogenization between full block and one sixth block are discussed.

42 ENGINEERING↗
DOE OSTI2021

Strong Enhancement of Nanoconfined Water Mobility by a Structure Breaking Salt

For the majority of the water present on earth, the two most important factors influencing its behavior are confinement, in either inorganic or organic matrixes, and the presence of solutes. Here, we investigate the effect of confinement in 3 nm pores on water diffusivity in aqueous solutions with archetypical solutes, a structure making (kosmotrope) NaCl and a structure breaking (chaotrope) KCl, up to 1.0 M in concentration. The water diffusivity in bulk aqueous solutions in such a concentration range is known to decrease very slightly in the presence of NaCl and increase very slightly in the presence of KCl. However, here we observe the water diffusivity in confined H 2 O–KCl increases by a factor of 2 compared to the pure water diffusivity in the same confinement. This unusually strong cumulative effect of confinement and a structure breaking additive may have profound implications for the mobility and transport of aqueous species in nature.

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

Single and Dual-Solute Transport in Quaternary Ammonium-Functionalized Anion Exchange Membranes: Isolating the Impact of Water Volume Fraction

Membranes are essential components of photoelectrochemical CO 2 reduction cells (PEC-CRC), as they regulate solute transport and impact device efficiency. These cells are a promising approach for converting carbon dioxide into valuable chemicals. However, electrochemical reactions involve not just a single electrolyte; membrane transport typically includes multiple electrolytes. When a mixture of solutes is present, the presence of one solute in the membrane can impact the diffusion and sorption of other solutes. In complex systems where multiple solutes are diffusing simultaneously through a hydrated, dense polymer membrane, the resulting transport behavior is poorly understood. This study highlights the individual and two solute transport behavior of potassium formate, acetate, and bicarbonate ions in neutral membranes and quaternary ammonium-functionalized anion exchange membranes (AEMs) based on phenyl acrylate polymer backbones. By isolating the effect of water content, we emphasize how fixed charge density and solute-membrane interactions influence transport behavior. In neutral membranes, the diffusivity follows the order of potassium formate > potassium acetate > potassium bicarbonate. However, in AEMs, the order shifts due to changes in ion-membrane interactions, with potassium bicarbonate diffusivity surpassing that of potassium acetate. These results underscore the importance of ion exchange membranes and ion-membrane interaction in modulating transport properties, providing valuable insights for optimizing membranes in electrochemical applications.

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

Chemical trends favoring interstitial cluster formation in bcc high-entropy alloys from first-principles calculations

Achieving high strength and ductility is a common goal in structural alloy design. Body-centered cubic high-entropy alloys (HEAs) commonly highlight the conflict between these properties, with stronger alloys being brittle and vice versa. Recent reports suggest interstitial solutes can be used to overcome this trade-off, in some cases providing both strength and ductility enhancements. This effect has been correlated with interstitial cluster formation, although the conditions favoring their formation remain incompletely understood. Using first-principles calculations of solution energies and diffusivities, we provide insights into thermodynamic and kinetic factors favoring interstitial solute clusters. Among C, N and O solutes, O interstitials display most desirable diffusion kinetics. Further, the results highlight the importance of local composition fluctuations in the HEAs to enable the formation of clusters of appreciable size. The results are explained in terms of bonding and distortion trends across solutes and HEA compositions to provide guidelines for alloy design.

Borges, Pedro P P O↗
DOE OSTI2023

Perivascular network segmentations derived from high-field MRI and their implications for perivascular and parenchymal mass transport in the rat brain

A custom segmentation workflow was applied to ex vivo high-field MR images of rat brains acquired following in vivo intraventricular contrast agent infusion to generate maps of the perivascular spaces (PVS). The resulting perivascular network segmentations enabled analysis of perivascular connections to the ventricles, parenchymal solute clearance, and dispersive solute transport within PVS. Numerous perivascular connections between the brain surface and the ventricles suggest the ventricles integrate into a PVS-mediated clearance system and raise the possibility of cerebrospinal fluid (CSF) return from the subarachnoid space to the ventricles via PVS. Assuming rapid solute exchange between the PVS and CSF spaces primarily by advection, the extensive perivascular network decreased the mean clearance distance from parenchyma to the nearest CSF compartment resulting in an over 21-fold reduction in the estimated diffusive clearance time scale, irrespective of solute diffusivity. This corresponds to an estimated diffusive clearance time scale under 10 min for amyloid-beta which suggests that the widespread distribution of PVS may render diffusion an effective parenchymal clearance mechanism. Additional analysis of oscillatory solute dispersion within PVS indicates that advection rather than dispersion is likely the primary transport mechanism for dissolved compounds greater than 66 kDa in the long (> 2 mm) perivascular segments identified here, although dispersion may be significant for smaller compounds in shorter perivascular segments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Impact of hydrophobic pendant phenyl groups on transport and co-transport of methanol and acetate in PEGDA-SPMAK cation exchange membranes

CO 2 reduction cells are innovative devices that reduce CO 2 into valuable chemicals (i.e. methanol (MeOH) and acetate (OAc – )) on the cathode. A major challenge with such devices is to develop ion exchange membranes that allow ion-selective transport (i.e., protons for cation exchange membranes, CEM) and suppress the crossover of CO 2 reduction products. To design such membranes, it is important to understand multi-solute transport behavior of these solutes. Previously, our group reported acetate diffusivities in co-transport with MeOH increase in sulfonated CEMs, where we speculated charge screening of the electrostatic interactions by co-diffusing MeOH has an impact. Here, crosslinked membranes fabricated by photopolymerization of poly(ethylene glycol) diacrylate (PEGDA), 3-sulfopropyl methacrylate potassium (SPMAK, SO 3 – -containing ionomer), and a phenyl-containing comonomer either phenoxyethyl acrylate (PEA) or poly(ethylene glycol) phenyl ether acrylate (PEGPEA)) are investigated. We observe OAc – diffusivities to both (1) PEGDA-SPMAK and (2) PEGDA-PEGPEA increase in co-diffusion with MeOH and those to (3) PEGDA-SPMAK/PEGPEA decrease. To rationalize this emergent co-transport behavior, we speculate (1) electrostatic interactions are interfered with by co-diffusing MeOH, (2) chain mobility (segmental dynamics) increases in the presence of MeOH and (3) for films with both hydrophobic (PEGPEA) and hydrophilic (SPMAK) comonomers chain mobility (segmental dynamics) decrease due to interactions between the comonomer sidechains depressing overall solute diffusivities. Finally, while further investigations are needed, this work contributes to improving our fundamental understanding of the relationships between polymer film chemistry, solute chemistry, and emergent cotransport behavior observed and described in this work (and others) towards enabling the design of improved membrane materials.

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

MPACT Theory Manual (V.4.3)

This manual presents the theory underlying the three-dimensional (3D) whole-core, pin-resolved neutron transport calculation methodologies employed in the MPACT code. MPACT’s primary goal is to provide accurate sub-pin power distributions in a computationally efficient manner. To accomplish this, MPACT offers several different transport method options. The 2D/1D method, in which 3D problems are decomposed into an axial stack of radial “slices,” is currently the most commonly used. In this option, two-dimensional (2D) planar solutions are provided by the method of characteristics (MOC), and axial solutions are provided via one-dimensional (1D) approximate diffusion, or P3 solutions. The radial and axial solutions are coupled by (i) axial and radial transverse leakages, and (ii) a global 3D coarse mesh finite difference (CMFD) solve, which provides both acceleration and stability to the solution iteration scheme. The subsequent chapters of this manual present a range of topics, including MOC, CMFD, axial nodal transport solvers, 2D/1D, self-shielding, depletion, thermal-hydraulics, and transient methods. Development of the underlying theory of the 2D/1D method is presented, diagrams are included to highlight important algorithmic flow, and important concepts are discussed as appropriate. This manual is intended to be self-sufficient, but references to published articles and other materials are included for further reading. MPACT is a relatively new code, with new capabilities and many computational methods that did not exist until recently. This manual is an even newer document, with chapters written by several different code contributors working under time constraints. For these reasons, the manual is not yet complete and finalized. Future versions of this manual will address current deficiencies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

High-accuracy method for modeling nucleation and growth of particles

State-of-the-art numerical models describing the kinetics of aerosol particle nucleation and growth from a cooling vapor primarily use a nodal method, in which particles that are smaller than the critical size are omitted from consideration because they are thermodynamically unfavorable. This omission is based on the assumption that most newly formed particles are above the critical size, so that subcritical-size particles are not important to take into account. Due to the nature of the nodal method, it suffers from numerical diffusion, which can cause an artificial broadening of the cluster size distribution leading to a significant overestimation of the number of large-size particles. To address these issues, we propose a more accurate numerical method that explicitly models particles of all sizes, and uses a special numerical scheme that substantially reduces the numerical diffusion and provides high solution accuracy and numerical stability. We extensively compare this novel method to the commonly used nodal solver of the general dynamic equation (GDE) for particle growth and demonstrate that it offers GDE solutions with higher accuracy with low numerical diffusion. Incorporating small subcritical clusters into the solution is crucial for: 1) more precise determination of the entire particle size distribution function and 2) wider applicability of the model to experimental studies with non-monotonic temperature variations leading to particle evaporation. The computational code implementing this numerical method in Python is available upon request.

42 ENGINEERING↗
DOE OSTI2023

An Efficient High-Order Solver for Diffusion Equations with Strong Anisotropy on Non-Anisotropy-Aligned Meshes

This paper concerns numerical solution of the diffusion equation with strong anisotropy on meshes not aligned with the anisotropic vector field. In order to resolve the numerical pollution for simulations on a non-anisotropy-aligned mesh and reduce the associated high computational cost we propose an effective preconditioner, extending our previous work. Similar to the anisotropy-aligned mesh case, we apply the auxiliary space preconditioning framework to design a preconditioner where a continuous finite element space is used as the auxiliary space for the discontinuous finite element space. The key component is an effective line smoother that can mitigate the high-frequency errors perpendicular to the magnetic field. We design a graph-based approach to find such a line smoother that is approximately perpendicular to the vector fields when the mesh does not align with the anisotropy. Finally, numerical experiments for several benchmark problems are presented, demonstrating the effectiveness and robustness of the proposed preconditioner when applied to Krylov iterative methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fast estimation of reaction rates in spherical and non-spherical porous catalysts

Here we present a methodology for modeling multi-step reaction rates in porous catalyst particles for use in CFD-DEM and two fluid models. Single-step effectiveness factors based on a Thiele modulus, while useful, cannot accurately capture the cascading reaction systems common in high temperature vapor-phase chemical reactors like fluidized catalytic cracking units and catalytic biomass fast pyrolysis systems. Instead, multi-step effectiveness vectors derived from steady-state solutions to the governing reaction-diffusion equations are needed. Solutions for various catalyst shapes are presented, including spheres, cylinders, and prisms. Computational challenges inherent in repeated evaluation of reaction rates with diffusion limitations are discussed, and an efficient implementation based on pre-computed lookup tables is proposed and demonstrated on a simulation of a fluidized bed reactor. Open-source code is provided for the compilation of reaction rate tables for use in ODE, DEM, and two-fluid models.

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

Heat conduction in an irregular magnetic field. Part 2. Heat transport as a measure of the effective non-integrable volume

Given the large anisotropy of transport processes in magnetized plasmas, the magnetic field structure can strongly impact heat diffusion: magnetic surfaces and cantori form barriers to transport while chaotic layers and island structures can degrade confinement. When a small but non-zero amount of perpendicular diffusion is included, the structure of the magnetic field becomes less important, allowing pressure gradients to be supported across chaotic regions and island chains. We introduce a metric for the effective volume over which the local parallel diffusion dominates based on the solution to the anisotropic heat diffusion equation. To validate this metric, we consider model fields with a single island chain and a strongly chaotic layer for which analytic predictions of the relative parallel and perpendicular transport can be made. We also analyse critically chaotic fields produced from different sets of perturbations, highlighting the impact of the mode number spectrum on the heat transport. Our results indicate that this metric coincides with the effective volume of non-integrability in the limit $\kappa _{\perp } \rightarrow 0$ , where $\kappa_{\perp}$ is the perpendicular diffusion coefficient. We propose that this metric be used to assess the impact of non-integrability on the heat transport in stellarator equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Vacancy-mediated transport and segregation tendencies of solutes in FCC Ni under diffusional creep

Nickel-based alloys are known for their excellent corrosion resistance and high-temperature strength. However, at high temperatures, polycrystals experience grain boundary-stress-driven vacancy diffusion leading to diffusion-driven creep. Understanding the fundamentals of defect migration is essential to predict diffusion creep behavior. In this talk, we explore the Nabarro-Herring diffusion creep and the associated solute segregation that is likely to result from the different diffusivities of the solute and solvent elements. Density functional theory calculations are performed to determine the formation, migration and binding energies of the solute and vacancy defects in stress-free and stressed states. The dilute-limit transport coefficients are then calculated using the self-consistent mean field theory as implemented in the Kinetic Cluster Expansion (KineCluE) code. Predictions of diffusion creep strain rates and solute segregation are then presented for various solutes, stress and temperatures. This work opens research avenues for exploration of diffusion creep via characterization of elemental segregation to grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Particle-Based Conditional Sampling Scheme for the Simulation of Transport in Fractured Rock With Diffusion Into Stagnant Water and Rock Matrix

In situ experiments and field-scale characterization studies have pointed out that, in fractured crystalline media, groundwater flow is highly channelized. This implies that, at the scale of a single fracture, only part of the fracture surface area is in contact with flowing water, while the rest of in-plane water is essentially stagnant and can be accessed by solutes via molecular diffusion. Despite their importance for contaminant retention, to date, there are no numerical or analytical approaches that could be used to assess the implication of stagnant water zones on solute transport in realistic large-scale Discrete Fracture Network-based models. In this study, we present an efficient and flexible algorithm for the simulation of transport in fractured rock with diffusion into stagnant water and rock matrix. The algorithm is a generalization of a previously developed numerical framework for time domain particle tracking in sparsely fractured rock. The key of the generalization is that total time in fracture (τ f ) is first evaluated using a Monte Carlo sampling and then a second sampling is performed conditioned on τ f . The algorithm has been successfully validated against existing independent solutions and the implication of diffusion into stagnant water and secondary diffusion into the matrix has been assessed for a realistic modeling scenario. The results show that, due to diffusion into stagnant water, contaminants are more strongly retarded. This increased retention is more significant for sorbing species, as a larger number of sorption sites is accessible. A high sensitivity to the flowing channel/stagnant water zone geometry has also been observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Atomic-level mechanisms of short-circuit diffusion in materials

This paper reviews the recent progress in understanding the atomic mechanisms of short-circuit diffusion along materials interfaces, such as grain and interphase boundaries, as well as lattice and interfacial dislocations/disconnections. Recent atomistic computer simulations have shown that short-circuit diffusion is dominated by collective atomic rearrangements in the form of strings and rings of mobile atoms. The process is dynamically heterogeneous in space and time and has many features in common with atomic dynamics in supercooled glass-forming liquids. We discuss examples of grain boundary, interphase boundary, and dislocation diffusion in metals and alloys, including the solute effect on the diffusion rates and mechanisms. Interphase boundaries are exemplified by Al–Si interfaces with diverse orientation relationships and atomic structures. The hierarchy of short-circuit diffusion paths in materials is reviewed by comparing the rates of grain boundary, interphase boundary, and dislocation diffusion. Future directions in the field of short-circuit diffusion in defect core regions are discussed.

36 MATERIALS SCIENCE↗