Electroosmotic Flow in Heterogeneous Microchannels

Electroosmosis in capillary zone electrophoresis with non-uniform zeta potential

Chapter 4 - Single-Phase Electrokinetic Flow in Microchannels

In this paper, a hybrid lattice-Boltzmann and finite-difference (LB-FD) model is applied to simulate the effects of three-dimensional surface roughness and electrokinetic heterogeneity on electroosmotic flow (EOF) in a microchannel. The lattice-Boltzmann (LB) method has been employed to obtain the flow field and a finite-difference (FD) method is used to solve the Poisson-Boltzmann (PB) equation for the electrostatic potential distribution. Numerical simulation of flow through a square cross-section microchannel with designed roughness is conducted and the results are critically analysed. The effects of surface heterogeneity on the electroosmotic transport are investigated for different roughness height, width, roughness interval spacing, and roughness surface potential. Numerical simulations reveal that the presence of surface roughness changes the nature of electroosmotic transport through the microchannel. It is found that the electroosmotic velocity decreases with the increase in roughness height and the velocity profile becomes asymmetric. For the same height of the roughness elements, the EOF velocity rises with the increase in roughness width. For the heterogeneously charged rough channel, the velocity profile shows a distinct deviation from the conventional plug-like flow pattern. The simulation results also indicate locally induced flow vortices which can be utilized to enhance the flow and mixing within the microchannel. The present study has important implications towards electrokinetic flow control in the microchannel, and can provide an efficient way to design a microfluidic system of practical interest.

The present study is concerned with unsteady electroosmotic flow (EOF) in a microchannel with the electric charge distribution described by the Poisson-Boltzmann (PB) equation. The nonlinear PB equation is solved by a systematic perturbation with respect to the parameter λ which measures the strength of the wall zeta potential relative to the thermal potential. In the small λ limits (λ<<1), we recover the linearized PB equation - the Debye-Hückel approximation. The solutions obtained by using only three terms in the perturbation series are shown to be accurate with errors <1% for λ up to 2. The accurate solution to the PB equation is then used to solve the electrokinetic fluid transport equation for two types of unsteady flow: transient flow driven by a suddenly applied voltage and oscillatory flow driven by a time-harmonic voltage. The solution for the transient flow has important implications on EOF as an effective means for transporting electrolytes in microchannels with various electrokinetic widths. On the other hand, the solution for the oscillatory flow is shown to have important physical implications on EOF in mixing electrolytes in terms of the amplitude and phase of the resulting time-harmonic EOF rate, which depends on the applied frequency and the electrokinetic width of the microchannel as well as on the parameter λ.

Electroosmotic flow in microchannels with arbitrary geometry and arbitrary distribution of wall charge

Toward tuning flow charecteristics in microchannel by nanotechnology and electrokinetic: Numerical simulation of heterogenous electroosmotic flow

This study analyzes the impact of slip-dependent zeta potential on the heat transfer characteristics of nanofluids in cylindrical microchannels with consideration of thermal radiation effects. An analytical model is developed, accounting for the coupling between surface potential and interfacial slip. The linearized Poisson-Boltzmann equation, along with the momentum and energy conservation equations, is solved analytically to obtain the electrical potential field, velocity field, temperature distribution, and Nusselt number for both slip-dependent (SD) and slip-independent (SI) zeta potentials. Subsequently, the effects of key parameters, including electric double-layer (EDL) thickness, slip length, nanoparticle volume fraction, thermal radiation parameters, and Brinkman number, on the velocity field, temperature field, and Nusselt number are discussed. The results show that the velocity is consistently higher for the SD zeta potential compared to the SI zeta potential. Meanwhile, the temperature for the SD case is higher than that for the SI case at lower Brinkman numbers, particularly for a thinner EDL. However, an inverse trend is observed at higher Brinkman numbers. Similar trends are observed for the Nusselt number under both SD and SI zeta potential conditions at different Brinkman numbers. Furthermore, for a thinner EDL, the differences in flow velocity, temperature, and Nusselt number between the SD and SI conditions are more pronounced.

Purpose The characteristics of fluid motions in micro-channel are strong fluid-wall surface interactions, high surface to volume ratio, extremely low Reynolds number laminar flow, surface roughness and wall surface or zeta potential. Due to zeta potential, an electrical double layer (EDL) is formed in the vicinity of the wall surface, namely, the stern layer (layer of immobile ions) and diffuse layer (layer of mobile ions). Hence, its competent designs demand more efficient micro-scale mixing mechanisms. This paper aims to therefore carry out numerical investigations of electro osmotic flow and mixing in a constricted microchannel by modifying the existing immersed boundary method. Design/methodology/approach The numerical solution of electro-osmotic flow is obtained by linking Navier–Stokes equation with Poisson and Nernst–Planck equation for electric field and transportation of ion, respectively. Fluids with different concentrations enter the microchannel and its mixing along its way is simulated by solving the governing equation specified for the concentration field. Both the electro-osmotic effects and channel constriction constitute a hybrid mixing technique, a combination of passive and active methods. In microchannels, the chief factors affecting the mixing efficiency were studied efficiently from results obtained numerically. Findings The results indicate that the mixing efficiency is influenced with a change in zeta potential (ζ), number of triangular obstacles, EDL thickness (λ). Mixing efficiency decreases with an increment in external electric field strength (Ex), Peclet number (Pe) and Reynolds number (Re). Mixing efficiency is increased from 28.2 to 50.2% with an increase in the number of triangular obstacles from 1 to 5. As the value of Re and Pe is decreased, the overall percentage increase in the mixing efficiency is 56.4% for the case of a mixing micro-channel constricted with five triangular obstacles. It is also vivid that as the EDL overlaps in the micro-channel, the mixing efficiency is 52.7% for the given zeta potential, Re and Pe values. The findings of this study may be useful in biomedical, biotechnological, drug delivery applications, cooling of microchips and deoxyribonucleic acid hybridization. Originality/value The process of mixing in microchannels is widely studied due to its application in various microfluidic devices like micro electromechanical systems and lab-on-a-chip devices. Hence, its competent designs demand more efficient micro-scale mixing mechanisms. The present study carries out numerical investigations by modifying the existing immersed boundary method, on pressure-driven electro osmotic flow and mixing in a constricted microchannel using the varied number of triangular obstacles by using a modified immersed boundary method. In microchannels, the theory of EDL combined with pressure-driven flow elucidates the electro-osmotic flow.

Non-isothermal effects in the slippage condition and absolute viscosity for an electroosmotic flow

In this study, we investigate the electro-osmotic flow pattern of two immiscible fluids through a microchannel with nonuniform wall potential, considering the impact of the surface tension along the liquid-liquid interface. The main objective of this work is to identify and predict the effects of the heterogeneous wall characteristics on the flow field and the deformation of the liquid-liquid interface. Both the Maxwell stress and the hydrodynamic stress are considered along the liquid-liquid interface with no charge accumulation, to predict the flow dynamics and interface deformation along the core due to the external electric field. For the flow analysis, the height of the channel is assumed to be much larger than the thickness of the electric double layers; hence Smoluchowski slip velocity conditions are applied along the channel walls to predict the electro-osmotic flow variation. Additionally, the transient deformation of the interface is considered when the flow field is considered to be quasisteady, predicting that the effects of inertia are minimal. The governing equations for the potential field and the flow field are solved analytically, using domain perturbation analysis for small deformations of the interface and a small amplitude of the wavy wall. Our results show that the heterogeneous wall topology generates swirling flow dynamics within the microchannel, and depending on the asymmetry in the zeta potential, the disturbances along the interface are enhanced or suppressed with the wavelength of the wavy microchannel.

Modulation of electroosmotic flows in electron-conducting microchannels by coupled quasi-reversible faradaic and adsorption-mediated depolarization

Influence of slip wall effect on a non-isothermal electro-osmotic flow of a viscoelastic fluid

In the present work, we realize a study about the influence of viscoelectric effect on the electroosmotic flow of Newtonian fluids in nanochannels formed by two parallel flat plates. In the problem, the channel walls have heterogeneous zeta potentials which follow a sinusoidal distribution; moreover, viscoelectric effects appear into the electric double layers when high zeta potentials are considered at the channel walls, modifying the fluid viscosity and the fluid velocity. To find the solution of flow field, the modified Poisson-Boltzmann, mass conservation and momentum governing equations, are solved numerically. In the results, combined effects from the zeta potential heterogeneities and viscosity changes yields different kind of flow recirculations controlled by the dephasing angle, amplitude and number of waves of the heterogeneities at the walls. The viscoelectric effect produces a decrease in the magnitude of velocity profiles and volumetric flow rate when the high zeta potentials are magnified. Additionally, the heterogeneous zeta potentials at the walls generate an induced pressure on the flow. This investigation extend the knowledge of electroosmotic flows under field effects for future mixing applications.

In this work the purely electro-osmotic flow of a viscoelastic liquid, which obeys the simplified Phan-Thien–Tanner (sPTT) constitutive equation, is solved numerically and asymptotically by using the lubrication approximation. The analysis includes Joule heating effects caused by an imposed electric field, where the viscosity function, relaxation time and electrical conductivity of the liquid are assumed to be temperature-dependent. Owing to Joule heating effects, temperature gradients in the liquid make the fluid properties change within the microchannel, altering the electric potential and flow fields. A consequence of the above is the appearance of an induced pressure gradient along the microchannel, which in turn modifies the normal plug-like electro-osmotic velocity profiles. In addition, it is pointed out that, depending on the fluid rheology and the used values of the dimensionless parameters, the velocity, temperature and pressure profiles in the fluid are substantially modified. Also, the finite thermal conductivity of the microchannel wall was considered in the analysis. The dimensionless temperature profiles in the fluid and the microchannel wall are obtained as function of the dimensionless parameters involved in the analysis, and the interactions between the coupled momentum, thermal energy and potential electric equations are examined in detail. A comparison between the numerical predictions and the asymptotic solutions was made, and reasonable agreement was found.

In this work, the effects of polarization of confining charged planar dielectric surfaces on induced electroosmotic flow are investigated. To this end, we use dissipative particle dynamics to model solvent and ionic particles together with a modified Ewald sum method to model electrostatic interactions and surfaces polarization. A relevant difference between counterions number density profiles, velocity profiles, and volumetric flow rates obtained with and without surface polarization for moderate and high electrostatic coupling parameters is observed. For low coupling parameters, the effect is negligible. An increase of almost 500% in volumetric flow rate for moderate/high electrostatic coupling and surface separation is found when polarizable surfaces are considered. The most important result is that the increase in electrostatic coupling substantially increases the electroosmotic flow in all studied range of separations when the dielectric constant of electrolytes is much higher than the dielectric constant of confining walls. For the higher separation simulated, an increase of around 340% in volumetric flow rate when the electrostatic coupling is increased by a factor of two orders of magnitude is obtained.