Stabilized lattice Boltzmann-Enskog method for compressible flows and its application to one- and two-component fluids in nanochannels

Electrokinetic Lattice Boltzmann Solver Coupled to Molecular Dynamics: Application to Polymer Translocation

We have developed a theoretical and computational approach to deal with systems that involve a disparate range of spatiotemporal scales, such as those composed of colloidal particles or polymers moving in a fluidic molecular environment. Our approach is based on a multiscale modeling that combines the slow dynamics of the large particles with the fast dynamics of the solvent into a unique framework. The former is numerically solved via Molecular Dynamics and the latter via a multicomponent Lattice Boltzmann. The two techniques are coupled together to allow for a seamless exchange of information between the descriptions. Being based on a kinetic multicomponent description of the fluid species, the scheme is flexible in modeling charge flow within complex geometries and ranging from large to vanishing salt concentration. The details of the scheme are presented and the method is applied to the problem of translocation of a charged polymer through a nanopores. Lastly, we discuss the advantages and complexities of the approach.

Charge Relaxation Dynamics of an Electrolytic Nanocapacitor

Understanding ion relaxation dynamics in overlapping electric double layers (EDLs) is critical for the development of efficient nanotechnology-based electrochemical energy storage, electrochemomechanical energy conversion, and bioelectrochemical sensing devices as well as the controlled synthesis of nanostructured materials. Here, a lattice Boltzmann (LB) method is employed to simulate an electrolytic nanocapacitor subjected to a step potential at t = 0 for various degrees of EDL overlap, solvent viscosities, ratios of cation-to-anion diffusivity, and electrode separations. The use of a novel continuously varying and Galilean-invariant molecular-speed-dependent relaxation time (MSDRT) with the LB equation recovers a correct microscopic description of the molecular-collision phenomena and enhances the stability of the LB algorithm. Results for large EDL overlaps indicated oscillatory behavior for the ionic current density, in contrast to monotonic relaxation to equilibrium for low EDL overlaps. Further, at low solvent viscosities and large EDL overlaps, anomalous plasmalike spatial oscillations of the electric field were observed that appeared to be purely an effect of nanoscale confinement. Employing MSDRT in our simulations enabled modeling of the fundamental physics of the transient charge relaxation dynamics in electrochemical systems operating away from equilibrium wherein Nernst-Einstein relation is known to be violated.

Charge transport in nanochannels: a molecular theory

We introduce a theoretical and numerical method to investigate the flow of charged fluid mixtures under extreme confinement. We model the electrolyte solution as a ternary mixture comprising two ionic species of opposite charge and a third uncharged component. The microscopy approach is based on kinetic theory and is fully self-consistent. It allows us to determine configurational properties, such as layering near the confining walls and the flow properties. We show that, under the appropriate assumptions, the approach reproduces the phenomenological equations used to describe electrokinetic phenomena, without requiring the introduction of constitutive equations to determine the fluxes. Moreover, we model channels of arbitrary shape and nanometric roughness, features that have important repercussions on the transport properties of these systems. Numerical simulations are obtained by solving the evolution dynamics of the one-particle phase-space distributions of each species by means of a lattice Boltzmann method for flows in straight and wedged channels. Results are presented for the microscopic density, the velocity profiles, and the volumetric and charge flow rates. Strong departures from electroneutrality are shown to appear at the molecular level.