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 [5–7], Tan [8–12], 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 [15–18]. 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 [20–23] 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, 28–30], 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).