Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nernst–Planck equation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Physics-based stabilized finite element approximations of the Poisson–Nernst–Planck equations

We present and analyze two stabilized finite element methods for solving numerically the Poisson–Nernst–Planck equations. The stabilization we consider is carried out by using a shock detector and a discrete graph Laplacian operator for the ion equations, whereas the discrete equation for the electric potential need not be stabilized. Discrete solutions stemmed from the first algorithm preserve both maximum and minimum discrete principles. For the second algorithm, its discrete solutions are conceived so that they hold discrete principles and obey an entropy law provided that an acuteness condition is imposed for meshes. Remarkably the latter is found to be unconditionally stable. We validate our methodology through transient numerical experiments that show convergence toward steady-state solutions.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2024

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗
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 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 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↗
NASA NTRS2013

Conductance of Ion Channels - Theory vs. Experiment

Transmembrane ion channels mediate a number of essential physiological processes in a cell ranging from regulating osmotic pressure to transmission of neural signals. Kinetics and selectivity of ion transport is of critical importance to a cell and, not surprisingly, it is a subject of numerous experimental and theoretical studies. In this presentation we will analyze in detail computer simulations of two simple channels from fungi - antiamoebin and trichotoxin. Each of these channels is made of an alpha-helical bundle of small, nongenomically synthesized peptides containing a number of rare amino acids and exhibits strong antimicrobial activity. We will focus on calculating ionic conductance defined as the ratio of ionic current through the channel to applied voltage. From molecular dynamics simulations, conductance can be calculated in at least two ways, each involving different approximations. Specifically, the current, given as the number of charges transferred through the channel per unit of time, can be obtained from the number of events in which ions cross the channel during the simulation. This method works well for large currents (high conductance values and/or applied voltages). If the number of crossing events is small, reliable estimates of current are difficult to achieve. Alternatively, conductance can be estimated assuming that ion transport can be well approximated as diffusion in the external potential given by the free energy profile. Then, the current can be calculated by solving the one-dimensional diffusion equation in this external potential and applied voltage (the generalized Nernst-Planck equation). To do so three ingredients are needed: the free energy profile, the position-dependent diffusion coefficient and the diffusive flux of ions into the channel. All these quantities can be obtained from molecular dynamics simulations. An important advantage of this method is that it can be used equally well to estimating large and small currents. In addition, once the free energy profile becomes available the full current-voltage dependence can be readily obtained. For both channels we carried out calculations using both approaches. We also tested the main assumptions underlying the diffusive model, such as uncorrelated nature of individual crossing events and Fickian diffusion. The accuracy and consistency of different methods will be discussed. Finally we will discuss how comparisons between calculated and measured ionic conductance and selectivity of transport can be used for determining structural models of the channels.

Pohorille, Andrew↗
DOE OSTI2026

Integral Kernel Methods for Nonlinear Parabolic-Elliptic Systems

Nonlinear parabolic-elliptic systems arise in many physical, biological, and chemical phenomena such as chemotaxis, ion transport, self-gravitating particles, and Brownian vortices. Existing methods struggle with the strong coupling and high nonlinearity and nonlocality of some of these systems, especially the ill-conditioned, convection-dominated problems. To overcome numerical difficulties, current approaches rely on initial guesses, preconditioning, or iterative techniques with no convergence guarantees. They might suffer from poor scalability, large memory usage, and difficulty to parallelize. Inspired by the connection of parabolic-elliptic systems to stochastic processes, we introduce a novel meshless, monolithic, and fully explicit method that naturally encapsulates the elliptic and parabolic operators into a single step which updates each node deterministically with global information. By being fully quadrature-based, it avoids solving systems of discretized equations and does not utilize initial guesses or preconditioning, while requiring little memory and being easy to parallelize. We first derive the method in an integral kernel formulation with quadratic complexity in the number of integration nodes and then leverage kernel-independent fast multipole methods (FMM) to present a scalable algorithm with linear complexity. We provide numerical examples for the Poisson-Nernst-Planck equations in one, two, and three dimensions, together with the derivation of the integral kernel for each case. Furthermore, the examples demonstrate the fast convergence and scalability of the FMM-accelerated algorithm, as well as its suitability for convection-dominated problems, making it competitive against traditional PDE solvers.

PDE systems↗
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 OSTI2026

Ion Transport in Charged Membranes: Linking Electric-Field-Driven Mechanisms to Pore Size via Perturbation Analysis

Ion-exchange membranes are a critical component in electrochemical systems. Nevertheless, the understanding and modeling of ion transport within these porous structures have been limited by particular complexity reductions, either ignoring the dimensionality of their porous network architecture or imposing geometric assumptions (i.e., overlapping double layers). Before addressing this morphology-transport gap, a framework that relates the driving forces of transport to the geometry of a single pore is required. In this work, our modeling domain consists of a two-dimensional single pore with charged walls, connecting two identical electrolyte reservoirs. Using the Poisson-Nernst-Planck equations and regular perturbation theory, we decouple the electric fields and analyze the driving forces of ion transport, specifically electromigration and induced electroosmosis within the pore. These processes are described as analytical functions of the interaction aspect ratio, ?, defined as the ratio of the pore radius to the Debye length. Using this parameter, our study (i) describes the interplay between electromigrative and electroosmotic mechanisms that set ionic conductivity, (ii) identifies a dimensionless group of intrinsic electrolyte properties that indicates the predominant driving force, and (iii) provides a qualitative, confinement-dependent perspective on selectivity in ion-conducting membranes.

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

Large deviations of ionic currents in dilute electrolytes

Here, we evaluate the exponentially rare fluctuations of the ionic current for a dilute electrolyte by means of macroscopic fluctuation theory. We consider the fluctuating hydrodynamics of a fluid electrolyte described by a stochastic Poisson–Nernst–Planck equation. We derive the Euler–Lagrange equations that dictate the optimal concentration profiles of ions conditioned on exhibiting a given current, whose form determines the likelihood of that current in the long-time limit. For a symmetric electrolyte under small applied voltages, number density fluctuations are small, and ionic current fluctuations are Gaussian with a variance determined by the Nernst–Einstein conductivity. Under large applied potentials, the ionic current distribution is generically non-Gaussian. Its structure is constrained thermodynamically by Gallavotti–Cohen symmetry and the thermodynamic uncertainty principle.

Farhadi, Jafar [University of California, Berkeley↗
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 OSTI2022

MOSCATO Development and Integration in Fiscal Year 2022

During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. We built upon the work for FY22 and improved the code from several different perspectives. First, we improved the user interface by adding a new component to the official Nek5000 input file (.par). This new part contains documents parameters like, salt properties (density, viscosity, Cp, thermal conductivity), diffusion coefficients, standard potential, etc. Second, we built a conversion script to extract salt properties from the MSTDB-TP salt database and write to MOSCATO input file. Third, we migrated the code to NekRS, which is the GPU branch of Nek5000 and suitable for next generation supercomputers. Verification and Validation (V&V) work was also continued in FY22. Two tasks were performed. The first V&V task involved the validation of the Poisson-Nernst-Planck equation solver and Butler-Volmer electrode kinetics, by comparing with numerical and experimental data about thermoelectric cells. The second task involved the comparisons to corrosion results from a thermal convection loop run during the MSRE era. Satisfactory agreement was obtained from both tasks.

Yuan, Haomin↗
NASA NTRS1995

A universal steady state I-V relationship for membrane current

A purely electrical mechanism for the gating of membrane ionic channel gives rise to a simple I-V relationship for membrane current. Our approach is based on the known presence of gating charge, which is an established property of the membrane channel gating. The gating charge is systematically treated as a polarization of the channel protein which varies with the external electric field and modifies the effective potential through which the ions migrate in the channel. Two polarization effects have been considered: 1) the up or down shift of the whole potential function, and 2) the change in the effective electric field inside the channel which is due to familiar effect of the effective reduction of the electric field inside a dielectric body because of the presence of surface charges on its surface. Both effects are linear in the channel polarization. The ionic current is described by a steady state solution of the Nernst-Planck equation with the potential directly controlled by the gating charge system. The solution describes reasonably well the steady state and peak-current I-V relationships for different channels, and when applied adiabatically, explains the time lag between the gating charge current and the rise of the ionic current. The approach developed can be useful as an effective way to model the ionic currents in axons, cardiac cells and other excitable tissues.

NASA Discipline Regulatory Physiology↗
NASA NTRS2005

Removal of organic contaminants by RO and NF membranes

Rejection characteristics of organic and inorganic compounds were examined for six reverse osmosis (RO) membranes and two nanofiltration (NF) membranes that are commercially available. A batch stirred-cell was employed to determine the membrane flux and the solute rejection for solutions at various concentrations and different pH conditions. The results show that for ionic solutes the degree of separation is influenced mainly by electrostatic exclusion, while for organic solutes the removal depends mainly upon the solute radius and molecular structure. In order to provide a better understanding of rejection mechanisms for the RO and NF membranes, the ratio of solute radius (r(i,s)) to effective membrane pore radius (r(p)) was employed to compare rejections. An empirical relation for the dependence of the rejection of organic compounds on the ratio r(i,s)/r(p) is presented. The rejection for organic compounds is over 75% when r(i,s)/r(p) is greater than 0.8. In addition, the rejection of organic compounds is examined using the extended Nernst-Planck equation coupled with a steric hindrance model. The transport of organic solutes is controlled mainly by diffusion for the compounds that have a high r(i,s)/r(p) ratio, while convection is dominant for compounds that have a small r(i,s)/r(p) ratio. c2005 Elsevier B.V. All rights reserved.

NASA Discipline Life Sciences Technologies↗
NASA NTRS2011

Molecular Dynamics Simulation of the Antiamoebin Ion Channel: Linking Structure and Conductance

Molecular dynamics simulations were carried out in order to ascertain which of the potential multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistant with its measured conductance. Estimates of the conductance obtained through counting ions that cross the channel and by solving the Nernst-Planck equation yield consistent results, indicating that the motion of ions inside the channel can be satisfactorily described as diffusive.The calculated conductance of octameric channels is markedly higher than the conductance measured in single channel recordings, whereas the tetramer appears to be non-conducting. The conductance of the hexamer was estimated to be 115+/-34 pS and 74+/-20 pS, at 150 mV and 75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. On this basis we propose that the antiamoebin channel consists of six monomers. Its pore is large enough to accommodate K(+) and Cl(-) with their first solvation shells intact. The free energy barrier encountered by K(+) is only 2.2 kcal/mol whereas Cl(-) encounters a substantially higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. Ion crossing events are shown to be uncorrelated and follow Poisson statistics. keywords: ion channels, peptaibols, channel conductance, molecular dynamics

Wilson, Michael A.↗
NASA NTRS2023

Multiphysics Modeling of Electrode Dissolution under Hydrodynamic Conditions in Low Conductivity Water

Efficient plug-flow electrochemical reactors are characterized by their ability to effectively convert species with significant ionic strength and/or with the assistance of a well-behaved electrolyte. Because of the classic electroconductivity requirement, electrochemical reactors are not particularly attractive options for applications with low-conductivity conditions, especially for low-pressure feeds. However, these reactors can still find application in dedicated processes such as ultra-clean water conditioning and biological-coupled separation systems. Consequently, a multi-physics model was developed using COMSOL Multiphysics to characterize different reactor-design conditions based on the Tertiary Current Distribution which accounts for the effects of expected variations in electrolyte composition and ionic strength on the electrochemical process, as well as solution resistance and electrode kinetics. The modelling framework presented here employs single-phase laminar fluid flow, the Nernst-Planck equation, the water-based electroneutrality condition, and concentration-dependent overpotentials. This investigation particularly explores the electrolysis process of “anodic dissolution,” in which elemental species dissociate from the solid matrix of the electrode and enter the liquid phase of the electrolyte as soluble ions. The reactor configuration of interest has a rectangular parallel-electrode design and is operated galvanostatically. A series of parametric studies are carried out to inspect the response of the system when hydrodynamic, kinetic, and geometric conditions are changed. The parametric sweep resulted in variable effluent concentrations and system voltage, which reveal underlaying optimization patterns for electrode dissolution in low conductivity water.

Multiphysics↗
DOE OSTI2025

Unraveling membrane electrode assembly design for electrochemical conversion of carbon dioxide to formate/formic acid

This work presents a one-dimensional continuum modeling approach to investigate various cell architectures used for electrochemical conversion of CO 2 to formate/formic acid. Ion transport is simulated by a system of generalized modified Poisson–Nernst–Planck (GMPNP) equations that reflect the reactive transport phenomena including steric effects as the electrolyte solutions become concentrated. In the cathode catalyst layer, ionic current contributions from both the supporting electrolyte and solid-state ionomer are considered. Voltage and CO 2 utilization breakdowns are utilized to deconvolute the impacts of the cell architecture. The origins of (bi)carbonate formation in the cathode are explored, as the subsequent decrease in CO 2 availability is a key reason for low faradaic efficiencies to formate/formic acid. In addition, the role of a supporting electrolyte (KOH) is investigated to understand its tradeoffs: while the K + ions can improve both conductivity and electrochemically active surface area in the cathode, the presence of OH − ions raises the pH and leads to deleterious formation of (bi)carbonates. To this end, we also present parametric studies on the concentration and flow rate of supplied KOH to the cell, to establish a path towards eliminating the need for a supporting electrolyte.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗