Electrokinetic instability due to streamwise conductivity gradients in microchip electrophoresis

We present an experimental and numerical investigation of electrokinetic instability (EKI) in microchannel flow with streamwise conductivity gradients, such as those observed during sample stacking in capillary electrophoresis. A plug of a low-conductivity electrolyte solution is initially sandwiched between two high-conductivity zones in a microchannel. This spatial conductivity gradient is subjected to an external electric field applied along the microchannel axis, and for sufficiently strong electric fields an instability sets in. We have explored the physics of this EKI through experiments and numerical simulations, and supplemented the results using scaling analysis. We performed EKI experiments at different electric field values and visualised the flow using a passive fluorescent tracer. The experimental data were analysed using the proper orthogonal decomposition technique to obtain a quantitative measure of the threshold electric field for the onset of instability, along with the corresponding coherent structures. To elucidate the physical mechanism underlying the instability, we performed high-resolution numerical simulations of ion transport coupled with fluid flow driven by the electric body force. Simulations reveal that the non-uniform electroosmotic flow due to axially varying conductivity field causes a recirculating flow within the low-conductivity region, and creates a new configuration wherein the local conductivity gradients are orthogonal to the applied electric field. This configuration leads to EKI above a threshold electric field. The spatial features of the instability predicted by the simulations and the threshold electric field are in good agreement with the experimental observations and provide useful insight into the underlying mechanism of instability.

Transitions and Instabilities in Imperfect Ion-Selective Membranes

Numerical investigation of the underlimiting, limiting, and overlimiting current modes and their transitions in imperfect ion-selective membranes with fluid flow through permitted through the membrane is presented. The system is treated as a three layer composite system of electrolyte-porous membrane-electrolyte where the Nernst–Planck–Poisson–Stokes system of equations is used in the electrolyte, and the Darcy–Brinkman approach is employed in the nanoporous membrane. In order to resolve thin Debye and Darcy layers, quasi-spectral methods are applied using Chebyshev polynomials for their accumulation of zeros and, hence, best resolution in the layers. The boundary between underlimiting and overlimiting current regimes is subject of linear stability analysis, where the transition to overlimiting current is assumed due to the electrokinetic instability of the one-dimensional quiescent state. However, the well-developed overlimiting current is inherently a problem of nonlinear stability and is subject of the direct numerical simulation of the full system of equations. Both high and low fixed charge density membranes (low- and high concentration electrolyte solutions), acting respectively as (nearly) perfect or imperfect membranes, are considered. The perfect membrane is adequately described by a one-layer model while the imperfect membrane has a more sophisticated response. In particular, the direct transition from underlimiting to overlimiting currents, bypassing the limiting currents, is found to be possible for imperfect membranes (high-concentration electrolyte). The transition to the overlimiting currents for the low-concentration electrolyte solutions is monotonic, while for the high-concentration solutions it is oscillatory. Despite the fact that velocities in the porous membrane are much smaller than in the electrolyte region, it is further demonstrated that they can dramatically influence the nature and transition to the overlimiting regimes. A map of the bifurcations, transitions, and regimes is constructed in coordinates of the fixed membrane charge and the Darcy number.