Numerical Solution for a Similar Flow between Two Disks in the Presence of a Magnetic Field


Numerical solutions are obtained for non-steady, incompressible fluid flow between two parallel disks which at time t are separated by a distance H(1-αt)1/2 and a magnetic field proportional to B0(1-αt) -1/2 is applied perpendicular to the disks where H denotes a representative length, BO denotes a representative magnetic field and α-1 denotes a representative time. Similarity transformations are used to convert the governing partial differential equations of motion in to ordinary differential form. The resulting ordinary differential equations are solved numerically using SOR method, Richardson extrapolation and Simpson’s (1/3) Rule. Our numerical scheme is straightforward, efficient and easy to program.

Share and Cite:

Hussain, S. , Kamal, M. , Ahmad, F. , Ali, M. , Shafique, M. and Hussain, S. (2013) Numerical Solution for a Similar Flow between Two Disks in the Presence of a Magnetic Field. Applied Mathematics, 4, 1163-1167. doi: 10.4236/am.2013.48155.

1. Introduction

The quest for similar solutions is particularly important with respect to the mathematical character of the solution. In cases where similar solutions exist, it is possible, to reduce the system of partial differential equations to one involving ordinary differential equations, which evidently constitute, a considerable mathematical simplifycation of the problem. Wang [1] studied a viscous fluid between two parallel plates, which are being squeezed or separated with normal velocity proportional to and found similarity solutions of the unsteady Navier-Stocks equations. Ishizawa [2] derived a similarity solution to the case of the unsteady laminar flow between two parallel disks. Tichy and Bourgin [3] found that a similarity solution does exist for the steady flow in a narrow channel of a gap width varying as, where and are constants. Bhupendra et al. [4] considered the problem of forced flow of an electrically conducting viscous incompressible fluid due to an infinite rotating disk under the influence of uniform magnetic field, applied normal to the flow. Pavlov [5] found an exact similarity solution of MHD boundary layer equations for the steady two dimensional flow of an electrically conducting incompressible fluid due to rotation of a plane elastic surface in the presence of a uniform transverse magnetic field. Guria et al. [6] obtained exact solution of hyderomagnetic flow between two porous disks rotating with same angular velocity about two non coincident axes in the presence of a uniform transverse magnetic field. Attia [7] studied the problem of steady flow and heat transfer of a conducting fluid due to the rotation of an infinite, non conducting porous disk in the presence of an external magnetic field. Sajid et al. [8] examined the MHD rotating flow of a viscous fluid over a shrinking sheet. Asir et al. [9] gave a new hybrid analytical algorithm to study the effects of uniform suction of a laminar, steady, incompressible magnetohyderodynamic electrically conducting fluid over a rotating disk. The purpose of present study is to obtain numerical solution for similar flows of a Newtonian fluid between two disks in the presence of a magnetic field. Usha and Vasudevan [10] studied a similar flow between two rotating disks in the presence of a magnetic field and obtained rather expensive solution of the problem to observe the effect of flow parameters on the velocity fluid.

2. Mathematical Analysis

It has been assumed the flow is axisymmetric, incomepressible and non-steady. The flow is between two parallel infinite disks, which are separated a distance apart, where denotes time. A magnetic field of strength is applied perpendicular to the two disks. The upper disk is moving with velocity towards the lower fixed disk. Cylindrical polar coordinates are used. The lower disk is at and the upper one is at where.

The equations of motion in component form become as follows:





where the subscripts denote the partial differentiation with respect to space coordinates is the density, the pressure and the coefficient of kinematics viscosity.

The boundary conditions are:


The following similarity transformations are used:


is the dimensionless variable. The equation of continuity is identically satisfied. Equations (2) and (3) take the forms below respectively.



where, , here denotes viscosity and denotes fluid electrical conductivity. Whence the differentiation of Equation (6) with respect and that of Equation (7) with respect to r yield:


Equation (8) is integrated to get:


where is constant of integration.

The boundary conditions in dimensionless form become;


3. Finite-Difference Equations

In order to solve Equation (9) numerically, we let

, (11)

and we obtain


The boundary conditions (10) become as:


The derivatives involved in Equation (12) are approximated by central difference approximation at a typical point of interval [0,1], we get.


where denotes a grid size and Equation (11) is integrated numerically. Also the symbols used denote and

4. Computational Procedure

Finite difference Equation (14) and the first order ordinary differential Equation (11) are solved simultaneously by using SOR method Smith [11, p.262] and Simpson’s (1/3) rule Gerald [12, p.293] with the formula given in Milne [13, p.48] respectively subject to the appropriate boundary conditions.

The order of the sequence of iterations is as follows:

1) The Equation (14) for the solution of is solved subject to the boundary conditions (13).

2) For the solution of we use the computed values of from above step in to Equation (11) and integrate by Simpson’s (1/3) rule.

3) The optimum value of the relaxation parameter is estimated to accelerate the convergence of the SOR method.

4) The SOR procedure is terminated when the following criterion is satisfied for q:

where denotes the number of iterations and stands for q.

The above steps 1 to 4 are repeated for higher grid levels and. The SOR procedure gives the solution of of order of accuracy due to second order finite differences used to approximate the derivatives while Simpson’s (1/3) rule gives the order of accuracy in the solution of. Higher order accuracy in the solution of on the basis of above solutions is achieved by using Richardson’s Extrapolation Burden [14, p.168]. The solution of order of accuracy in the following tables for computation of is the most accurate and accepted solution.

5. Results and Discussion

The numerical solutions examine the way in which the flow pattern changes with the squeeze Reynolds number Re and Hartmann number M. Hamza [15] investigated this problem for the range, 0.0 < M ≤ 30.0. We extended the previous work and analyzed the problem for the parameters involved in the range and.

Table 1 shows the values of different parameters used in the numerical procedure. The numerical results for have been computed for different values of flow parameters namely Re and. The accuracy of numerical results is checked by comparing the results on three different grid sizes namely h = 0.025, 0.012 and 0.006. The comparison of is shown in the Tables 2- 4 using Richardson extrapolation method. Graphically, the results have been demonstrated in Figures 1-6. It is found that for fixed Re, there is a slight increase in the

Table 1. Optimum value of relaxation parameter used in SOR method.

Table 2. M = 40.0, Re = 0.01, M = 50.0, Re = 0.01.

Table 3. M = 10.0, Re = 10.0, M = 0.0, Re = 15.0.

Table 4. M = 10.0, Re = 15.0, M = 5.0, Re = 20.0.

Figure 1. Graph of for different values of M when Re = 0.01.

Figure 2. Graph of for different values of M when Re = 1.0.

value of near the disks and a slight decrease in the region of the mid plane with increase in. This in-

Figure 3. Graph of for different values of M when Re = 5.0.

Figure 4. Graph of for different values of M when Re = 10.0.

crease and decrease become more prominent with more increase in, also the radial velocity profiles become more flat in the interior region for all values of Re. On

Figure 5. Graph of for different values of M when Re = 15.0.

Figure 6. Graph of for different values of M when Re = 20.0.

the other hand, for fixed values of, the magnitude of radial flow decreases with increase in Re. The value of is determined by hit and trial to satisfy the approximate zero-order perturbation results given in [15] as.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] C. Y. Wang, “The Squeezing of a Fluid between Two Plates,” Journal of Applied Mechanics, Vol. 43, No. 4, 1976, pp. 579-583. doi:10.1115/1.3423935
[2] S. Ishizawa, “The Unsteady Laminar Flow between Two Parallel Disks with Arbitrarily Varying Gap Width,” Bul letin of JSME, Vol. 9, No. 35, 1966, pp. 533-550. doi:10.1299/jsme1958.9.533
[3] J. A. Tichy and P. Bourgin, “A Similarity Solution for Flow in a Narrow Channel of Varying Gap,” ASME Jour nal of Applied Mechanics, Vol. 53, No. 4, 1986, pp. 943-946. doi:10.1115/1.3171885
[4] K. S. Bhupendera, K. J. Abhay and R. C. Chaudhary, “MHD Forced Flow of a Conducting Viscous Fluid through a Porous Medium Induced by an Impervious Ro tating Disk Ram,” Journal of Physics, Vol. 52, No. 1, 2007, pp. 73-84.
[5] K. B. Pavlov, “Magnetohydrodynamic Flow of an In compressible, Viscous Fluid Caused by the Deformation of a Plane Surface,” Magnitnaya Gidrodinamika, Vol. 4, 1974, pp. 146-147.
[6] M. Guria, B. K. Das, R. N. Jana and C. E. Imark, “Hy deromagnetic Flow between Two Porous Disks Rotating about Non Coincident Axes,” Acta Mechanica Sinica, Vol. 24, No. 5, 2008, pp. 489-496. doi:10.1007/s10409-008-0158-x
[7] H. A. Attia, “On the Effectiveness of the Ion Slip on the Steady Flow of a Conducting Fluid Due to a Porous Ro tating Disk with Heat Transfer,” Tamkang Journal of Science and Engineering, Vol. 9, No. 3, 2006, pp. 185-213.
[8] M. Sajid, T. Javed and T. Hayat, “MHD Rotating Flow of a Viscous Fluid over a Shrinking Surface,” Nonlinear Dynamics, Vol. 51, No. 1-2, 2008, pp. 259-265. doi:10.1007/s11071-007-9208-3
[9] B. El-Asir, A. Mansoor, S. Bataineh and R. Arar, “A New Hybrid Analytical Analysis of the Magnetohydrodynamic Flow over a Rotating Disk under Uniform Suction,” Jour nal of Applied Sciences, Vol. 6, No. 5, 2006, pp. 1059-1065. doi:10.3923/jas.2006.1059.1065
[10] R. Usha and S. Vasudevan, “A Similar Flow between Two Rotating Disks in the Presence of a Magnetic Field,” Journal of Applied Mechanics, Vol. 60, No. 3, 1993, pp. 707-714. doi:10.1115/1.2900862
[11] G. D. Smith, “Numerical Solution of Partial Differential Equation,” Clarendon Press, Oxford, 1979.
[12] C. F. Gerald, “Applied Numerical Analysis,” Addison Wesley Pub., New York, 1989.
[13] W. E. Milne, “Numerical Solution of Differential Equa tion,” Dover Pub., Dover, 1970.
[14] R. L. Burden, “Numerical Analysis,” Prindle, Weber & Schmidt, Boston, 1985.
[15] E. A. Hamza, “A Similar Flow between Two Disks in the Presence of a Magnetic Field,” IC/87/287, Trieste, 1987.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.