Abstract and Applied Analysis

Volume 2014 (2014), Article ID 789147, 10 pages

http://dx.doi.org/10.1155/2014/789147

Research Article

Solution of Time Periodic Electroosmosis Flow with Slip Boundary

Qian Sun,¹ Yonghong Wu,¹ Lishan Liu,² and B. Wiwatanapataphee³

¹Department of Mathematics and Statistics, Curtin University, Perth, WA 6845, Australia

²School of Mathematical Sciences, Qufu Normal University, Qufu, Shandong 273165, China

³Department of Mathematics, Faculty of Science, Mahidol University, Bangkok 10400, Thailand

Received 27 March 2014; Revised 28 April 2014; Accepted 22 May 2014; Published 16 June 2014

Academic Editor: Xinguang Zhang

Copyright © 2014 Qian Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Recent research confirms that slip of a fluid on the solid surface occurs at micrometer scale. Slip on solid surface may cause the change of interior material deformation which consequently leads to the change of velocity profile and stress field. This paper concerns the time periodic electroosmotic flow in a channel with slip boundary driven by an alternating electric field, which arises from the study of particle manipulation and separation such as flow pumping and mixing enhancement. Although exact solutions to various flow problems of electroosmotic flows under the no-slip condition have been obtained, exact solutions for problems under slip boundary conditions have seldom been addressed. In this paper, an exact solution is derived for the time periodic electroosmotic flow in two-dimensional straight channels under slip boundary conditions.

1. Introduction

One of the major focuses in scientific research during the past decades has been the study of material behaviour at micro- and nanoscale. The subject of micro- and nanofluidics concerns mass and momentum transfer in micro- and nanoscale systems or around micro- or nanosized objects, as well as the design and application of biological and engineering devices and systems involving fluid flow through microchannels [1–3]. The study of microflows has attracted great attention from the science and engineering communities in the effort of developing a better understanding of the mechanism of microflows and consequently developing better models and control. Traditionally, in the model for fluid flows, the no-slip boundary condition is used; that is, the fluid velocity relative to the solid is assumed to be zero on the fluid-solid interface. However, it has been established that the no-slip condition is a hypothesis instead of a condition deduced from any principle. Evidence of slip of a fluid on the solid surface was reported by many researchers [4, 5].

In recent years, time periodic electroosmosis has been studied intensively due to the recent development of microfluidic devices. Time periodic electroosmosis is known as AC electroosmosis, and the driving force is generated by an alternative electric field. Dose and Guiochon published numerical results for impulsive started electroosmotic flow [6]. Green et al. studied the electroosmosis on a planner microelectrode under both steady and unsteady fields [7]. Soderman and Jonsson investigated the electroosmotic flows for various geometries and derived various analytical results for the flow in two-dimensional microchannels and for the flow over a flat plate [8]. By using the green function method, Kang and coworkers [9] analysed the instantaneous velocity of the flow through a microchannel with time-dependent electric field. Jian et al. [10] studied the time periodic electroosmotic flow through circular microchannels. They discussed two different kinds of situations, which are the electroosmotic flow through the parallel plates and circular tubes.

Although analytical and numerical solutions to many flow problems of Newtonian fluids have been established under the no-slip assumption [11–13], very few exact solutions for time-dependent problems have been obtained for the slip case [14, 15]. Hence, we investigate, in this paper, the time periodic electroosmotic flows in a microchannel under the Navier slip assumption. The rest of the paper is organized as follows. In Section 2, the mathematical equations for the channel flow are presented, which include the partial differential equations and the boundary conditions based on the slip assumption. In Section 3, the general solution for the velocity of the fluid in the channel is derived with some arbitrary constants. Then the slip boundary condition and symmetric condition are used to determine the arbitrary constants. Numerical investigations are then given in Section 4, followed by a concluding remark in Section 5.

2. Problem Description and Mathematical Formation

The problem considered is the electrically driven flow of fluid through channels with height . The geometry of the channels and the coordinates used are shown in Figure 1.

Figure 1: The flow channel and coordinate system used.

Based on the principle of continuous mechanics, the governing equations for the flow of an incompressible ionized fluid are the standard Navier-Stokes equations:

where is the pressure, and are, respectively, the fluid density and the electric charge density, and is the velocity vector with being the fluid viscosity. The last term on the right side denotes the electroosmotic body forces and denotes the externally applied electric field.

We assume fully developed symmetric channel flow, and thus the swirling flow and consequently the velocity components in and directions can be neglected. This allows us to use the superposition principle because the resulted unsteady Navier-Stokes equations for this case are linear in velocity. This means that we can solve exactly the pressure and electroosmotically driven flows separately and then combine the two solutions to get the solution of the problem with the combined effects. The form of the electroosmotic body force, in general, depends on the form of the electric field applied on the system. In this paper, we will focus on sinusoidally driven flows.

Taking into account the above assumptions, (1) reduces to

where and are, respectively, the magnitude and the frequency of the external electric field . With (2), we can then determine the velocity of the flow driven by the time periodic electric field.

Based on the theory of electric field [16], the in (2) can be determined by

where denotes the average number of positive or negative ions in the buffer, is the valance, is the electron charge, is the absolute temperature, is the electric potential, and is the Boltzmann constant. Let be the normalized electroosmotic potential in which is the zeta potential and let be a parameter; then we have

Based on the work in [16], is governed by the Poisson-Boltzmann equation:

in which is the permittivity, and we have assumed that the effect of time fluctuation of the external electric field on the electric potential is negligible.

It has been established in the literature that the effective thickness of the electric double layer (EDL) depends on the ionic energy parameter . For example, for , the effective electric double layer thickness is about 4.217 , where is related to other system parameters by

In most microfluidic applications, the value of is about 10–100 nm, which is generally much smaller than the height of a typical microchannel. Hence, the electroosmotic potential generally reaches zero in the channel center, and its distribution along the channel height can be expressed by [16, 17]

where denotes the distance from the channel wall. From the definition of the Debye-Huckel parameter , we can write the electroosmotic body force as

where is the Helmholtz-Smoluchowski velocity. Then, (2) can be written as

To completely determine the fluid flow in the channel, the differential equation (2) must be supplemented by a set of boundary and initial conditions. To describe the slip characteristics of a fluid on the solid surface, the Navier slip condition is used; namely, the velocity of fluid, relative to the solid surface, is assumed to be proportional to the sheer stress on the fluid-solid surface; that is,

where is the so-called slip length. We should address here that, for , condition (10) reduces to the no-slip boundary condition.

In addition, due to the symmetry of velocity about the symmetric plane , we have

Now the partial differential equation (9) together with conditions (10) and (11) constitutes a complete boundary value problem for the fluid flow in the channel.

3. Solution of the Boundary Value Problems

By introducing the characteristic time and the characteristic length , (9) can be expressed in normalized form

in which denotes the nondimensional time, denotes the nondimensional distance, and denotes the nondimensional velocity. From (12), we have

where is a normalized parameter expressed as a function of the Debye length , the electric field excitation frequency (), and the kinematic viscosity . The parameter represents the ratio of the Debye length to a diffusion length scale () given by

Let be the solution of the following equation:

Then, from the linearity and the superposition principle, the solution of (13) will be . It is noted that (15) admits solution of the form where we have written as the imaginary part of , and thus we proceed to solve (15) first.

The derivative of with respect to time can then be written as

and the second space derivatives can be written as

Therefore (15) becomes

For the above to be true for all values, we require

The general solution of (19) includes a complimentary solution and a particular solution. The complimentary solution can be derived by solving the associated homogenous equation, while the particular solution can be obtained by using the variation of parameter methods. Through some calculation, we have

and consequently we have

Now we proceed to determine the arbitrary constants in (21). In terms of the nondimensional time , the nondimensional distance , and the nondimensional velocity , the boundary conditions (10) and (11) become

As , (22) will be satisfied if satisfies the following boundary conditions:

Noticing that , the above boundary conditions become

To apply the boundary condition (24), we first need to calculate the spatial derivative of . From (7) and (20), we have

At the centre of the channel, and and so from (24) and (25) we have

From (20), (25), and the first equation of (24), we get

and hence

Substituting (28) into (26) yields

4. Numerical Investigation

In this section, we investigate the characteristic of the electric driven flow. In particular, we will investigate the influence of various model parameters on the flow behaviour, including the value, the energy parameter , the channel height , and the slip parameter . Table 1 shows the solutions of and for various sets of the model parameters.

Table 1: and values for different sets of model parameters.

0.1
1

3

5

0.01
1

3

5

0.001
1

3

5

In the first example of investigation, we set the channel height to be where the typically has value of hundreds of nanometres. Other parameters are chosen as and . Figures 2 and 3 show the distribution of velocity along the channel height obtained from the model with various and values at , as functions of for different values of the parameters and . Figure 4 presents the velocity distribution and , at for and . The results show that as increases, the thickness of the layer with significant electrically driven flow decreases and essentially the velocity at the centre of the channel is zero for as shown in Figures 3(d) and 4(a). It is also found that the slip parameter has significant influence on the flow behaviour. At , as increases from 0 to 0.5, the flow velocity increases from 1 to about 1.685.

(a)

(b)

(c)

(d)

Figure 2: Effect of slip length on velocity profile obtained from the model with , , and for various values: (a) ; (b) ; (c) ; (d) .

(a)

(b)

(c)

(d)

Figure 3: Effect of on velocity profile obtained from the model with , , and for various slip lengths: (a) ; (b) ; (c) ; (d) .

(a) when

(b) when

Figure 4: 3D plot of and obtained from the model with , , and .

To investigate the effect of the energy parameter on the flow behaviour, we choose four values of including 1, 3, 5, and 7. Figure 5 shows the influence of on the velocity distribution for , , and two different values of 0.001 and 0.03. For a channel with fixed dimension, filled with a given buffer solution, the value of depends directly on the frequency of the external electric field. Hence, in practice, it is possible to obtain the desired velocity field by controlling the frequency by electric switch [7]. We found that as increases, the flow velocity increases. At , as increases from 1 to 7, the flow velocity at the centre increases from 1 to 1.4 as shown in Figure 5. At , as increases from 1 to 7, the flow at the centre approaches zero.

(a)

(b)

Figure 5: Effect of the energy parameter on the velocity profile obtained from the model with , , and for two different values: (a) and (b) .

At the typical values of , , and , the distributions of velocity along the height of the channel at various instants of time in an oscillation cycle of electric field are shown in Figure 6. The magnitude of the velocity varies with time, which is because of the fluctuation of the external electric field in time.

Figure 6: 2D and 3D plots of velocity profile at various instants of time obtained from the model with , , and .

Figure 7 shows the influence of channel height on the magnitude of velocity. It can be noted that the magnitude of velocity is influenced by the channel height , while the pattern of velocity profile along the channel is the same. Figure 8 shows the influence of on the first derivative of the velocity.

(a)

(b)

Figure 7: Effect of the channel half-height on the velocity profile obtained from the model with , , and for two different values: (a) and (b) .

(a) for

(b) for

Figure 8: 3D plot of for two values: (a) and (b) .

5. Concluding Results

In this paper, we study the time periodic electroosmotic flow in a two-dimensional microchannel with slip boundary. Some analytical and numerical solutions for the slip velocity field under various different conditions have been obtained.

Our results show that the flow behaviour depends on the parameters and as well as the slip parameter and channel height. The parameter represents the ratio of the electric double layer thickness to a characteristic diffusion length scale , while is an energy parameter.

We found that, depending on the value of the parameter , the flow pattern may change from oscillatory plug flow to uniform flat flow. Generally, as increases, the velocity at the centre drops to zero dramatically. We also found that the slip parameter has significant influence on the flow behaviour. As increases, the flow rate increases.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Qian Sun, Yonghong Wu, and Lishan Liu were supported financially by the National Natural Science Foundation of China (11371221), the Program for Scientific Research Innovation Team in Colleges and Universities of Shandong Province, and the Specialized Research Foundation for the Doctoral Program of Higher Education of China (20123705110001).

References

B. Bourlon, J. Wong, C. Mikó, L. Forró, and M. Bockrath, “A nanoscale probe for fluidic and ionic transport,” Nature Nanotechnology, vol. 2, no. 2, pp. 104–107, 2007.

H. Herwig and O. Hausner, “Critical view on new results in micro-fluid mechanics : an example,” International Journal of Heat and Mass Transfer, vol. 46, no. 5, pp. 935–937, 2003.

H. Huang, T. S. Lee, and C. Shu, “Lattice Boltzmann method simulation gas slip flow in long microtubes,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 17, no. 6, pp. 587–607, 2007.

G. Chauveteau, “Rodlike polymer solution flow through fine pores: influence of pore size on rheological behavior,” Journal of Rheology, vol. 26, no. 2, pp. 111–142, 1982.

R. Pit, H. Hervet, and L. Leger, “A sequential quadratic programming algorithm for nonconvex, nonsmooth constrained optimization,” Physics Review Letters, vol. 2, no. 85, p. 980, 2000.

E. V. Dose and G. J. Guiochon, “Time scales of transient processes in capillary electrophoresis,” Journal of Chromatography, vol. 652, no. 1, pp. 263–275, 1993.

N. G. Green, A. Ramos, A. Gonzalez, H. Morgan, and A. Castellanos, “Fluid flow induced by nonuniform ac electric fields in electrolytes on microelectrodes. III. Observation of streamlines and numerical simulation,” Physical Review E, vol. 66, no. 2, Article ID 026305, 2002.

O. Soderman and B. Jonsson, “Electro-osmosis: velocity profiles in different geometries with both temporal and spatial resolution,” The Journal of Chemical Physics, vol. 105, no. 23, pp. 10300–10311, 1996.

Y. J. Kang, C. Yang, and X. Y. Huang, “Dynamic aspects of electroosmotic flow in a cylindrical microcapillary,” International Journal of Engineering Science, vol. 40, no. 20, pp. 2203–2221, 2002.

Y. J. Jian, L. G. Yang, and Q. S. Liu, “Time periodic electro-osmotic flow through a microannulus,” Physics of Fluids, vol. 22, no. 4, Article ID 042001, 2010.

B. Wiwatanapataphee, Y. H. Wu, J. Archapitak, P. F. Siew, and B. Unyong, “A numerical study of the turbulent flow of molten steel in a domain with a phase-change boundary,” Journal of Computational and Applied Mathematics, vol. 166, no. 1, pp. 307–319, 2004.

B. Wiwatanapataphee, D. Poltem, Y. H. Wu, and Y. Lenbury, “Simulation of pulsatile flow of blood in stenosed coronary artery bypass with graft,” Mathematical Biosciences and Engineering, vol. 3, no. 2, pp. 371–383, 2006.

J. M. Hill, Y. H. Wu, and B. Wiwatanapataphee, “Analysis of flux flow and the formation of oscillation marks in the continuous caster,” Journal of Engineering Mathematics, vol. 36, no. 4, pp. 311–326, 1999.

Y. H. Wu, B. Wiwatanapataphee, and M. Hu, “Pressure-driven transient flows of Newtonian fluids through microtubes with slip boundary,” Physica A: Statistical Mechanics and Its Applications, vol. 387, no. 24, pp. 5979–5990, 2008.

B. Wiwatanapataphee, Y. H. Wu, M. Hu, and K. Chayantrakom, “A study of transient flows of Newtonian fluids through micro-annuals with a slip boundary,” Journal of Physics A: Mathematical and Theoretical, vol. 42, no. 6, Article ID 065206, 2009.

P. Dutta and A. Beskok, “Analytical solution of time periodic electroosmotic flows: analogies to stokes' second problem,” Analytical Chemistry, vol. 73, no. 21, pp. 5097–5102, 2001.

R. J. Hunter, Zeta Potential in Colloid Science: Principles and Applications, Academic Press, New York, NY, USA, 1981.