SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research

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

Muhammad Afzal Rana1, Sadia Siddiqa2* and Saima Noor3

Author Affiliations

1 Department of Basic Sciences, Riphah International University, Sector I-14, Islamabad 44000, Pakistan

2 Department of Mathematics, COMSATS Institute of Information Technology, Attock, Pakistan

3 Department of Mathematics, COMSATS Institute of Information Technology, Abbottabad, Pakistan

For all author emails, please log on.

Boundary Value Problems 2014, 2014:146  doi:10.1186/s13661-014-0146-y

The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2014/1/146

Received:15 November 2013
Accepted:28 May 2014
Published:9 August 2014

© 2014 Rana et al.; licensee Springer

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


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.

viscous fluid; MHD flow; orthogonal rheometer; eccentric rotation; exact solutions


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 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 and each point of the plane is performing non-torsional oscillations. A uniform magnetic field 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:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M1">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M2">View MathML</a>

where V is the velocity field, T the Cauchy stress tensor, J the electric current density, ρ the density and the material time derivative. In addition, B () is the total magnetic field in which b and 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M3">View MathML</a>

in which E is the electric field, 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 . 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 takes the form

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M4">View MathML</a>

The velocity field is chosen as follows:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M5">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M6">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M7">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M8">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M9">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M10">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M11">View MathML</a>

where 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M12">View MathML</a>

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 plane gives the inequality:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M13">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M14">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M15">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M16">View MathML</a>

This angular velocity is same for all planes π which are parallel to the plane . 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M17">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M18">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M19">View MathML</a>

Eq. (17) takes the form

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M20">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M21">View MathML</a>

The general solutions of the differential equation (20) are

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M22">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M23">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M24">View MathML</a>

The solutions for non-torsional oscillation of the type 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M25">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M26">View MathML</a>

where and denote solutions due to steady rotation and

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M27">View MathML</a>

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 .

Now we will establish that by making a suitable change of axes, that is, the change that satisfies the condition that the plane be parallel to the plane , 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 .

(i) A translation parallel to the plane and defined by the equations:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M28">View MathML</a>

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

(ii) If and , 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:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M29">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M30">View MathML</a>

with , , . Equation (29) contains only two arbitrary real constants a and n.

(iii) If and , 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:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M31">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M32">View MathML</a>

with , , . Equation (31) involves only one arbitrary real constant a.

(iv) If and , 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:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M33">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M34">View MathML</a>

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

(v) Lastly, when the curve Γ by a proper choice of axes can be written as

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M35">View MathML</a>

with ; n is any arbitrary constant. This transformation has no counterpart in the corresponding steady problem. Likewise, curves and , similar to the curves and , can be obtained for .

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M36">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M37">View MathML</a>

Let be the projection of the curve on the plane. The equation of in polar coordinates is

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M38">View MathML</a>

Therefore, the curve 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 -. It can be further noted that dominates over the other terms, the curve is wrapped around the surface of revolution 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 and rotating with the same angular velocity ω about non-coincident axes and 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 and , respectively. Let the lower disk rotate about the point and the upper disk about . Taking the middle point O of as the origin and perpendicular axes Ox and Oy lying in the plane 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 and , respectively, where , , and , are given complex constants. The governing equation for this motion is given by Eq. (17) whereas the boundary conditions are

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M39">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M40">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M41">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M42">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M43">View MathML</a>

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 and . 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 and be the coordinates of the point of intersection of the curve Γ and the plane and also the non-torsional oscillations are prescribed on the plane . Then the boundary conditions will be of the form

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M44">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M45">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M46">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M47">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M48">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M49">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M50">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M51">View MathML</a>

The minus sign is taken when . Thus we have

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M52">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M53">View MathML</a>

where . 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M54">View MathML</a>

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 , , , , and , 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 .

Results and discussion

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M55">View MathML</a>

Then the components of the vorticity are given by

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M56">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M57">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M58">View MathML</a>

where is an arbitrary constant, , are given by

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M59">View MathML</a>


<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M60">View MathML</a>

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 applied by the fluid at any point on the upper disk in directions are given by

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M61">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M62">View MathML</a>

where f and g are defined in Eqs. (48) and (49). Similarly, the stress applied on the lower plate at any point 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:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M63">View MathML</a>

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 . The velocity field given by Eq. (54) satisfies the condition (62), provided

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M64">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M65">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M66">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M67">View MathML</a>

This solution includes the symmetric solution of Kasiviswanathan and Rao [13] as a special case. Taking and 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 in Eq. (45), we obtain solutions given by Rao and Kasiviswanathan [10]. By taking , (together with the amplitudes of the oscillations) in Eq. (45), we get

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M68">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M69">View MathML</a>

and the corresponding boundary conditions are

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M70">View MathML</a>

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 and . Since in the limit , 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M71">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M72">View MathML</a>

Eliminating the pressure gradient we obtain

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M73">View MathML</a>

The boundary conditions are given by

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M74">View MathML</a>

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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M75">View MathML</a>

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 is chosen arbitrarily. For and , 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2014/1/146/mathml/M76">View MathML</a>

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


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

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.


  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)

  2. Berker, R: Intégration des équations du mouvement d’un fluide visqueux incompressible. Handbuch der Physik, pp. 1–384. Springer, Berlin (1963)

  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.. 31, 265–280 (1979)

  4. Berker, R: An exact solution of the Navier-Stokes equation. The vortex curvilinear axis. Int. J. Eng. Sci.. 20, 217–230 (1982). Publisher Full Text OpenURL

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

  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.. 19, 1401–1409 (1981). Publisher Full Text OpenURL

  7. Rajagopal, KR: On the flow of a simple fluid in an orthogonal rheometer. Arch. Ration. Mech. Anal.. 79, 39–47 (1982). Publisher Full Text OpenURL

  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.. 13, 1554–1570 (2008). Publisher Full Text OpenURL

  9. Siddiqui, AM, Rana, MA, Ahmed, N: Magnetohydrodynamics flow of a Burgers’ fluid in an orthogonal rheometer. Appl. Math. Model.. 34, 2881–2892 (2010). Publisher Full Text OpenURL

  10. Rao, AR, Kasiviswanathan, SR: On exact solutions of the unsteady Navier-Stokes equations-the vortex with instantaneous curvilinear axis. Int. J. Eng. Sci.. 25, 337–349 (1987). Publisher Full Text OpenURL

  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.. 40, 205–213 (1970). Publisher Full Text OpenURL

  12. Thornley, C: On Stokes and Rayleigh layers in a rotating system. Q. J. Mech. Appl. Math.. 21, 451–461 (1968). Publisher Full Text OpenURL

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

  14. Mohanty, HK: Hydromagnetic flow between two rotating disks with noncoincident parallel axes of rotation. Phys. Fluids. 15, 1456–1458 (1972). Publisher Full Text OpenURL

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

  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.. 16, 943–949 (1978). Publisher Full Text OpenURL

  17. Rao, AR, Rao, PR: On the magnetohydrodynamic flow between eccentrically rotating disks. Int. J. Eng. Sci.. 21, 359–372 (1983). Publisher Full Text OpenURL

  18. Ersoy, HV: MHD flow of an Oldroyd-B fluid between eccentric rotating disks. Int. J. Eng. Sci.. 37, 1973–1984 (1999). Publisher Full Text OpenURL

  19. Rajagopal, KR: Flow of viscoelastic fluids between rotating plates. Theor. Comput. Fluid Dyn.. 3, 185–206 (1992). Publisher Full Text OpenURL

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

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

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