`Mathematical Problems in EngineeringVolume 2008, Article ID 641431, 10 pageshttp://dx.doi.org/10.1155/2008/641431`
Research Article

## Exact Solutions of Rayleigh-Stokes Problem for Heated Generalized Maxwell Fluid in a Porous Half-Space

1Department of Fundamental Sciences, Yancheng Institute of Technology, Yancheng, Jiangsu 224003, China
2Department of Mathematics, School of Science, Beijing Jiaotong University, Beijing 100044, China

Received 8 October 2007; Accepted 24 April 2008

Academic Editor: Katica Hedrih

Copyright © 2008 Changfeng Xue and Junxiang Nie. 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

The Rayleigh-Stokes problem for a generalized Maxwell fluid in a porous half-space with a heated flat plate is investigated. For the description of such a viscoelastic fluid, a fractional calculus approach in the constitutive relationship model is used. By using the Fourier sine transform and the fractional Laplace transform, exact solutions of the velocity and the temperature are obtained. Some classical results can be regarded as particular cases of our results, such as the classical solutions of the first problem of Stokes for Newtonian viscous fluids, Maxwell fluids, and Maxwell fluids in a porous half-space.

#### 1. Introduction

Many materials such as clay coatings, drilling mud, suspensions, certain oils and greases, polymer melts, elastomers, and many emulsions have been treated as non-Newtonian fluids. It is difficult to suggest a single model which can exhibit all properties of non-Newtonian fluids as is done for Newtonian fluids. They cannot be described in a simple model as for the Newtonian fluids, and there has been much confusion over the classification of the non-Newtonian fluids. For this reason, many models of constitutive equations have been proposed. Recently, fractional calculus has encountered much success in the description of constitutive relations of viscoelastic fluids. The starting point of the fractional derivative model of viscoelastic fluids is usually a classical differential equation which is modified by replacing the time derivative of an integer order with the so-called Riemann-Liouville fractional calculus operator. This generalization allows one to define precisely noninteger order integrals or derivatives. Bagley and Torvik [1], Friedrich [2], Huang et al. [3], He et al. [4], Xu and Tan [57], Tan [812], and Shen et al. [13, 14] have sequentially introduced the fractional calculus approach into various rheology problems. Fractional derivatives have been found to be quite flexible in describing viscoelastic behaviors.

The first problem of Stokes for the flat plate, like the Rayleigh-Stokes problem for an edge, has received much attention because of its practical importance [1518]. This unsteady flow problem examines the diffusion of vorticity in a half-space filled with a viscous incompressible fluid that is set to motion when an infinite flat plate suddenly assumes a constant velocity parallel to itself from the rest. By using similarity transformation of variables, the exact solution corresponding to a Newtonian fluid was obtained in an elegant form by Stokes. But for a Maxwell fluid, it is difficult to obtain such a strict similarity solution. Further, the order of the equation of motion for such a fluid is higher than that of the Navier-Stokes equation and thus, in general, one needs conditions in addition to the usual adherence boundary condition. Friedrich first investigated this problem and gave a few exact solutions [19].

On the other hand, the problems of the fluid through porous media play an important role in various practical applications. But very little efforts have so far been made to discuss these problems. Recently, based on the local volume averaging technique and the balance of forces acting on a volume element of viscoelastic fluids in porous media, Tan and Masuoka [2023] developed a modified Darcy-Brinkman model for viscoelastic fluid flows in porous media. Based on their works, in this paper, we extend Stokes' first problem to that of a heated generalized Maxwell fluid with fractional derivative in a porous half-space. Exact solutions of the velocity and the temperature are obtained by using the Fourier sine transform and the fractional Laplace transform. Some previous classical results can be recovered from our results, such as the classical solutions of the first problem of Stokes for Newtonian viscous fluids and those for Maxwell fluids.

#### 2. Basic Equations

The well-known Darcy law states that, in the flow of a Newtonian fluid through a porous medium, the pressure gradient caused by the friction drag is directly proportional to the velocity, namely,where is the Darcian velocity, is the dynamic viscosity of the fluid, and is the permeability of the porous medium, respectively. The Darcy law, however, does not consider the boundary effect on the flow, which unfortunately cannot be ignored in many applications. In 1947, Brinkman proposed an equation to describe the locally averaged flow in a porous medium [24]. In the Brinkman equation, the pressure gradient, the divergence of the viscous stress tensor, and the friction force exerted by the porous medium satisfywhere is the effective dynamic viscosity of the medium. Although the Brinkman equation takes the boundary effect on the flow into account, it is justified for and can thus be applied to only steady-state flows of a Newtonian fluid in porous media. For viscoelastic fluid flows in porous media, on the other hand, only a few mathematical macroscopic filtration models have been proposed. In particular, by analogy with the constitutive equation of a Maxwell fluid [25, 26], the following phenomenological model has been introduced: where is the relaxation time of the generalized Maxwell fluid, is fractional coefficient, and is Riemann-Liouville fractional derivative, which is defined as in [27]where is the gamma function.

Let us now follow the same procedure as used in [22] to investigate the flow of a generalized Maxwell fluid in a porous half-space. With the positive -axis of the rectangular coordinate system pointing up, a generalized Maxwell fluid flows through the porous half-space above and in contact with a flat plate occupying the -plane. At time , both the fluid and the plate are at rest. At time , the flat plate suddenly starts to slide in its plane with a constant speed . Under these conditions, no flow occurs in the - and the -directions, and the flow velocity at a given point in the porous half-space depends only on its -coordinate and time , that is, the intrinsic velocity field takes the formwhere represents the unit vector in the -direction. The intrinsic velocity V, averaged over the pore space, is related to the Darcian velocity by , where is the porosity of the porous medium. In this work, employing a local volume averaging technique [20, 2830], we will consider a modified, more rigorous Darcy law for Stokes' first problem with a generalized Maxwell fluid in the porous half-space. In fact, in terms of the pressure gradient, the local volume-averaged balance of linear momentum can be given bywhere is the fluid density, is the viscous stress tensor, and is the Darcy resistance, namely, the viscous damping force caused by the micropore structure of a porous medium, respectively. Under the assumption thatas in [23, 31], next we discuss the viscoelastic properties of the fluid using the upper convected derivative generalized Maxwell fluid.

It should be pointed out that the classical viscous Newtonian fluid is just a special case of this model at . Inserting (2.5) into (2.7) and using the initial condition , we get and

Note that the pressure gradient given in (2.3) can be taken as a measure of the same resistance to flow in the bulk of the porous medium [26]. Since the Darcy resistance in (2.6) is also in effect a measure of the flow resistance offered by the solid matrix, from (2.3) we get

Using (2.9) and (2.10) to eliminate and in (2.6) and recalling that , we have the momentum equation of Stokes' first problem for a generalized Maxwell fluid in the porous half-space

Ignoring the pressure gradient in the -direction finally leads to the equation of motion as follows:

In addition, the initial and the boundary conditions arerespectively, and the natural condition is

We will assume that there is a local thermal equilibrium. Then, , where and are the temperatures of the solid and the fluid phases in porous media, respectively. The thermal energy equation is then [32, 33]where is the specific heat capacity of the system formed by both the fluid and the solid phases. More precisely, , where is the density of the solid phase, is the specific heat capacity of the solid phase, and is the specific heat capacity of the fluid at a constant pressure, respectively. is the equivalent thermal conductivity of the fluid in the porous half-space, and is the contribution due to viscous dissipation [34, 35]

Equation (2.16) is a unified form compatible with an expression derived from the Navier-Stokes equation for a fluid clear of solid [36, 37].

The corresponding initial and boundary conditions are

Moreover, the natural conditions are

#### 3. Exact Solutions

Employing the nondimensional quantitiesthe dimensionless motion and energy equations can be written as (for brevity the dimensionless marks “” are omitted here)where , , , , . is the kinematic viscosity. The dimensionless viscous dissipation, namely, , has the following form:

##### 3.1. Solution of the Velocity Field

Multiplying both sides of (3.2) and (3.3) by integrating with respect to from to , and then using the boundary conditions of (3.4) and (3.5), in fact we havewhere denotes the Fourier sine transform of with respect to . In order to obtain an exact solution of (3.11) subject to the initial condition of (3.12), the fractional Laplace transform is used. Letbe the image function of , where is the transform parameter. Using Laplace transform principle of sequential fractional derivatives [27], we get

In order to avoid the burdensome calculations of residues and contour integrals, we apply the discrete inverse Laplace transform method here. First, we rewrite (3.14) in a series form

Then, applying the inversion formulae term by term for the Laplace transform, (3.15) becomesin whichdenotes the generalized Mittag-Leffler function. Here, we use the following property of the generalized Mittag-Leffler function's inverse Laplace transform [27]

Inverting (3.18) by the Fourier sine inverse transform, we obtain an exact solution to the velocity field as follows:

In a special case, when , corresponding to Stokes' first problem for a Maxwell fluid in a porous half-space, (3.14) can be written as

Correspondingly, (3.19) can be written aswhich is just the dimensionless velocity solution of a Maxwell fluid in a porous half-space obtained in another form by Tan [22].

In a special case, when and , we havewhich is just the dimensionless velocity solution of a clear Maxwell fluid obtained in another form by Jordan et al. [38].

In another special case, when , and , corresponding to the classical version of Stokes' first problem for a Newtonian fluid, (3.14) can be written as

Correspondingly, (3.19) can be written as

Equation (3.24) can be simplified further as follows:where is the error function of Gauss. This is just the well-known result of Stokes' first problem for a clear Newtonian fluid.

##### 3.2. Solution of Temperature Field

In order to determine the temperature distribution in the porous half-space, we multiply both sides of (3.6) and (3.7) by and then integrate with respect to from to . Using (3.8) and (3.9), we getwhere and denote the Fourier sine transforms of and with respect to , respectively. Since the solution of velocity is obtained by (3.21), is a known function. The solution of the ordinary differential equation (3.26) subject to the initial condition (3.27) can be easily obtained as

Inverting (3.28) by the Fourier sine inverse transform, and integrating by parts the term containing , we get

In a special case, when , the fluid is at rest, and the temperature distribution is the same for both a Maxwell fluid and a Newtonian one. Then, (3.29) takes the form

From (3.30), we easily obtain that

If, moreover, the plate is kept at a constant temperature, (3.31) reduces to a well-known form,

#### 4. Conclusion

In this work, we have presented some results about the generalized Maxwell fluid in a porous half-space on a heated flat plate. An exact solution of the velocity field was obtained by using the Fourier sine transform and the fractional Laplace transform. The temperature distribution in the generalized Maxwell fluid in the porous half-space subject to a linear flow on a heated flat plate was determined. Some classical results can be considered as particular cases of our results, such as the solutions corresponding to a Newtonian viscous fluid, a clear Maxwell fluid, and a Maxwell fluid in a porous half-space. They can be easily obtained from the results presented in this paper by letting , , ; , and , respectively. It has also been shown that the fractional constitutive relationship model is more flexible than the conventional model in describing the properties of viscoelastic fluids.

#### Acknowledgments

This work was supported by National Natural Science Foundation of China (Grants no. 10572006, 10632010, 40628001) and National Key Basic Research Program (2006CB705803).

#### References

1. R. L. Bagley and P. J. Torvik, “A theoretical basis for the application of fractional calculus to viscoelasticity,” Journal of Rheology, vol. 27, no. 3, pp. 201–210, 1983.
2. C. Friedrich, “Relaxation and retardation functions of the Maxwell model with fractional derivatives,” Rheologica Acta, vol. 30, no. 2, pp. 151–158, 1991.
3. J. Huang, G. He, and C. Liu, “Analysis of general second-order fluid flow in double cylinder rheometer,” Science in China Series A, vol. 40, no. 2, pp. 183–190, 1997.
4. G. He, J. Huang, and C. Liu, “General second order fluid flow in a pipe,” Applied Mathematics and Mechanics, vol. 16, no. 9, pp. 825–831, 1995.
5. M.-Y. Xu and W.-C. Tan, “Theoretical analysis of the velocity field, stress field and vortex sheet of generalized second order fluid with fractional anomalous diffusion,” Science in China Series A, vol. 44, no. 11, pp. 1387–1399, 2001.
6. M.-Y. Xu and W.-C. Tan, “The representation of the constitutive equation of viscoelastic materials by the generalized fractional element networks and its generalized solutions,” Science in China Series A, vol. 32, no. 8, pp. 673–681, 2002.
7. M.-Y. Xu and W.-C. Tan, “Intermediate processes and critical phenomena: theory, method and progress of fractional operators and their applications to modern mechanics,” Science in China Series G, vol. 49, no. 3, pp. 257–272, 2006.
8. W.-C. Tan and M.-Y. Xu, “The impulsive motion of flat plate in a general second grade fluid,” Mechanics Research Communication, vol. 29, no. 1, pp. 3–9, 2002.
9. W.-C. Tan and M.-Y. Xu, “Plane surface suddenly set in motion in a viscoelastic fluid with fractional Maxwell model,” Acta Mechanica Sinica, vol. 18, no. 4, pp. 342–349, 2002.
10. W.-C. Tan, F. Xian, and L. Wei, “An exact solution of unsteady Couette flow of generalized second grade fluid,” Chinese Science Bulletin, vol. 47, no. 21, pp. 1783–1785, 2002.
11. W.-C. Tan, W. X. Pan, and M.-Y. Xu, “A note on unsteady flows of a viscoelastic fluid with the fractional Maxwell model between two parallel plates,” International Journal of Non-Linear Mechanics, vol. 38, no. 5, pp. 645–650, 2003.
12. W.-C. Tan and M.-Y. Xu, “Unsteady flows of a generalized second grade fluid with the fractional derivative model between two parallel plates,” Acta Mechanica Sinica, vol. 20, no. 5, pp. 471–476, 2004.
13. F. Shen, W.-C. Tan, Y.-H. Zhao, and T. Masuoka, “Decay of vortex velocity and diffusion of temperature in a generalized second grade fluid,” Applied Mathematics and Mechanics, vol. 25, no. 10, pp. 1151–1159, 2004.
14. F. Shen, W.-C. Tan, Y. Zhao, and T. Masuoka, “The Rayleigh-Stokes problem for a heated generalized second grade fluid with fractional derivative model,” Nonlinear Analysis: Real World Applications, vol. 7, no. 5, pp. 1072–1080, 2006.
15. I. Teipel, “The impulsive motion of a flat plate in a viscoelastic fluid,” Acta Mechanica, vol. 39, no. 3-4, pp. 277–279, 1981.
16. K. R. Rajagopal, “A note on unsteady unidirectional flows of a non-Newtonian fluid,” International Journal of Non-Linear Mechanics, vol. 17, no. 5-6, pp. 369–373, 1982.
17. K. R. Rajagopal and A. S. Gupta, “On a class of exact solutions to the equations of motion of a second grade fluid,” International Journal of Engineering Science, vol. 19, no. 7, pp. 1009–1014, 1981.
18. K. R. Rajagopal and A. S. Gupta, “An exact solution for the flow of a non-Newtonian fluid past an infinite porous plate,” Meccanica, vol. 19, no. 2, pp. 158–160, 1984.
19. C. Friedrich, “Relaxation and retardation functions of the Maxwell model with fractional derivatives,” Rheologica Acta, vol. 30, no. 2, pp. 151–158, 1991.
20. W.-C. Tan and T. Masuoka, “Stokes' first problem for a second grade fluid in a porous half-space with heated boundary,” International Journal of Non-Linear Mechanics, vol. 40, no. 4, pp. 515–522, 2005.
21. W.-C. Tan and T. Masuoka, “Stokes' first problem for an Oldroyd-B fluid in a porous half space,” Physics of Fluids, vol. 17, no. 2, Article ID 023101, 7 pages, 2005.
22. W.-C. Tan, “Velocity overshoot of start-up flow for a Maxwell fluid in a porous half-space,” Chinese Physics, vol. 15, no. 11, pp. 2644–2650, 2006.
23. W.-C. Tan and T. Masuoka, “Stability analysis of a Maxwell fluid in a porous medium heated from below,” Physics Letters A, vol. 360, no. 3, pp. 454–460, 2007.
24. F. Capuani, D. Frenkel, and C. P. Lowe, “Velocity fluctuations and dispersion in a simple porous medium,” Physical Review E, vol. 67, no. 52, Article ID 056306, 8 pages, 2003.
25. G. M. Alisaev and A. K. Mirzadjanzade, “For the calculation of delay phenomenon in filtration theory,” Izvestiya Vysshikh Uchebnykh Zavedeniy. Neft' i Gaz, vol. 6, pp. 71–78, 1975.
26. B. Khuzhayorov, J.-L. Auriault, and P. Royer, “Derivation of macroscopic filtration law for transient linear viscoelastic fluid flow in porous media,” International Journal of Engineering Science, vol. 38, no. 5, pp. 487–504, 2000.
27. I. Podlubny, Fractional Differential Equations, Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
28. K. Vafai and C. L. Tien, “Boundary and inertia effects on flow and heat transfer in porous media,” International Journal of Heat and Mass Transfer, vol. 24, no. 2, pp. 195–203, 1981.
29. T. Masuoka and Y. Takatsu, “Turbulence characteristics in porous media,” in Transport Phenomena in Porous Media Vol. II, D. B. Ingham and I. Pop, Eds., pp. 231–256, Pergamon Press, Oxford, UK, 2002.
30. D. B. Ingham and I. Pop, Eds., Transport Phenomena in Porous Media II, Pergamon Press, Oxford, UK, 2002.
31. P. N. Kaloni and J. X. Lou, “Stability of Hadley circulations in a Maxwell fluid,” Journal of Non-Newtonian Fluid Mechanics, vol. 103, no. 2-3, pp. 167–186, 2002.
32. A. Younes, “On modelling the multidimensional coupled fluid flow and heat or mass transport in porous media,” International Journal of Heat and Mass Transfer, vol. 46, no. 2, pp. 367–379, 2003.
33. D. A. Nield and B. Adrian, Convection in Porous Media, Springer, Berlin, Germany, 2nd edition, 1999.
34. A. K. Al-Hadhrami, L. Elliott, and D. B. Ingham, “Combined free and forced convection in vertical channels of porous media,” Transport in Porous Media, vol. 49, no. 3, pp. 265–289, 2002.
35. A. K. Al-Hadhrami, L. Elliott, and D. B. Ingham, “A new model for viscous dissipation in porous media across a range of permeability values,” Transport in Porous Media, vol. 53, no. 1, pp. 117–122, 2003.
36. E. R. G. Eckert and R. M. Drake, Analysis of Heat and Mass Transfer, McGraw-Hill, New York, 1977.
37. D. A. Nield, A. V. Kuznetsov, and M. Xiong, “Thermally developing forced convection in a porous medium: parallel plate channel with walls at uniform temperature, with axial conduction and viscous dissipation effects,” International Journal of Heat and Mass Transfer, vol. 46, no. 4, pp. 643–651, 2003.
38. P. M. Jordan, A. Puri, and G. Boros, “On a new exact solution to Stokes' first problem for Maxwell fluids,” International Journal of Non-Linear Mechanics, vol. 39, no. 8, pp. 1371–1377, 2004.