Abstract
The velocity field and the adequate shear stress corresponding to the first problem of Stokes for generalized Burgers’ fluids are determined in simple forms by means of integral transforms. The solutions that have been obtained, presented as a sum of steady and transient solutions, satisfy all imposed initial and boundary conditions. They can be easily reduced to the similar solutions for Burgers, Oldroyd-B, Maxwell, and second-grade and Newtonian fluids. Furthermore, as a check of our calculi, for small values of the corresponding material parameters, their diagrams are almost identical to those corresponding to the known solutions for Newtonian and Oldroyd-B fluids. Finally, the influence of the rheological parameters on the fluid motions, as well as a comparison between models, is graphically illustrated. The non-Newtonian effects disappear in time, and the required time to reach steady-state is the lowest for Newtonian fluids.
1. Introduction
There is evidence that the interest of the workers in non-Newtonian fluids is on the leading edge during the last few years. Many researchers have the opinion that flows of such fluids are important in industry and technology. Several investigations in the field cite a wide variety of applications in rheological problems in biological sciences, geophysics, and chemical and petroleum industries [1]. It is an established fact that unlike the Newtonian fluids, the flows of non-Newtonian fluids cannot be analyzed by a single constitutive equation. This is due to the rheological properties of non-Newtonian fluids. The understanding of flows of such fluids has progressed via a number of theoretical, computational, and experimental efforts. The resulting equations of such fluids are in general of higher order than the Navier-Stokes equation and one needs additional conditions for a unique solution [2, 3]. Specifically to obtain an analytic solution for such flows is not an easy task. In spite of several challenges, many investigations regarding the analytic solutions for flows of non-Newtonian fluids have been performed [4–19].
Many models are accorded to describe the rheological behavior of non-Newtonian fluids [20, 21]. They are usually classified as fluids of differential, rate and integral type. Amongst the non-Newtonian fluids, the rate-type fluids are those which take into account the elastic and memory effects. The simplest subclasses of rate type fluids are those of Maxwell and Oldroyd-B fluids. But these fluid models do not exhibit rheological properties of many real fluids such as asphalt in geomechanics and cheese in food products. Recently, a thermodynamic framework has been put into place to develop the one-dimensional rate type model known as Burgers’ model [22] to the frame-indifferent three-dimensional form by Murali Krishnan and Rajagopal [23]. This model has been successfully used to describe the motion of the earth’s mantle. The Burgers’ model is the preferred model to describe the response of asphalt and asphalt concrete [24]. This model is mostly used to model other geological structures, such as Olivine rocks [25] and the propagation of seismic waves in the interior of the earth [26]. In the literature, the vast majority of the flows of the rate-type models has been discussed using Maxwell and Oldroyd-B models. However, the Burgers’ model has not received much attention in spite of its diverse applications. We here mention some of the studies [27–33] made by using Burgers’ model.
The purpose of this work is to established exact solutions corresponding to the first problem of Stokes for generalized Burgers’ fluids. Actually, we determine the velocity and the adequate shear stress corresponding to the motion of such a fluid over a plane wall, which initially is at rest and is suddenly moved in its own plane withe a constant velocity. The general solutions, obtained by means of Fourier sine and Laplace transforms, are presented under integral form in terms of the elementary functions and can be reduced to the similar solutions for Burgers fluids. As a check of their correctness, we also showed that for small values of the rheological parameters , , , and or and only, the diagrams of the general solutions are almost identical to those corresponding to the known solutions for Newtonian, respectively, Oldroyd-B fluids. The influence of the material parameters on the fluid motion, as well as a comparison between some models, is also underlined by graphical illustrations. The non-Newtonian effects disappear in time and the Newtonian fluid flows faster.
2. Basic Governing Equations
The Cauchy stress tensor for an incompressible generalized Burgers’ fluid is characterized by the following constitutive equations [30–33]: where denotes the indeterminate spherical stress, is the extra-stress tensor, is the first Rivlin-Ericksen tensor ( being the velocity gradient), is the dynamic viscosity, and (<) are relaxation and retardation times, and are new material parameters of the generalized Burgers’ fluid (having the dimension of ), and denotes the upper convected derivative defined in [30–33].
This model includes as special cases the Burgers’ model , Oldroyd-B model , Maxwell model , and the Newtonian fluid model when . In some special flows, like those to be here considered, the governing equations corresponding to generalized Burgers’ fluids resemble those for second-grade fluids.
For the problem under consideration, we assume a velocity field and an extra-stress tensor of the form where i is the unit vector along the -coordinate direction. For these flows, the constraint of incompressibility is automatically satisfied. If the fluid is at rest up to the moment , then Equations (2.1) and (2.3) imply , and the meaningful equation where is the nonzero shear stress. In the absence of body forces, the balance of linear momentum reduces to Eliminating between (2.4) and (2.5) and assuming that there is no pressure gradient in the flow direction, we find the governing equation under the form
Consider an incompressible generalized Burgers’ fluid occupying the space above a flat plate perpendicular to the -axis. Initially, the fluid is at rest and at the moment the plate is impulsively brought to the constant velocity in its plane. Due to the shear, the fluid above the plate is gradually moved. Its velocity is of the form (2.3)1, while the governing equations are given by (2.6) and (2.4). The relevant problem under initial and boundary conditions [34–36] is where is the Heaviside function. Moreover, the natural conditions have to be also satisfied. They are consequences of the fact that the fluid is at rest at infinity and there is no shear in the free stream.
3. Solution of the Problem
In order to determine the exact solution, we shall use the Fourier sine transforms [37]. Multiplying both sides of (2.6) by , integrating the result with respect to from 0 to infinity, and taking into account the boundary conditions (2.8) and (2.9), we find that where and are delta function and its derivative and the Fourier sine transform of defined by has to satisfy the initial conditions By applying the Laplace transform to (3.1) and having in mind the initial conditions (2.7), we find that Now, for a more suitable presentation of the final results, we rewrite (3.4) in the following equivalent form: Inverting (3.5) by means of the Fourier sine formula, we can write as Finally, in order to obtain the velocity field , we apply the inverse Laplace transform to (3.6) and use (A.1) from the Appendix. As a result, we find for the velocity field, the following simple expression: where are the roots of the algebraic equation . In above relations (see the Cardano’s formulae [38]), are the roots of the algebraic equation , where From Routh-Hurwitz’s principle [39], we get if , provided . The corresponding shear stress (see also (A.2)) is obtained in the same way from (2.4).
4. Special Cases
4.1. Burgers’ Fluid
Making into (3.7) and (3.11), we obtain the velocity field and the associated shear stress corresponding to a Burgers’ fluid performing the same motion.
4.2. Oldroyd-B Fluid
Making and into (3.6) and following the same way as before, we get the velocity field (see also (A.3)) and the shear stress corresponding to an Oldroyd-B fluid. In the above relations, and (4.1) is identical to (15) from [40].
4.3. Maxwell Fluid
Making the limit of (4.1) and (4.2) as , we obtain the solutions corresponding to a Maxwell fluid. The new roots and are given by
4.4. Second-Grade Fluid
It is worthwhile pointing out that the similar solutions for second-grade fluids can be also obtained as limiting case of our solutions. Indeed, if we do not take into consideration the restriction and make into (4.1) and (4.2), we recover the expressions corresponding to a second-grade fluid. The solution (4.6) is identical to that from [40, equation (16)] or [36, equation (14)].
4.5. Newtonian Fluid
Finally, making into (4.4) or into (4.6) and (4.7), the solutions for a Newtonian fluid are achieved. The above equations for and can be written under classical forms corresponding to the first problem of Stokes.
5. Numerical Results and Discussion
In order to reveal some relevant physical aspects of the obtained results, several graphs are sketched in this section. A series of diagrams of the velocity and the shear stress against were performed for different situations with typical values. For example, we chose , , and for simplicity, and different values for , , , and were chosen to illustrate their effects on the fluid motion. From Figure 1, it is clear that the velocity is an increasing function with respect to , while the shear stress in absolute value decreases with regard to . Both are decreasing functions with respect to . Figure 2 shows the variations of the two physical entities with respect to the kinematic viscosity . As it was to be expected, both the velocity and the shear stress (of course, in absolute value) are increasing functions with respect to .
(a)
(b)
(a)
(b)
The influence of the relaxation and retardation times and on the fluid motion is underlined by Figures 3 and 4. Their effects, as expected, are opposite. More exactly, both the velocity and the shear stress are decreasing functions with respect to and increasing ones with regard to . Figures 5 and 6 show the influence of the other two material parameters on the fluid motion. From these figures, it is clear that and have opposite effects upon velocity on the whole domain and shear stress on a part only. More exactly, the velocity of the fluid is everywhere a decreasing function with respect to and an increasing one with regard to exist. The shear stress is a decreasing function of on the whole domain and of near the plate. The effects of and on the fluid motion and of and upon velocity are qualitatively the same.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
Finally, for comparison, the profiles of the velocity and the shear stress corresponding to three models (Newtonian, Oldroyd-B, and generalized Burgers’) are together depicted in Figure 7 for the same values of and the common material constants. It is clearly seen from these figures that the Newtonian fluid is the swiftest and the generalized Burgers’ fluid is the slowest. Furthermore, the non-Newtonian effects disappear in time and the behavior of Oldroyd-B and generalized Burgers’ fluids, as it results from Figure 8, can be well enough approximated by that of the Newtonian fluids. From the expressions (4.6) and (4.8) of the velocity field , it results that the required time to reach the steady state is lower for Newtonian fluids in comparison with second-grade fluids. A comparison with other types of fluids has been also realized by graphical illustrations.
(a)
(b)
(a)
(b)
6. Concluding Remarks
In this paper, the velocity field and the adequate shear stress corresponding to the first problem of Stokes for generalized Burgers’ fluids are determined using the Fourier sine and Laplace transforms. The solutions that have been obtained are presented under integral form in terms of the elementary functions , , and and satisfy all imposed initial and boundary conditions. They are written as a sum of steady and transient solutions and can be easily reduced to give the similar solutions for Burgers’ fluids. The steady solutions are the same for both types of fluids if the conditions and are satisfied. Furthermore, they are also identical to the steady solutions corresponding to Oldroyd-B, Maxwell, and second-grade and Newtonian fluids. The required time to reach the steady-state can be easily determined by graphical illustrations. It depends of the material constants and differs from a fluid to another one.
The general solutions (3.7) and (3.11) presented in the simplest forms, and their correctness has been graphically verified by comparison with the known solutions for Oldroyd-B and Newtonian fluids. More exactly, from Figures 9 and 10, it clearly results that for small values of the material constants and or , , , and , as expected, the diagrams of these solutions are almost identical to those corresponding to Oldroyd-B and Newtonian fluids, respectively. Finally, in order to bring light on some relevant physical aspects of the obtained results, the influence of the material parameters on the fluid motion is underlined by graphical illustrations. A comparison between the Newtonian, Oldrotd-B, and generalized Burgers’ fluid is also realized. The main outcomes of this study are as follows.(i)The general solutions (3.7) and (3.11) are presented under simple forms as a sum of steady and transient solutions. They have been immediately particularized to give the similar solutions for Burgers’ fluids.(ii)As a check of our calculi, we showed that for small values of the material constants and or , , , and the diagrams of these solutions as it was to be expected, are almost identical to those corresponding to Oldroyd-B and Newtonian fluids.(iii)The velocity and the shear stress (in absolute value) are increasing functions with respect to .(iv)The relaxation and retardation times, and , as expected, have opposite effects on the fluid motion. Both the velocity and the shear stress (in absolute value) are decreasing functions with respect to and increasing ones with regard to .(v)The other two material constants and have opposite effects on the velocity on the whole flow domain. Their effect on the shear stress is qualitatively the same near the plate and different in rest. Roughly speaking, the effects of and on the fluid velocity are qualitatively the same as those of and .(vi)The Newtonian fluid is the swiftest, and the generalized Burgers’ fluid is the slowest. The non-Newtonian effects disappear in time, and the required time to reach the steady-state is the lowest for Newtonian fluid.
(a)
(b)
(a)
(b)
Appendix
Acknowledgments
The author Muhammad Jamil is highly thankful and grateful to the Abdus Salam School of Mathematical Sciences, GC University, Lahore, Pakistan; Department of Mathematics, NED University of Engineering and Technology, Karachi-75270, Pakistan; also Higher Education Commission of Pakistan for generous support, facilitating this research work.