Modeling Electrokinetic Flows with the Discrete Ion Stochastic Continuum Overdamped Solvent Algorithm
arXiv:2204.14167 · doi:10.1103/PhysRevE.106.035104
Abstract
In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the Discrete Ion Stochastic Continuum Overdamped Solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-H{ü}ckel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a `dry diffusion' approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.
arXiv comment27 pages, 15 figures
References in corpus (16)
- Induced-Charge Electro-Osmosis
- A Stochastic Immersed Boundary Method for Fluid-Structure Dynamics at Microscopic Length Scales
- Stochastic Eulerian Lagrangian Methods for Fluid-Structure Interactions with Thermal Fluctuations
- An Immersed Boundary Method for Rigid Bodies
- Brownian Dynamics without Green's Functions
- Gaussian-Like Immersed Boundary Kernels with Three Continuous Derivatives and Improved Translational Invariance
- Fluctuating force-coupling method for simulations of colloidal suspensions
- Efficient Variable-Coefficient Finite-Volume Stokes Solvers
- Low Mach Number Fluctuating Hydrodynamics for Electrolytes
- Simulating Brownian suspensions with fluctuating hydrodynamics
- Coupling a nano-particle with isothermal fluctuating hydrodynamics: Coarse-graining from microscopic to mesoscopic dynamics
- A fast spectral method for electrostatics in doubly-periodic slit channels
- Low Mach Number Fluctuating Hydrodynamics of Binary Liquid Mixtures
- Spatially Adaptive Stochastic Methods for Fluid-Structure Interactions Subject to Thermal Fluctuations in Domains with Complex Geometries
- Fluctuating hydrodynamics of electrolytes at electroneutral scales
- A Discrete Ion Stochastic Continuum Overdamped Solvent Algorithm for Modeling Electrolytes