Fluctuating hydrodynamics of electrolytes at electroneutral scales

Aleksandar Donev, Andrew J. Nonaka, Changho Kim, Alejandro L. Garcia, John B. Bell

Research output: Contribution to journalArticle

Abstract

At mesoscopic scales electrolyte solutions are modeled by the fluctuating generalized Poisson-Nernst-Planck (PNP) equations [J.-P. Péraud, Phys. Rev. Fluids 1, 074103 (2016)2469-990X10.1103/PhysRevFluids.1.074103]. However, at length and time scales larger than the Debye scales, electrolytes are effectively electroneutral and the charged-fluid PNP equations become too stiff to solve numerically. Here we formulate the isothermal incompressible equations of fluctuating hydrodynamics for reactive multispecies mixtures involving charged species in the electroneutral limit and design a numerical algorithm to solve these equations. Our model does not assume a dilute electrolyte solution but rather treats all species on an equal footing, accounting for cross diffusion and nonideality using Maxwell-Stefan theory. By enforcing local electroneutrality as a constraint, we obtain an elliptic equation for the electric potential that replaces the Poisson equation in the fluctuating PNP equations. We develop a second-order midpoint predictor-corrector algorithm to solve either the charged-fluid or electroneutral equations with only a change of the elliptic solver. We use the electroneutral algorithm to study a gravitational fingering instability, triggered by thermal fluctuations, at an interface where an acid and base react to neutralize each other. Our results demonstrate that, because the four ions diffuse with very different coefficients, one must treat each ion as an individual species and cannot treat the acid, base, and salt as neutral species. This emphasizes the differences between electrodiffusion and classical Fickian diffusion, even at electroneutral scales.

Original languageEnglish (US)
Article number043701
JournalPhysical Review Fluids
Volume4
Issue number4
DOIs
StatePublished - Apr 1 2019

Fingerprint

Fluctuating Hydrodynamics
Electrolyte
Electrolytes
Hydrodynamics
Fluids
Ions
Siméon Denis Poisson
Acids
Poisson equation
Fluid
Salts
Cross-diffusion
Predictor-corrector
Midpoint
Electric Potential
Electric potential
Poisson's equation
Salt
Numerical Algorithms
Length Scale

ASJC Scopus subject areas

  • Computational Mechanics
  • Modeling and Simulation
  • Fluid Flow and Transfer Processes

Cite this

Fluctuating hydrodynamics of electrolytes at electroneutral scales. / Donev, Aleksandar; Nonaka, Andrew J.; Kim, Changho; Garcia, Alejandro L.; Bell, John B.

In: Physical Review Fluids, Vol. 4, No. 4, 043701, 01.04.2019.

Research output: Contribution to journalArticle

Donev, Aleksandar ; Nonaka, Andrew J. ; Kim, Changho ; Garcia, Alejandro L. ; Bell, John B. / Fluctuating hydrodynamics of electrolytes at electroneutral scales. In: Physical Review Fluids. 2019 ; Vol. 4, No. 4.
@article{4985fde3a5a54a1e89fd4ae4049e50ae,
title = "Fluctuating hydrodynamics of electrolytes at electroneutral scales",
abstract = "At mesoscopic scales electrolyte solutions are modeled by the fluctuating generalized Poisson-Nernst-Planck (PNP) equations [J.-P. P{\'e}raud, Phys. Rev. Fluids 1, 074103 (2016)2469-990X10.1103/PhysRevFluids.1.074103]. However, at length and time scales larger than the Debye scales, electrolytes are effectively electroneutral and the charged-fluid PNP equations become too stiff to solve numerically. Here we formulate the isothermal incompressible equations of fluctuating hydrodynamics for reactive multispecies mixtures involving charged species in the electroneutral limit and design a numerical algorithm to solve these equations. Our model does not assume a dilute electrolyte solution but rather treats all species on an equal footing, accounting for cross diffusion and nonideality using Maxwell-Stefan theory. By enforcing local electroneutrality as a constraint, we obtain an elliptic equation for the electric potential that replaces the Poisson equation in the fluctuating PNP equations. We develop a second-order midpoint predictor-corrector algorithm to solve either the charged-fluid or electroneutral equations with only a change of the elliptic solver. We use the electroneutral algorithm to study a gravitational fingering instability, triggered by thermal fluctuations, at an interface where an acid and base react to neutralize each other. Our results demonstrate that, because the four ions diffuse with very different coefficients, one must treat each ion as an individual species and cannot treat the acid, base, and salt as neutral species. This emphasizes the differences between electrodiffusion and classical Fickian diffusion, even at electroneutral scales.",
author = "Aleksandar Donev and Nonaka, {Andrew J.} and Changho Kim and Garcia, {Alejandro L.} and Bell, {John B.}",
year = "2019",
month = "4",
day = "1",
doi = "10.1103/PhysRevFluids.4.043701",
language = "English (US)",
volume = "4",
journal = "Physical Review Fluids",
issn = "2469-990X",
publisher = "American Physical Society",
number = "4",

}

TY - JOUR

T1 - Fluctuating hydrodynamics of electrolytes at electroneutral scales

AU - Donev, Aleksandar

AU - Nonaka, Andrew J.

AU - Kim, Changho

AU - Garcia, Alejandro L.

AU - Bell, John B.

PY - 2019/4/1

Y1 - 2019/4/1

N2 - At mesoscopic scales electrolyte solutions are modeled by the fluctuating generalized Poisson-Nernst-Planck (PNP) equations [J.-P. Péraud, Phys. Rev. Fluids 1, 074103 (2016)2469-990X10.1103/PhysRevFluids.1.074103]. However, at length and time scales larger than the Debye scales, electrolytes are effectively electroneutral and the charged-fluid PNP equations become too stiff to solve numerically. Here we formulate the isothermal incompressible equations of fluctuating hydrodynamics for reactive multispecies mixtures involving charged species in the electroneutral limit and design a numerical algorithm to solve these equations. Our model does not assume a dilute electrolyte solution but rather treats all species on an equal footing, accounting for cross diffusion and nonideality using Maxwell-Stefan theory. By enforcing local electroneutrality as a constraint, we obtain an elliptic equation for the electric potential that replaces the Poisson equation in the fluctuating PNP equations. We develop a second-order midpoint predictor-corrector algorithm to solve either the charged-fluid or electroneutral equations with only a change of the elliptic solver. We use the electroneutral algorithm to study a gravitational fingering instability, triggered by thermal fluctuations, at an interface where an acid and base react to neutralize each other. Our results demonstrate that, because the four ions diffuse with very different coefficients, one must treat each ion as an individual species and cannot treat the acid, base, and salt as neutral species. This emphasizes the differences between electrodiffusion and classical Fickian diffusion, even at electroneutral scales.

AB - At mesoscopic scales electrolyte solutions are modeled by the fluctuating generalized Poisson-Nernst-Planck (PNP) equations [J.-P. Péraud, Phys. Rev. Fluids 1, 074103 (2016)2469-990X10.1103/PhysRevFluids.1.074103]. However, at length and time scales larger than the Debye scales, electrolytes are effectively electroneutral and the charged-fluid PNP equations become too stiff to solve numerically. Here we formulate the isothermal incompressible equations of fluctuating hydrodynamics for reactive multispecies mixtures involving charged species in the electroneutral limit and design a numerical algorithm to solve these equations. Our model does not assume a dilute electrolyte solution but rather treats all species on an equal footing, accounting for cross diffusion and nonideality using Maxwell-Stefan theory. By enforcing local electroneutrality as a constraint, we obtain an elliptic equation for the electric potential that replaces the Poisson equation in the fluctuating PNP equations. We develop a second-order midpoint predictor-corrector algorithm to solve either the charged-fluid or electroneutral equations with only a change of the elliptic solver. We use the electroneutral algorithm to study a gravitational fingering instability, triggered by thermal fluctuations, at an interface where an acid and base react to neutralize each other. Our results demonstrate that, because the four ions diffuse with very different coefficients, one must treat each ion as an individual species and cannot treat the acid, base, and salt as neutral species. This emphasizes the differences between electrodiffusion and classical Fickian diffusion, even at electroneutral scales.

UR - http://www.scopus.com/inward/record.url?scp=85065045569&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85065045569&partnerID=8YFLogxK

U2 - 10.1103/PhysRevFluids.4.043701

DO - 10.1103/PhysRevFluids.4.043701

M3 - Article

AN - SCOPUS:85065045569

VL - 4

JO - Physical Review Fluids

JF - Physical Review Fluids

SN - 2469-990X

IS - 4

M1 - 043701

ER -