Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nernst-Planck”

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.

12 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A scalable DG solver for the electroneutral Nernst-Planck equations

The robust, scalable simulation of flowing electrochemical systems is increasingly important due to the synergy between intermittent renewable energy and electrochemical technologies such as energy storage and chemical manufacturing. The high Péclet regime of many such applications prevents the use of off-the-shelf discretization methods. In this work, we present a high-order Discontinuous Galerkin scheme for the electroneutral Nernst-Planck equations. The chosen charge conservation formulation allows for the specific treatment of the different physics: upwinding for advection and migration, and interior penalty for diffusion of ionic species as well the electric potential. Similarly, the formulation enables different treatments in the preconditioner: AMG for the potential blocks and ILU-based methods for the advection-dominated concentration blocks. Here we evaluate the convergence rate of the discretization scheme through numerical tests. Strong scaling results for two preconditioning approaches are shown for a large 3D flow-plate reactor example.

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

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗
DOE OSTI2022

Finite Element Method for Electrochemical Transport

This code provides Finite Element solvers for Electrochemical Transport. A few examples from the literature are reproduced, with a focus on CO2 electrolysis. The Discontinuous Galerkin scheme for the electroneutral Nernst-Planck equations is from Roy, T., Andrej, J. and Beck, V.A., 2021. A scalable DG solver for the electroneutral Nernst-Planck equations. arXiv preprint arXiv:2112.09271. This work also includes a scalable preconditioner.

Beck, VictorA↗
DOE OSTI2024

An integrated approach to optimizing concentration shock wave electrodialysis using 2D multicell simulation and response surface models

Shock wave electrodialysis (SWED) is a highly promising technique for energy-efficient ion separation in the context of a circular economy. This paper presents a approach way of modeling and improving SWED using a two-dimensional multicell model combined with the COMSOL program and response surface methodology. The model integrates the Nernst-Planck equation, Darcy's law, and first-order electroosmosis to examine the local concentration, flux of ionic species, distribution of current, and velocity of flow in SWED cells under various operating conditions. We first illustrate the clear depiction of concentration, velocity, and electric potential distribution through contours which aids in identifying optimal operating conditions and designing scalable SWED systems. The results emphasize the significance of surface charge density and voltage in influencing the features of shock waves for obtaining effective ion separation while optimizing energy consumption and improving current efficiency by controlling the retention time of feed flow. Here, this study defines two crucial characteristics of shock waves, namely the length of the flat depletion zone of a fully developed shock wave (shock wave height) and the distance of shock wave propagation (shock wave length). These properties significantly impact separation performance, as determined by the simulation results. Additionally, the response surface methodology is incorporated with the COMSOL models to develop predictive models and graph responses, enabling a more comprehensive understanding of the interactions between parameters and performance indicators, such as removal ratio, energy consumption, and water recovery. Finally, this work suggests design tactics for expanding SWED processes and outlines potential areas for further research. This research provides valuable insights into the prospective applications, design optimization, and scalability of SWED in the field of electrokinetic separation technologies for green chemistry and a circular economy.

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

Validating corrosion models: A comparison of governing equations

Experimental validation of Finite Element Method (FEM) models varying electrochemical governing equations, inclusion of chemical reactions, and time on the resultant damage profile for two galvanic couples is explored. Two anode materials (Magnesium AZ31 and Carbon Steel) in contact with a cathode (Stainless Steel 304 L) were modeled in/exposed to NaCl (1 and 0.1 M respectively for the anode materials) for up to one week. The physics approach, inclusion of chemical reactions, and the boundary conditions required to accurately represent the damage profile in FEM models depended on the galvanic couple materials and, ultimately, the corrosion rate. For high rates of corrosion (i.e., magnesium anode), the Nernst-Planck equation with Electroneutrality was sufficient to describe the damage, while, for low rates of corrosion (i.e., carbon steel anode), the Laplace equation was sufficient. In all cases, the most complete governing equation (Nernst-Planck-Poisson Equation) was not necessary to accurately describe the damage. Precipitation reactions in solution also played a critical role in the predicted damage profile, especially for high corrosion rate systems. Finally, for short time periods (< 6 h), the choice of governing equations does not significantly influence damage profile results. Overall, the choice of physics to reduce error in simulations relies on the boundary conditions, geometry, conductivity of the solution, electrochemical potential differences, and time of exposure. The above results are discussed with regard to accuracy and computational savings.

Carbon steel↗
DOE OSTI2022

Numerical modeling of ion transport and adsorption in porous media: A pore-scale study for capacitive deionization desalination

In this report a pore-scale model is presented to simulate the dynamic ion transport and adsorption processes in porous electrodes used for capacitive deionization (CDI). The Stokes equation governing water flow is solved using the lattice Boltzmann method and Nernst-Planck equation describing ion transport is solved using the finite volume method. The ion adsorption process is considered at the surface of carbon electrodes. After validation against analytical solutions and published results, the model is used to study the coupled water flow, ion transport and adsorption in both two-dimensional and three-dimensional porous CDI electrodes at the pore scale. The effect of electrode microstructure, electrical potential and flow velocity on the adsorption processes is quantitatively investigated, and the relative importance of various parameters is determined. The presented model can be a powerful numerical tool to quantitatively analyze ion transport and adsorption in porous electrodes, and may provide useful information for the design and optimization of CDI electrodes and operating conditions for desired desalination efficiency, water throughput and cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Dimensionality effects on multicomponent ionic transport and surface complexation in porous media

Coulombic interactions between charged species in pore water and at surface/solution interfaces are of pivotal importance for multicomponent reactive transport in porous media. In this study, we investigate the impact of domain dimensionality on electrostatically coupled dispersion and surface-solution reactions during transport of acidic plumes and major ions in porous media. Column and quasi two-dimensional flow-through experiments were performed, with identical silica porous media and under the same advection-dominated conditions. Equal mass fluxes of different electrolyte solutions (i.e., HCl - pH ~ 2.8, NaBr - 100 mM, HCl - pH ~ 2.8 plus NaBr - 100 mM) were continuously injected in the 1-D and 2-D experiments and breakthrough curves of pH and major ions were measured at the outlet of the domains. The presence of pronounced ionic strength gradients in the transverse direction in the 2-D setup caused distinct retardation and transport behaviors of protons and major ions which were not observed in the one-dimensional column experiments. Furthermore, in the cases of salt electrolytes injection, considerably enhanced release of H + (>61%) from the quartz surface was observed in the multidimensional system compared to the one-dimensional setup. Reactive transport modeling was performed to reproduce the experimental outcomes and to analyse the coupling between transport processes, based on the Nernst-Planck formulation of diffusive/dispersive fluxes and on surface complexation reactions at the solid-solution interface. Electrostatic interactions between Na + , B – , and H + , and deprotonation of the quartz surface upon the formation of sodium outer-sphere complexes, are the primary controllers of the spatial and temporal features displayed by the pH and major ions measurements. The reactive transport simulations allowed us to interpret the experimental observations, to visualize the distribution and spatio-temporal evolution of dissolved and solid species, to identify a spatially heterogeneous zonation of Coulombic interactions with distinct behavior at the fringe and core of the injected plumes in the multidimensional setup, and to quantify the different components of the Nerst-Planck fluxes of the charged solutes. This study demonstrates that the domain dimensionality directly affects electrostatic interactions between charged aqueous species in the pore water and surface complexation reactions at the solid-solution interface. The non-trivial effects of dimensionality on multicomponent ionic transport result in a significantly different behavior in 1-D and 2-D systems.

58 GEOSCIENCES↗
DOE OSTI2021

Highly selective ion separations based on counter-flow electromigration in nanoporous membranes

Highly selective ion separations are vital for recovering important salts, and membrane-based techniques may enable environmentally friendly ion separations that operate continuously. This study demonstrates that simply opposing convective flow with electromigration in track-etched membranes leads to remarkable selectivities among monovalent ions. Because the membrane-pore diameters are orders of magnitude larger than the ion diameters, advection moves all ions at the same velocity, whereas electrophoretic velocities are proportional to ion mobilities. Thus during counter-flow cation electromigration, less mobile cations have lower electromigration velocity components and pass through porous membranes more rapidly than more mobile cations. Using membranes with 400 nm pores, counter-flow electromigration gives Li + /K + selectivities of 100 and Li + /Na + selectivities around 30 in mixed-salt studies. In this work, numerical simulations based on the extended Nernst-Planck equation agree with trends in experimental data and highlight the importance of high Péclet numbers, high current-to-flow ratios, and uniform current distributions for achieving high selectivities. Importantly, experimental separations give high selectivities at 0.3 M ionic strength and low Li + /K + ratios, which suggest these methods may also work in brine solutions. Nevertheless, energy costs for such separations are high.

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

Salt and ion transport in a series of crosslinked AMPS/PEGDA hydrogel membranes

Salt permeability and ionic conductivity are critical properties for membranes in water purification and energy applications. Both properties depend on individual ion sorption and diffusion coefficients, which are significantly influenced by polymer chemical and physical parameters such as fixed charge concentration and membrane water content. However, systematic studies connecting polymer structure to ion transport properties are still lacking. In this study, a series of uncharged and charged membranes were synthesized using poly(ethylene glycol) diacrylate (PEGDA) as a cross-linker and 2-acrylamido-2-methyl-1-propanesulfonic acid (AMPS) as a charged monomer. Membrane fixed charge concentration and water uptake were systematically varied by adjusting AMPS content in the pre-polymerization mixture. Salt sorption and permeability coefficients and ionic conductivity of these membranes were measured as a function of NaCl solution concentration (0.01–1 M). Combining the solution-diffusion model and Nernst-Planck equation, individual ion diffusion coefficients were calculated. Experimental Na + diffusion coefficients for all materials were well described by the Mackie and Meares tortuosity model, highlighting the strong influence of water content on ion diffusivity in both uncharged and charged polymers. Model predictions for Cl – diffusion coefficients agree reasonably well with experimental values, with some deviations occurring in more highly charged membranes. Finally, this discrepancy might result from interactions not captured by the Mackie and Meares model (e.g., fixed charge-ion interactions).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2023

Capturing membrane structure and function in lattice Boltzmann models

Here, we develop a mesoscopic approach to model the nonequilibrium behavior of membranes at the cellular scale. Relying on lattice Boltzmann methods, we develop a solution procedure to recover the Nernst-Planck equations and Gauss's law. A general closure rule is developed to describe mass transport across the membrane, which is able to account for protein-mediated diffusion based on a coarse-grained representation. We demonstrate that our model is able to recover the Goldman equation from first principles and show that hyperpolarization occurs when membrane charging dynamics are controlled by multiple relaxation timescales. The approach provides a promising way to characterize non-equilibrium behaviors that arise due to the role of membranes in mediating transport based on realistic three-dimensional cell geometries.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Modeling the performance and faradaic efficiency of solid oxide electrolysis cells using doped barium zirconate perovskite electrolytes

Y-doped BaZrO 3 (BaZr 1 x Y x O 3 δ , or “BZY”), a proton-conducting ceramic featuring high bulk conductivity and good chemical stability, is a promising electrolyte material for solid oxide electrolysis cells. Further doping with Ce and/or Yb (creating materials “BCZY” and “BCZYYb”) can improve conductivity and sintering properties, but at significant penalty to cells’ faradaic efficiency (FE). Studies have proposed that reduction of lattice Ce can occur in the hydrogen electrode, which consumes some hydrogen produced by the hydrogen evolution reaction, leading to decreased FE. Despite studies suggesting this phenomenon, the mechanism is largely unknown. In this work, we developed a multiphysics model to study the transport of multiple defect species and the performance of BZY, BCZY, and BCZYYb, capturing the tradeoff between enhanced performance at the cost of FE for BCZY and BCZYYb electrolytes compared to BZY. We also found that increasing the water content of the anode gas supply lowers the current output of the cell but results in better FE. The model, which uses several parameters previously unavailable in the literature, was validated to experiments varying temperature, steam water content, and electrolyte material, as well as two performance metrics (performance curves and FE). Results verify and explain observed trends, informing future work on Ce-doped BZY electrolytes.

08 HYDROGEN↗