Effect of viscoelasticity on the flow pattern and the volumetric flow rate in electroosmotic flows through a microchannel

Combined electromechanically driven pulsating flow of nonlinear viscoelastic fluids in narrow confinements

Controlled microscale transport is at the core of many scientific and technological advancements, including medical diagnostics, separation of biomolecules, etc., and often involves complex fluids. One of the challenges in this regard is to actuate flows at small scales in an energy efficient manner, given the strong viscous forces opposing fluid motion. We try to address this issue here by probing a combined time-periodic pressure and electrokinetically driven flow of a viscoelastic fluid obeying the simplified linear Phan-Thien-Tanner model, using numerical as well as asymptotic tools, in view of the fact that oscillatory fields are less energy intensive. We establish that the interplay between oscillatory electrical and mechanical forces can lead to complex temporal mass flow rate variations with short-term bursts and peaks in the flow rate. We further demonstrate that an oscillatory pressure gradient or an electric field, in tandem with another steady actuating force can indeed change the net throughput significantly-a paradigm that is not realized in Newtonian or other simpler polymeric liquids. Our results reveal that the extent of augmentation in the flow rate strongly depends on the frequency of the imposed actuating forces along with their waveforms. We also evaluate the streaming potential resulting from an oscillatory pressure-driven flow and illustrate that akin to the volume throughput, the streaming potential also shows complex temporal variations, while its time average gets augmented in the presence of a time-periodic pressure gradient in a nonlinear viscoelastic medium.

Ion transport and current rectification in a charged conical nanopore filled with viscoelastic fluids

The ionic current rectification (ICR) is a non-linear current-voltage response upon switching the polarity of the potential across nanopore which is similar to the I–V response in the semiconductor diode. The ICR phenomenon finds several potential applications in micro/nano-fluidics (e.g., Bio-sensors and Lab-on-Chip applications). From a biological application viewpoint, most biological fluids (e.g., blood, saliva, mucus, etc.) exhibit non-Newtonian visco-elastic behavior; their rheological properties differ from Newtonian fluids. Therefore, the resultant flow-field should show an additional dependence on the rheological material properties of viscoelastic fluids such as fluid relaxation time \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\lambda )$$\end{document}(λ) and fluid extensibility \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\varepsilon )$$\end{document}(ε). Despite numerous potential applications, the comprehensive investigation of the viscoelastic behavior of the fluid on ionic concentration profile and ICR phenomena has not been attempted. ICR phenomena occur when the length scale and Debye layer thickness approaches to the same order. Therefore, this work extensively investigates the effect of visco-elasticity on the flow and ionic mass transfer along with the ICR phenomena in a single conical nanopore. The Poisson–Nernst–Planck (P–N–P) model coupled with momentum equations have been solved for a wide range of conditions such as, Deborah number, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\le De \le 100$$\end{document}1≤De≤100, Debye length parameter, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\le \kappa R_t \le 50$$\end{document}1≤κRt≤50, fluid extensibility parameter, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.05\le \varepsilon \le 0.25$$\end{document}0.05≤ε≤0.25, applied electric potential, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-40\le V \le 40$$\end{document}-40≤V≤40, and surface charge density \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma = -10$$\end{document}σ=-10 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-50$$\end{document}-50. Limited results for Newtonian fluid (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$De = 0$$\end{document}De=0, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon = 0$$\end{document}ε=0) have also been shown in order to demonstrate the effectiveness of non-Newtonian fluid behaviour over the Newtonian fluid behaviour. Four distinct novel characteristics of electro-osmotic flow (EOF) in a conical nanopore have been investigated here, namely (1) detailed structure of flow field and velocity distribution in viscoelastic fluids (2) influence of Deborah number and fluid extensibility parameter on ionic current rectification (ICR) (3) volumetric flow rate calculation as a function of Deborah number and fluid extensibility parameter (4) effect of viscoelastic parameters on concentration distribution of ions in the nanopore. At high applied voltage, both the extensibility parameter and Deborah number facilitate the ICR phenomena. In addition, the ICR phenomena are observed to be more pronounced at low values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa R_t$$\end{document}κRt than the high values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa R_t$$\end{document}κRt. This effect is due to the overlapping of the electric double layer at low values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa R_t$$\end{document}κRt.

Electroosmotic Flow of Viscoelastic Fluid through a Constriction Microchannel

Electroosmotic flow (EOF) has been widely used in various biochemical microfluidic applications, many of which use viscoelastic non-Newtonian fluid. This study numerically investigates the EOF of viscoelastic fluid through a 10:1 constriction microfluidic channel connecting two reservoirs on either side. The flow is modelled by the Oldroyd-B (OB) model coupled with the Poisson-Boltzmann model. EOF of polyacrylamide (PAA) solution is studied as a function of the PAA concentration and the applied electric field. In contrast to steady EOF of Newtonian fluid, the EOF of PAA solution becomes unstable when the applied electric field (PAA concentration) exceeds a critical value for a fixed PAA concentration (electric field), and vortices form at the upstream of the constriction. EOF velocity of viscoelastic fluid becomes spatially and temporally dependent, and the velocity at the exit of the constriction microchannel is much higher than that at its entrance, which is in qualitative agreement with experimental observation from the literature. Under the same apparent viscosity, the time-averaged velocity of the viscoelastic fluid is lower than that of the Newtonian fluid.

Electroosmotic flow of generalized Burgers’ fluid with Caputo–Fabrizio derivatives through a vertical annulus with heat transfer

Applying the electric field to a fluid confined between capillary surfaces is the most recent technique used for studying the fluid movement. This technique is known as electroosmotic flow (EOF). The problem of the Caputo–Fabrizio time-fractional derivative of the electroosmotic generalized Burgers’ fluid through a vertical annulus with free convection heat transfer has been investigated. The annulus walls kept at constant values of a zeta potential. The dimensionless governing equations have been solved with the help of the Laplace and finite Hankel transforms. The inversion of Laplace for the transformed functions is calculated numerically. The essential features of the electroosmotic flow (EOF) and the related thermal characteristics are clearly illustrated by the variations in the velocity, the flow rate, the temperature and the Nusselt number. Moreover, comparisons with other non-Newtonian fluids have been discussed. It was found that the presence of the electric double layer (EDL) accelerates the fluid near the walls of the annulus, while in the core region, the reverse flow is noticed.

Electroosmotic Flow Behavior of Viscoelastic LPTT Fluid in a Microchannel

In many research works, the fluid medium in electroosmosis is considered to be a Newtonian fluid, while the polymer solutions and biological fluids used in biomedical fields mostly belong to the non-Newtonian category. Based on the finite volume method (FVM), the electroosmotic flow (EOF) of viscoelastic fluids in near-neutral (pH = 7.5) solution considering four ions (K+, Cl-, H+, OH-) is numerically studied, as well as the viscoelastic fluids' flow characteristics in a microchannel described by the Linear Phan-Thien-Tanner (LPTT) constitutive model under different conditions, including the electrical double layer (EDL) thickness, the Weissenberg number (Wi), the viscosity ratio and the polymer extensibility parameters. When the EDL does not overlap, the velocity profiles for both Newtonian and viscoelastic fluids are plug-like and increase sharply near the charged wall. Compared with Newtonian fluid at Wi = 3, the viscoelastic fluid velocity increases by 5 times and 9 times, respectively, under the EDL conditions of kH = 15 and kH = 250, indicating the shear thinning behavior of LPTT fluid. Shear stress obviously depends on the viscosity ratio and different Wi number conditions. The EOF is also enhanced by the increase (decrease) in polymer extensibility parameters (viscosity ratio). When the extensibility parameters are large, the contribution to velocity is gradually weakened.

Electroosmotic Flow of Viscoelastic Fluid in a Nanoslit

The electroosmotic flow (EOF) of viscoelastic fluid in a long nanoslit is numerically studied to investigate the rheological property effect of Linear Phan-Thien-Tanner (LPTT) fluid on the fully developed EOF. The non-linear Poisson-Nernst-Planck equations governing the electric potential and the ionic concentration distribution within the channel are adopted to take into account the effect of the electrical double layer (EDL), including the EDL overlap. When the EDL is not overlapped, the velocity profiles for both Newtonian and viscoelastic fluids are plug-like and increase sharply near the charged wall. The velocity profile resembles that of pressure-driven flow when the EDL is overlapped. Regardless of the EDL thickness, apparent increase of velocity is obtained for viscoelastic fluid of larger Weissenberg number compared to the Newtonian fluid, indicating the shear thinning behavior of the LPTT fluid. The effect of the Weissenberg number on the velocity distribution is less significant as the degree of EDL overlapping increases, due to the overall decrease of the shear rate. The increase (decrease) of polymer extensibility (viscosity ratio) also enhances the EOF of viscoelastic fluid.

Analytical Solution of Electro-Osmotic Peristalsis of Fractional Jeffreys Fluid in a Micro-Channel

The electro-osmotic peristaltic flow of a viscoelastic fluid through a cylindrical micro-channel is studied in this paper. The fractional Jeffreys constitutive model, including the relaxation time and retardation time, is utilized to describe the viscoelasticity of the fluid. Under the assumptions of long wavelength, low Reynolds number, and Debye-Hückel linearization, the analytical solutions of pressure gradient, stream function and axial velocity are explored in terms of Mittag-Leffler function by Laplace transform method. The corresponding solutions of fractional Maxwell fluid and generalized second grade fluid are also obtained as special cases. The numerical analysis of the results are depicted graphically, and the effects of electro-osmotic parameter, external electric field, fractional parameters and viscoelastic parameters on the peristaltic flow are discussed.

Dynamic Electroosmotic Flows of Power-Law Fluids in Rectangular Microchannels

Dynamic characteristics of electroosmosis of a typical non-Newtonian liquid in a rectangular microchannel are investigated by using numerical simulations. The non-Newtonian behavior of liquids is assumed to obey the famous power-law model and then the mathematical model is solved numerically by using the finite element method. The results indicate that the non-Newtonian effect produces some noticeable dynamic responses in electroosmotic flow. Under a direct current (DC) driving electric field, it is found that the fluid responds more inertly to an external electric field and the steady-state velocity profile becomes more plug-like as the flow behavior index decreases. Under an alternating current (AC) driving electric field, the fluid is observed to experience more significant acceleration and the amplitude of oscillating velocity becomes larger as the fluid behavior index decreases. Furthermore, our investigation also shows that electroosmotic flow of power-law fluids under an AC/DC combined driving field is enhanced as compared with that under a pure DC electric field. These dynamic predictions are of practical use for the design of electroosmotically-driven microfluidic devices that analyze and process non-Newtonian fluids such as biofluids and polymeric solutions.

Electrokinetics of non-Newtonian fluids: a review

This work presents a comprehensive review of electrokinetics pertaining to non-Newtonian fluids. The topic covers a broad range of non-Newtonian effects in electrokinetics, including electroosmosis of non-Newtonian fluids, electrophoresis of particles in non-Newtonian fluids, streaming potential effect of non-Newtonian fluids and other related non-Newtonian effects in electrokinetics. Generally, the coupling between non-Newtonian hydrodynamics and electrostatics not only complicates the electrokinetics but also causes the fluid/particle velocity to be nonlinearly dependent on the strength of external electric field and/or the zeta potential. Shear-thinning nature of liquids tends to enhance electrokinetic phenomena, while shear-thickening nature of liquids leads to the reduction of electrokinetic effects. In addition, directions for the future studies are suggested and several theoretical issues in non-Newtonian electrokinetics are highlighted.

Effects of non-Newtonian power law rheology on mass transport of a neutral solute for electro-osmotic flow in a porous microtube

Mass transport of a neutral solute for a power law fluid in a porous microtube under electro-osmotic flow regime is characterized in this study. Combined electro-osmotic and pressure driven flow is conducted herein. An analytical solution of concentration profile within mass transfer boundary layer is derived from the first principle. The solute transport through the porous wall is also coupled with the electro-osmotic flow to predict the solute concentration in the permeate stream. The effects of non-Newtonian rheology and the operating conditions on the permeation rate and permeate solute concentration are analyzed in detail. Both cases of assisting (electro-osmotic and poiseulle flow are in same direction) and opposing flow (the individual flows are in opposite direction) cases are taken care of. Enhancement of Sherwood due to electro-osmotic flow for a non-porous conduit is also quantified. Effects if non-Newtonian rheology on Sherwood number enhancement are observed.

Electrokinetically induced alterations in dynamic response of viscoelastic fluids in narrow confinements

We investigate a dynamical interplay between interfacial electrokinetics and a combined dissipative and elastic behavior of flow through narrow confinements, in analogy with spatiotemporal hydrodynamics of porous media. In particular, we investigate the effects of streaming potential on the pertinent dynamic responses, by choosing a Maxwell fluid model for representing the consequent electro-hydrodynamic characteristics. We transform the pertinent governing equation to the frequency domain, so that a dynamic generalization of Darcy's law in the presence of streaming potential effects can be effectively realized. We show that the frequencies corresponding to local maxima in the dynamic permeability also correspond to local maxima in the induced streaming potential. We also bring out the effects of Stern layer conductivity on the dynamic permeability. Our analytical estimates do reveal that serious overestimations in the commonly portrayed notion of massive amplifications of dynamic permeability at resonating frequencies may be possible, if interactions between spontaneous electrochemical interfacial phenomena and pulsating pressure-gradient-driven viscoelastic transport are trivially ignored.

Coupling electrokinetics and rheology: Electrophoresis in non-Newtonian fluids

We present a theoretical scheme to calculate the electrophoretic motion of charged colloidal particles immersed in complex (non-Newtonian) fluids possessing shear-rate-dependent viscosities. We demonstrate that this non-Newtonian rheology leads to an explicit shape and size dependence of the electrophoretic velocity of a uniformly charged particle in the thin-Debye-layer regime, in contrast to electrophoresis in Newtonian fluids. This dependence is caused by non-Newtonian stresses in the bulk (electroneutral) fluid outside the Debye layer, whose magnitude is naturally characterized in an electrophoretic Deborah number.

Steric-effect induced alterations in streaming potential and energy transfer efficiency of non-newtonian fluids in narrow confinements

In this work, we explore the possibilities of utilizing the combined consequences of interfacial electrokinetics and rheology toward augmenting the energy transfer efficiencies in narrow fluidic confinements. In particular, we consider the exploitation of steric effects (i.e., effect of finite size of the ionic species) in non-Newtonian fluids over small scales, to report dramatic augmentations in the streaming potential, for shear-thickening fluids. We first derive an expression for the streaming potential considering strong electrical double layer interactions in the confined flow passage and the consequences of the finite conductance of the Stern layer, going beyond the Debye-Hückel limit. With a detailed accounting for the excluded volume effects of the ionic species and their interaction with pertinent interfacial phenomena of special type of rheological fluids such as the power law fluids in the above-mentioned formalism, we demonstrate that a confluence of the steric interactions with the non-Newtonian transport characteristics may result in giant augmentations in the energy transfer efficiency for shear-thickening fluids under appropriate conditions.

Abatement of mixing in shear-free elongationally unstable viscoelastic microflows

The addition of minute amounts of chemically inert polyacrylamide polymer to liquids results in large instabilities under steady electro-osmotic pumping through 2 : 1 constrictions, demonstrating that laminar flow conditions can be broken in electro-osmotic flow of viscoelastic material. By excluding shear and imposing symmetry we create a platform where only elongational viscoelastic instabilities, and diffusion, affect mixing. In contrast to earlier studies with significant shear that found up to orders of magnitude increase in mixing we find that inclusion of polymers excites large viscoelastic instabilities yet mixing is reduced relative to polymer-free liquids. The absolute decrease in mixing we find is consistent with the understanding that adding polymer increases viscosity while viscoelastic flows progress towards elastic turbulence, a type of mild (Batchelor) turbulence, and indicates that electro-osmotic pumped devices are an ideal platform for studying viscoelastic instabilities without supplementary factors.

The EOF of polymer solutions

The EOF of polymer solutions is analysed in the framework of continuum fluid mechanics and the standard electrokinetic model. Two key aspects are taken into consideration: the non-Newtonian character of the fluid and the polymer concentration near the interface, which greatly modify the fluid viscosity in the region where electroosmosis takes place. A satisfactory mathematical model is derived for the electroosmotic mobility of solutions that present polymer depletion at the wall. The case of solutions containing polymers that adsorb onto the wall is briefly reviewed, and a preliminary approach is discussed for the limit of strong polymer adsorption. In order to illustrate the theoretical discussions, experimental data obtained from aqueous solutions of carboxymethyl cellulose in fused-silica capillaries are presented. Relevant results are observed, which are appropriately captured by the modelling proposed. The fundamental phenomena discussed in this work are of interest in microfluidics and electrophoresis.