Oscillating electroosmotic flow in channels and capillaries with modulated wall charge distribution

Electrolyte-filled channels with modulated wall charge distribution subjected to an applied DC electric field, form time-independent vortices whose sense of circulation is determined by the field direction [Physical Review Letters $ \mathbf{75}, 755, (1995)$]. In this paper we show that an electrolyte in a channel or cylindrical capillary subjected to an external \emph{alternating} (AC) electric field gives rise to various laminar flow structures, including vortices whose sense of circulation changes with the period of oscillation of the applied AC field. The introduction of a period of oscillation lifts certain degeneracies associated with its time-independent counterpart. Although, in general, the mass flux vanishes, the charge flux is nonzero. The flow is accompanied by a longitudinal (oscillating) advective current that displays hysteresis accompanied by a diverging and negative self-similar conductance that depends on the applied voltage [Nano Letters $\mathbf{10}, 2674, (2010)$]. We show that this behavior can be interpreted with respect to a ``memory retention time’’, that depends on frequency, viscosity and the Debye length and could thus form the impetus for investigating control protocols of signal carriers.

💡 Research Summary

In this paper the authors investigate the electro‑osmotic flow that arises when a 1:1 electrolyte fills a rectangular micro‑channel or a cylindrical capillary whose walls carry a spatially modulated surface charge, and an alternating (AC) electric field is applied. Building on earlier work that showed static vortices for a DC field and uniform wall charge, the authors demonstrate that an AC field introduces a time‑dependent body force that generates oscillating vortex patterns whose sense of rotation follows the instantaneous sign of the applied field.

The theoretical framework starts from the Navier‑Stokes equations with an electric body force  F = ρ_e(E − ∇ϕ), where the charge density ρ_e = −εκ²ϕ follows from the linearized Poisson–Boltzmann (Debye‑Hückel) approximation. The external field is taken as E = E_∥ e^{−iωt} x̂ (or along the axis for the capillary) and the wall charge is prescribed as σ_±(x)=±σ₀ cos(qx) (or σ = σ₀ cos(qz) for the capillary). By introducing a streamfunction ψ(x,z,t)=ψ(z) e^{−iωt} cos(qx) the governing equations reduce to a fourth‑order ordinary differential equation in z.

A key novelty is the appearance of a complex wave number k = (1 + i)/δ, where the penetration depth δ = √(2ν/ω) depends on the fluid kinematic viscosity ν and the driving frequency ω. Together with the charge‑modulation wave number q, this yields a second complex wave number K = √(k² − q²). These two distinct length scales (δ and the Debye length κ⁻¹) lift the degeneracy that exists in the DC case, where only the purely imaginary modes ±iq appear. Consequently, the velocity field contains both oscillatory (∝ cos(Kz), sin(Kz)) and exponential (∝ cosh(qz), sinh(qz)) contributions, and the vortex sense reverses each half‑cycle of the AC field.

For the rectangular channel the authors solve the boundary‑value problem analytically, obtaining explicit expressions for the streamfunction (Eq. 16) and the velocity components u and w (Eqs. 18‑19). The solutions reduce smoothly to the known DC results as ω → 0. Two charge configurations are examined: symmetric (σ_+=σ_−) and antisymmetric (σ_+=−σ_−). The former yields a single vortex row, while the latter produces a pair of counter‑rotating vortices whose relative phase can be tuned continuously by shifting the wall‑charge phase.

The cylindrical geometry is treated analogously. With a wall charge σ = σ₀ cos(qz) and an axial AC field, the flow consists of a stack of toroidal vortices (Fig. 3). The analytical solution again displays the combined exponential‑oscillatory structure, confirming that the phenomenon is robust to geometry.

Beyond the flow field, the paper analyses the ionic current. The axial charge flux J_z = ∫ρ_e v_z dA is non‑zero even though the net mass flux averages to zero over a period. The current–voltage (I‑V) characteristic exhibits a hysteresis loop whose area A_loop ∝ ∮V dI represents the energy dissipated per cycle. By varying the driving frequency, the loop area reaches a maximum at a particular ω*; the corresponding time τ_mem ≈ 1/ω* is interpreted as a “memory retention time” that depends on ν, κ, and q. Notably, for certain parameter regimes the differential conductance dI/dV becomes negative and diverges, reminiscent of memristive behavior. The authors discuss this negative, self‑similar conductance in the context of recent literature on “memristive electro‑osmotic” systems.

The discussion emphasizes the practical implications: the ability to switch vortex direction and to modulate the ionic current by simply adjusting the AC frequency or amplitude provides a route to reconfigurable micro‑fluidic devices. Potential applications include rapid mixing, targeted solute transport, on‑chip pumping without moving parts, and even fluidic logic elements that exploit the hysteretic I‑V response as a memory element. The work also suggests that engineered wall‑charge patterns could be used to tailor the effective penetration depth and thus the frequency response, opening a design space for fluidic devices that combine electro‑kinetic actuation with information‑processing capabilities.

In summary, the paper presents a comprehensive analytical treatment of AC‑driven electro‑osmotic flow in channels with spatially modulated wall charge, uncovers new dynamical features such as vortex reversal, hysteretic ionic currents, and negative conductance, and points toward novel micro‑fluidic functionalities based on these memory‑like effects.