Skip to main content

On exact solution of unsteady MHD flow of a viscous fluid in an orthogonal rheometer

Abstract

This paper studies the unsteady MHD flow of a viscous fluid in which each point of the parallel planes are subject to the non-torsional oscillations in their own planes. The streamlines at any given time are concentric circles. Exact solutions are obtained and the loci Γ of the centres of these concentric circles are discussed. It is shown that the motion so obtained gives three infinite sets of exact solutions in the geometry of an orthogonal rheometer in which the above non-torsional oscillations are superposed on the disks. These solutions reduce to a single unique solution when symmetric solutions are looked for. Some interesting special cases are also obtained from these solutions.

Introduction

Berker [1] has defined the ‘pseudo plane motions’ of the first kind that: if the streamlines in a plane flow are contained in parallel planes but the velocity components are dependent on the coordinate normal to the planes. Berker [2] has obtained a class of exact solutions to the Navier-Stokes equations belonging to the above type of flows. The ingenious analysis of Berker [3], [4] of the flow of an incompressible viscous fluid confined between two coaxially or non-coaxially rotating infinite parallel disks has shown the existence of an infinite number of exact solutions characterizing the asymmetric and symmetric motions. The symmetric solutions presented by Berker [3] belong to a new class and cannot be obtained as limiting cases of any of the known solutions. It was observed that these results are relevant in understanding the motions of the fluid in a device called an orthogonal rheometer.

Owing to the importance of exact solutions, many authors [5]–[9] have analyzed the flows in this geometry for different types of fluids. Rao and Kasiviswanathan [10] have investigated the unsteady viscous flow in the orthogonal rheometer and established the existence of three infinite sets of exact solutions. The flow of an elastico-viscous fluid in the above geometry has been investigated by Abbott and Walters [11]. They obtained the unique exact solution of the problem and it was considered as a particular case of the infinite set of solutions which satisfies the condition for symmetric solutions. Thornley [12] has considered an exact solution for the flow of an incompressible fluid confined between two infinite disks when both the fluid and the disk were in rigid-body rotation with constant angular velocity and one of the disks is used to perform non-torsional oscillations in its own plane in a rotating frame of reference.

The MHD flow of Newtonian and non-Newtonian fluids between two non-coaxial rotating disks has been one of the subjects which has attracted attention of many workers. Kasiviswanathan and Rao [13] investigated the unsteady MHD flow of a viscous, homogeneous, incompressible, electrically conducting fluid between two infinite, parallel, insulated disks rotating with the same angular velocity about two non-coincident axes in which the disks were subjected to non-torsional oscillations. Using the Laplace transform two sets of infinite solutions are presented and these solutions reduce to unique solution when symmetric solutions are looked for. Mohanty [14] obtained an exact solution for the viscous, homogeneous, incompressible, conducting fluid between eccentrically rotating disks having the same angular velocity in the presence of uniform transverse magnetic field. Further, Erkman [15] studied the same problem but for an induced magnetic field. Murthy and Ram [16] examined an exact solution for the case of flow due to non-coaxial rotations of a porous disk and a fluid at infinity under the application of a uniform transverse magnetic field. Rao and Rao [17] also considered the same physical situation as discussed by [14] but they considered the second grade fluid. Ersoy [18] obtained an exact solution of an Oldroyd-B fluid between two infinite, parallel, insulated disks rotating about non-coincident axes in the presence of a uniform transverse magnetic field. Moreover, Rajagopal [19] reviewed problems related to both Newtonian and non-Newtonian fluids between two parallel disks rotating about the same or distinct axes with the same and different angular velocities.

In this paper we have generalised the work of Rao and Kasiviswanathan [10] to the case of unsteady MHD flow of a conducting, incompressible viscous fluid between two infinite, parallel, insulated disks rotating with the same angular velocity about non-coincident axes normal to the disks. Each point of the parallel planes is subjected to non-torsional oscillations in their own planes. The streamlines at a given instant are concentric circles in each plane parallel to a fixed plane π 0 and each point of the plane is performing non-torsional oscillations. Exact solutions are obtained and the loci of the centres of concentric circles are discussed. In the geometry of an orthogonal rheometer where non-torsional oscillations are imposed on the disks relative to the basic rotations, the existence of three infinite sets of exact solutions is established. Three cases arise i.e. when: (i) ω>δ, (ii) ω<δ and (iii) ω=δ. These solutions reduce to a single unique solution when symmetric solutions are looked for.

Governing equations

Consider the flow of an incompressible viscous fluid filling the entire three-dimensional space in which the streamlines at a given instant are concentric circles in each plane parallel to a fixed plane π 0 and each point of the plane is performing non-torsional oscillations. A uniform magnetic field B 0 is applied along the z-axis. The unsteady flow of an incompressible fluid in the presence of a magnetic field is governed by the following equations:

ρ D V D t =divT+J×B,
(1)
V=0,
(2)

where V is the velocity field, T the Cauchy stress tensor, J the electric current density, ρ the density and D/Dt the material time derivative. In addition, B (= B 0 +b) is the total magnetic field in which b and B 0 denote the induced magnetic field and imposed magnetic field, respectively. In the absence of displacement currents, the Maxwell equations [20], [21] and the modified Ohm law [22] can be written as

B = 0 , × B = μ m J , × E = B t , J = σ [ E + V × B ] , }
(3)

in which E is the electric field, μ m the magnetic permeability and σ the electric field conductivity. It is assumed that all the thermo-physical fluid properties are constant throughout the flow field region. Also, the induced magnetic field b is negligible compared with the applied magnetic field B 0 . It is further supposed that no energy is added or extracted from the fluid by the electric field which implies that there is no electric field present in the fluid flow region. Under these assumptions the electromagnetic force J×B takes the form

J×B=σ B 0 2 V.
(4)

The velocity field is chosen as follows:

u=u(x,y,z,t),v=v(x,y,z,t),w=0,
(5)

where u, v, w are the components of velocity in the Cartesian coordinate system (x,y,z). For this motion, the streamlines for all time t are contained in planes parallel to the fixed plane π 0 and the plane Oxy is parallel to this fixed plane. The constitutive equation of the viscous fluid is

T=pI+μ A 1 +ρf,
(6)

with p being the pressure, I the unit tensor, A 1 the Rivilin-Ericksen tensor, μ the dynamic viscosity and f the body force. The Rivilin-Ericksen tensor A 1 is defined by

A 1 =L+ L T ,L=V,
(7)

in which T denotes the transpose. The Navier-Stokes equations governing the flow in view of Eqs. (1) and (4)-(7) are given by

ψ y t + ψ y ψ y x ψ x ψ y y = ϕ x +v ( ψ x x + ψ y y + ψ z z ) y σ B 0 2 ψ y ρ ,
(8)
ψ x t + ψ y ψ x x ψ x ψ y x = ϕ y +v ( ψ x x + ψ y y + ψ z z ) x σ B 0 2 ψ x ρ ,
(9)
u= ψ y ,v= ψ x ,
(10)
ϕ= p ρ +Ω,
(11)

where ψ(x,y,z,t) is the stream function, Ω the potential due to external force and v the kinematic coefficient of viscosity. The subscript denotes the partial derivatives with respect to corresponding variables throughout the paper. From the equation of motion in the z-direction one finds ϕ to be independent of z. In the present analysis the authors are interested in getting solutions of an unsteady MHD fluid motion in which the streamlines at any instant are concentric circles in each plane with normal line parallel to the Oz-axis. Therefore, consider

x=f(z,t),y=g(z,t),
(12)

and let these be the loci Γ of the centres of these circles which are fixed points at a given instant t and belong to the instantaneous axis of the vortex. The condition that Γ must not be a straight line perpendicular to the Oxy plane gives the inequality:

f z 2 + g z 2 0.
(13)

Following an analysis similar to Berker [4], the general form of the stream function for this type of unsteady MHD flow is given by

ψ=K [ ( x f ( z , t ) ) 2 + ( y g ( z , t ) ) 2 ] .
(14)

This stream function ψ has to satisfy Eqs. (8) and (9). According to Eq. (10), the velocity components are given by the equations:

u=2K ( y g ( z , t ) ) ,v=2K ( x f ( z , t ) ) .
(15)

Let c be any constant and let π be the plane z=c. Equation (15) shows that at a given time  t 0 , the fluid particles belonging to the plane π move like the points of a solid plate which rotates around a fixed point with coordinates (f(c, t 0 ),g(c, t 0 ),c) and with constant angular velocity ω, given by

ω=2K.
(16)

This angular velocity is same for all planes π which are parallel to the plane Oxy. Thus (15) should be interpreted as a motion due to rigid rotation superposed by some unsteady motion.

Exact solutions

Eliminating ϕ from (8) and (9) by differentiating with respect to z and making use of (14) and (16) we obtain

ν W z z z ( σ B 0 2 ρ + i ω ) W z = W z t ,
(17)

where

W(z,t)=f(z,t)+ig(z,t).
(18)

We investigate the unsteady MHD motion generated by subjecting each of the parallel planes to a non-torsional oscillation of the type e i δ t or e i δ t , where δ is the frequency of oscillation. Letting

W=F(z) e i δ t ,
(19)

Eq. (17) takes the form

F z z z ( N + i ω δ v ) F z =0,
(20)

where

N= σ B 0 2 ρ ν .
(21)

The general solutions of the differential equation (20) are

F ( z ) = A + B exp ( N + 2 i β 1 2 z ) + C exp ( N + 2 i β 1 2 z ) , ω > δ , F ( z ) = A + B exp ( N 2 i β 2 2 z ) + C exp ( N 2 i β 2 2 z ) , ω < δ , F ( z ) = A + B exp ( N z ) + C exp ( N z ) , ω = δ , }
(22)

where

β 1 2 = ( ω δ ) / 2 v , β 2 2 = ( δ ω ) / 2 v , }
(23)

and the arbitrary complex constants A, B and C are assumed to be

A= A 1 +i A 2 ,B= B 1 +i B 2 ,C= C 1 +i C 2 .
(24)

The solutions for non-torsional oscillation of the type e i δ t are obtained by replacing δ by −δ at appropriate places in the above equations. The general solutions of the differential equation (17) in view of Eqs. (18), (19), (22) and (24) are given by

f = f 0 + ( A 1 cos δ t + A 2 sin δ t ) + e β r 1 z [ B 1 cos ( β i 1 z δ t ) B 2 sin ( β i 1 z δ t ) ] f = + e β r 1 z [ C 1 cos ( β i 1 z + δ t ) + C 2 sin ( β i 1 z + δ t ) ] ,
(25)
g = g 0 + ( A 2 cos δ t A 1 sin δ t ) + e β r 1 z [ B 1 sin ( β i 1 z δ t ) + B 2 cos ( β i 1 z δ t ) ] g = + e β r 1 z [ C 2 cos ( β i 1 z + δ t ) C 1 sin ( β i 1 z + δ t ) ] ,
(26)

where f 0 and g 0 denote solutions due to steady rotation and

β r 1 +i β i 1 = N + 2 i β 1 2 .
(27)

Similarly, expressions for ω<δ and ω=δ can be obtained. It is worthy to mention that the solutions (25) and (26) reduce to the solutions presented by Rao and Kasiviswanathan [10] when B 0 =0.

Now we will establish that by making a suitable change of axes, that is, the change that satisfies the condition that the plane Oxy be parallel to the plane π 0 , it is possible to reduce the number of arbitrary constants contained in the curve Γ given in Eqs. (25) and (26). The change to be made is the result of transformations described below for ω>δ.

  1. (i)

    A translation parallel to the plane Oxy and defined by the equations:

    x 1 =x( A 1 cosδt+ A 2 sinδt), y 1 =y( A 2 cosδt A 1 sinδt).
    (28)

    This transformation gets rid of the constants A 1 and A 2 contained in Eqs. (25) and (26).

  2. (ii)

    If B 1 2 + B 2 2 >0 and C 1 2 + C 2 2 >0, a translation parallel to the axis Oz and a rotation about this axis. This translation and rotation is so chosen that the equation of the curve Γ in the new axes reduces to the following form:

    Γ 1 : x = a x 0 + a ( x B cos δ t + y B sin δ t ) , y = a y 0 + a ( y B cos δ t x B sin δ t ) , }
    (29)

    where

    x 0 = cosh ( β 0 z ) cos ( β 0 z + n ) , y 0 = sinh ( β 0 z ) sin ( β 0 z + n ) , x B = cosh ( β i 1 z ) cos ( β i 1 z + n ) , y B = sinh ( β i 1 z ) sin ( β i 1 z + n )
    (30)

    with β 1 , a>0, β 0 = ω / 2 υ >0. Equation (29) contains only two arbitrary real constants a and n.

  3. (iii)

    If B 1 2 + B 2 2 >0 and C 1 = C 2 =0, a translation parallel to the axis Oz. This translation is so chosen that the equations of the curve Γ with respect to the new axes reduce to the following form:

    Γ 2 : x = a x 0 + a exp ( β r 1 z ) cos ( β i 1 z δ t ) , y = a y 0 + a exp ( β r 1 z ) sin ( β i 1 z + δ t ) , }
    (31)

    where

    x 0 =exp( β 0 z)cos( β 0 z), y 0 =exp( β 0 z)sin( β 0 z),
    (32)

    with a>0, β 0 0, β 1 0. Equation (31) involves only one arbitrary real constant a.

  4. (iv)

    If C 1 2 + C 2 2 >0 and B 1 =0= B 2 , again a translation parallel to the axis Oz. This translation is so chosen that the equations of the curve Γ in the new axes reduce to the following form:

    Γ 3 : x = a x 0 a exp ( β r 1 z ) cos ( β i 1 z + δ t ) , y = a y 0 a exp ( β r 1 z ) sin ( β i 1 z δ t ) , }
    (33)

    where

    x 0 =exp( β 0 z)cos( β 0 z), y 0 =exp( β 0 z)sin( β 0 z),
    (34)

    with a>0, β 0 0, β 1 0. Equation (33) again contains only one constant a. It is worth mentioning that the curves Γ 1 - Γ 3 reduce to the curves given by Rao and Kasiviswanathan [10] when B 0 =0.

  5. (v)

    Lastly, when ω=δ the curve Γ by a proper choice of axes can be written as

    Γ 4 : x = ( a z + b z 2 ) cos ( ω t + n ) , y = ( a z + b z 2 ) sin ( ω t + n ) , }
    (35)

    with a,b>0; n is any arbitrary constant. This transformation has no counterpart in the corresponding steady problem. Likewise, curves Γ 5 and Γ 6 , similar to the curves Γ 1 and Γ 2 , can be obtained for ω<δ.

Equations (29) of the curve Γ 1 in cylindrical coordinates can be written in the form

r= a 2 e ( β 0 z ) [ 1 + f 1 ( β 0 , β 1 , z , t ) ] 1 / 2 ,
(36)

where f 1 is negligible for large z or for very small δ. So for large z (36) reduces to

r= a 2 e β 0 z .
(37)

Let C 1 be the projection of the curve Γ 1 on the Oxy plane. The equation of C 1 in polar coordinates is

r= a 2 e θ ,θ= β 0 z.
(38)

Therefore, the curve C 1 tends towards the logarithmic spiral (38) for large value of z, having the properties given by Berker [4]. A similar discussion holds for the curves Γ 2 - Γ 6 . It can be further noted that e β 0 z dominates over the other terms, the curve Γ 1 is wrapped around the surface of revolution S 1 the axis of which is Oz and the meridian (in cylindrical coordinates) of which is given in Eq. (36).

Unsteady MHD flow in an orthogonal rheometer

Consider the unsteady MHD flow of a conducting, incompressible viscous fluid between two infinite, parallel, insulated disks D 1 and D 2 rotating with the same angular velocity ω about non-coincident axes z and z normal to the planes of the disks and performing non-torsional oscillations in their own planes. This model represents the motion of the fluid in the instrument called an orthogonal rheometer.

We choose the Cartesian coordinate system with the z-axis perpendicular to the disks and the lower and the upper disks are located at z=h and z=h, respectively. Let the lower disk rotate about the point P 1 ( x 1 , y 1 ,h) and the upper disk about P 2 ( x 1 , y 1 ,h). Taking the middle point O of P 1 P 2 as the origin and perpendicular axes Ox and Oy lying in the plane z=0 so that the z-axis is perpendicular to these axes. The superposed unsteady motion is created by subjecting each point of the lower and upper disk to non-torsional oscillations of the type a ¯ e i δ t + b ¯ e i δ t and c ¯ e i δ t + d ¯ e i δ t , respectively, where a ¯ , b ¯ , c ¯ and d ¯ , are given complex constants. The governing equation for this motion is given by Eq. (17) whereas the boundary conditions are

W=( x 1 +i y 1 )+ a ¯ e i δ t + b ¯ e i δ t ,on z=h,
(39)
W=( x 1 +i y 1 )+ c ¯ e i δ t + d ¯ e i δ t ,on z=h.
(40)

Assume the solution for Eq. (17) to be of the form

W= a ¯ F 1 (z) e i δ t + b ¯ F 2 (z) e i δ t + F 3 (z),
(41)

where the first two terms correspond to the oscillatory motion and the third term gives the solution due to the basic rigid rotation. The corresponding boundary conditions are

F 1 (h)=1, F 2 (h)=1, F 3 (h)=( x 1 +i y 1 ),
(42)

and

F 1 (h)= c ¯ a ¯ , F 2 (h)= d ¯ b ¯ , F 3 (h)= x 1 +i y 1 .
(43)

The governing differential equation (17) is of the third order, giving three sets of infinite solutions; one corresponding to the steady part and the other two corresponding to the oscillatory parts, proportional to e i δ t and e i δ t . As there are only two boundary conditions (42) and (43) to determine each complex function, therefore its solution is not unique. We introduce one more arbitrary boundary condition in addition to the boundary conditions (42) and (43) in order to determine the solution. Let x p and y p be the coordinates of the point of intersection of the curve Γ and the plane Oxy and also the non-torsional oscillations are prescribed on the plane Oxy. Then the boundary conditions will be of the form

W(0)= x P +i y P + a 0 e i δ t + b 0 e i δ t
(44)

giving

F 1 (0)= a 0 a ¯ , F 2 (0)= b 0 b ¯ , F 3 (0)= x p +i y p .
(45)

The solutions of the differential equation (17), satisfying the boundary conditions (42)-(44), are given by

W = { c ¯ + a ¯ 2 a 0 2 cosh N + λ 3 2 z 1 cosh N + λ 3 2 h 1 + c ¯ a ¯ 2 sinh N + λ 3 2 z sinh N + λ 3 2 h + a 0 } e i δ t + { d ¯ + b ¯ 2 b 0 2 cosh N + λ 1 2 z 1 cosh N + λ 1 2 h 1 + d ¯ b ¯ 2 sinh N + λ 1 2 z sinh N + λ 1 2 h + b 0 } e i δ t + ( x p + i y p ) { 1 cosh N + λ 0 2 z 1 cosh N + λ 0 2 h 1 } + ( x 1 + i y 1 ) sinh N + λ 0 2 z sinh N + λ 0 h ,
(46)

where

λ 0 2 = 2 i β 0 2 , λ 1 2 = 2 i β 1 2 = 2 i β 2 2 , λ 3 2 = 2 i β 3 2 , β 0 2 = ω 2 ν , β 1 2 = ( ω δ ) / 2 ν , β 2 2 = ( δ ω ) / 2 ν , β 3 2 = ( ω + δ ) / 2 ν .
(47)

Setting

c ¯ a ¯ 2 = a 1 + i a 2 , c ¯ + a ¯ 2 a 0 2 = b 1 + i b 2 , a 0 = a 1 + i a 2 , d ¯ b ¯ 2 = c 1 + i c 2 , d ¯ + b ¯ 2 b 0 2 = d 1 + i d 2 , b 0 = b 1 + i b 2 ,
(48)

in Eq. (45) and equating the real and imaginary parts, we get for ω>δ

f = [ ( c 2 R 1 + c 1 Q 1 + d 1 T 1 + d 2 S 1 ) ( a 2 R 3 + a 1 Q 3 + b 2 S 3 + b 1 T 3 ) a 2 + b 2 ] sin δ t + [ ( a 1 R 3 a 2 Q 3 + b 1 S 3 b 2 T 3 ) + ( c 1 R 1 c 2 Q 1 + d 1 S 1 d 2 T 1 ) + a 1 + b 1 ] cos δ t + x 1 R 0 y 1 Q 0 + x p ( 1 S 0 ) + y p T 0
(49)

and

g = [ ( a 1 R 3 a 2 Q 3 + b 1 S 3 b 2 T 3 ) ( c 1 R 1 c 2 Q 1 + d 1 S 1 d 2 T 1 ) + a 1 b 1 ] sin δ t + [ ( a 2 R 3 + a 1 Q 3 + b 2 S 3 + b 1 T 3 ) + ( c 2 R 1 + c 1 Q 1 + d 1 T 1 + d 2 S 1 ) + a 2 + b 2 ] cos δ t + x 1 Q 0 + y 1 R 0 x p T 0 + y p ( 1 S 0 ) ,
(50)

where

sinh N + λ j 2 z sinh N + λ j 2 h = R j ±i Q j , cosh N + λ j 2 z 1 cosh N + λ j 2 h 1 = S j ±i T j ,j=0,1,3.
(51)

The minus sign is taken when j=2. Thus we have

R j = [ ξ j ( z ) ξ j ( h ) + η j ( z ) η j ( h ) ] / Δ 1 j , Q j = [ η j ( z ) ξ j ( h ) ξ j ( z ) η j ( h ) ] / Δ 1 j , S j = [ ϕ j ( z ) ϕ j ( h ) + χ j ( z ) χ j ( h ) ] / Δ 2 j , T j = [ χ j ( z ) ϕ j ( h ) ϕ j ( z ) χ j ( h ) ] / Δ 2 j ,
(52)
ξ j ( z ) = sinh ( β r j z ) cos ( β ι j z ) , η j ( z ) = cosh ( β r j z ) sin ( β ι j z ) , ϕ j ( z ) = cosh ( β r j z ) cos β i j z 1 , χ j ( z ) = sinh ( β r j z ) sin ( β i j z ) , Δ 1 j = ξ j 2 ( h ) + η j 2 ( h ) , Δ 2 j = ϕ j 2 ( h ) + χ j 2 ( h ) ,
(53)

where j=0,1,2,3. Similarly, we can get the solutions for ω<δ. The solution for the case ω=δ is obtained from Eqs. (48) and (49) by taking the limit δω. In this case we obtain

R 1 = R 2 = sinh N z sinh N h , S 1 = S 2 = cosh N z 1 cosh N h 1 , Q 1 = Q 2 = T 1 = T 2 =0.
(54)

We have thus obtained the unsteady MHD motion of viscous fluid in the orthogonal rheometer for the three different cases of ω, i.e., (i) ω>δ, (ii) ω<δ, (iii) ω=δ. The solutions for f and g contain six arbitrary constants, namely a 1 , a 2 , b 1 , b 2 , x p and y p , giving three sets of infinite solutions. Hence the existence of three sets of infinite solutions for the unsteady MHD motion of the fluid in an orthogonal rheometer in which the non-torsional oscillations are superposed on the disks are established. These solutions for f and g reduce to the solutions of Rao and Kasiviswanathan [10] for B 0 0.

Results and discussion

In view of Eqs. (15) and (16) the velocity components are given by

u=ω ( y g ( z , t ) ) ,v=ω ( x f ( z , t ) ) ,w=0.
(55)

Then the components (ξ,η,ζ) of the vorticity are given by

ξ=ω f z ,η=ω g z ,ζ=2ω,
(56)

where f and g are given by Eqs. (25) and (26). Substituting Eqs. (14) and (16) in Eqs. (8) and (9), we get

ϕ x = ω 2 x + ω ( v g z z ω f g t σ B 0 2 ρ g ) + σ B 0 2 ρ ω y , ϕ y = ω 2 y ω ( v f z z + ω g f t σ B 0 2 ρ f ) σ B 0 2 ρ ω x .
(57)

Using Eqs. (11), (25) and (26) in (56) we obtain

p + ρ Ω = ρ ω 2 2 ( 1 σ 2 B 0 4 ρ 2 ω 2 ) [ ( x X 0 ) 2 ( 1 + σ B 0 2 ρ ω ) + ( y Y 0 ) 2 ( 1 σ B 0 2 ρ ω ) ] + p 0 , B 0 2 < ρ ω σ , p + ρ Ω = ρ ω 2 2 ( σ 2 B 0 4 ρ 2 ω 2 1 ) [ ( y Y 0 ) 2 ( σ B 0 2 ρ ω 1 ) ( x + X 0 ) 2 ( σ B 0 2 ρ ω + 1 ) ] + p 0 , B 0 2 > ρ ω σ , p + ρ Ω = ρ ω 2 [ ( x X 0 ) 2 4 Y 0 y ] + p 0 , σ 2 B 0 4 ρ 2 ω 2 = 1 ,
(58)

where p 0 is an arbitrary constant, X 0 (h,t), Y 0 (h,t) are given by

X 0 = x p ( 1 + ϕ 0 ( h ) / Δ 20 ) + y p χ 0 ( h ) / Δ 20 X 0 = [ ( 1 + δ ω ) { ( b 2 χ 3 ( h ) + b 1 ϕ 3 ( h ) ) / Δ 23 a 1 } X 0 = + ( 1 δ ω ) { ( d 1 ϕ 1 ( h ) + d 2 χ 1 ( h ) ) / Δ 21 b 1 } ] cos δ t X 0 = [ ( 1 + δ ω ) { ( b 1 χ 3 ( h ) b 2 ϕ 3 ( h ) ) / Δ 23 + a 2 } X 0 = + ( 1 δ ω ) { ( d 2 ϕ 1 ( h ) d 1 χ 1 ( h ) ) / Δ 21 b 2 } ] sin δ t X 0 = N ν ω [ y p ( 1 ϕ 0 ( h ) / Δ 20 ) x p χ 0 ( h ) X 0 = + { ( b 2 ϕ 3 ( h ) b 1 χ 3 ( h ) ) / Δ 23 + ( d 1 χ 1 ( h ) d 2 ϕ 1 ( h ) ) / Δ 21 a 2 b 2 } cos δ t X 0 = + { ( b 1 ϕ 3 ( h ) + b 2 χ 3 ( h ) ) / Δ 23 ( d 1 ϕ 1 ( h ) + d 2 χ 1 ( h ) ) / Δ 21 a 1 + b 1 } sin δ t ] ,
(59)
Y 0 = y p ( 1 + ϕ 0 ( h ) / Δ 20 ) x p χ 0 ( h ) / Δ 20 Y 0 = + [ ( 1 + δ ω ) { ( b 1 χ 3 ( h ) b 2 ϕ 3 ( h ) ) / Δ 23 + a 2 } Y 0 = + ( 1 δ ω ) { ( d 1 χ 1 ( h ) d 2 ϕ 1 ( h ) ) / Δ 21 + b 2 } ] cos δ t Y 0 = + [ ( 1 + δ ω ) { a 1 ( b 2 χ 3 ( h ) + b 1 ϕ 3 ( h ) ) / Δ 23 } Y 0 = + ( 1 δ ω ) { ( d 1 ϕ 1 ( h ) + d 2 χ 1 ( h ) ) / Δ 21 b 1 } ] sin δ t Y 0 = + N ν ω [ x p ( 1 + ϕ 0 ( h ) / Δ 20 ) y p χ 0 ( h ) / Δ 20 Y 0 = + { ( b 1 ϕ 3 ( h ) + b 2 χ 3 ( h ) ) / Δ 23 + ( d 1 ϕ 1 ( h ) + d 2 χ 1 ( h ) ) / Δ 21 a 1 b 1 } cos δ t Y 0 = + { ( d 2 ϕ 1 ( h ) d 1 χ 1 ( h ) ) / Δ 21 ( b 2 ϕ 3 ( h ) b 1 χ 3 ( h ) ) / Δ 23 + a 2 b 2 } sin δ t ] .
(60)

The results (57)-(59) reduce to the results of Rao and Kasiviswanathan [10] in the absence of MHD and Berker [4] in the absence of MHD and oscillations.

The components of stress t 1 applied by the fluid at any point (x,y,h) on the upper disk in (x,y,z) directions are given by

t 1 x =μ( u z + w x ), t 1 y =μ( ν z + w y ), t 1 z =p2μ w z ,
(61)

where p is the pressure and μ the coefficient of viscosity. Using the values of (u,v,w) given in Eq. (54), we obtain from Eq. (60)

t 1 x =μω g z (t,h), t 1 y =μω f z (t,h), t 1 z =p,
(62)

where f and g are defined in Eqs. (48) and (49). Similarly, the stress applied on the lower plate at any point (x,y,h) can be calculated.

A number of interesting special cases can be derived from the general solutions given by Eq. (45). A solution is said to be symmetric solution with respect to a given point O, if the velocity field is a symmetric vector field with respect to the point O. Therefore, the symmetric solution has to satisfy the following condition:

V(x,y,z)=V(x,y,z).
(63)

Among the solutions obtained for the motion in the orthogonal rheometer, we need to find whether there exist solutions which are symmetric with respect to: (i) the point O, (ii) the mid point of the line segment P 1 P 2 . The velocity field given by Eq. (54) satisfies the condition (62), provided

f(z,t)=f(z,t)andg(z,t)=g(z,t),
(64)

which implies that f and g are odd functions of z and

f(0)=0,g(0)=0.
(65)

The latter equation shows that the curve Γ must pass through the origin O of the axes. Therefore, for a symmetric solution the constants in Eq. (45) should satisfy

c ¯ + a ¯ 2 a 0 = d ¯ + b ¯ 2 b 0 =0, a 0 = b 0 = x p = y p =0.
(66)

Thus, for a symmetric solution all the arbitrary constants given in the boundary condition (44) vanish and the solution thus obtained is unique. This symmetric solution is given by

W = { c ¯ a ¯ 2 sinh N + λ 3 2 z sinh N + λ 3 2 h } e i δ t + { d ¯ b ¯ 2 sinh N + λ 1 2 z sinh N + λ 1 2 h } e i δ t + ( x 1 + i y 1 ) sinh N + λ 0 2 z sinh N + λ 0 h .
(67)

This solution includes the symmetric solution of Kasiviswanathan and Rao [13] as a special case. Taking y 1 =0 and y P =0 in Eq. (45), we obtain the solutions of the MHD flow due to non-torsionally oscillating eccentric rotating disks given by Kasiviswanathan and Rao [13]. In the limit B 0 0 in Eq. (45), we obtain solutions given by Rao and Kasiviswanathan [10]. By taking B 0 0, δ0 (together with the amplitudes of the oscillations) in Eq. (45), we get

W(z)=( x 1 +i y 1 ) sinh λ 0 z sinh λ 0 h +( x P +i y P ) [ 1 cosh λ 0 z 1 cosh λ 0 z 1 ] ,
(68)

which corresponds to the solutions given by Berker [4]. As δ0, B 0 0, ω0 in Eq. (45), we get

W ( z ) = [ c ¯ + a ¯ 2 a 0 2 + d ¯ + b ¯ 2 b 0 2 ( x P + i y P ) ] z 2 h 2 + [ c ¯ a ¯ 2 + d ¯ b ¯ 2 + x 1 + i y 1 ] z h + ( x P + i y P + a 0 + b 0 ) ,
(69)

and the corresponding boundary conditions are

W ( h ) = x 1 + i y 1 + c ¯ + d ¯ , W ( h ) = ( x 1 + i y 1 ) + a ¯ + b ¯ , W ( 0 ) = x P + i y P + a 0 + b 0 .
(70)

This case corresponds to the flow of a fluid confined between two plates moving with velocities of different magnitudes and direction. It reduces to the solutions given by Kasiviswanathan and Rao [13] for y P =0 and y 1 =0. Since in the limit ω0, all the velocity components given in Eq. (15) become zero; therefore we cannot take this velocity field for the above problem. However, in order to understand the nature and physical significance of this case, we choose the velocity field given by Kasiviswanathan and Rao [13]. We have

u=f(z),v=g(z),w=0.
(71)

In the absence of a magnetic field, Eq. (1) in view of Eq. (70) yields

μ f z z p x =0,μ g z z p y =0, p z =0.
(72)

Eliminating the pressure gradient we obtain

W z z z =0.
(73)

The boundary conditions are given by

W ( h ) = u 1 + i u 1 = U 1 , W ( h ) = u 2 + i u 2 = U 2 , W ( 0 ) = u 3 + i u 3 = U 3 ,
(74)

where U 3 is chosen arbitrarily. The solution of Eq. (72) subject to the boundary conditions (73) yields

W(z)= 1 2 [ U 1 + U 2 2 U 3 ] z 2 h 2 + 1 2 [ U 1 U 2 ] z h + U 3 ,
(75)

which is similar to the solution given in Eq. (68). The solution given by the above equation yields a set of infinite number of solutions as U 3 is chosen arbitrarily. For U 1 = U 2 and U 3 =0, Eq. (74) reduces to a simple shearing motion and the solution so obtained is a symmetric solution. The stream function for this flow is given by

ψ(z)=yf(z)xg(z)+C,
(76)

and the streamlines are given by a family of parallel straight lines in each of the plane where z=constant. In Eq. (75), C denotes an arbitrary constant.

Conclusion

In this study we have considered the unsteady MHD flow of a viscous fluid in which each point of the parallel planes is subjected to the non-torsional oscillations in their own planes. Exact solutions for the described model are presented. The major findings of the present work are summarised as follows:

The existence of three sets of infinite numbers of solutions for the unsteady MHD motion of a viscous fluid in an orthogonal rheometer in which the non-torsional oscillations are superposed on the disks is established. These solutions reduce to a single unique solution when symmetric solutions are looked for.

The results of Rao and Kasiviswanathan [10] are recovered for B 0 =0 and Berker [4] in the absence of MHD and oscillations.

The present work also includes the solutions of Kasiviswanathan and Rao [13] as a special case.

References

  1. Berker, RA: Sur Quelques cas d’ intéegration des équations du mouvement d’un fluide visqueux incompressible. Théses françaises de l’entre-deux-guerres (1936)

    Google Scholar 

  2. Berker R: Intégration des équations du mouvement d’un fluide visqueux incompressible. In Handbuch der Physik. Springer, Berlin; 1963:1-384.

    Google Scholar 

  3. Berker R: A new solution of the Navier-Stokes equation for the motion of a fluid contained between two parallel plates rotating about the same axis. Arch. Mech. Stosow. 1979, 31: 265-280.

    MathSciNet  Google Scholar 

  4. Berker R: An exact solution of the Navier-Stokes equation. The vortex curvilinear axis. Int. J. Eng. Sci. 1982, 20: 217-230. 10.1016/0020-7225(82)90017-9

    Article  MathSciNet  Google Scholar 

  5. Rajagopal KR, Gupta AS: Flow and stability of second grade fluid between rotating planes. Arch. Mech. 1981, 33: 663-674.

    MathSciNet  Google Scholar 

  6. Rajagopal KR, Gupta AS: Flow and stability of a second grade fluid between two parallel plates rotating about noncoincident axes. Int. J. Eng. Sci. 1981, 19: 1401-1409. 10.1016/0020-7225(81)90037-9

    Article  MathSciNet  Google Scholar 

  7. Rajagopal KR: On the flow of a simple fluid in an orthogonal rheometer. Arch. Ration. Mech. Anal. 1982, 79: 39-47. 10.1007/BF02416565

    Article  Google Scholar 

  8. Siddiqui AM, Rana MA, Ahmed N: Effects of Hall current and heat transfer on MHD flow of a Burgers’ fluid due to a pull of eccentric rotating disks. Commun. Nonlinear Sci. Numer. Simul. 2008, 13: 1554-1570. 10.1016/j.cnsns.2006.10.005

    Article  MathSciNet  Google Scholar 

  9. Siddiqui AM, Rana MA, Ahmed N: Magnetohydrodynamics flow of a Burgers’ fluid in an orthogonal rheometer. Appl. Math. Model. 2010, 34: 2881-2892. 10.1016/j.apm.2009.12.022

    Article  MathSciNet  Google Scholar 

  10. Rao AR, Kasiviswanathan SR: On exact solutions of the unsteady Navier-Stokes equations-the vortex with instantaneous curvilinear axis. Int. J. Eng. Sci. 1987, 25: 337-349. 10.1016/0020-7225(87)90040-1

    Article  MathSciNet  Google Scholar 

  11. Abbott TNG, Walters K: Rheometrical flow systems, part 2. Theory for the orthogonal rheometer, including an exact solution of the Navier-Stokes equations. J. Fluid Mech. 1970, 40: 205-213. 10.1017/S0022112070000125

    Article  Google Scholar 

  12. Thornley C: On Stokes and Rayleigh layers in a rotating system. Q. J. Mech. Appl. Math. 1968, 21: 451-461. 10.1093/qjmam/21.4.451

    Article  Google Scholar 

  13. Kasiviswanathan SR, Rao AR: On exact solutions of unsteady MHD flow between eccentrically rotating disks. Arch. Mech. 1987, 39: 411-418.

    Google Scholar 

  14. Mohanty HK: Hydromagnetic flow between two rotating disks with noncoincident parallel axes of rotation. Phys. Fluids 1972, 15: 1456-1458. 10.1063/1.1694107

    Article  Google Scholar 

  15. Erkman CS: Magnetohydrodynamic flow of a viscous, incompressible fluid between rotating parallel plates. Lett. Appl. Eng. Sci. 1975, 3: 51-59.

    Google Scholar 

  16. Murthy SN, Ram RKP: MHD flow and heat transfer due to eccentric rotations of a porous disk and a fluid at infinity. Int. J. Eng. Sci. 1978, 16: 943-949. 10.1016/0020-7225(78)90053-8

    Article  Google Scholar 

  17. Rao AR, Rao PR: On the magnetohydrodynamic flow between eccentrically rotating disks. Int. J. Eng. Sci. 1983, 21: 359-372. 10.1016/0020-7225(83)90120-9

    Article  Google Scholar 

  18. Ersoy HV: MHD flow of an Oldroyd-B fluid between eccentric rotating disks. Int. J. Eng. Sci. 1999, 37: 1973-1984. 10.1016/S0020-7225(99)00010-5

    Article  Google Scholar 

  19. Rajagopal KR: Flow of viscoelastic fluids between rotating plates. Theor. Comput. Fluid Dyn. 1992, 3: 185-206. 10.1007/BF00417912

    Article  Google Scholar 

  20. Sherchiff J: A Textbook of Magnetohydrodynamics. Pergamon, Oxford; 1965.

    Google Scholar 

  21. Davidson PA: An Introduction to Magnetohydrodynamics. Cambridge University Press, Cambridge; 2001.

    Book  Google Scholar 

  22. Sutton GW, Sherman A: Engineering Magnetohydrodynamics. McGraw-Hill, New York; 1965.

    Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sadia Siddiqa.

Additional information

Competing interests

The authors have no competing interests.

Authors’ contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Authors’ original file for figure 7

Authors’ original file for figure 8

Authors’ original file for figure 9

Authors’ original file for figure 10

Authors’ original file for figure 11

Authors’ original file for figure 12

Authors’ original file for figure 13

Authors’ original file for figure 14

Authors’ original file for figure 15

Authors’ original file for figure 16

Authors’ original file for figure 17

Authors’ original file for figure 18

Authors’ original file for figure 19

Authors’ original file for figure 20

Authors’ original file for figure 21

Authors’ original file for figure 22

Authors’ original file for figure 23

Authors’ original file for figure 24

Authors’ original file for figure 25

Authors’ original file for figure 26

Authors’ original file for figure 27

Authors’ original file for figure 28

Authors’ original file for figure 29

Authors’ original file for figure 30

Authors’ original file for figure 31

Authors’ original file for figure 32

Authors’ original file for figure 33

Authors’ original file for figure 34

Authors’ original file for figure 35

Authors’ original file for figure 36

Authors’ original file for figure 37

Authors’ original file for figure 38

Authors’ original file for figure 39

Authors’ original file for figure 40

Authors’ original file for figure 41

Authors’ original file for figure 42

Authors’ original file for figure 43

Authors’ original file for figure 44

Authors’ original file for figure 45

Authors’ original file for figure 46

Authors’ original file for figure 47

Authors’ original file for figure 48

Authors’ original file for figure 49

Authors’ original file for figure 50

Authors’ original file for figure 51

Authors’ original file for figure 52

Authors’ original file for figure 53

Authors’ original file for figure 54

Authors’ original file for figure 55

Authors’ original file for figure 56

Authors’ original file for figure 57

Authors’ original file for figure 58

Authors’ original file for figure 59

Authors’ original file for figure 60

Authors’ original file for figure 61

Authors’ original file for figure 62

Authors’ original file for figure 63

Authors’ original file for figure 64

Authors’ original file for figure 65

Authors’ original file for figure 66

Authors’ original file for figure 67

Authors’ original file for figure 68

Authors’ original file for figure 69

Authors’ original file for figure 70

Authors’ original file for figure 71

Authors’ original file for figure 72

Authors’ original file for figure 73

Authors’ original file for figure 74

Authors’ original file for figure 75

Authors’ original file for figure 76

Authors’ original file for figure 77

Authors’ original file for figure 78

Authors’ original file for figure 79

Authors’ original file for figure 80

Authors’ original file for figure 81

Authors’ original file for figure 82

Authors’ original file for figure 83

Authors’ original file for figure 84

Authors’ original file for figure 85

Authors’ original file for figure 86

Authors’ original file for figure 87

Authors’ original file for figure 88

Authors’ original file for figure 89

Authors’ original file for figure 90

Authors’ original file for figure 91

Authors’ original file for figure 92

Authors’ original file for figure 93

Authors’ original file for figure 94

Authors’ original file for figure 95

Authors’ original file for figure 96

Authors’ original file for figure 97

Authors’ original file for figure 98

Authors’ original file for figure 99

Authors’ original file for figure 100

Authors’ original file for figure 101

Authors’ original file for figure 102

Authors’ original file for figure 103

Authors’ original file for figure 104

Authors’ original file for figure 105

Authors’ original file for figure 106

Authors’ original file for figure 107

Authors’ original file for figure 108

Authors’ original file for figure 109

Authors’ original file for figure 110

Authors’ original file for figure 111

Authors’ original file for figure 112

Authors’ original file for figure 113

Authors’ original file for figure 114

Authors’ original file for figure 115

Authors’ original file for figure 116

Authors’ original file for figure 117

Authors’ original file for figure 118

Authors’ original file for figure 119

Authors’ original file for figure 120

Authors’ original file for figure 121

Authors’ original file for figure 122

Authors’ original file for figure 123

Authors’ original file for figure 124

Authors’ original file for figure 125

Authors’ original file for figure 126

Authors’ original file for figure 127

Authors’ original file for figure 128

Authors’ original file for figure 129

Authors’ original file for figure 130

Authors’ original file for figure 131

Authors’ original file for figure 132

Authors’ original file for figure 133

Authors’ original file for figure 134

Authors’ original file for figure 135

Authors’ original file for figure 136

Authors’ original file for figure 137

Authors’ original file for figure 138

Authors’ original file for figure 139

Authors’ original file for figure 140

Authors’ original file for figure 141

Authors’ original file for figure 142

Authors’ original file for figure 143

Authors’ original file for figure 144

Authors’ original file for figure 145

Authors’ original file for figure 146

Authors’ original file for figure 147

Authors’ original file for figure 148

Authors’ original file for figure 149

Authors’ original file for figure 150

Authors’ original file for figure 151

Authors’ original file for figure 152

Authors’ original file for figure 153

Authors’ original file for figure 154

Authors’ original file for figure 155

Authors’ original file for figure 156

Authors’ original file for figure 157

Authors’ original file for figure 158

Authors’ original file for figure 159

Authors’ original file for figure 160

Authors’ original file for figure 161

Authors’ original file for figure 162

Authors’ original file for figure 163

Authors’ original file for figure 164

Authors’ original file for figure 165

Authors’ original file for figure 166

Authors’ original file for figure 167

Authors’ original file for figure 168

Authors’ original file for figure 169

Authors’ original file for figure 170

Authors’ original file for figure 171

Authors’ original file for figure 172

Authors’ original file for figure 173

Authors’ original file for figure 174

Authors’ original file for figure 175

Authors’ original file for figure 176

Authors’ original file for figure 177

Authors’ original file for figure 178

Authors’ original file for figure 179

Authors’ original file for figure 180

Authors’ original file for figure 181

Authors’ original file for figure 182

Authors’ original file for figure 183

Authors’ original file for figure 184

Authors’ original file for figure 185

Authors’ original file for figure 186

Authors’ original file for figure 187

Authors’ original file for figure 188

Authors’ original file for figure 189

Authors’ original file for figure 190

Authors’ original file for figure 191

Authors’ original file for figure 192

Authors’ original file for figure 193

Authors’ original file for figure 194

Authors’ original file for figure 195

Authors’ original file for figure 196

Authors’ original file for figure 197

Authors’ original file for figure 198

Authors’ original file for figure 199

Authors’ original file for figure 200

Authors’ original file for figure 201

Authors’ original file for figure 202

Authors’ original file for figure 203

Authors’ original file for figure 204

Authors’ original file for figure 205

Authors’ original file for figure 206

Authors’ original file for figure 207

Authors’ original file for figure 208

Authors’ original file for figure 209

Authors’ original file for figure 210

Authors’ original file for figure 211

Authors’ original file for figure 212

Authors’ original file for figure 213

Authors’ original file for figure 214

Authors’ original file for figure 215

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rana, M.A., Siddiqa, S. & Noor, S. On exact solution of unsteady MHD flow of a viscous fluid in an orthogonal rheometer. Bound Value Probl 2014, 146 (2014). https://doi.org/10.1186/s13661-014-0146-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-014-0146-y

Keywords