Streaming potential and electroosmotic flow in heterogeneous circular microchannels with nonuniform zeta potentials: requirements of flow rate and current continuities

Real surfaces are typically heterogeneous, and microchannels with heterogeneous surfaces are commonly found due to fabrication defects, material impurities, and chemical adsorption from solution. Such surface heterogeneity causes a nonuniform surface potential along the microchannel. Other than surface heterogeneity, one could also pattern the various surface potentials along the microchannels. To understand how such variations affect electrokinetic flow, we proposed a model to describe its behavior in circular microchannels with nonuniform surface potentials. Unlike other models, we considered the continuities of flow rate and electric current simultaneously. These requirements cause a nonuniform electric field distribution and pressure gradient along the channel for both pressure-driven flow (streaming potential) and electric-field-driven flow (electroosmosis). The induced nonuniform pressure and electric field influence the electrokinetic flow in terms of the velocity profile, the flow rate, and the streaming potential.     

3-D transient electrophoretic motion of a spherical particle in a T-shaped rectangular microchannel

This paper considers the electrophoretic motion of a spherical particle in an aqueous electrolyte solution in a T-shaped rectangular microchannel, where the size of the channel is close to that of the particle. This is a complicated transient process where the electric field, the flow field, and the particle motion are coupled together. A theoretical model was developed to investigate the influences of the applied electric potentials, the zeta potentials of the channel and the particle, and the size of the particle on the particle motion. A direct numerical simulation method using the finite element method is employed. This method employs a generalized Galerkin finite element formulation that incorporates both equations of the fluid flow and equations of the particle motion into a single variational equation where the hydrodynamic interactions are eliminated. The ALE method is used to track the surface of the particle at each time step. The numerical results show that the electric field in the T-shaped microchannel is influenced by the presence of the particle, and that the particle motion is influenced by the applied electric potentials and the zeta potentials of the channel and the particle. The path of the particle motion is dominated by the local electric field and the ratio of the zeta potential of the channel to that of the particle. The particle's velocity is also dependent on its size in a small channel.     

The Electrophoretic Mobility and Electric Conductivity of a Concentrated Suspension of Colloidal Spheres with Arbitrary Double-Layer Thickness

The electrophoresis in a monodisperse suspension of dielectric spheres with an arbitrary thickness of the electric double layers is analytically studied. The effects of particle interactions are taken into account by employing a unit cell model, and the overlap of the double layers of adjacent particles is allowed. The electrokinetic equations, which govern the ionic concentration distributions, the electric potential profile, and the fluid flow field in the electrolyte solution surrounding the charged sphere in a unit cell, are linearized assuming that the system is only slightly distorted from equilibrium. Using a perturbation method, these linearized equations are solved with the surface charge density (or zeta potential) of the particle as the small perturbation parameter. Analytical expressions for the electrophoretic mobility of the colloidal sphere in closed form correct to O(zeta) are obtained. Based on the solution of the electrokinetic equations in a cell, a closed-form formula for the electric conductivity of the suspension up to O(zeta(2)) is derived from the average electric current density. Comparisons of the results of the cell model with different conditions at the outer boundary of the cell are made for both the electrophoretic mobility and the electric conductivity. Copyright 2001 Academic Press.     

Temperature gradient focusing of bio-analyte in a microfluidic channel dealing with non-Newtonian electrolyte considering temperature-dependent zeta potential

Temperature gradient focusing (TGF) relies on establishing a precise balance between the electrophoretic motility of a target analyte and the advective flow of the background electrolyte (BGE) to locally concentrate the analyte in a microfluidic configuration. This paper presents a finite-element-based numerical analysis where the coupled electric field and the transport equations are solved to describe the effects of the shear-dependent apparent viscosity of a non-Newtonian BGE on the localized concentration buildup of a charged bio-sample inside a microchannel by TGF via Joule heating. Effects of the temperature-dependent nature of the wall zeta potential and the flow behavior index (n) of BGE on the flow, thermal, and species concentration profiles inside the microchannel have been investigated. Study using a fluorescein-Na analyte sample shows that the maximum normalized analyte concentration (C_max/C₀) reduces as the zeta potential increases linearly with temperature. The maximum concentration enhancement is achieved when the BGE displays the Newtonian rheology. For example, C_max/C₀ increases 134- to 280-fold when n is increased from 0.8 to 1 (pseudoplastic regime) and again reduces to 190-fold when n increases further from 1 to 1.2 (dilatant regime).     

Sample transport and electrokinetic injection in a microchip device for chemical cytometry

Sample transport and electrokinetic injection bias are well characterized in capillary electrophoresis and simple microchips, but a thorough understanding of sample transport on devices combining electroosmosis, electrophoresis, and pressure-driven flow is lacking. In this work, we evaluate the effects of electric fields from 0 to 300 V/cm, electrophoretic mobilities from 10(-4) to 10(-6) cm(2)/Vs, and pressure-driven fluid velocities from 50 to 250 μm/s on sample injection in a microfluidic chemical cytometry device. By studying a continuous sample stream, we find that increasing electric field strength and electrophoretic mobility result in improved injection and that COMSOL simulations accurately predict sample transport. The effects of pressure-driven fluid velocity on injection are complex, and relative concentration values lie on a surface defined by pressure-driven flow rates. For high-mobility analytes, this surface is flat, and sample injection is robust despite fluctuations in flow rate. For lower mobility analytes, the surface becomes steeper, and injection depends strongly on pressure-driven flow. These results indicate generally that device design must account for analyte characteristics and specifically that this device is suited to high-mobility analytes. We demonstrate that for a suitable pair of peptides fluctuations in injection volume are correlated; electrokinetic injection bias is minimized; and electrophoretic separation is achieved.     

Characterization of electrokinetic mobility of microparticles in order to improve dielectrophoretic concentration

Insulator-based dielectrophoresis (iDEP), an efficient technique with great potential for miniaturization, has been successfully applied for the manipulation of a wide variety of bioparticles. When iDEP is applied employing direct current (DC) electric fields, other electrokinetic transport mechanisms are present: electrophoresis and electroosmotic flow. In order to concentrate particles, iDEP has to overcome electrokinetics. This study presents the characterization of electrokinetic flow under the operating conditions employed with iDEP; in order to identify the optimal conditions for particle concentration employing DC-iDEP, microparticle image velocimetry (μPIV) was employed to measure the velocity of 1-μm-diameter inert polystyrene particles suspended inside a microchannel made from glass. Experiments were carried out by varying the properties of the suspending medium (conductivity from 25 to 100 μS/cm and pH from 6 to 9) and the strength of the applied electric field (50–300 V/cm); the velocities values obtained ranged from 100 to 700 μm/s. These showed that higher conductivity and lower pH values for the suspending medium produced the lowest electrokinetic flow, improving iDEP concentration of particles, which decreases voltage requirements. These ideal conditions for iDEP trapping (pH = 6 and σ _m = 100 μS/cm) were tested experimentally and with the aid of mathematical modeling. The μPIV measurements allowed obtaining values for the electrokinetic mobilities of the particles and the zeta potential of the glass surface; these values were used with a mathematical model built with COMSOL Multiphysics software in order to predict the dielectrophoretic and electrokinetic forces exerted on the particles; the modeling results confirmed the μPIV findings. Experiments with iDEP were carried out employing the same microparticles and a glass microchannel that contained an array of cylindrical insulating structures. By applying DC electric fields across the insulating structures array, it was seen that the dielectrophoretic trapping was improved when the electrokinetic force was the lowest; as predicted by μPIV measurements and the mathematical model. The results of this study provide guidelines for the selection of optimal operating conditions for improving insulator-based dielectrophoretic separations and have the potential to be extended to bioparticle applications. Figure Comparison of experimental measurements and mathematical modeling of electrokinetic and dielectrophoretic effects on microparticles

Theoretical analysis of electrokinetic flow and heat transfer in a microchannel under asymmetric boundary conditions

The main theme of the present work is to investigate the electrokinetic effects on liquid flow and heat transfer in a flat microchannel of two parallel plates under asymmetric boundary conditions including wall-sliding motion, unequal zeta potentials, and unequal heat fluxes on two walls. Based on the Debye-Huckel approximation, an electrical potential solution to the linearized Poisson-Boltzmann equation is obtained and employed in the analysis. The analytic solutions of the electrical potential, velocity distributions, streaming potential, friction coefficient, temperature distribution, and heat transfer rate are obtained, and thereby the effects of electrokinetic separation distance (K), zeta-potential level (zeta;(1)), ratio of two zeta potentials (r(zeta) identical with zeta;(2)/zeta;(1)), wall-sliding velocity (u(w)), and heat flux ratio (r(q) identical with q"(2)/q"(1)) are investigated. The present results reveal the effects of wall-sliding and zeta-potential ratio on the hydrodynamic nature of microchannel flow, and they are used to provide physical interpretations for the resultant electrokinetic effects and the underlying electro-hydrodynamic interaction mechanisms. In the final part the results of potential and velocity fields are applied in solving the energy equation. The temperature distributions and heat transfer characteristics under the asymmetrical kinematic, electric, and thermal boundary conditions considered presently are dealt with.     

Continuous microfluidic DNA and protein trapping and concentration by balancing transverse electrokinetic forces

Sample pre-concentration can be a critical element to improve sensitivity of integrated microchip assays. In this work a converging Y-inlet microfluidic channel with integrated coplanar electrodes was used to investigate transverse DNA and protein migration under uniform direct current (DC) electric fields to assess the ability to concentrate a sample prior to other enzymatic modifications or capillary electrophoretic separations. Employing a pressure-driven flow to perfuse the microchannel, negatively charged samples diluted in low and high ionic strength buffers were co-infused with a receiving buffer of the same ionic strength into a main daughter channel. Experimental results demonstrated that, depending of the buffer selection, different DNA migration and accumulation dynamics were seen. Charged analytes could traverse the channel width and accumulate at the positive bias electrode in a low electroosmotic mobility, high electrophoretic mobility, high ionic strength buffer or migrated towards an equilibrium position within the channel in a high electroosmotic mobility, high electrophoretic mobility, low ionic strength buffer. The various migration behaviours are the result of a balance between the electrophoretic force and a drag force induced by a recirculating electroosmotic flow generated across the channel width due to the bounding walls. Under continuous flow conditions, DNA samples were concentrated several-fold by balancing these transverse electrokinetic forces. The electrokinetic trapping technique presented here is a simple technique which could be expanded to concentrate or separate other analytes as a preconditioning step for downstream processes.     

Nonlinear electrophoresis of dielectric particles in Newtonian fluids

In classical electrokinetics, the electrophoretic velocity of a dielectric particle is a linear function of the applied electric field. Theoretical studies have predicted the onset of nonlinear electrophoresis at high electric fields because of the nonuniform surface conduction over the curved particle. However, experimental studies have been left behind and are insufficient for a fundamental understanding of the parametric effects on nonlinear electrophoresis. We present in this work a systematic experimental study of the effects of buffer concentration, particle size, and particle zeta potential on the electrophoretic velocity of polystyrene particles in a straight rectangular microchannel for electric fields of up to 3 kV/cm. The measured nonlinear electrophoretic particle velocity is found to exhibit a 2(±0.5)-order dependence on the applied electric field, which appears to be within the theoretically predicted 3- and 3/2-order dependences for low and high electric fields, respectively. Moreover, the obtained nonlinear electrophoretic particle mobility increases with decreasing buffer concentration (for the same particle) and particle size (for particles with similar zeta potentials) or increasing particle zeta potential (for particles with similar sizes). These observations are all consistent with the theoretical predictions for high electric fields.     

Revisit of wall‐induced lateral migration in particle electrophoresis through a straight rectangular microchannel: Effects of particle zeta potential

Previous studies have reported a lateral migration in particle electrophoresis through a straight rectangular microchannel. This phenomenon arises from the inherent wall‐induced electrical lift that can be exploited to focus and separate particles for microfluidic applications. Such a dielectrophoretic‐like force has been recently found to vary with the buffer concentration. We demonstrate in this work that the particle zeta potential also has a significant effect on the wall‐induced electrical lift. We perform an experimental study of the lateral migration of equal‐sized polystyrene particles with varying surface charges under identical electrokinetic flow conditions. Surprisingly, an enhanced focusing is observed for particles with a faster electrokinetic motion, which indicates a substantially larger electrical lift for particles with a smaller zeta potential. We speculate this phenomenon may be correlated with the particle surface conduction that is a strong function of particle and fluid properties.     

Near-wall velocities of particles suspended in shear flow and a streamwise electric field

Particles with a diameter of ∼0.5 µm in a dilute (volume fractions φ_∞ < 4 × 10⁻³) suspension assemble into highly elongated structures called “bands” under certain conditions in combined Poiseuille and electroosmotic flows in opposite directions through microchannels at particle-based Reynolds numbers Re_p < < 1. The particles are first concentrated near, then form “bands” within ∼6 µm of, the channel wall. The experiments described here examine the near-wall dynamics of individual “tracer” particles during the initial concentration, or accumulation, of particles, and the steady-state stage when the particles have formed relatively stable bands at different near-wall shear rates and electric field magnitudes. Surprisingly, the near-wall upstream particle velocities are found to be consistently greater in magnitude than the expected values based on the particles being convected by the superposition of both flows and subject to electrophoresis, which is in the same direction as the Poiseuille flow. However, the particle velocities scale linearly with the change in electric field magnitude, suggesting that the particle dynamics are dominated by linear electrokinetic phenomena. If this discrepancy with theory is only due to changes in particle electrophoresis, electrophoresis is significantly reduced to values as small as 20%–50% of the Smoluchowski relation, or well below previous model predictions, even for high particle potentials.     

Transport of charged samples in fluidic channels with large zeta potentials

In this article, we present an analysis on the transport of charged samples through micro- and nanofluidic channels with large zeta potentials (|zeta| > (kBT)/e). Using the Method of Moments formulation, the diffusion-convection equation has been solved to evaluate the mean velocity and the dispersion of analyte bands in a parallel-plate device under electrokinetically- and pressure-driven flow conditions. The effect of electromigration induced by the lateral electric field within the Debye layer has been quantified in our work using a Peclet number (Pe t) based on the characteristic electrophoretic velocity of the solute molecules in the transverse direction. It has been shown that while the effects of transverse electromigration on analyte transport only depends on the product Pe t zeta* for |zeta*| = (ezeta)/kBT < 1, both these parameters independently affect the flow of charged species in large zeta potential systems. For a given value of Pe t zeta*, the mean velocity and the slug dispersivity can vary by as much as an order of magnitude in going from a small zeta potential system (|zeta*| < 1) to a channel with |zeta*| = 4.     

Ultrafast electrokinetics

The influence of a high electric field applied on both fluid flow and particle velocities is quantified at large Peclet numbers. The experiments involved simultaneous particle image velocimetry and flow rate measurements. These are conducted in polydimethylsiloxane channels with spherical nonconducting polystyrene particles and DI water as the background flow. The high electric field tests produced up to three orders of magnitude higher electrokinetic velocities than any previous reports. The maximum electroosmotic velocity and electrophoretic velocity measured were 3.55 and 2.3 m/s. Electrophoretic velocities are measured over the range of 100 V/cm < E < 250 000 V/cm. The results are separated according to the different nonlinear theoretical models, including low and high Peclet numbers, and weak and strong concentration polarization. They show good agreement with the models. Such fast velocities could be used for flow separation, mixing, transport, control, and manipulation of suspended particles as well as microthrust generation among other applications.     

Dielectrophoretic choking phenomenon of a deformable particle in a converging‐diverging microchannel

The translational motion of small particles in an electrokinetic fluid flow through a constriction can be enhanced by an increase of the applied electric potential. Beyond a critical potential, however, the negative dielectrophoresis (DEP) can overpower other forces to prevent particles that are even smaller than the constriction from passing through the constriction. This DEP choking phenomenon was studied previously for rigid particles. Here, the DEP choking phenomenon is revisited for deformable particles, which are ubiquitous in many biomedical applications. Particle deformability is measured by the particle shear modulus, and the choking conditions are reported through a parametric study that includes the channel geometry, external electric potential, and particle zeta potential. The study was carried out using a numerical model based on an arbitrary Lagrangian‐Eulerican (ALE) finite‐element method.     

Dielectrophoretic focusing of particles in a microchannel constriction using DC‐biased AC flectric fields

This paper presents a fundamental study of particle electrokinetic focusing in a single microchannel constriction. Through both experiments and simulations, we demonstrate that such dielectrophoresis‐induced particle focusing can be implemented in a much smaller magnitude of DC‐biased AC electric fields (10 kV/m in total) as compared to pure DC electric fields (up to 100 kV/m). This is attributed to the increase in the ratio of cross‐stream particle dielectrophoretic velocity to streamwise electrokinetic velocity as only the DC field component contributes to the latter. The effects of the 1 kHz frequency AC to DC electric field ratio on particle trajectories and velocity variations through the microchannel constriction are also examined, which are found to agree with the simulation results.     

A velocity program using the Kanade–Lucas–Tomasi feature-tracking algorithm with demonstration for pressure and electroosmosis conditions

Investigating microfluidic flow profiles is of interest in the microfluidics field for the determination of various characteristics of a lab-on-a-chip system. Microparticle tracking velocimetry uses computational methods upon recording video footage of microfluidic flow to ultimately visualize motion within a microfluidic system across all frames of a video. Current methods are computationally expensive or require extensive instrumentation. A computational method suited to microparticle tracking applications is the robust Kanade–Lucas–Tomasi (KLT) feature-tracking algorithm. This work explores a microparticle tracking velocimetry program using the KLT feature-tracking algorithm. The developed program is demonstrated using pressure-driven and EOF and compared with the respective mathematical fluid flow models. An electrostatics analysis of EOF conditions is performed in the development of the mathematical using a Poisson's Equation solver. This analysis is used to quantify the zeta potential of the electroosmotic system. Overall, the KLT feature-tracking algorithm presented in this work proved to be highly reliable and computationally efficient for investigations of pressure-driven and EOF in a microfluidic system.     

Electrophoresis of concentrically and eccentrically positioned cylindrical particles in a long tube

We study analytically and numerically the electrophoretic motion of cylindrical particles translating slowly in long tubes filled with an electrolyte solution and subjected to axial electric fields. Both thin and thick double layers are considered. Of particular interest is the case when the particle's and tube's radii are of the same order of magnitude. The model accounts for the flow induced by the particle's motion (the particle acts as a leaky piston) and the electroosmotic flow in the tube. The electrophoretic velocity of the particle and the forces and torques acting on it are determined as functions of the particle's radius, length, and position, the particle's and tube's zeta potentials, the tube's length, and the externally imposed pressures. When the particle is positioned off center, it rotates and its trajectory traces an oscillatory path.     

Reduction in microparticle adsorption using a lateral interconnection method in a PDMS‐based microfluidic device

Microparticle adsorption on microchannel walls occurs frequently due to nonspecific interactions, decreasing operational performance in pressure‐driven microfluidic systems. However, it is essential for delicate manipulation of microparticles or cells to maintain smooth fluid traffic. Here, we report a novel microparticle injection technique, which prevents particle loss, assisted by sample injection along the direction of fluid flow. Sample fluids, including microparticles, mammalian (U937), and green algae (Chlorella vulgaris) cells, were injected directly via a through hole drilled in the lateral direction, resulting in a significant reduction in microparticle attachment. For digital microfluidic application, the proposed regime achieved a twofold enhancement of single‐cell encapsulation compared to the conventional encapsulation rate, based on a Poisson distribution, by reducing the number of empty droplets. This novel interconnection method can be straightforwardly integrated as a microparticle or cell injection component in integrated microfluidic systems.