Skip to main content

Exact solution for motion of an Oldroyd-B fluid over an infinite flat plate that applies an oscillating shear stress to the fluid

Abstract

The unsteady motion of an Oldroyd-B 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 steady-state 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 steady-state solutions that are periodic in time and independent of the initial conditions. Finally, the required time to reach the steady-state 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.

1 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 non-Newtonian fluids, has considerably grown due to their multiple applications. Among the many models that have been used to describe the behavior of non-Newtonian 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 Oldroyd-B 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 Oldroyd-B 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 steady-state 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 non-Newtonian 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 closed-form expression for the starting solutions corresponding to the motion induced by a flat plate that applies an oscillating shear to an Oldroyd-B fluid has not been given before.

Consequently, the aim of this article is to determine starting solutions for the unsteady motion of an incompressible Oldroyd-B 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 steady-state 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 steady-state 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 Oldroyd-B fluid is characterized by the following constitutive equations [9, 14, 17]

T = - p I + S , S + λ S . - L S - S L T = μ A + λ r A . - L A - A L T ,
(1)

where T is the Cauchy stress tensor, -p I denotes the indeterminate spherical stress, S is the extra-stress tensor, L is the velocity gradient, A = L + LT is the first Rivlin-Ericksen 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

V = V ( y , t ) = u ( y , t ) i ,
(2)

where i denotes the unit vector along the x-direction 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 extra-stress 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]

u ( y , t ) t + λ 2 u ( y , t ) t 2 = ν 1 + λ r t 2 u ( y , t ) y 2 ; y , t > 0 ,
(3)

where v = μ /ρ is the kinematic viscosity and ρ is the constant density of the fluid. The non-trivial shear stress τ(y, t) = S xy (y, t) satisfies the partial differential equation [[24], Equation (4)]

1 + λ t τ ( y , t ) = μ 1 + λ r t u ( y , t ) y ; y , t > 0 .
(4)

3 Statement of the problem

Let us consider an incompressible Oldroyd-B 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

u ( y , 0 ) = 0 , u ( y , 0 ) t = 0 , τ ( y , 0 ) = 0 , y > 0 ,
(5)
1 + λ t τ ( y , t ) y = 0 = μ 1 + λ r t u ( y , t ) y y = 0 = f sin ω t o r f cos ω t t > 0 .
(6)

Moreover, the natural condition

u ( y , t ) 0 a s y ,
(7)

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

V ( y , t ) = u c ( y , t ) + i u s ( y , t ) , T ( y , t ) = τ c ( y , t ) + i τ s ( y , t ) ,
(8)

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

V ( y , t ) t + λ 2 V ( y , t ) t 2 = ν 1 + λ r t 2 V ( y , t ) y 2 y , t > 0 ,
(9)
1 + λ t T ( y , t ) = μ 1 + λ r t V ( y , t ) y y , t > 0 ,
(10)
V ( y , 0 ) = 0 , V ( y , 0 ) t = 0 , T ( y , 0 ) = 0 y > 0 ,
(11)
1 + λ t T ( y , t ) y = 0 = μ 1 + λ r t V ( y , t ) y y = 0 = f e i ω t t > 0 ,
(12)
V ( y , t ) 0 a s y t > 0 .
(13)

4.1 Calculation of the velocity field

In order to determine the solution of initial-boundary values problem (9), (11)1,2, (12)2 and (13), we first take the Laplace transform [25] of Equation (9) and obtain

q V ̄ ( y , q ) + λ q 2 V ̄ ( y , q ) = ν 1 + λ r q 2 V ̄ ( y , q ) y 2 ,
(14)

where the Laplace transform V ̄ ( y , q ) of function V (y, t) has to satisfy the conditions

V ̄ ( y , q ) y y = 0 = f μ ( q - i ω ) ( 1 + λ r q ) ,
(15)
V ̄ ( y , q ) 0 a s y .
(16)

Multiplying Equation (14) 2 π cos ( y ξ ) , integrating the result with respect to y from 0 to infinity and using Equations (15) and (16), we obtain

V ̄ c ( ξ , q ) = - 2 π f ρ 1 ( q - i ω ) [ λ q 2 + ( 1 + λ r ν ξ 2 ) q + ν ξ 2 ] ,
(17)

where

V ̄ c ( ξ , q ) = 2 π 0 V ̄ ( y , q ) cos ( y ξ ) d y ,
(18)

denotes the Fourier cosine transform [26] of function V ̄ ( y , q ) . Equation (17) can be written as

V ̄ c ( ξ , q ) = V ̄ c 1 ( ξ , q ) + V ̄ c 2 ( ξ , q ) ,
(19)

where

V ̄ c 1 ( ξ , q ) = - 2 π f ρ ( ν ξ 2 - λ ω 2 ) - i ω ( 1 + λ r ν ξ 2 ) ( ν ξ 2 - λ ω 2 ) 2 + ω 2 ( 1 + λ r ν ξ 2 ) 2 1 q - i ω ,
(20)
V ̄ c 2 ( ξ , q ) = 2 π f ρ ( ν ξ 2 - λ ω 2 ) - i ω ( 1 + λ r ν ξ 2 ) ( ν ξ 2 - λ ω 2 ) 2 + ω 2 ( 1 + λ r ν ξ 2 ) 2 λ q + [ i λ ω + ( 1 + λ r ν ξ 2 ) ] λ q 2 + ( 1 + λ r ν ξ 2 ) q + ν ξ 2 .
(21)

Applying the inverse Laplace transform and then the inverse Fourier cosine transform to Equation (20), we get the following expression

V 1 ( y , t ) = - 2 π f e i ω t ρ ν 2 α 0 ( ν ξ 2 - λ ω 2 ) - i ω ( 1 + λ r ν ξ 2 ) ( ξ 2 - β 2 ) 2 + γ 2 cos ( y ξ ) d ξ ,
(22)

where

α = 1 + λ r 2 ω 2 , β 2 = ω 2 ( λ - λ r ) ν ( 1 + λ r 2 ω 2 ) , γ = ω ( 1 + λ λ r ω 2 ) ν ( 1 + λ r 2 ω 2 ) .

Using (A1) and (A2) from Appendix 1, we obtain the following simplified expression

V 1 ( y , t ) = f μ ν ω e - y B ( 1 + λ r 2 ω 2 ) ( 1 + λ 2 ω 2 ) 4 cos ω t - y A + φ + π 2 + i sin ω t - y A + φ + π 2
(23)

where

2 A 2 = β 4 + γ 2 + β 2 , 2 B 2 = β 4 + γ 2 - β 2 a n d tan φ = 1 + λ 2 ω 2 - λ ω 1 + λ r 2 ω 2 1 + λ r 2 ω 2 + λ r ω 1 + λ 2 ω 2 .

Now, for Equation (21), we introduce the function

F ( ξ , q ) = λ q + i λ ω + ( 1 + λ r ν ξ 2 ) λ q 2 + ( 1 + λ r ν ξ 2 ) q + ν ξ 2 ,
(24)

which can be written in the following equivalent form

F ( ξ , q ) = q + b ( ξ ) 2 λ q + b ( ξ ) 2 λ 2 - c ( ξ ) 2 λ 2 + b ( ξ ) + i 2 λ ω c ( ξ ) c ( ξ ) 2 λ q + b ( ξ ) 2 λ 2 - c ( ξ ) 2 λ 2 ,
(25)

where

b ( ξ ) = 1 + λ r ν ξ 2 , c ( ξ ) = ( 1 + λ r ν ξ 2 ) 2 - 4 ν λ ξ 2 .

Applying the inverse Laplace transform and then the inverse Fourier transform to Equation (21) and using Equation (25), we obtain the following expression

V 2 ( y , t ) = 2 π f ν μ α 0 cos ( y ξ ) ( ξ 2 - β 2 ) 2 + γ 2 ( ν ξ 2 - λ ω 2 ) c h c ( ξ ) t 2 λ + ( ν ξ 2 + λ ω 2 ) b ( ξ ) c ( ξ ) s h c ( ξ ) t 2 λ e - b ( ξ ) t 2 λ d ξ - i 2 π f ω ν μ α 0 cos ( y ξ ) ( ξ 2 - β 2 ) 2 + γ 2 b ( ξ ) c h c ( ξ ) t 2 λ + b 2 ( ξ ) - 2 λ ( ν ξ 2 - λ ω 2 ) c ( ξ ) s h c ( ξ ) t 2 λ e - b ( ξ ) t 2 λ d ξ .
(26)

Finally, the velocity corresponding to the cosine oscillations of the shear is given by

u c ( y , t ) = f μ ν ω e - y B ( 1 + λ r 2 ω 2 ) ( 1 + λ 2 ω 2 ) 4 cos ω t - y A + φ + π 2 + 2 π f ν μ α × 0 cos ( y ξ ) ( ξ 2 - β 2 ) 2 + γ 2 ( ν ξ 2 - λ ω 2 ) c h c ( ξ ) t 2 λ + ( ν ξ 2 + λ ω 2 ) b ( ξ ) c ( ξ ) s h c ( ξ ) t 2 λ e - b ( ξ ) t 2 λ d ξ ,
(27)

while that corresponding to sine oscillations has the form

u s ( y , t ) = f μ ν ω e - y B ( 1 + λ r 2 ω 2 ) ( 1 + λ 2 ω 2 ) 4 sin ω t - y A + φ + π 2 - 2 π f ω ν μ α × 0 cos ( y ξ ) ( ξ 2 - β 2 ) 2 + γ 2 b ( ξ ) c h ( c ( ξ ) t 2 λ ) + b 2 ( ξ ) - 2 λ ( ν ξ 2 - λ ω 2 ) c ( ξ ) s h ( c ( ξ ) t 2 λ ) e - b ( ξ ) t 2 λ d ξ .
(28)

The starting solutions (27) and (28) corresponding to cosine and sine oscillations of the shear, are presented as a sum between the steady-state 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 steady-state solutions

u c s ( y , t ) = f μ ν ω e - y B ( 1 + λ r 2 ω 2 ) ( 1 + λ 2 ω 2 ) 4 cos ω t - y A + φ + π 2 = f μ ν ω e - y B ( 1 + λ r 2 ω 2 ) ( 1 + λ 2 ω 2 ) 4 sin ( ω t - y A + φ + π ) ,
(29)

respectively,

u s s ( y , t ) = f μ ν ω e - y B ( 1 + λ r 2 ω 2 ) ( 1 + λ 2 ω 2 ) 4 sin ω t - y A + φ + π 2 .
(30)

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):

T ̄ ( y , q ) = 2 π f μ ρ 0 ξ sin ( y ξ ) 1 1 + λ q 1 q - i ω λ r q + 1 λ q 2 + ( 1 + λ r ν ξ 2 ) q + ν ξ 2 d ξ .
(31)

Equation (31) can be written as

T ̄ ( y , q ) = T ̄ 1 ( y , q ) + T ̄ 2 ( y , q ) + T ̄ 3 ( y , q ) ,
(32)

where

T ̄ 1 ( y , q ) = - 2 f π ( 1 + λ 2 ω 2 ) 1 - i λ ω q + 1 λ 0 sin ( y ξ ) ξ d ξ ,
(33)
T ̄ 2 ( y , q ) = 2 f π ( 1 + λ 2 ω 2 ) 1 q - i ω 0 ξ sin ( y ξ ) ξ 2 - β 2 - λ ω γ ( ξ 2 - β 2 ) 2 + γ 2 - i λ ω ξ 2 - λ ω β 2 + γ ( ξ 2 - β 2 ) 2 + γ 2 d ξ ,
(34)
T ̄ 3 ( y , q ) = 2 π f μ ρ 0 ξ sin ( y ξ ) M ( ξ ) q + N ( ξ ) λ q 2 + ( 1 + λ r ν ξ 2 ) q + ν ξ 2 d ξ ,
(35)

with

M ( ξ ) = λ ν ( 1 + λ 2 ω 2 ) ξ 2 ( λ ω γ - β 2 ) + ( γ 2 + β 4 ) + i [ ξ 2 ( γ + λ ω β 2 ) - λ ω ( γ 2 + β 4 ) ] ξ 2 ( ( ξ 2 - β 2 ) 2 + γ 2 ) , N ( ξ ) = 1 ω ( 1 + λ 2 ω 2 ) - ξ 2 ( γ + λ ω β 2 ) + λ ω ( γ 2 + β 4 ) + i [ ξ 2 ( λ ω γ - β 2 ) + ( γ 2 + β 4 ) ] ( ξ 2 - β 2 ) 2 + γ 2 ) .

Applying the inverse Laplace transform to the Equation (33) and using (A3) from Appendix 1, we get

T 1 ( y , t ) = - f 1 + λ 2 ω 2 + i f λ ω 1 + λ 2 ω 2 e - t λ .
(36)

Similarly, using (A4) and (A5) from Appendix 1, from Equation (34), we obtain the following suitable form for T2(y, t)

T 2 ( y , t ) = f e - y B 1 + λ 2 ω 2 cos ( ω t - y A ) + λ ω sin ( ω t - y A ) + i ( sin ( ω t - y A ) - λ ω cos ( ω t - y A ) ) .
(37)

A straightforward calculation leads to the following simplified form

T 2 ( y , t ) = f e - y B 1 + λ 2 ω 2 cos ( ω t - y A - ψ ) + i sin ( ω t - y A - ψ )
(38)

where tan ψ = λω.

Now, we consider the function

G ( ξ , q ) = M ( ξ ) q + N ( ξ ) λ q 2 + ( 1 + λ r ν ξ 2 ) q + ν ξ 2 ,
(39)

which can be written in the following equivalent form:

G ( ξ , q ) = M ( ξ ) q + b ξ 2 λ λ ( q + b ξ 2 λ ) 2 - ( b ξ 2 λ ) 2 + - b ( ξ ) M ( ξ ) λ c ( ξ ) + 2 N ( ξ ) c ( ξ ) c ξ 2 λ ( q + b ξ 2 λ ) 2 - ( c ξ 2 λ ) 2 .
(40)

Applying the inverse Laplace transform to Equation (35) and using Equation (40), we obtain

T 3 ( y , t ) = 2 f π ( 1 + λ 2 ω 2 ) 0 sin ( y ξ ) ξ [ ( ξ 2 - β 2 ) 2 + ( γ ) 2 ] ( ξ 2 ( λ ω γ - β 2 ) + ( γ 2 + β 4 ) ) c h c ( ξ ) t 2 λ - p ( ξ ) ω c ( ξ ) s h c ( ξ ) t 2 λ + i [ ξ 2 ( γ + λ ω β 2 ) - λ ω ( γ 2 + β 4 ) ] c h c ( ξ ) t 2 λ + r ( ξ ) ω c ( ξ ) s h c ( ξ ) t 2 λ e - b ( ξ ) t 2 λ d ξ ,
(41)

where

p ( ξ ) = ν ξ 4 [ λ r ω ( λ ω γ - β 2 ) + 2 ( γ + λ ω β 2 ) ] + ω ( γ 2 + β 4 ) + ω ξ 2 [ ( λ ω γ - β 2 ) - ν ( γ 2 + β 4 ) ( 2 λ - λ r ) ] , r ( ξ ) = ν ξ 4 [ 2 ( λ ω γ - β 2 ) - λ r ω ( γ + λ ω β 2 ) ] + ξ 2 [ ν ( γ 2 + β 4 ) ( 2 + λ λ r ω 2 ) - ω ( γ + λ ω β 2 ) ] + λ ω 2 ( γ 2 + β 4 ) .

Using Equations (32), (36), (38), and (41), the shear stress corresponding to cosine oscillations of the shear can be written in the form

τ c ( y , t ) = - f 1 + λ 2 ω 2 e - t λ + f e - y B 1 + λ 2 ω 2 cos ( ω t - y A - ψ ) + 2 f π ( 1 + λ 2 ω 2 ) 0 sin ( y ξ ) ξ [ ( ξ 2 - β 2 ) 2 + ( γ ) 2 ] × ( ξ 2 ( λ ω γ - β 2 ) + ( γ 2 + β 4 ) ) c h c ( ξ ) t 2 λ - p ( ξ ) ω c ( ξ ) s h c ( ξ ) t 2 λ e - b ( ξ ) t 2 λ d ξ .
(42)

Also, the shear stress corresponding to sine oscillations is given by

τ s ( y , t ) = f λ ω 1 + λ 2 ω 2 e - t λ + f e - y B 1 + λ 2 ω 2 sin ( ω t - y A - ψ ) + 2 f π ( 1 + λ 2 ω 2 ) 0 sin ( y ξ ) ξ [ ( ξ 2 - β 2 ) 2 + ( γ ) 2 ] × ξ 2 ( γ + λ ω β 2 ) - λ ω ( γ 2 + β 4 ) c h c ( ξ ) t 2 λ + r ( ξ ) ω c ( ξ ) s h c ( ξ ) t 2 λ e - b ( ξ ) t 2 λ d ξ .
(43)

Of course, the shear stresses given by Equations (42) and (43) are also presented as a sum of steady-state and transient solutions. The steady-state solutions

τ c s ( y , t ) = f e - y B 1 + λ 2 ω 2 cos ( ω t - y A - ψ ) a n d τ s s ( y , t ) = f e - y B 1 + λ 2 ω 2 sin ( ω t - y A - ψ ) ,

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

u c M ( y , t ) = f μ ν ω e - y B ( 1 + λ 2 ω 2 ) 4 cos ω t - y A + φ + π 2 + 2 f π μ ν 0 cos ( y ξ ) ξ 2 - λ ω 2 ν 2 + ω ν 2 × ( ν ξ 2 - λ ω 2 ) c h 1 - 4 λ ν ξ 2 2 λ t + ν ξ 2 + λ ω 2 1 - 4 λ ν ξ 2 s h 1 - 4 λ ν ξ 2 2 λ t e - t 2 λ d ξ ,
(44)
u s M ( y , t ) = f μ ν ω e - y B ( 1 + λ 2 ω 2 ) 4 sin ω t - y A + φ + π 2 - 2 π f ω μ ν 0 cos ( y ξ ) ξ 2 - λ ω 2 ν 2 + ω ν 2 × c h 1 - 4 λ ν ξ 2 2 λ t + 1 - 2 λ ( ν ξ 2 - λ ω 2 ) 1 - 4 λ ν ξ 2 s h 1 - 4 λ ν ξ 2 2 λ t e - t 2 λ d ξ ,
(45)

corresponding to a Maxwell fluid performing the same motion. Similarly, from Equations (42) and (43), we obtain the corresponding shear stresses

τ c M ( y , t ) = - f 1 + λ 2 ω 2 e - t λ + f e - y B 1 + λ 2 ω 2 cos ( ω t - y A - ψ ) + 2 f π 0 sin ( y ξ ) e - t 2 λ ξ ( ν ξ 2 - λ ω 2 ) 2 + ω 2 × ω 2 c h 1 - 4 λ ν ξ 2 2 λ t - ω 2 + 2 ν ξ 2 ( ν ξ 2 - λ ω 2 ) 1 - 4 λ ν ξ 2 s h 1 - 4 λ ν ξ 2 2 λ t d ξ .
(46)
τ s M ( y , t ) = f λ ω 1 + λ 2 ω 2 e - t λ + f e - y B 1 + λ 2 ω 2 sin ( ω t - y A - ψ ) + 2 f ω π 0 sin ( y ξ ) e - t 2 λ ξ ( ν ξ 2 - λ ω 2 ) 2 + ω 2 × ( ν ξ 2 - λ ω 2 ) c h 1 - 4 λ ν ξ 2 2 λ t + ν ξ 2 + λ ω 2 1 - 4 λ ν ξ 2 s h 1 - 4 λ ν ξ 2 2 λ t d ξ .
(47)

5.2 λ→ 0 (second grade fluid)

Using Equations (27), (28), (42), (43) and the limits

lim λ 0 e - b ξ t 2 λ c h c ξ t 2 λ = lim λ 0 e - b ξ t 2 λ s h c ξ t 2 λ = 1 2 exp - γ ξ 2 t 1 + ν λ r t ,

we obtain solutions corresponding to a fluid of second grade:

u c S G ( y , t ) = f μ ν ω e - y B 1 + δ 2 4 cos ω t - y A + φ + π 2 + 2 f π μ 1 + δ 2 0 ξ 2 cos ( y ξ ) ξ 2 + β 2 2 + γ 2 exp - ν ξ 2 t 1 + α ξ 2 d ξ ,
(48)
u s S G ( y , t ) = f μ ν ω e - y B 1 + δ 2 4 sin ω t - y A + φ + π 2 - 2 f π μ ν 1 + δ 2 0 1 + α ξ 2 cos ( y ξ ) ξ 2 + β 2 2 + γ 2 exp - ν ξ 2 t 1 + α ξ 2 d ξ ,
(49)
τ c S G ( y , t ) = f e - y B cos ω t - y A - 2 f ν γ π ω 0 ξ 3 sin ( y ξ ) ξ 2 + β 2 2 + γ 2 1 + α ξ 2 exp - ν ξ 2 t 1 + α ξ 2 d ξ ,
(50)
τ s S G ( y , t ) = f e - y B sin ω t - y A - 2 f π γ 0 ξ 3 sin ( y ξ ) ξ 2 + β 2 2 + γ 2 exp - ν ξ 2 t 1 + α ξ 2 d ξ ,
(51)

where

α = ν λ r , δ = α ω ν , β 2 = α ω 2 ν 2 ( 1 + δ 2 ) , γ = ω ν ( 1 + δ 2 ) , 2 A 2 = γ ( 1 + δ 2 - δ ) , 2 B 2 = γ ( 1 + δ 2 + δ ) , t a n φ = ( 1 + δ 2 - δ ) .
(52)

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)]

u c N ( y , t ) = f μ ν ω e - y ω 2 ν cos ω t - y ω 2 ν + 3 π 4 + 2 f μ π 0 ξ 2 cos ( y ξ ) ξ 4 + ω ν 2 e - ν ξ 2 t dξ,
(53)
u s N ( y , t ) = f μ ν ω e - y ω 2 ν sin ω t - y ω 2 ν + 3 π 4 - 2 f μ π ω ν 0 cos ( y ξ ) ξ 4 + ω ν 2 e - ν ξ 2 t dξ,
(54)
τ c N ( y , t ) =f e - y ω 2 ν cos ω t - y ω 2 ν - 2 f π 0 ξ 3 sin ( y ξ ) ξ 4 + ( ω ν ) 2 e - ν ξ 2 t dξ,
(55)
τ s N ( y , t ) =f e - y ω 2 ν sin ω t - y ω 2 ν + 2 f π ω ν 0 ξ sin ( y ξ ) ξ 4 + ( ω ν ) 2 e - ν ξ 2 t dξ,
(56)

corresponding to the flow of a Newtonian fluid.

6 Numerical results and conclusions

In this article, the unsteady motion of an incompressible Oldroyd-B 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 steady-state 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 steady-state 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 steady-state 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)]

v 0 , t = V cos ω t o r v 0 , t = V sin ω t ; t > 0 .
(57)

The velocity field (see [[13], Equations (3.11) and (6.2)])

v s N ( y , t ) =V e - y ω 2 ν sin ω t - y ω 2 ν + 2 V π ω ν 0 ξ sin ( y ξ ) ξ 4 + ω ν 2 e - ν ξ 2 t dξ,
(58)

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 steady-state. More exactly, in practice it is necessary to know the approximate time after which the fluid is moving according to the steady-state solutions. For this, the variations of the starting and steady-state 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 steady-state solutions is meaningful. This difference decreases in time and it is clearly seen from figures that the required time to reach the steady-state 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
figure 1

Comparison between the starting velocity u s ( y , t ) and the steady-state 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
figure 2

Comparison between the starting velocity u s ( y , t ) and the steady-state 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
figure 3

Comparison between the starting velocity u s ( y , t ) and the steady-state 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
figure 4

Comparison between the starting velocity u c ( y , t ) and the steady-state 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
figure 5

Comparison between the starting velocity u c ( y , t ) and the steady-state 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
figure 6

Comparison between the starting velocity u c ( y , t ) and the steady-state 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
figure 7

Comparison between the starting velocity and steady-state 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
figure 8

Comparison between the starting velocity and steady-state 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 steady-state 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 steady-state 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 state-state decreases with respect to λ and increases with regard to λ r for both types of oscillating shear. The steady-state is rather obtained for a Maxwell fluid in comparison with an Oldroyd-B 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 steady-state is lower for second grade fluids in comparison with Newtonian fluids.

Appendix

0 x 2 cos ( m x ) ( x 2 - b 2 ) 2 + c 2 dx= π e - m B 2 c A cos ( m A ) - B sin ( m A ) ,
(A1)
0 cos ( m x ) ( x 2 - b 2 ) 2 + c 2 dx= π e - m B 2 c ( A 2 + B 2 ) A cos ( m A ) + B sin ( m A ) ,
(A2)
0 sin ( y ξ ) ξ = π 2 ,
(A3)
0 x sin ( m x ) ( x 2 - b 2 ) ( x 2 - b 2 ) 2 + c 2 dx= π 2 e - m B cos ( m A ) ,
(A4)
0 x sin ( m x ) x ( x 2 - b 2 ) 2 + c 2 dx= π 2 c e - m B sin ( m A ) ,
(A5)
0 sin ( m x ) x ( x 2 - b 2 ) 2 + c 2 dx= π 2 c ( b 4 + c 2 ) c + b 2 sin ( m A ) - c cos ( m A ) . exp ( - m B ) ,
(A6)

where

b 2 = ω 2 ( λ - λ r ) ν ( 1 + λ r 2 ω 2 ) , c = ω ( 1 + λ λ r ω 2 ) ν ( 1 + λ r 2 ω 2 ) , 2 A 2 = b 4 + c 2 + b 2 , 2 B 2 = b 4 + c 2 - b 2 .

References

  1. Rajagopal KR, Srinivasa AR: A thermodynamic frame work for rate type fluid models. J Non-Newtonian Fluid Mech 2000, 88: 207-227. 10.1016/S0377-0257(99)00023-3

    Article  Google Scholar 

  2. Tanner RY: Note of the Rayleigh problem for a visco-elastic fluid. Z Angew Math Phys 1962, 13: 573-580. 10.1007/BF01595580

    Article  MathSciNet  Google Scholar 

  3. Waters ND, King MJ: Unsteady flow of an elastico-viscous liquid. Rheol Acta 1970, 9: 345-355. 10.1007/BF01975401

    Article  Google Scholar 

  4. Rajagopal KR: A note on unsteady unidirectional flows of a non-Newtonian fluid. Int J Non-linear Mech 1982, 17: 369-373. 10.1016/0020-7462(82)90006-3

    Article  Google Scholar 

  5. Erdogan ME: A note on an unsteady flow of a viscous fluid due to an oscillating plane wall. Int J Non-linear Mech 2000, 35: 1-6. 10.1016/S0020-7462(99)00019-0

    Article  MathSciNet  Google Scholar 

  6. Penton R: The transient for Stokes' oscillating plane: a solution in terms of tabulated functions. J Fluid Mech 1968, 31: 819-825. 10.1017/S0022112068000509

    Article  Google Scholar 

  7. Rajagopal KR: Longitudinal and torsional oscillations of a rod in a non-Newtonian fluid. Acta Mech 1983, 49: 281-285. 10.1007/BF01236358

    Article  Google Scholar 

  8. Puri P, Kythe PK: Stokes' first and second problems for Rivlin-Ericksen fluids with nonclassical heat conduction. ASME J Heat Trans 1998, 120: 44-50. 10.1115/1.2830061

    Article  Google Scholar 

  9. Rajagopal KR, Bhatnagar RK: Exact solutions for some simple flows of an Oldroyd-B fluid. Acta Mech 1995, 113: 233-239. 10.1007/BF01212645

    Article  MathSciNet  Google Scholar 

  10. Hayat T, Asghar S, Siddiqui AM: Periodic unsteady flows of a non-Newtonian fluid. Acta Mech 1998, 131: 169-175. 10.1007/BF01177223

    Article  MathSciNet  Google Scholar 

  11. Siddiqui AM, Hayat T, Asghar S: Periodic flows of a non-Newtonian fluid between two parallel plates. Int J Non-linear Mech 1999, 34: 895-899. 10.1016/S0020-7462(98)00063-8

    Article  Google Scholar 

  12. Fetecau C, Fetecau C: Starting solutions for some unsteady unidirectional flows of a second grade fluid. Int J Eng Sci 2005, 43: 781-789. 10.1016/j.ijengsci.2004.12.009

    Article  MathSciNet  Google Scholar 

  13. 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 2006, 44: 788-796. 10.1016/j.ijengsci.2006.04.010

    Article  MathSciNet  Google Scholar 

  14. Aksel N, Fetecau C, Scholle M: Starting solutions for some unsteady unidirectional flows of Oldroyd-B fluids. Z Angew Math Phys 2006, 57: 815-831. 10.1007/s00033-006-0063-8

    Article  MathSciNet  Google Scholar 

  15. Vieru D, Akhtar W, Fetecau C, Fetecau C: Starting solutions for the oscillating motion of a Maxwell fluid in cylindrical domains. Mecca 2007, 42: 573-583. 10.1007/s11012-007-9081-7

    Article  Google Scholar 

  16. 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 2008, 11: 584-591.

    Article  MathSciNet  Google Scholar 

  17. Fetecau C, Hayat T, Fetecau C: Starting solutions for oscillating motions of Oldryd-B fluids in cylindrical domains. Int J Non-Newtonian Fluid Mech 2008, 153: 191-201. 10.1016/j.jnnfm.2008.02.005

    Article  Google Scholar 

  18. 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 2010, 11: 3744-3751. 10.1016/j.nonrwa.2010.02.004

    Article  MathSciNet  Google Scholar 

  19. Anjum A, Ayub M, Khan M: Starting solutions for oscillating motions of an Oldroyd-B fluid over a plane wall. Commun Nonlinear Sci Numer Simul 2012, 17: 472-482. 10.1016/j.cnsns.2011.05.004

    Article  MathSciNet  Google Scholar 

  20. Khan M, Arshad M, Anjum A: On exact solutions of Stokes second problem for MHD Oldroyd-B fluid. Nucl Eng Des 2012, 243: 20-32.

    Article  Google Scholar 

  21. Waters ND, King MJ: Unsteady flow of an elastico-viscous liquid. Rheol Acta 1970, 9: 345-355. 10.1007/BF01975401

    Article  Google Scholar 

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

    MathSciNet  Google Scholar 

  23. 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 2011, 9: 816-825. 10.2478/s11534-010-0073-1

    Google Scholar 

  24. Fetecau C, Prasad SC, Rajagopal KR: A note on the flow induced by a constantly accelerating plate in an Oldryod-B fluid. Appl Math Model 2007, 31: 647-654. 10.1016/j.apm.2005.11.032

    Article  Google Scholar 

  25. Debnath L, Bhatta D (Eds): Integral Transforms and Their Applications. Chapman and Hall, Boca-Raton; 2007.

    Google Scholar 

  26. Sneddon IN: Fourier Transforms. McGraw Hill, New York; 1951.

    Google Scholar 

Download references

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.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Nazish Shahid.

Additional information

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 steady-state 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.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Shahid, N., Rana, M. & Siddique, I. Exact solution for motion of an Oldroyd-B fluid over an infinite flat plate that applies an oscillating shear stress to the fluid. Bound Value Probl 2012, 48 (2012). https://doi.org/10.1186/1687-2770-2012-48

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2012-48

Keywords