Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electroneutrality”

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.

15 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Electrochemical Processes Breaking Strict Phase Electroneutrality in Microemulsions

Phase electroneutrality is a near-omnipresent assumption in electrochemical reactions. Even media such as biphasic oil-water interfaces and kinetically stable emulsions follow this “rule.” This implies all electron transfer reactions must be accompanied by ion transfer to maintain charge balance within a phase. However, due to the small domain size, electroneutrality is known to break down in the electric double layer adjacent to the electrode. Therefore, the appropriateness of the electroneutrality approximation for media with nanoscale domains, such as microemulsions, is questionable. Here, in this work, we show definitively, for the first time, that phase electroneutrality can be violated during an electron transfer reaction in microemulsions. Cyclic voltammograms show that electrochemical oxidation of rubrene is possible in microemulsions when sodium hydroxide is the only supporting electrolyte. Rubrene is a hydrophobic redox compound that is insoluble in water in neutral and radical cation forms, and hydroxide anions are too hydrophilic to transfer to the oil (toluene) phase. Contrary to observations made in emulsion systems, rubrene oxidation in microemulsions seemingly occurs by a mechanism where coupled ion transfer required for rigorous phase electroneutrality is negated. Our results demonstrate that electron transfer reactions within nanometer-sized domains are not subject to restrictions that govern systems with larger domains.

Barth, Brian A. [Univ. of Tennessee, Knoxville, TN↗
DOE OSTI2020

Breakdown of electroneutrality in nanopores

Ion transport in extremely narrow nanochannels has gained increasing interest in recent years due to unique physical properties at the nanoscale and the technological advances that allow us to study them. It is tempting to approach this confined regime with the theoretical tools and knowledge developed for membranes and microfluidic devices, and naively apply continuum models, such as the Poisson-Nernst-Planck and Navier-Stokes equations. However, it turns out that some of the most basic principles we take for granted in larger systems, such as the complete screening of surface charge by counter-ions, can break down under extreme confinement. Additionally, we show that in a truly one-dimensional system of ions interacting with three-dimensional electrostatic interactions, the screening length is exponentially large, and can easily exceed the macroscopic length of a nanotube. Without screening, electroneutrality breaks down within the nanotube, with fundamental consequences for ion transport and electrokinetic phenomena. In this work, we build a general theoretical framework for electroneutrality breakdown in nanopores, focusing on the most interesting case of a one-dimensional nanotube, and show how it provides an elegant interpretation for the peculiar scaling observed in experimental measurements of ionic conductance in carbon nanotubes.

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

Interfacial electroneutrality controls transport of asymmetric salts through charge-patterned mosaic membranes

Membranes that selectively enhance target solute permeation while rejecting competing species are essential for precision separations. This study introduces charge-patterned mosaic membranes (CMMs) that selectively transport divalent asymmetric salts by leveraging a net-neutral membrane–solution interface. This mechanism, dictated by the charge ratio of positive and negative domains on the membrane surface and the balance of cations and anions in the salt, is supported by analytical, numerical, and experimental results. Analytical solutions identified cationic domain coverages (f + ) of 33%, 50%, and 66% as optimal for the selective transport of +2:−1 salts, +1:−1 salts, and +1:−2 salts, respectively, under conditions where the pattern size (L) is significantly larger than the Debye length. Numerical simulations and experiments using CMMs with alternating charged-stripes inkjet-printed onto nanostructure copolymer substrates confirmed these findings. By varying stripe widths to control f + , pressure-driven filtration experiments demonstrated selective enrichment of MgCl 2 and K 2 SO 4 at the predicted f + values, with deviations from these values leading to salt rejection. These results highlight the pivotal role of a net-neutral interface in enabling asymmetric salt enrichment. This study positions CMMs as a versatile platform for tuning ion selectivity, addressing challenges in resource recovery, water treatment, and precision separations.

additive manufacturing↗
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 OSTI2020

Molecular dynamics and continuum analyses of the electrokinetic zeta potential in nanostructured slit channels

This work presents a theoretical and numerical study of electrokinetic flow and the zeta potential for the case of slit channels with nanoscale surface features of dimensions comparable to the Debye length, by employing molecular dynamics simulations and continuum-level analyses. Additionally, a simple analytical model for considering the average effect of such surface features producing nanoscale roughness is proposed by employing matched asymptotic solutions for the charge density and fluid flow field, and matching conditions that satisfy electroneutrality and the Onsager reciprocal relation between the electroosmotic flow rate and streaming current. The proposed analytical model quantitatively accounts for results from molecular dynamics simulations that consider the presence of ion solvation shells and surface hydration layers. Our analysis indicates that a simultaneous knowledge of the electroosmotic and pressure-driven flow rate or streaming current can be instrumental to unambiguously determine the zeta potential in the presence of nanoscale surface features.

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

Phase-field modeling of solid-state metathesis reactions with the charge neutrality constraint

In this work, we present a phase-field model that captures the evolution of ionic concentrations and phase fractions during solid-state metathesis (SSM) reactions where diffusion limits the rate of transformation. The evolution of the mole fraction of each ion is obtained via governing equations that describe the reduction of a free energy, which includes an energy landscape with local minima located at compositions corresponding to stable products. We utilized two Lagrange multipliers to impose constraints of electroneutrality as well as on the sum of mole fractions, which were then eliminated to derive set of two partial differential equations that describe the dynamics of the mole fraction evolution. From these governing equations, the expressions for effective mobilities for the cations and the anions were obtained. We first study the effect of mobilities of ions on the reaction kinetics, using a simple model considering the ions with an identical absolute value of charge numbers. The simulation results show that the overall characteristic mobility, defined as the sum of the two effective ionic mobilities, provides an excellent measure of the rate at which reaction progresses and that the ratio of the effective mobilities of the anions and the cations signifies the manner by which the reaction progresses. We then generalize the model to consider ions with different charge numbers and tuned the mobility of ions based on their diffusion coefficients reported in the literature and experimental data from a thin-film experiment for the synthesis of FeS 2 to demonstrate the capability of the model to predict the phase evolution during SSM reactions. In particular, the simulation predicts nonplanar phase evolution, which is recently observed in thin-film reactions for the synthesis of FeS 2 via transmission electron microscopy. The approach can serve as a basis for models for phase transformations in other multiphase ionic mixtures, such as in all-solid-state batteries and in ionic liquids.

36 MATERIALS SCIENCE↗
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 OSTI2025

Selective concentration of multivalent cations via electrodialysis through tandem membranes

Selective electrodialysis is attractive for recovering and recycling valuable ions. However, most studies in this field focus on membranes that are selective for monovalent ions. Multivalent cations preferentially partition into cation-exchange membranes, but low mobilities and concentration polarization limit electrodialysis selectivities for these cations. This work employs electrodialysis through an “in-series” combination of a Nafion membrane and a multilayer polyelectrolyte film (MPF)-coated Nafion membrane to concentrate and isolate La 3+ or Mg 2+ from monovalent ions. The bare Nafion membrane gives a La 3+ /Na + selectivity up to 4 (depending on conditions), which enables concentration of La 3+ with a high current efficiency. The MPF-coated membrane allows selective passage of Na + while rejecting La 3+ to increase the purity of the concentrated multivalent cation. For example, the concentration of La 3+ in the accumulation solution can reach 9 times its original concentration in the feed solution with a purity of 93 mol% (99 wt%) when starting with equimolar mixtures of NaNO 3 and La(NO 3 ) 3 . This study examines how the separation varies with the type of polyelectrolyte, the applied voltage, and the composition of feed, accumulation, and cathode solutions. Here, a high nitric acid concentration (0.5 M) in the cathode compartment allows monovalent anions to pass through the MPF-coated Nafion membrane to maintain the electroneutrality of the accumulation solution. With a low acid concentration (0.04 M) in the cathode chamber, water splitting occurs at the MPF-coated membrane and La(OH) 3 precipitates.

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

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N 1888 ][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the shortand long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ~ 0.44 Å –1 ), charge (k ~ 0.75 Å –1 ), and neighboring or adjacent (k ~ 1.46 Å –1 ) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N + ···N – and N + ···O – , compete for optimal packing around the polar head of the [N 1888 ] + cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, S ZZ , and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N 1888 ][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

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

Structural basis for inhibition of the Cation-chloride cotransporter $\mathrm{NKCC1}$ by the diuretic drug bumetanide

Cation-chloride cotransporters (CCCs) NKCC1 and NKCC2 catalyze electroneutral symport of 1 Na + , 1K + , and 2 Cl - across cell membranes. NKCC1 mediates trans-epithelial Cl - secretion and regulates excitability of some neurons and NKCC2 is critical to renal salt reabsorption. Both transporters are inhibited by the so-called loop diuretics including bumetanide, and these drugs are a mainstay for treating edema and hypertension. Here, our single-particle electron cryo-microscopy structures supported by functional studies reveal an outward-facing conformation of NKCC1, showing bumetanide wedged into a pocket in the extracellular ion translocation pathway. Based on these and the previously published inward-facing structures, we define the translocation pathway and the conformational changes necessary for ion translocation. We also identify an NKCC1 dimer with separated transmembrane domains and extensive transmembrane and C-terminal domain interactions. We further define an N-terminal phosphoregulatory domain that interacts with the C-terminal domain, suggesting a mechanism whereby (de)phosphorylation regulates NKCC1 by tuning the strength of this domain association.

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

Impact of surface hydrophilicity on the ordering and transport properties of bicontinuous microemulsions

Microemulsions (MEs) have many industrial applications, where recent developments have shown that MEs can be utilized for electrochemical applications, including potentially in redox flow batteries. However, understanding the structure and dynamics of these systems, including at a surface, is needed to direct and rationally control their electrochemical behavior. While bulk solution measurements have provided insight into their structure, their assembly at an interface also impacts the electron (to the electrode) and ion (across the surfactant) charge transfer processes in the system. To address this shortcoming, neutron reflectivity experiments and molecular simulations have been completed that document the near surface structure of a series of deuterated water (D 2 O)/polysorbate-20/toluene MEs on hydrophilic and amphiphilic surfaces. These results show that the microemulsions form complex layered structures near a hard electrode surface, where most layers are not purely one component. Decreasing the D 2 O concentration in the ME increases the number of and purity of the layers established on the solid surface. These lamellar-type layers transition from the surface to the bulk microemulsion as a series of mixed layers (i.e., containing oil, water, and surfactant) that are consistent with perforated lamellae. Additionally, these mixed lamellae appear to become more perforated with oil and water pathways on an amphiphilic surface. Furthermore, the purity and thickness of these layers will influence the accessibility of an electrode by redox active species, as well as ion transport required to satisfy the electroneutrality condition.

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

Structure of the human cation–chloride cotransport KCC1 in an outward-open state

Cation–chloride cotransporters (CCCs) catalyze electroneutral symport of Cl - with Na + and/or K + across membranes. CCCs are fundamental in cell volume homeostasis, transepithelia ion movement, maintenance of intracellular Cl - concentration, and neuronal excitability. Here, we present a cryoelectron microscopy structure of human K + –Cl - cotransporter (KCC)1 bound with the VU0463271 inhibitor in an outward-open state. In contrast to many other amino acid–polyamine–organocation transporter cousins, our first outward-open CCC structure reveals that opening the KCC1 extracellular ion permeation path does not involve hinge-bending motions of the transmembrane (TM) 1 and TM6 half-helices. Instead, rocking of TM3 and TM8, together with displacements of TM4, TM9, and a conserved intracellular loop 1 helix, underlie alternate opening and closing of extracellular and cytoplasmic vestibules. We show that KCC1 intriguingly exists in one of two distinct dimeric states via different intersubunit interfaces. Our studies provide a blueprint for understanding the mechanisms of CCCs and their inhibition by small molecule compounds.

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

Suppression of electroconvective and morphological instabilities by an imposed cross flow of the electrolyte

Electroconvective (EC) instability and its influence on surface morphological perturbations (Morph) are important in many applications, including electrodialysis, batteries, and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the EC and Morph instabilities using two approaches. In the bulk analysis, we derive the asymptotic solutions for small and large wave numbers by neglecting the space charge layer and imposing a second kind electroosmosis slip velocity boundary condition on the electroneutral bulk region. In the full analysis, the instability of the entire region of the liquid electrolyte is analyzed using the ultraspherical spectral method. Both studies show that the flow significantly affects the EC instability. The imposed flow distorts the concentration field, causing a sheltering effect which hinders the ion transport from low- to high-concentration regions, and therefore suppresses the EC instability below a certain wave number. In combination with the viscous stabilization of the high wave number modes by the space charge layer, a sufficiently strong imposed flow can fully suppress the EC instability. The increment in the critical voltage for the instability onset is roughly proportional to the square root of the product of imposed velocity and double layer thickness. The imposed flow has a smaller effect on the Morph instability, except that it may remove the morphological modes resulting from the EC instability and thereby change the wave number of the most unstable mode.

42 ENGINEERING↗