Abstract and Applied Analysis

Volume 2014 (2014), Article ID 530605, 13 pages

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

Research Article

Transient Flows of Newtonian Fluid through a Rectangular Microchannel with Slip Boundary

Benchawan Wiwatanapataphee,¹ Yong Hong Wu,² and Suharsono Suharsono³

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

²Department of Mathematics & Statistics, Curtin University of Technology, Perth, WA 6845, Australia

³Department of Mathematics, Lampung University, Lampung 35145, Indonesia

Received 12 March 2014; Accepted 22 May 2014; Published 29 June 2014

Academic Editor: Xinguang Zhang

Copyright © 2014 Benchawan Wiwatanapataphee 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.

We study the transient flow of a Newtonian fluid in rectangular microchannels taking into account boundary slip. An exact solution is derived by using the separation of variables in space and Fourier series expansion in time. It is found that, for different forms of driving pressure field, the effect of boundary slip on the flow behavior is qualitatively different. If the pressure gradient is constant, the flow rate is almost linearly proportional to the slip parameter when is large; if the pressure gradient is in a waveform, as the slip parameter increases, the amplitude of the flow rate increases until approaching a constant value when becomes sufficiently large.

1. Introduction

In recent years, many researchers worldwide focus on the study of behavior of materials at micro- and nanoscales [1, 2], leading to the development of many biological and engineering systems and devices. Most of these systems and devices involve fluid flow in microchannels, called microflows [1, 3–7]. Typical examples include drug delivery systems [8], fuel cell devices, energy conversion, and biological sensing devices [9]. As the functional characteristics of these systems depend on the behavior of fluid flow in the systems, the study of microflows is important and has attracted more and more attention from the engineering and science communities in order to derive a better understanding of the mechanism of microflows and consequently better design and control of the devices and systems [1, 6, 10].

The field equations governing the flow of Newtonian fluids are the continuity equation and the Navier-Stokes equations. These equations are subject to a set of boundary conditions. The no-slip boundary condition is used traditionally; namely, the tangential fluid velocity relative to the solid is zero on the fluid-solid interface [11]. However, recent molecular dynamic simulations and experiments in micrometer scale have shown that the fluid flow in microsystems is granular and slip may occur on the fluid-solid interface [10, 12–17]. Hence, for the study of microflow, it is important to take into account the boundary slip of fluids on the fluid-solid interfaces.

Over the last few decades, intensive research has been carried out to study various problems of fluid flow with the no-slip assumption or a slip boundary condition [1, 5, 17–33]. For flows under the no-slip assumption, exact solutions to many problems have been obtained and are available in the literature [11, 34–37]. But for slip flow very few exact solutions are available in the literature. Exact solutions for the fluid flow in circular microtubes and circular microannuals with boundary slip have been derived and discussed in the papers [17, 25, 28].

For microchannels with rectangular cross-section, a no-slip solution has been obtained [38–43]. For the slip case steady state solution has also been obtained [15, 35, 44–48]. However, so far no exact solution has been derived for the transient flow of fluids through rectangular channel under pulsatile pressure. As many microsystems and devices have microchannels of rectangular cross-section, it is important to derive exact analytical solutions for the behavior of transient flow through rectangular microchannels with slip boundary.

Based on the current development in the field, in this paper, we study the time-dependent flow of incompressible Newtonian fluids through a rectangular microchannel taking into account boundary slip. The rest of the paper is organized as follows. In Section 2, we give the underlying initial boundary value problem for the transient slip flow. In Section 3, we derive the exact solution for the velocity field. In Section 4, an analysis is conducted to study the effect of the slip length on the flow behavior. Finally, a conclusion is presented in Section 5.

2. Governing Boundary Value Problem

Consider the unsteady flow of an incompressible Newtonian fluid in a rectangular channel of cross-section dimension with the -axes being in the axial direction as shown in Figure 1. The field equations governing the flow include the Navier-Stokes equations and the continuity equation. As the flow is symmetric about the -plane and the -plane and is fully developed, there is no cross-sectional flow and thus the velocity components in the and directions vanish; that is, .

Figure 1: Rectangular microchannel.

Hence, from the continuity equation

and the Navier-Stokes equation

the axial velocity, , is governed by the following equation:

As a large class of functions may be expressed by Fourier series, we consider, in this work, the flow of a fluid driven by the pressure gradient that may be expressed in the form of Fourier series given by

To define the problem completely, we supplement the field equation by the boundary condition. To take into account the boundary slip, the so-called Navier-slip boundary condition is used. On the fluid-solid interface and , the axial fluid velocity, relative to the solid surface, is assumed to be proportional to the shear stress on the interface. Let be the unit normal vector of the surface of the fluid, and let be the positive tangential direction. Also let the fluid velocity on the wall direction be , and let the velocity of the solid in the tangential direction of the surface be . Then, as shown in our previous work [17], the Navier-slip boundary condition can be expressed by

where the negative sign indicates that the surface traction force which acted on the fluid by the solid is opposite to the tangential velocity of fluid relative to the solid. For our problems in the coordinate system, and . On the surface , and , and so and and consequently (5), on the surface , becomes

On the surface , we have and , and hence and and consequently (5), on the surface , becomes

Similarly, the boundary condition on the surface is

We will remark here that, for , the slip boundary conditions (6)–(8) reduce to the no-slip condition; on the other extreme, where , (6)–(8) become the surface traction conditions for perfectly smooth surfaces; that is, .

3. Exact Solutions for Transient Velocity and Stress Fields

To solve the partial differential equation (3), complex numbers are used to express the Fourier series for the pressure gradient; namely,

where and .

From the symmetry of the problem and the linearity of (3), we get , where is defined by

Through some mathematical analysis, we get that, for , the boundary value problem (10) has solution of the form

where is determined by

with and being governed, respectively, by the following boundary value problems:

Thus, the remaining work for finding and consequently and then is to solve the two BVPs: BVP1 and BVP2. We first solve BVP1 by the separation of variables. From the PDE and the homogeneous boundary conditions of BVP1, we obtain

where are the roots of the nonlinear equation

are respectively eigenvalues and eigenfunctions of BVP1 and

The eigenvalues are the solution of the following nonlinear equation:

which has infinite number of solutions as shown by the intersection of the graphs and where in Figure 2.

Figure 2: Infinite number of solutions of the two functions and .

It can also be proved that the eigenfunctions are orthogonal and thus the coefficients can be determined from the nonhomogeneous boundary condition of BVP1 by

Similarly, the solution of the BVP2 is

where are the roots of the nonlinear equation

Substituting (14) and (20) into (12) yields the solution

For , proceeding as for finding , we obtain

where and are as defined in (17) and (22) with ; that is, and ; and are as follows:

Hence collecting all solutions of the subproblems, we have

where and are determined, respectively, by (15) and (21), and are defined by (17) and (22), respectively, and are defined by (26) and (27), and and are defined by (19) and (23), respectively.

Now we determine the exact solutions of the flow rate and the stresses in the fluid. From the velocity solution (28), we obtain the flow rate as follows:

where and denote, respectively, the flow rate corresponding to the constant component and the nth harmonic component of the pressure gradient and

The stress tensor in the fluid can then be determined by the constitutive equation

where is an identity matrix and the deformation rate depends on the velocity by

As , we get

From (28) and (32)–(34), we obtain and

Thus, from the constitutive equation (32), we get

where denotes the pressure gradient while is arbitrary and may be chosen to meet certain pressure condition.

4. Investigation of the Effect of Boundary Slip on the Flow Behavior

Based on the exact solutions obtained, we investigate the flow behaviour and the influence of the slip length in this section. As a general pressure field can be expressed by a Fourier series in the form of (4), from the principle of superposition, the exact solution of the problem is the superposition of the solutions corresponding to the constant pressure gradient plus the solutions corresponding to the sine or cosine waveform pressure gradients. In this work, without loss of generality, we consider flow problems under two different cases of driving pressure fields including the case with a constant pressure gradient and the case with a sine waveform pressure gradient. For simplicity, we introduce four dimensionless variables as follows:

Case 1 (). For this case, and for . The constant pressure gradient means that the pressure gradient does not depend on time. From (28), (30), (35), and (36), we get the following normalized velocity, flow rate, and shear stresses:

As the pressure gradient does not depend on time, the normalized velocity and flow rate as well as stresses are influenced by the slip length only, which is implicitly contained in and .

Case 2 (). For this case, the pressure gradient is sinusoidal with amplitude , and , , for . From (27), we have

Let

where

Then

Using (19) and (23), through a lengthy derivation, we obtain

where

For convenience in discussion, transform (47) into the following form:

where and denote the amplitude and phase angles of the normalized velocity defined, respectively, by

The flow rate is

Using real arithmetic, through lengthy calculation, we obtain the normalized flow rate

where

In this study, we analyze the flow pattern through the rectangular microchannel having the same size of the cross-sectional area of for ,  Pa·s, , , and .

For Case 1, as the pressure gradient is constant, the normalized velocity and flow rate as well as shear stress vary with the slip length and the geometry of the cross-section, namely, the ratio of the cross-sectional area. The influences of the slip length and the ratio on the flow behavior are demonstrated by analyzing the solutions graphically. Figures 3 and 4 show, respectively, the three-dimensional and two-dimensional velocity profiles on a cross-section of the channel obtained from (39) for four different values of the slip length, . The results indicate that the axial velocity increases significantly when the slip length increases. Figure 5 shows the effects of the slip length and ratio on the flow rate . It is found that the flow rate is linear in slip length. For the same size of the cross-sectional area with the variation of the ratio , the flow rate increases significantly as the ratio increases.

(a)

(b)

(c)

(d)

Figure 3: Axial velocity profile on the cross-section of the channel with the area of and ratio of 3/4 for constant pressure gradient under various values of slip lengths : (a) ; (b) ; (c) ; (d) .

(a)

(b)

Figure 4: 2D velocity profile along the -axis and -axis on the cross-section of the channel with the area of and ratio of 3/4 for constant pressure gradient under different slip lengths values.

(a)

(b)

Figure 5: Influence of slip length and the ratio on the flow rate on the cross-section of the channel with the constant area of for a constant gradient pressure.

For Case 2, as the pressure gradient depends on time, the normalized velocity and flow rate as well as shear stress vary with time and the slip length. We investigate the axial velocity on a cross-section of the channel having the ratio of 1 and the slip length of 0.001 under the frequency with at various instants of time. Figures 6 and 7 show the transient velocity obtained from (43) for the slip length . For , the axial velocity has similar profile with , as is constant on the velocity equation (43) so that the curves coincide. The influences of the slip length, the frequency , and the ratio of the cross-section of the channel on the flow behavior are illustrated by analyzing the solutions graphically. Figure 8 shows the transient flow rate on a cross-section of the channel for four different values of the slip length, . Figure 9 presents the influence of on the amplitude of the flow rate for various frequencies with four different values: . The result indicates that the dependence of on is different for different . At high frequency (high ), the amplitude of the flow rate increases initially as increases but it then tends to a constant value once the slip length becomes sufficiently large. Figure 10 shows the influence of the ratio under the same cross-sectional area size on the flow rate for different slip lengths . It is noted that as the ratio increases the flow rate increases but tends to a constant value as becomes sufficiently large.

(a)

(b) .25

(c) .5

(d)

(e)

Figure 6: Axial velocity profile on the cross-section of the channel with the area of and ratio under the frequency with and slip length of 0.001 for various instants of time.

(a)

(b)

Figure 7: Axial velocity along the -axis and -axis on the cross-section of the channel with the area of and ratio = 1 under the frequency with and slip length of 0.001 at various instants of time.

Figure 8: Influence of slip length on the flow rate on the cross-section of the channel with the area of and ratio of 1 under the frequency with at various instants of time.

Figure 9: Influence of slip length on the amplitude of the flow rate on the cross-section of the channel with the area of and ratio of 1 under different frequencies with different .

Figure 10: Influence of ratio of and slip length on the amplitude of the flow rate on the cross-section of the channel with constant area of under the frequency with .

5. Conclusions

In this paper, we derive an exact solution for the unsteady flow of an incompressible Newtonian fluid in a rectangular microchannel with a Navier-slip boundary. From the explicit analytical solutions of the velocity and flow rate, we investigate the effect of the slip length and the geometry of the cross-section on the flow of the fluid through the channel. The investigation shows the following.

(1)For the flow through rectangular microchannels with constant pressure gradient, the axial velocity increases faster in the center of the cross-section than in other areas as the slip length increases, while for the flows, due to the waveform pressure gradient, the velocity changes significantly as the slip length increases.(2)For flow driven by a constant pressure gradient, the flow rate is linear with respect to the slip length and the ratio (), and, for different values of with constant (), the flow rate increases when increases and also when the slip length increases.(3)The amplitudes of flow rate initially increase significantly as slip length increases but tend to a constant value when becomes sufficiently large for various values. This profile is similar to the case shown in [39].

Conflict of Interests

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

Acknowledgments

The authors wish to acknowledge the support of the Faculty of Science, Mahidol University. The second author gratefully acknowledges the support of the Australian Research Council through a discovery project grant. The authors are also grateful to both referees for their comments and suggestions which have led to significant improvement of this paper.

References

R. B. Bird, W. E. Stewart, and E. N. Lightfoot, Transport Phenomena, John Wiley & Sons, New York, NY, USA, 1960.

P. Huang and K. S. Breuer, “Direct measurement of slip length in electrolyte solutions,” Physics of Fluids, vol. 19, no. 2, Article ID 028104, 2007.

D. You and P. Moin, “Effects of hydrophobic surfaces on the drag and lift of a circular cylinder,” Physics of Fluids, vol. 19, no. 8, Article ID 081701, 2007.

M. Gad-el-Hak, “The fluid mechanics of microdevices—the Freeman scholar lecture,” Journal of Fluids Engineering, vol. 121, no. 1, pp. 5–33, 1999.

C. Ho and Y. Tai, “Micro-electro-mechanical-systems (MEMS) and fluid flows,” Annual Review of Fluid Mechanics, vol. 30, pp. 579–612, 1998.

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 & Fluid Flow, vol. 17, no. 6, pp. 587–607, 2007.

F. Saidi, “Non-Newtonian flow in a thin film with boundary conditions of Coulomb's type,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 86, no. 9, pp. 702–721, 2006.

S. J. Liao, Beyond Perturbation: Introduction to Homotopy Analysis Method, Chapman & Hall/CRC Press, Boca Raton, Fla, USA, 2003.

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.

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.

T. C. Papanastasiou, “Flows of materials with yield,” Journal of Rheology, vol. 31, no. 5, pp. 385–404, 1987.

J. P. Pascal, “Instability of power-law fluid flow down a porous incline,” Journal of Non-Newtonian Fluid Mechanics, vol. 133, no. 2-3, pp. 109–120, 2006.

R. Pit, H. Hervet, and L. Léger, “Direct experimental evidence of slip in hexadecane: solid interfaces,” Physical Review Letters, vol. 85, no. 5, pp. 980–983, 2000.

P. A. Thompson and S. M. Troian, “A general boundary condition for liquid flow at solid surfaces,” Nature, vol. 389, no. 6649, pp. 360–362, 1997.

G. Tunc and Y. Bayazitoglu, “Heat transfer in rectangular microchannels,” International Journal of Heat and Mass Transfer, vol. 45, no. 4, pp. 765–773, 2002.

C. Y. Wang, “Exact solutions of the unsteady Navier-Stokes equations,” Applied Mechanics Reviews, vol. 42, no. 11, pp. S269–S282, 1989.

B. Wiwatanapataphee, Y. H. Wu, M. Hu, and K. Chayantrakon, “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.

A. M. Afonso, L. L. Ferrás, J. M. Nóbrega, M. A. Alves, and F. T. Pinho, “Pressuredriven electrokinetic slip flows of viscoelastic fluids in hydrophobic microchannels,” Microfluidics and Nanofluidics, vol. 26, pp. 1–12, 2013.

S. R. Deshmukh and D. G. Vlachos, “CFD simulations of coupled, countercurrent combustor/reformer microdevices for hydrogen production,” Industrial and Engineering Chemistry Research, vol. 44, no. 14, pp. 4982–4992, 2005.

Z. Duan and B. He, “Discussion: slip-flow in microchannels of non-circular cross sections,” Journal of Fluids Engineering, vol. 136, no. 5, 2014.

M. Fichman and G. Hetsroni, “Viscosity and slip velocity in gas flow in microchannels,” Physics of Fluids, vol. 17, no. 12, Article ID 123102, pp. 1–5, 2005.

A. He, “Interfacial instability of compressible slip flows in a microchannel,” Physical Review E, vol. 87, no. 5, Article ID 053006, 2013.

T. Kuo, D. M. Cannon Jr., M. A. Shannon, P. W. Bohn, and J. V. Sweedler, “Hybrid three-dimensional nanofluidic/microfluidic devices using molecular gates,” Sensors and Actuators, A: Physical, vol. 102, no. 3, pp. 223–233, 2003.

H. B. Lee, I. W. Yeo, and K. K. Lee, “Water flow and slip on NAPL-wetted surfaces of a parallel-walled fracture,” Geophysical Research Letters, vol. 34, no. 19, Article ID L19401, 2007.

M. T. Matthews and J. M. Hill, “Newtonian flow with nonlinear Navier boundary condition,” Acta Mechanica, vol. 191, no. 3-4, pp. 195–217, 2007.

S. Palle and S. Aliabadi, “Slip flow and heat transfer in rectangular and circular microchannels using hybrid {FE}/{FV} method,” International Journal for Numerical Methods in Engineering, vol. 89, no. 1, pp. 53–70, 2012.

L. Szalmas, “Slip-flow boundary condition for straight walls in the lattice Boltzmann model,” Physical Review E, vol. 73, no. 6, Article ID 066710, 2006.

Q. Sun, Y. H. Wu, L. Liu, and B. Wiwatanapataphee, “Study of a Newtonian fluid through circular channels with slip boundary taking into account electrokinetic effect,” Abstract and Applied Analysis, vol. 2013, Article ID 718603, 9 pages, 2013.

F. Tahmouresi and S. K. Das, “Analytical modeling of gaseous slip flow in parabolic microchannels,” Journal of Fluids Engineering, vol. 136, no. 7, pp. 10–1115, 2014.

A. Tamayol and K. Hooman, “Slip-flow in microchannels of non-circular cross sections,” Journal of Fluids Engineering, vol. 133, no. 9, Article ID 091202, 2011.

G. H. Tang, X. F. Li, Y. L. He, and W. Q. Tao, “Electroosmotic flow of non-Newtonian fluid in microchannels,” Journal of Non-Newtonian Fluid Mechanics, vol. 157, no. 1-2, pp. 133–137, 2009.

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.

C. Zhang, Y. Chen, Z. Deng, and M. Shi, “Role of rough surface topography on gas slip flow in microchannels,” Physical Review E, vol. 86, no. 1, Article ID 016319, 2012.

Z. Duan and Y. S. Muzychka, “Slip flow in elliptic microchannels,” International Journal of Thermal Sciences, vol. 46, no. 11, pp. 1104–1111, 2007.

Z. Duan and Y. S. Muzychka, “Slip flow in non-circular microchannels,” Microfluidics and Nanofluidics, vol. 3, no. 4, pp. 473–484, 2007.

S. C. Hunter, Mechanics of Continuous Media, John Wiley & Sons, New York, NY, USA, 1983.

T. R. Salamon, D. E. Bornside, R. C. Armstrong, and R. A. Brown, “Local similarity solutions in the presence of a slip boundary condition,” Physics of Fluids, vol. 9, no. 5, pp. 1235–1247, 1997.

E. B. Arkilic, “Slip flow in microchannel,” in Proceeding of Rarefield Gas Dynamic Symposium, Oxford, UK, July 1994.

M. Bahrami, M. M. Yovanovich, and J. R. Culham, “Pressure drop of fully-developed, laminar flow in microchannel of arbitrary cross-section,” Journal of Fluids Engineering, vol. 128, no. 5, pp. 1036–1044, 2006.

R. Baviere, F. Ayela, S. le Person, and M. Favre-Marinet, “Experimental characterization of water flow through smooth rectangular microchannels,” Physics of Fluids, vol. 17, no. 9, Article ID 098105, 2005.

G. Chakraborty, “A note on methods for analysis of flow through microchannels,” International Journal of Heat and Mass Transfer, vol. 51, no. 17-18, pp. 4583–4588, 2008.

W. A. Ebert and E. M. Sparrow, “Slip flow in rctangular and annular ducts,” Journal of Fluids Engineering, vol. 87, no. 4, pp. 1018–1024, 1965.

V. O'Brien, “Pulsatile fully developed flow in rectangular channels,” Journal of the Franklin Institute, vol. 800, no. 3, pp. 225–230, 1975.

C. Aubert and S. Colin, “High-order boundary conditions for gaseous flows in rectangular microducts,” Microscale Thermophysical Engineering, vol. 5, no. 1, pp. 41–54, 2001.

E. Lauga and T. M. Squires, “Brownian motion near a partial-slip boundary: a local probe of the no-slip condition,” Physics of Fluids, vol. 17, no. 10, Article ID 103102, 2005.

E. Lauga and C. Cossu, “A note on the stability of slip channel flows,” Physics of Fluids, vol. 17, no. 8, pp. 1–4, 2005.

G. L. Morini and M. Spiga, “Slip flow in rectangular microtubes,” Microscale Thermophysical Engineering, vol. 2, no. 4, pp. 273–282, 1998.

M. Spiga and G. L. Morini, “Laminar heat transfer in rectangular ducts,” Trends in Heat, Mass & Momentum Transfer, vol. 3, pp. 19–30, 1997.