Abstract
An analysis is made of the unsteady flow of a third-grade fluid in a porous medium. A modified Darcy's law is considered in the flow modelling. Reduction and solutions are obtained by employing similarity and numerical methods. The effects of pertinent parameters on the flow velocity are studied through graphs.
1. Introduction
Recently, it has been recognised in industrial and technological applications that non-Newtonian fluids are more appropriate than viscous fluids. However, there is no model which can alone predict the behaviour of all non-Newtonian fluids. The governing equations of non-Newtonian fluids are of higher order than the Navier-Stokes equations. Therefore, the adhering boundary conditions are not sufficient, and one needs additional boundary conditions for a unique solution. Excellent critical reviews in this direction have been given by Rajagopal [1, 2], Rajagopal et al. [3], and Rajagopal and Kaloni [4]. Amongst the several non-Newtonian fluid models, much attention has been paid to the simplest subclass of viscoelastic fluids known as the second grade. However, this model is not capable of describing the shear thinning and thickening phenomena for steady flow over a rigid boundary. The third-grade fluid model represents a further, although inconclusive, attempt towards a more comprehensive description of the behaviour of viscoelastic fluids. Also, the flows of such fluids in porous medium are quite prevalent in many engineering fields such as enhanced oil recovery, paper and textile coating, and composite manufacturing processes. Some important studies dealing with the flows of non-Newtonian fluids are made by Rajagopal and Na [5, 6], Rajagopal and Gupta [7], Hayat et al. [8–11], Hayat and Ali [12], Ariel et al. [13], Hayat and Kara [14], Abdel-Malek et al. [15], Wafo Soh [16] and Chen et al. [17], and Fetecau and Fetecau [18–20].
Recently, Tan and Masuoka [21] analysed the Stokes' first problem for a second-grade fluid in a porous medium. In another paper, Tan and Masuoka [22] studied the Stokes' first problem for an Oldroyd-B fluid in a porous medium. In these investigations, the authors have used the modified Darcy's law.
The main goal of this paper is to determine analytical solutions for an unsteady flow of a third-grade fluid over a moved plate. The relevant problem is formulated using modified Darcy's law of a third-grade fluid. Two types of analytical solutions are presented and discussed. A numerical solution is also presented.
2. Problem Formulation
Let us introduce a Cartesian coordinate system OXYZ with -axis in the upward direction. The third-grade fluid fills the porous space and is in contact with an infinite moved plate at . For unidirectional flow, the velocity field is where the above definition of velocity automatically satisfies the incompressibility condition. The equation of motion in a porous medium without body forces is where is the fluid density, is the material time differentiation, is the Cauchy stress tensor, and is the Darcy's resistance in a porous space. The Cauchy stress tensor of an incompressible third-grade fluid has the form [23] in which is the pressure, is the identity tensor, () and () are the material constants, and () are the first Rivlin-Ericksen tensors [24] which may be defined through the following equations: In studying fluid dynamics, it is assumed that the flow meets the Clausius-Duhem inequality and that the specific Helmholtz free energy of the fluid is a minimum at equilibrium when [25] On the basis of constitutive equation in an Oldroyd-B fluid, the following expression in a porous medium has been proposed [26]: where and are the relaxation and retardation times, and and are the porosity and permeability of the porous medium, respectively. It should be pointed out that for , (2.6) reduces to the expression which holds for a Maxwell fluid [26] and when , it reduces to that of second-grade fluid [22].
Keeping the analogy of (2.6) with the constitutive equation of an extra stress tensor in an Oldroyd-B fluid, the following expression in the present problem has been suggested: Since the pressure gradient in (2.7) can also be interpreted as a measure of the flow resistance in the bulk of the porous medium, and is the measure of the flow resistance offered by the solid matrix in -direction, then From (2.1) to (2.5) and (2.8), we have The relevant boundary and initial conditions are in which is the reference velocity.
3. Solutions of the Problem
We rewrite (2.9) as where
3.1. Lie Symmetry Analysis
The Lie symmetry analysis reveals that (3.1) admits two sets of symmetry generators depending on the value of . The Appendix provides details of the symmetry analysis of (3.1).
Case 1 (). We obtain a two-dimensional Lie algebra generated by
Case 2 (). We find a three-dimensional Lie algebra generated by
3.2. Travelling Wave Solutions
We now look for invariant solutions under the operator , which represents wave-front-type travelling wave solutions with constant wave speed . The invariant is given by Substituting (3.5) into (3.1) yields a third-order ordinary differential equation for , It can be seen that this equation admits the solution provided that and hence (3.1) subject to (2.10)–(2.12) admits the solution This solution is plotted in Figures 1 and 4 for various values of the emerging parameters.
On the other hand, we find group-invariant solutions corresponding to operators which give meaningful physical solutions of the initial and boundary value problems (2.9) to (2.12). This means and .
3.3. Group-Invariant Solutions Corresponding to
The invariant solution admitted by is the steady-state solution The substitution of (3.10) into (3.1) yields the second-order ordinary differential equation for : subject to boundary conditions where is sufficiently large, and is a constant ( is taken to be a constant). Let then (3.11) transforms to The integration of (3.15) gives where is a constant. Equation (3.16) is equivalent to the following first-order ODE in : One can solve this numerically subject to the boundary conditions (3.12) and (3.13). This solution is plotted in Figure 5.
3.4. Group-Invariant Solutions Corresponding to
The invariant solution admitted by is where as yet is an undetermined function of . Substituting (3.18) into (3.1) yields the linear second-order ordinary differential equation From (2.10) to (2.11), the appropriate boundary conditions for (3.19) are where
We solve (3.19) subject to the boundary conditions (3.20) for positive and obtain
The solutions (3.18) are plotted for positive in Figure 6. This solution is similar to (3.9) except that we do not have a condition like (3.8) here.
3.5. Numerical Solution
We present the numerical solution of (3.1) subject to the initial and boundary conditions: where is an arbitrary function of .
This solution is plotted using Mathematica's solver NDSolve.
4. Results and Discussion
In order to see the variation of various physical parameters on the velocity, Figures 1–7 have been plotted.
The effect of unsteadiness on the velocity profile is shown in Figure 1. This figure depicts that velocity decreases for large values of time. Clearly, the variation of velocity is observed for . For , the velocity profile remains the same. In other words, one can say that steady-state behaviour is achieved for .
The influence of the wave speed on the velocity profile has been presented in Figure 2. It is revealed that velocity decreases by increasing . Moreover, the effects of the fluid parameters and are given in Figures 3 and 4, respectively. These figures depict that and have opposite roles on the velocity. These figures show that velocity increases for large values of whereas it decreases for increasing . In Figure 5, the steady-state solution is plotted, and the velocity profile is the same as observed in the case of travelling wave solution.
Further, the analytical solutions (3.19) for is plotted in Figure 6. Here as indicated in Figure 6, the velocity profile decreases for large values of . Ultimately when , there is almost no variation in velocity.
Finally in Figure 7, we have plotted numerically the velocity profile for small variations of time, and it is observed that the velocity decreases as time increases, which is the the same observation made previously for the analytical solutions.
Appendix
The operator is a generator of Lie point symmetry of (3.1) if where in whichand the total derivative operators areSubstituting the expansion of (A.4) into the symmetry condition (A.2) and separating them by powers of the derivatives of , since , and are independent of the derivatives of , lead to the overdetermined system of linear partial differential equations (note that and are not zero):The solution of this linear system (A.6) gives rise to two cases and . In the former, we obtain and for the second case In both (A.7) and (A.8), the are constants. Setting one of the constants equal to one and the rest of the constants to zero results in the generators given in Section 3.1.