Abstract
The unsteady motion of an OldroydB fluid over an infinite flat plate is studied by means of the Laplace and Fourier transforms. After time t = 0, the plate applies cosine/sine oscillating shear stress to the fluid. The solutions that have been obtained are presented as a sum of steadystate and transient solutions and can be easily reduced to the similar solutions corresponding to Newtonian or Maxwell fluids. They describe the motion of the fluid some time after its initiation. After that time when the transients disappear, the motion is described by the steadystate solutions that are periodic in time and independent of the initial conditions. Finally, the required time to reach the steadystate is established by graphical illustrations. It is lower for cosine oscillations in comparison with sine oscillations of the shear, decreases with respect to ω and λ and increases with regard to λ_{r}.
Mathematical Subject Classification (2010): 76A05; 76A10.
PACS: 47.50.d; 47.85.g.
Keywords:
OldroydB fluid; oscillating shear stress; exact solutions1 Introduction
The laminar flow of a great number of fluids such as polymeric liquids, food products, paints and so forth cannot be adequately described by means of the classical linearly viscous Newtonian model. The interest into motion problems of such fluids, also called nonNewtonian fluids, has considerably grown due to their multiple applications. Among the many models that have been used to describe the behavior of nonNewtonian fluids, the rate type models have received much attention. The first systematic thermodynamic study of such models is that of Rajagopal and Srinivasa [1], within which models for a variety of rate type viscoelastic fluids can be obtained. They showed that the OldroydB fluid is one which stores energy like a linearized elastic solid, its dissipation however being due to two dissipative mechanisms that implies that they arise from a mixture of two viscous fluids. The first exact solutions corresponding to some motions of OldroydB fluids seem to be those of Tanner [2] and Waters and King [3].
Over the past few decades, the unsteady flows of viscoelastic fluids caused by the oscillations of the boundary are of considerable interest. Rajagopal [4] found steadystate solutions for some oscillating motions of second grade fluids and Erdogan [5] provided two starting solutions for the motion of a viscous fluid due to cosine and sine oscillations of a flat plate. Flows of the fluids due to oscillating boundary for different constitutive models can be found into [620]. However, all these articles deal with motion problems in which the velocity is given on the boundary. To the best of our knowledge, the first exact solutions for motions of nonNewtonian fluids due to an infinite plate that applies a shear stress to the fluid seem to be those of Waters and King [21] and Bandelli et al. [22]. Recently [23], interesting solutions have been obtained for the motion of Newtonian fluids induced by an infinite plate that applies oscillating shear stresses to the fluid. This is very important as in some problems, what is specified is the force applied on the boundary. The "no slip" boundary condition may not be necessarily applicable to flows of polymeric fluids that can slip or slide on the boundary. Thus, the shear stress boundary condition is particularly meaningful and a closedform expression for the starting solutions corresponding to the motion induced by a flat plate that applies an oscillating shear to an OldroydB fluid has not been given before.
Consequently, the aim of this article is to determine starting solutions for the unsteady motion of an incompressible OldroydB fluid due to an infinite plate that applies an oscillating shear to the fluid. Such exact solutions, which are not common in the literature, provide an important check for numerical methods that are used to study flows of such fluids in a complex domain. They are presented as a sum of steadystate and transient solutions and satisfy both the governing equations and all imposed initial and boundary conditions. Furthermore, the similar solutions for Maxwell and Newtonian fluids can easily be obtained as limiting cases of general solutions. Finally, the influence of the material parameters on the fluid motion and the required time to reach the steadystate are determined by graphical illustrations. This time is lower for the cosine oscillations in comparison with the sine oscillations of the shear, decreases with respect to the relaxation time λ and the frequency ω of the shear and increases with respect to the retardation time λ_{r}.
2 Governing equations
An incompressible OldroydB fluid is characterized by the following constitutive equations [9,14,17]
where T is the Cauchy stress tensor, pI denotes the indeterminate spherical stress, S is the extrastress tensor, L is the velocity gradient, A = L + L^{T }is the first RivlinEricksen tensor, μ is the dynamic viscosity, λ and λ_{r }are the relaxation and retardation times, the superscript T indicates the transpose operation and the superposed dot denotes the material time derivative. In the following analysis, we will consider a unidirectional flow whose velocity field is given by
where i denotes the unit vector along the xdirection of the Cartesian coordinate system x, y, and z. For such a flow, the constraint of incompressibility is automatically satisfied. We also assume that the extrastress tensor S, as well as the velocity V, depends only on y and t. In the absence of a pressure gradient in the flow direction and neglecting body forces, the governing equation is given by [14]
where v = μ /ρ is the kinematic viscosity and ρ is the constant density of the fluid. The nontrivial shear stress τ(y, t) = S_{xy}(y, t) satisfies the partial differential equation [[24], Equation (4)]
3 Statement of the problem
Let us consider an incompressible OldroydB fluid at rest, over an infinite plate. After time t = 0 the plate applies an oscillating shear to the fluid (f sin ωt or f cos ωt, where f and ω are constants). Owing to the shear, the fluid is gradually moved. Its velocity has the form of Equation (2), the governing equation is given by Equation (3) and the appropriate initial and boundary conditions are
Moreover, the natural condition
also has to be satisfied.
4 Exact solutions
In the following, let us denote by u_{s}(y, t), τ_{s}(y, t) and u_{c}(y, t), τ_{c}(y, t), the solutions corresponding to the two problems and by
the complex velocity and the complex tension, respectively. In view of the above equations, the functions V (y, t) and T (y, t) have to be solutions of the next initial and boundary values problems
4.1 Calculation of the velocity field
In order to determine the solution of initialboundary values problem (9), (11)_{1,2}, (12)_{2 }and (13), we first take the Laplace transform [25] of Equation (9) and obtain
where the Laplace transform of function V (y, t) has to satisfy the conditions
Multiplying Equation (14) , integrating the result with respect to y from 0 to infinity and using Equations (15) and (16), we obtain
where
denotes the Fourier cosine transform [26] of function . Equation (17) can be written as
where
Applying the inverse Laplace transform and then the inverse Fourier cosine transform to Equation (20), we get the following expression
where
Using (A_{1}) and (A_{2}) from Appendix 1, we obtain the following simplified expression
where
Now, for Equation (21), we introduce the function
which can be written in the following equivalent form
where
Applying the inverse Laplace transform and then the inverse Fourier transform to Equation (21) and using Equation (25), we obtain the following expression
Finally, the velocity corresponding to the cosine oscillations of the shear is given by
while that corresponding to sine oscillations has the form
The starting solutions (27) and (28) corresponding to cosine and sine oscillations of the shear, are presented as a sum between the steadystate and transient solutions. They describe the motion of the fluid some time after its initiation. After that time, in which the transients disappear, the starting solutions tend to the steadystate solutions
respectively,
which are periodic in time and independent of the initial conditions. However, they satisfy the governing equation and boundary conditions. As expected, they differ with a phase shift.
4.2 Calculation of shear stress
In order to obtain the corresponding shear stresses, we apply the Laplace transform to Equation (10) and the inverse Fourier cosine transform to Equation (17). Combining the results, we get the following expression for the Laplace transform of the complex tension T (y, t):
Equation (31) can be written as
where
with
Applying the inverse Laplace transform to the Equation (33) and using (A_{3}) from Appendix 1, we get
Similarly, using (A_{4}) and (A_{5}) from Appendix 1, from Equation (34), we obtain the following suitable form for T_{2}(y, t)
A straightforward calculation leads to the following simplified form
where tan ψ = λω.
Now, we consider the function
which can be written in the following equivalent form:
Applying the inverse Laplace transform to Equation (35) and using Equation (40), we obtain
where
Using Equations (32), (36), (38), and (41), the shear stress corresponding to cosine oscillations of the shear can be written in the form
Also, the shear stress corresponding to sine oscillations is given by
Of course, the shear stresses given by Equations (42) and (43) are also presented as a sum of steadystate and transient solutions. The steadystate solutions
also differ by a phase shift. The property seems to be not true for the transient solutions. This is the reason that we separately gave the starting solutions for both cosine and sine oscillations of the shear stress on the boundary.
5 Particular cases
5.1 λ_{r }= 0 (Maxwell fluid)
Making λ_{r }→ 0 into Equations (27) and (28), we obtain the velocity fields
corresponding to a Maxwell fluid performing the same motion. Similarly, from Equations (42) and (43), we obtain the corresponding shear stresses
5.2 λ → 0 (second grade fluid)
Using Equations (27), (28), (42), (43) and the limits
we obtain solutions corresponding to a fluid of second grade:
where
5.3 λ → 0, λ_{r }→ 0 (Newtonian fluid)
Making, λ → 0 and, λ_{r }→ 0 into Equations (27), (28), (42), and (43) or, λ → 0 into Equations (44)(47), or α = 0 into Equations (48)(52), we recover the solutions [[23], Equations (20)(23)]
corresponding to the flow of a Newtonian fluid.
6 Numerical results and conclusions
In this article, the unsteady motion of an incompressible OldroydB fluid over an infinite plate that applies an oscillating shear stress to the fluid is studied by means of integral transforms. The starting solutions that have been obtained for velocity and shear stress are presented as a sum of steadystate and transient solutions. They describe the motion of the fluid some time after its initiation. After that time, when the transients disappear, the starting solutions tend to the steadystate solutions that are periodic in time and independent of the initial conditions. However, they satisfy the governing equations and boundary conditions. Furthermore, as it was to be expected, the steadystate solutions corresponding to the cosine oscillations of the shear differ with a phase shift from those due to the sine oscillations of the shear. This property is not true for the transient components of solutions. This is the reason that we separately gave the starting solutions for both cosine and sine oscillations of the shear stress on the boundary.
By making λ_{r }= 0 into general solutions (27), (28), (42), and (43), we obtain the similar solutions (44)(47) corresponding to a Maxwell fluid performing the same motion. Making λ → 0 into Equations (27), (28), (42), and (43) we obtain the similar solutions (48)(51) corresponding to a second grade fluid performing the same motion. Solutions (44)(47) with λ → 0, respectively (48)(51) with α = 0 give the similar solutions (53)(56) for Newtonian fluids. It is worth pointing out that the expressions of τ_{cN }(y, t) and τ_{sN }(y, t) are identical, as form, with those of v_{cN }(y, t) and v_{sN }(y, t) corresponding to a similar motion with the boundary conditions [[13], Equation (3.1)]
The velocity field (see [[13], Equations (3.11) and (6.2)])
for instance, has the same form as τ_{sN }(y, t) given by Equation (56). This is not a surprise because, for Newtonian fluids, Equations (3) and (4) together with the balance of linear momentum lead to a governing equation for shear stress of the same form as that for velocity.
Generally, the starting solutions for unsteady motions of fluids are important for those who need to eliminate the transients from their rheological measurements. Consequently, an important problem regarding the technical relevance of these solutions is to find the required time to get the steadystate. More exactly, in practice it is necessary to know the approximate time after which the fluid is moving according to the steadystate solutions. For this, the variations of the starting and steadystate velocities with the distance from the plate are depicted in Figures 1, 2, 3, 4, 5, 6, 7, and 8 for sine and cosine oscillations of the shear stress on the boundary. At small values of time, the difference between the starting and steadystate solutions is meaningful. This difference decreases in time and it is clearly seen from figures that the required time to reach the steadystate for the sine oscillations is higher in comparison with the cosine oscillations of the shear. This is obvious, because at t = 0 the shear stress on the boundary is zero for sine oscillations.
Figure 1. Comparison between the starting velocity u_{s}(y, t) and the steadystate velocity u_{ss}(y, t), for f = 6, ν = 0.004, μ = 3.92, λ = 8, λ_{r }= 4 and different values of time t and frequency ω.
Figure 2. Comparison between the starting velocity u_{s}(y, t) and the steadystate velocity u_{ss}(y, t), for f = 6, ν = 0.004, μ = 3.92, λ_{r }= 0.5, ω = 0.08 and different values of time t and the relaxation time λ.
Figure 3. Comparison between the starting velocity u_{s}(y, t) and the steadystate velocity u_{ss}(y, t), for f = 6, ν = 0.004, μ = 3.92, λ = 4, ω = 0.08 and different values of time t and the retardation time λ_{r}.
Figure 4. Comparison between the starting velocity u_{c}(y, t) and the steadystate velocity u_{cs}(y, t), for f = 6, ν = 0.004, μ = 3.92, λ = 8, λ_{r }= 4, and different values of time t and frequency ω.
Figure 5. Comparison between the starting velocity u_{c}(y, t) and the steadystate velocity u_{cs}(y, t), for f = 6, ν = 0.004, μ = 3.92, λ_{r }= 0.5, ω = 0.08 and different values of time t and the relaxation time λ.
Figure 6. Comparison between the starting velocity u_{c}(y, t) and the steadystate velocity u_{cs}(y, t), for f = 6, ν = 0.004, μ = 3.92, λ = 4, ω = 0.08 and different values of time t and the retardation time λ_{r}.
Figure 7. Comparison between the starting velocity and steadystate velocity for second grade and Newtonian fluids for f = 6, ν = 0.004, μ = 3.92, α = 0.016 and different values of time t and frequency ω (sine oscillations).
Figure 8. Comparison between the starting velocity and steadystate velocity for second grade and Newtonian fluids for f = 6, ν = 0.004, μ = 3.92, α = 0.016 and different values of time t and frequency ω (cosine oscillations).
Of course, the required time to reach the steadystate depends on the material constants and the frequency ω of the shear. Figures 1 and 4 show the influence of ω on the fluid motion. Furthermore, the required time to reach the steadystate decreases for increasing ω. The influence of the relaxation and retardation times λ and λ_{r }on the fluid motion is underlined by Figures 2, 3, 5, and 6. The two parameters, as expected, have opposite effects on the motion. The required time to reach the statestate decreases with respect to λ and increases with regard to λ_{r }for both types of oscillating shear. The steadystate is rather obtained for a Maxwell fluid in comparison with an OldroydB fluid having the same relaxation time λ. Figures 7 and 8 show the influence of ω on the motion of second grade and Newtonian fluids. As expected, the required time to reach the steadystate is lower for second grade fluids in comparison with Newtonian fluids.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
The authors' contribution for this article is as follows: NS 35%, MR 35%, and IS 30%. All authors participated in the design of the study and established exact solutions. NS and MR have had a special contribution to get the simplest form of steadystate solutions. IS prepared the most part of graphical illustrations while NS and MR participated in the sequence alignment and drafted the manuscript. All authors participated in drawing conclusions, reading and approving the manuscript.
Appendix
Acknowledgements
The authors would like to express their gratitude to reviewers for fruitful remarks and suggestions regarding the initial version of the manuscript. N. Shahid and M. Rana are highly thankful and grateful to Abdus Salam School of Mathematical Sciences GC University, Lahore, Pakistan and Higher Education Commission of Pakistan for generously supporting and facilitating this research work.
References

Rajagopal, KR, Srinivasa, AR: A thermodynamic frame work for rate type fluid models. J NonNewtonian Fluid Mech. 88, 207–227 (2000). Publisher Full Text

Tanner, RY: Note of the Rayleigh problem for a viscoelastic fluid. Z Angew Math Phys. 13, 573–580 (1962). Publisher Full Text

Waters, ND, King, MJ: Unsteady flow of an elasticoviscous liquid. Rheol Acta. 9, 345–355 (1970). Publisher Full Text

Rajagopal, KR: A note on unsteady unidirectional flows of a nonNewtonian fluid. Int J Nonlinear Mech. 17, 369–373 (1982). Publisher Full Text

Erdogan, ME: A note on an unsteady flow of a viscous fluid due to an oscillating plane wall. Int J Nonlinear Mech. 35, 1–6 (2000). Publisher Full Text

Penton, R: The transient for Stokes' oscillating plane: a solution in terms of tabulated functions. J Fluid Mech. 31, 819–825 (1968). Publisher Full Text

Rajagopal, KR: Longitudinal and torsional oscillations of a rod in a nonNewtonian fluid. Acta Mech. 49, 281–285 (1983). Publisher Full Text

Puri, P, Kythe, PK: Stokes' first and second problems for RivlinEricksen fluids with nonclassical heat conduction. ASME J Heat Trans. 120, 44–50 (1998). Publisher Full Text

Rajagopal, KR, Bhatnagar, RK: Exact solutions for some simple flows of an OldroydB fluid. Acta Mech. 113, 233–239 (1995). Publisher Full Text

Hayat, T, Asghar, S, Siddiqui, AM: Periodic unsteady flows of a nonNewtonian fluid. Acta Mech. 131, 169–175 (1998). Publisher Full Text

Siddiqui, AM, Hayat, T, Asghar, S: Periodic flows of a nonNewtonian fluid between two parallel plates. Int J Nonlinear Mech. 34, 895–899 (1999). Publisher Full Text

Fetecau, C, Fetecau, C: Starting solutions for some unsteady unidirectional flows of a second grade fluid. Int J Eng Sci. 43, 781–789 (2005). Publisher Full Text

Fetecau, C, Fetecau, C: Starting solutions for the motion of a second grade fluid due to longitudinal and torsional oscillations of a circular cylinder. Int J Eng Sci. 44, 788–796 (2006). Publisher Full Text

Aksel, N, Fetecau, C, Scholle, M: Starting solutions for some unsteady unidirectional flows of OldroydB fluids. Z Angew Math Phys. 57, 815–831 (2006). Publisher Full Text

Vieru, D, Akhtar, W, Fetecau, C, Fetecau, C: Starting solutions for the oscillating motion of a Maxwell fluid in cylindrical domains. Mecca. 42, 573–583 (2007). Publisher Full Text

Nazar, M, Fetecau, C, Vieru, D, Fetecau, C: New exact solutions corresponding to the second problem of Stokes for second grade fluids. Nonlinear Anal Real World Appl. 11, 584–591 (2008)

Fetecau, C, Hayat, T, Fetecau, C: Starting solutions for oscillating motions of OldrydB fluids in cylindrical domains. Int J NonNewtonian Fluid Mech. 153, 191–201 (2008). Publisher Full Text

Zheng, L, Zhao, F, Zhang, X: Exact solutions for generalized Maxwell fluid flow due to oscillatory and constantly accelerating plate. Nonlinear Anal Real World Appl. 11, 3744–3751 (2010). Publisher Full Text

Anjum, A, Ayub, M, Khan, M: Starting solutions for oscillating motions of an OldroydB fluid over a plane wall. Commun Nonlinear Sci Numer Simul. 17, 472–482 (2012). Publisher Full Text

Khan, M, Arshad, M, Anjum, A: On exact solutions of Stokes second problem for MHD OldroydB fluid. Nucl Eng Des. 243, 20–32 (2012)

Waters, ND, King, MJ: Unsteady flow of an elasticoviscous liquid. Rheol Acta. 9, 345–355 (1970). Publisher Full Text

Bandelli, R, Rajagopal, KR, Galdi, GP: On some unsteady motions of fluids of second grade. Arch Mech (Arch Mech Stos). 47, 661–676 (1995)

Fetecau, C, Vieru, D, Fetecau, C: E^{®}ect of the side walls on the motion of a viscous fluid induced by an infinite plate that applies an oscillating shear stress to the fluid. Cent Eur J Phys. 9, 816–825 (2011). Publisher Full Text

Fetecau, C, Prasad, SC, Rajagopal, KR: A note on the flow induced by a constantly accelerating plate in an OldryodB fluid. Appl Math Model. 31, 647–654 (2007). Publisher Full Text

Debnath L, Bhatta D (eds.): Integral Transforms and Their Applications, Chapman and Hall, BocaRaton (2007)

Sneddon, IN: Fourier Transforms, McGraw Hill, New York (1951)