- Research
- Open access
- Published:
Closed-form solutions of transient electro-osmotic flow driven by AC electric field in a microannulus
Boundary Value Problems volume 2014, Article number: 215 (2014)
Abstract
The time-periodic electro-osmotic flow of Newtonian fluids through a microannulus is studied in the Debye-Hückel approximation. Analytical series solutions for velocity and flow rate are presented with the help of an integral transform. The expression for the distribution of the velocity profile consists of a time-dependent oscillating part and a time-dependent generating or transient one, and the normalized velocity function is independent of the Reynolds number, which is very different from previous results. Then the effects of the electrokinetic width K, the wall zeta potential ratio β, and the frequency of applied external electric field ω on the distribution of the velocity profiles and flow rates are discussed numerically and theoretically. Some new physical and chemical phenomena are found theoretically. We point out that the electro-osmotic flow driven by an alternating electric field is not periodic in time, but quasi-periodic. There is a phase shift between voltage and flow, which is only dependent on the frequency of the external electric field.
1 Introduction
When an electric field is applied to the fluids in a channel, the walls of which are charged, the migration of the ions present in excess in the double layer induces the motion of the bulk solution due to viscous drag. This phenomenon provides an attractive means of manipulating liquids in microdevices, and it has been widely used in different microdevices and for various applications, such as microfractionation [1], [2], electrophoresis [3], and microspray generation systems [4].
Time-periodic electro-osmotic flow is also known as AC electro-osmosis, and it is driven by an alternating electric field. It is very important for biotechnology and separation science. Recently, various studies analyzed the time-periodic electro-osmotic flow theory and modeling in different geometry. Dutta and Beskok [5] were among the early researchers who analytically investigated the time-periodic electro-osmotic flow between two parallel plates, illustrating interesting similarities or dissimilarities with the Stokes second problem. Based on the method proposed by Dutta and Beskok, many researchers studied time-periodic electro-osmotic flows through microchannels, and some new results are given. General solutions were developed by Xuan and Li [6] for direct current and alternating current electro-osmotic flows in microfluidic channels with arbitrary cross-sectional geometry and arbitrary distribution of wall charge. Jian et al. and his colleagues investigated the flow behavior of time-periodic electro-osmosis in a cylindrical microannulus [7], [8].
Unfortunately, due to the incorrect critical assumption of the form of velocity distribution, the results given in these researches are not correct, and some very important physical phenomena have not been found theoretically. In their researches, these authors believed that the velocity profiles will be oscillatory, and they assumed that these oscillations are instantaneous responses of the externally applied electric fields, i.e., they have the same frequency. True, the electro-osmotic flows should really be generated by the applied time-periodic electric fields, and the flows may be time periodic. But, as we know, there is a phase difference between phase voltage and phase current, and the flow in the microchannel should need some time to start. In other words, there is a phase difference between the applied electric fields and the electro-osmotic flows. On the other hand, on the basis of the aforementioned ‘assumption’, the obtained analytical solutions of velocities are represented as complex functions, which is unreasonable in physics. So, as a result, the solutions given in these research papers are, generally speaking, incorrect.
In fact, the phase shift between the applied electric field and the flow response has been proved by Nayak [9], as well as some other researchers [10]. The steady/unsteady electro-osmotic flow in an infinitely extended cylindrical channel with diameters ranging from 10 to 100 nm has been investigated by Nayak [9], and the degree of the phase shift between the velocity field and the applied electric field is found numerically. Using the backwards-Euler time stepping numerical method, Luo [10] clarified the relationship between the changes in the axial-flow velocity and the intensity of the applied electric field. Erickson and Li [11] developed the analytical solution for the AC electro-osmotic flow through a rectangular microchannel for the case of a sinusoidal applied electric field. Shilov et al.[12] discussed the mechanisms for different times after the application of the electrical field according to the relationship between the dipole moment and the electrophoretic mobility.
The aim of the present paper is to present the analytical solutions for the time-periodic electro-osmotic flow of Newtonian fluids through a microannulus. Analytical solutions are rare. Not only do they represent electro-osmotic flows through fundamental cross-sectional shapes but they also serve as standards for asymptotic and fully numerical methods. Most important of all, some new physical and chemical phenomena can be found from the analytical solutions.
2 Governing equations
The motions of an ionized, incompressible Newtonian fluid with electro-osmotic body forces are governed by the following Navier-Stokes equation:
where P is the pressure, ρ is the fluid density, μ is the dynamic viscosity, and the tensor V is a divergence-free velocity field, i.e., subject to the non-slip boundary conditions on the walls, is the externally applied electric field, and is the electric charge density, which can be expressed by a potential distribution ψ; we have
and
here is the bulk electrolyte concentration of a binary electrolyte dissociating into cations and anions of valence , e is the electron charge, is the Boltzmann constant, and T is the absolute temperature.
In the present study, we assume the surface potential is small enough, then with the help of Debye-Hückel approximation and cylindrical coordinate system , Eq. (2) is linearized to
where is the Debye-Hückel parameter and means the Debye length.
Because of the effect of the electric field, the fluid in the capillary will flow along the axis direction. Neglecting the pressure gradient along the axis, the Cauchy momentum equation in cylindrical coordinate system with AC electric field can be expressed as
where is the axial velocity, t is time, and is AC electric field, is the magnitude, and ω is the frequency of the unsteady external electric field E.
In the present research, the geometric shape of the microchannel is considered, as shown in Figure 1. An electrolyte fluid is flowing unsteadily in the annular region between two uniform coaxial circular cylinders with inner radius () and outer radius . The chemical interaction of the electrolyte liquid and solid wall generates an electric double layer (EDL), a very thin charged liquid layer at the solid-liquid interface. The outer and inner wall zeta potentials are and , respectively. Here, and are small enough, so that the Debye-Hückel linearization approximation is available.
Define the following dimensionless variables:
here . Substituting the above dimensionless variables into (4) and (5) yields the governing equation for the potential distribution with boundary conditions
and the equations for the flow with boundary and initial conditions
Here , and are normalized wall potentials.
3 Analytical solutions
The general solution of (7) has the form
where and are undetermined constants, and are the modified Bessel functions of the first and second kind of order zero, respectively. Considering the boundary conditions (8) and (9), we have
and the constants A and B are
and
Here is defined as the ratio of the zeta potentials of the inner wall to that of the outer wall.
We now consider the solution of the governing equation (10); here it is not convenient to use the classical method of separation of variables because of the nonhomogeneities of the master governing equation (10). It is for this reason that we consider the integral-transform technique, and this method provides a systematic, efficient, and straightforward approach for the solution of both homogeneous and nonhomogeneous, steady-state, and time-dependent initial and boundary-value problems.
The integral-transform pair in the r variable for the function is defined as [13]
where
and is the m th positive root of .
Applying the above integral transform (16) to (7)-(13) yields
Equation (21) is an ordinary differential equation with initial condition (22), and its solution can be given directly as
where is the phase difference or the phase shift, and can be obtained from (2) with the help of the aforementioned integral transform (16). Then, substituting into Eq. (23) yields the distribution of the velocity in the capillary,
4 Results and discussion
4.1 The generation of the flow
From the expression (24), we find that the velocity field of electro-osmotic flow in the capillary generated by the external applied electric field is not time periodic. In particular, the distribution of the velocity can be written as a sum of a time-dependent oscillating part and a time-dependent generating part :
where
and
Here and are defined by (18) and (19), respectively.
The first five roots of for different α are listed in Table 1. It can be seen that the minimum of increases with increasing α. When , for the minimum of we have , and . As a result, we can draw the conclusion that the generating part of the solution (24) will tend to zero in a very short time, which results in the electro-osmotic flow reaching a steady ‘periodic’ state. Additionally, it is worth pointing out that the increasing frequency of the applied external electric field accelerates the generation of flow in the microannulus.
In the sense of the above discussion, the generating part can also be called the transient part. In other words, the electro-osmotic flow generated by the AC electric field is quasi-periodic. In spite of this, the generating part of the solution is very important for the researcher in this field, since it explains both the characteristics of electro-osmotic flow and the practical applications due to rapid development of the biochip technology [14]. Furthermore, in a study of the stability of a colloidal system, Overbeek [15] pointed out that the relaxation time for surface charges (about 10−6 to 104 s) and the time scale for Brownian coagulation (about 10−7 to 10−5 s) are very different; the aggregation of colloidal particles may occur earlier than the equilibrium of the electrical conditions near a surface. In these cases, the steady-state analysis on the electrical condition near a charged surface is unrealistic, and an extension of the conventional treatment to a temporal description is inevitable, and this is the significance of the present study.
4.2 Special cases
In particular, when , i.e., , where is the Heaviside step function,
which yields
Then we get the distribution of velocity profile when the applied external electric field remains constant from time , i.e., the electric field follows a step-change:
When , then we have , i.e., the inner wall and outer one have the same zeta potential,
When , , i.e., the inner wall is not charged, then the solution of velocity (24) reduces to
4.3 Effects of K on velocity profiles and flow rates
For given values , , , and , plots of the normalized velocity as a function of the non-dimensional radius for selected values of time t and for four different values of the electrokinetic width K are shown in Figure 2. It is clear from this figure that the flow in the microannulus is similar to general pipe flow for small K as shown in the figure when , because the EDL is thicker for small K. When the value of K is large enough, the velocity of the flow increases with increasing electrokinetic width K, and this phenomenon can also be proved analytically from the expression (24), where the term is included. The same conclusion is true for the flow rate of fluid, which can be obtained by integrating Eq. (24),
To demonstrate the effect of K on the flow rate, plots of the flow rate as a function of t using Eq. (33), for selected values of K and ω, are displayed in Figure 3. The physical mechanism for the above phenomena is that a large K indicates a thinner EDL, then the region of bulk flow (mass flow) becomes wider, which results in more mass of fluid through the pipe during the same period.
Additionally, some authors drew the conclusion that the flow rate is proportional to the cross-sectional area of the channel for large K, and the flow rate is quadratic as in the leading-order behavior for small K[16]–[18]; it is noteworthy that this conclusion is incorrect. In this research, the authors discussed the asymptotic expansion of for large and small K, respectively, with the help of series summation formulas. Unfortunately, they ignored the fact that is a monotonically increasing infinite subsequence. In other words, for a given K, no matter how large it is, there exists a natural number N such that if , then .
On the other hand, from Figure 3 and Eqs. (26)-(27), it is important to note that because the transient part of the velocity solution decays in a short time, if the researchers consider the electro-osmotic flow of fluids for time , the time-dependent oscillating part can be used as a good approximation for simplicity of computation.
4.4 Effects of ω on velocity profiles
Figure 4 indicates the axial velocity distributions along the axial direction during the period of transient response from to for different values of ω. For given , , , and , the oscillation of flow is enhanced by the increasing frequency of the applied external electric field. However, the mean velocity of flow decreases as ω increases from to . An interesting phenomenon is found: there is almost no flow in the areas far away from the EDL when ω is large enough, as in the case of . As the frequency of external electric field increases, the flow rate gets closer to zero, and it keeps its oscillation, as shown in Figure 5, which is obtained from (33) as a function of ω at fixed . The physical interpretation is that the large frequency AC electric field makes the ions in the fluid oscillate around an equilibrium position. Mathematically and analytically, the increasing ω decreases the term in Eq. (24),
as a result, the velocity reliably decreases with increasing frequency of the applied external electric field. Away from the walls of the microannulus, the EDL can be divided into the compact layer and the diffuse double layer. Within the diffuse layer, the motion of the ions is subject to the zeta potential of the EDL. Then we have the conclusion that the main contribution to the flow rate is the zeta potential when the frequency of the applied AC electric field is large enough.
Additionally, the effects of ω on the normalized flow rate and phase shift are shown in Figure 6 for given . With the increasing frequency of the applied external electrical field, there appears to be a decrease of the phase shift between voltage and flow. In fact, according to Eq. (23), mathematically tends to be zero for large enough ω. At the same time, Figure 6 shows the increasing frequency decreases the flow rate , which is consistent with the results shown in Figure 4 and Figure 5.
4.5 Effects of wall zeta potential ratio β on velocity profiles
Electro-osmotic flow is the result of the interaction between the applied external electric field and the electric double layer, so it is necessary to discuss the effect of the wall’s zeta potential ratio β of the inner to the outer cylinder on the velocity profiles. For each fixed time t, a sequence of curves for different values of β is shown in Figure 7. From the figure, it can be seen that the direction of the flow is directly correlated with the polarity of the charged channel wall. When , flow reversal is observed, which is caused by an adverse zeta potentials on the walls. When , the velocity of flow increases with increasing β, this can also be proved theoretically with the help of the average velocity in the microannulus, i.e., , and the following relationship is obtained:
5 Summary and conclusion
It has not been a accurate task to find the analytic solutions for time-periodic electro-osmotic flow (AC electro-osmotic). In present research, we pointed out and corrected the errors in the published research articles in this field, and we obtained an analytical solutions of distribution of velocity profiles and flow rates for time-periodic electro-osmotic flow in microannulus. With the help of numerical plots, our analysis has resulted some remarks. The velocity field of the electro-osmosis flow consists of two parts, a time-dependent oscillating part and a time-dependent generating or transient one. The transient part tends to zero very fast with the increasing time. The electro-osmosis flow driven by an alternating electric field is not periodic in time, but quasi-periodic. There is a phase shift between applied external electric field and the electro-osmotic flow, which is only dependent of the frequency of external electric field, and it is less than . The increasing values of K decrease the average velocity and flow rate, particularly, there is almost no flow in the region far from the walls when ω is large enough. Additionally, the increasing frequency of applied external electric field accelerates the generation of flow in the microannulus. The corresponding physical interpretation is given.
References
Effenhauser CS, Manz A, Widmer HM: Manipulation of sample fractions on a capillary electrophoresis chip. Anal. Chem. 1995, 67: 2284-2287. 10.1021/ac00109a057
Raymond DE, Manz A, Widmer HM: Continuous separation of high molecular weight compounds using a microliter volume free-flow electrophoresis microstructure. Anal. Chem. 1996, 68: 2515-2522. 10.1021/ac950766v
Harrison DJ, Manz A, Fan ZH, Ludi H, Widmer HM: Capillary electrophoresis and sample injection systems integrated on a planar glass chip. Anal. Chem. 1992, 64: 1926-1932. 10.1021/ac00041a030
Ramsey RS, Ramsey JM: Generating electrospray from microchip devices using electroosmotic pumping. Anal. Chem. 1997, 69: 1174-1178. 10.1021/ac9610671
Dutta P, Beskok A: Analytical solution of time periodic electroosmotic flows: analogies to Stokes’ second problem. Anal. Chem. 2001, 73: 5097-5102. 10.1021/ac015546y
Xuan X, Li D: Electroosmotic flow in microchannels with arbitrary geometry and arbitrary distribution of wall charge. J. Colloid Interface Sci. 2005, 289: 291-303. 10.1016/j.jcis.2005.03.069
Jian Y, Yang L, Liu Q: Time periodic electro-osmotic flow through a microannulus. Phys. Fluids 2010., 22: 10.1063/1.3358473
Bao LP, Jian YJ, Chang L, Su J, Zhang HY, Liu QS: Time periodic electroosmotic flow of the generalized Maxwell fluids in a semicircular microchannel. Commun. Theor. Phys. 2013, 59: 615-622. 10.1088/0253-6102/59/5/16
Nayak AK: An analysis of steady/unsteady electroosmotic flows through charged cylindrical nano-channels. Theor. Comput. Fluid Dyn. 2013.
Luo WJ: Transient electroosmotic flow induced by DC or AC electric fields in a curved microtube. J. Colloid Interface Sci. 2004, 278: 497-507. 10.1016/j.jcis.2004.06.017
Erickson D, Li D: Analysis of alternating current electroosmotic flows in a rectangular microchannel. Langmuir 2003, 19: 5421-5430. 10.1021/la027035s
Shilov VN, Delgado AV, González-Caballero F, Horno J, López-García JJ, Grosse C: Polarization of the electrical double layer. Time evolution after application of an electric field. J. Colloid Interface Sci. 2000, 232: 141-148. 10.1006/jcis.2000.7152
Özisik MN: Heat Conduction. Wiley, New York; 1993.
Wong PK, Chen C, Wang T, Ho C: Electrokinetic bioprocessor for concentrating cells and molecules. Anal. Chem. 2004, 76: 6908-6914. 10.1021/ac049479u
Overbeek JTG: Recent developments in the understanding of colloid stability. J. Colloid Interface Sci. 1977, 58: 408-422. 10.1016/0021-9797(77)90151-5
Wang CY, Liu YH, Chang CC: Analytical solution of electro-osmotic flow in a semicircular microchannel. Phys. Fluids 2008., 20: 10.1063/1.2939399
Wang CY, Chang CC: EOF using the Ritz method: application to superelliptic microchannels. Electrophoresis 2007, 28: 3296-3301. 10.1002/elps.200700001
Chang CC, Kuo CY, Wang CY: Unsteady electroosmosis in a microchannel with Poisson-Boltzmann charge distribution. Electrophoresis 2011, 32: 3341-3347. 10.1002/elps.201100181
Acknowledgements
This work is supported by the National Natural Science Foundation of China (Nos. 11002083, 51279093), the National Basic Research Program of China (2013CB0360000).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Wang, S., Zhao, M. Closed-form solutions of transient electro-osmotic flow driven by AC electric field in a microannulus. Bound Value Probl 2014, 215 (2014). https://doi.org/10.1186/s13661-014-0215-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-014-0215-2