Three-dimensional coherent structures of electrokinetic instability

Electrokinetic and Electroconvective Effects in Ternary Electrolyte Near Ion-Selective Microsphere

The paper presents theoretical and experimental investigations of the behavior of an electrolyte solution with three types of ions near an ion-selective microparticle with electrokinetically and pressure-driven flow. A special experimental cell has been developed for the investigations. An anion-selective spherical particle composed of ion-exchange resin is fixed in the center of the cell. An enriched region with a high salt concentration appears at the anode side of the particle when an electric field is turned on, according to the nonequilibrium electrosmosis behavior. A similar region exists near a flat anion-selective membrane. However, the enriched region near the particle produces a concentration jet that spreads downstream akin to a wake behind an axisymmetrical body. The fluorescent cations of Rhodamine-6G dye are chosen as the third species in the experiments. The ions of Rhodamine-6G have a 10-fold lower diffusion coefficient than the ions of potassium while bearing the same valency. This paper shows that the concentration jet behavior is described accurately enough with the mathematical model of a far axisymmetric wake behind a body in a fluid flow. The third species also forms an enriched jet, but its distribution turns out to be more complex. The concentration of the third species increases in the jet with an increase in pressure gradient. The pressure-driven flow stabilizes the jet, yet electroconvection has been observed near the microparticle for sufficiently strong electric fields. The electrokinetic instability and the electroconvection partially destroy the concentration jet of salt and the third species. The conducted experiments show good qualitative agreement with the numerical simulations. The presented results could be used in future for implementing microdevices based on membrane technology for solving problems of detection and preconcentration, and thus simplifying chemical and medical analyses utilizing the superconcentration phenomenon. Such devices are called membrane sensors, and are actively being studied.

Novel electroosmotic micromixer configuration based on ion-selective microsphere

In this paper, a micromixer of a new configuration is presented, consisting of a spherical chamber in the center of which an ion-selective microsphere is placed. Stratified liquid is introduced through the chamber via inlet and outlet holes under an external pressure gradient and an external electric field is directed in such a way that the resulting electroosmotic flow is directed against the pressure-driven flow, resulting in mixing. The investigation is carried out by direct numerical simulation on a super-computer. Optimal values of the applied electric field are determined to yield strong mixing. Above this optimal mixing regime, a number of instabilities and bifurcations are realized, which qualitatively coincide with those occurring during electrophoresis of an ion-selective microgranule. As shown by our calculation, these instabilities do not lead to an enhanced mixing. The resulting electroconvective vortices remain confined near the surface of the microgranule, and do not sufficiently perturb the stratified fluid flow further from the granule. On the other hand, another type of instability caused by the salt concentration gradient can generate sufficiently strong oscillations to enhance mixing. However, this only occurs when the external electric field is sufficiently high that the electroosmotic flow is comparable to the pressure-driven flow. This ultimately leads to creation of reverse flows of the liquid and cessation of the device operation. Thus, it was shown that the best mixing occurs in the absence of electrokinetic instability. Based on the data obtained, it is possible to select the necessary geometric characteristics of the micromixer to achieve the optimal mixing mode for a given set of liquids, which may be ten times more effective than passive mixers at the same flow rates. A comparison with the experimental data of the other authors confirms the effectiveness of this device and its other capabilities. Furthermore, the basic device design can be operated in other modes, for example, an electrohydrodynamic pump, a streaming current generator, or even a micro-reactor, depending on the system parameters and choice of an ion-selective granule.

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.

Scaling laws of electroconvective flow with finite vortex height near permselective membranes

In a steady state, the linear scaling laws are confirmed between the intensity characteristics of electroconvective (EC) vortex (including the vortex height and electroosmotic slip velocity) and the applied voltage for the nonshear EC flow with finite vortex height near permselective membranes. This finding in the nonshear EC flow is different from the shear EC flow [Kwak et al., Phys. Rev. Lett. 110, 114501 (2013)10.1103/PhysRevLett.110.114501] and indicates that the local concentration gradient has a significant improvement in the analysis of slip velocity. Further, our study reveals that the EC vortex is mainly driven by the second peak effect of the Coulomb thrust in the extended space-charge layer, and the linear scaling law exhibited by the Coulomb thrust is an essential reason for the linear scaling laws of vortex intensity. The scaling laws proposed in this paper are supported by our direct numerical simulation data and previous experimental observations [Rubinstein et al., Phys. Rev. Lett. 101, 236101 (2008)10.1103/PhysRevLett.101.236101].

Pattern Formation of Three-Dimensional Electroconvection on a Charge Selective Surface

When a charge selective surface consumes or transports only cations or anions in the electrolyte, biased ion rejection initiates hydrodynamic instability, resulting in vortical fluid motions called electroconvection. In this Letter, we describe the first laboratory observation of three-dimensional electroconvection on a charge selective surface. Combining experiment and scaling analysis, we successfully categorized three distinct patterns of 3D electroconvection according to [(Ra_{E})/(Re^{2}Sc)] [electric Rayleigh number (Ra_{E}), Reynolds number (Re), Schmidt number (Sc)] as (i) polygonal, (ii) transverse, or (iii) longitudinal rolls. If Re increases or Ra_{E} decreases, pure longitudinal rolls are presented. On the other hand, transverse rolls are formed between longitudinal rolls, and two rolls are transformed as polygonal one at higher Ra_{E} or lower Re. In this pattern selection scenario, Sc determines the critical electric Rayleigh number (Ra_{E}^{*}) for the onset of each roll, resulting in Ra_{E}^{*}∼Re^{2}Sc. We also verify that convective ion flux by electroconvection (represented by an electric Nusselt number Nu_{E}) is fitted to a power law, Nu_{E}∼[(Ra_{E}-Ra_{E}^{*})/(Re^{2}Sc)]^{α_{1}}Re^{α_{2}}Pe^{α_{3}} [Péclet number (Pe)], where each term represents the characteristics of electroconvection, shear flow, and ion transport.

Controllable pH Manipulations in Micro/Nanofluidic Device Using Nanoscale Electrokinetics

Recently introduced nanoscale electrokinetic phenomenon called ion concentration polarization (ICP) has been suffered from serious pH changes to the sample fluid. A number of studies have focused on the origin of pH changes and strategies for regulating it. Instead of avoiding pH changes, in this work, we tried to demonstrate new ways to utilize this inevitable pH change. First, one can obtain a well-defined pH gradient in proton-received microchannel by applying a fixed electric current through a proton exchange membrane. Furthermore, one can tune the pH gradient on demand by adjusting the proton mass transportation (i.e., adjusting electric current). Secondly, we demonstrated that the occurrence of ICP can be examined by sensing a surrounding pH of electrolyte solution. When pH > threshold pH, patterned pH-responsive hydrogel inside a straight microchannel acted as a nanojunction to block the microchannel, while it did as a microjunction when pH < threshold pH. In case of forming a nanojunction, electrical current significantly dropped compared to the case of a microjunction. The strategies that presented in this work would be a basis for useful engineering applications such as a localized pH stimulation to biomolecules using tunable pH gradient generation and portable pH sensor with pH-sensitive hydrogel.

Electro-Kinetic Instability in a Laminar Boundary Layer Next to an Ion Exchange Membrane

The electro-kinetic instability in a pressure driven shear flow near an ion exchange membrane is considered. The electrochemical system, through which an electrical potential drop is applied, consists in a polarization layer in contact with the membrane and a bulk. The numerical investigation contained two aspects: analysis of the instability modes and description of the Lagrangian transport of fluid and ions. Regarding the first aspect, the modes were analyzed as a function of the potential drop. The analysis revealed how the spatial distribution of forces controls the dynamics of vortex association and dissociation. In particular, the birth of a counter-clockwise vortex between two clockwise vortices, and the initiation of clusters constituting one or two envelopes wrapping a vortex group, were examined. In regards to the second aspect, the trajectories were computed with the fourth order Runge Kutta scheme for the time integration and with the biquadratric upstream scheme for the spatial and time interpolation of the fluid velocity and the ion flux. The results for the periodic mode showed two kinds of trajectories: the trochoidal motion and the longitudinal one coupled with a periodic transverse motion. For the aperiodic modes, other mechanisms appeared, such as ejection from the mixing layer, trapping by a growing vortex or merging vortices. The analysis of the local velocity field, the vortices’ shape, the spatial distribution of the forces and the ion flux components explained these trajectories.

Electroconvection in a Viscoelastic Electrolyte

Direct numerical simulations of a liquid electrolyte with polymer additives demonstrate that viscoelasticity promotes an earlier transition from steady to unsteady electroconvective flow. Viscoelasticity also decreases the overlimiting current resulting from convection by up to 40%. Both of these effects would reduce the time-averaged spatial variability of ion flux suggesting that polymeric fluids may inhibit dendrite growth. Polymer relaxation near a surface destabilizes the flow structures and decreases the time duration of high current fluxes. This mechanism of polymer-induced flux reduction is general to wall bounded flows with transfer of mass, heat or momentum.

Modes of electrokinetic instability for imperfect electric membranes

The direct transition to overlimiting current bypassing the stage of limiting currents is considered for imperfect membranes. Instability of the quiescent steady-state one-dimensional solution, which is the result of a balance of diffusion and electromigration, is investigated on the basis of the full Nernst-Planck-Poisson-Stokes system and a simplified quasielectroneutral system. A three-layer geometry, electrolyte-porous membrane-electrolyte, is considered. The usual assumption of a constant electrochemical potential along the membrane surface is removed from consideration. The effect of bulk and surface effects on the instability and transition to the overlimiting currents is evaluated for a different membrane selectivity. It becomes clear that for sufficiently small fixed charge concentration (large ion concentration in the electrolyte), the monotonic instability is replaced by an oscillatory one. The dependence of instability on the membrane porosity is found to be weak.

Numerical simulation of electrochemical desalination

We present an effective numerical approach to simulate electrochemically mediated desalination of seawater. This new membraneless, energy efficient desalination method relies on the oxidation of chloride ions, which generates an ion depletion zone and local electric field gradient near the junction of a microchannel branch to redirect sea salt into the brine stream, consequently producing desalted water. The proposed numerical model is based on resolution of the 3D coupled Navier-Stokes, Nernst-Planck, and Poisson equations at non-uniform spatial grids. The model is implemented as a parallel code and can be employed to simulate mass-charge transport coupled with surface or volume reactions in 3D systems showing an arbitrarily complex geometrical configuration.

Helical vortex formation in three-dimensional electrochemical systems with ion-selective membranes

The rate of electric-field-driven transport across ion-selective membranes can exceed the limit predicted by Nernst (the limiting current), and encouraging this "overlimiting" phenomenon can improve efficiency in many electrochemical systems. Overlimiting behavior is the result of electroconvectively induced vortex formation near membrane surfaces, a conclusion supported so far by two-dimensional (2D) theory and numerical simulation, as well as experiments. In this paper we show that the third dimension plays a critical role in overlimiting behavior. In particular, the vortex pattern in shear flow through wider channels is helical rather than planar, a surprising result first observed in three-dimensional (3D) simulation and then verified experimentally. We present a complete experimental and numerical characterization of a device exhibiting this recently discovered 3D electrokinetic instability, and show that the number of parallel helical vortices is a jump-discontinuous function of width, as is the overlimiting current and overlimiting conductance. In addition, we show that overlimiting occurs at lower fields in wider channels, because the associated helical vortices are more readily triggered than the planar vortices associated with narrow channels (effective 2D systems). These unexpected width dependencies arise in realistic electrochemical desalination systems, and have important ramifications for design optimization.

On the Dynamical Regimes of Pattern-Accelerated Electroconvection

Recent research has established that electroconvection can enhance ion transport at polarized surfaces such as membranes and electrodes where it would otherwise be limited by diffusion. The onset of such overlimiting transport can be influenced by the surface topology of the ion selective membranes as well as inhomogeneities in their electrochemical properties. However, there is little knowledge regarding the mechanisms through which these surface variations promote transport. We use high-resolution direct numerical simulations to develop a comprehensive analysis of electroconvective flows generated by geometric patterns of impermeable stripes and investigate their potential to regularize electrokinetic instabilities. Counterintuitively, we find that reducing the permeable area of an ion exchange membrane, with appropriate patterning, increases the overall ion transport rate by up to 80%. In addition, we present analysis of nonpatterned membranes, and find a novel regime of electroconvection where a multivalued current is possible due to the coexistence of multiple convective states.

Thermoelectroconvection near charge-selective surfaces

A new kind of instability caused by Joule heating near charge-selective surfaces (permselective membranes, electrodes, or systems of micro- and nanochannels) is investigated theoretically using a model based on the Rubinstein-Zaltzman approach. A simple relation is derived for the marginal stability curves: Joule heating can either destabilize or stabilize the steady state, depending on the location of the space charge region relative to the gravity vector. For the destabilizing case, the short-wave Rubinstein-Zaltzman instability is replaced by a long-wave thermal instability. The physical mechanism of the thermal instability is found to be very different from Rayleigh-Bénard convection, and is based on a nonuniform distribution of the electrical conductivity in the electrolyte. The study is complemented by numerical investigations both of linear and nonlinear instabilities near a charge-selective surface. There is a good qualitative agreement with the analytics. A possible explanation of the discrepancy between the experimental data and our previous theoretical voltage-current characteristics is highlighted.

Bridging the gap between an isolated nanochannel and a communicating multipore heterogeneous membrane

To bridge the gap between single and isolated pore systems to multipore systems, such as membranes and electrodes, we studied an array of nanochannels with varying interchannel spacing that controlled the degree of channel communication. Instead of treating them as individual channels connected in parallel or an assembly like a homogeneous membrane, this study resolves the pore-pore interaction. We found that increased channel isolation leads to current intensification, whereas at high voltages electroconvective effects control the degree of communication via suppression of the diffusion layer growth.