Approximate Analysis of MHD Squeeze Flow between Two Parallel Disks with Suction or Injection by Homotopy Perturbation Method
An analysis has been performed to study magneto-hydrodynamic (MHD) squeeze flow between two parallel infinite disks where one disk is impermeable and the other is porous with either suction or injection of the fluid. We investigate the combined effect of inertia, electromagnetic forces, and suction or injection. With the introduction of a similarity transformation, the continuity and momentum equations governing the squeeze flow are reduced to a single, nonlinear, ordinary differential equation. An approximate solution of the equation subject to the appropriate boundary conditions is derived using the homotopy perturbation method (HPM) and compared with the direct numerical solution (NS). Results showing the effect of squeeze Reynolds number, Hartmann number and the suction/injection parameter on the axial and radial velocity distributions are presented and discussed. The approximate solution is found to be highly accurate for the ranges of parameters investigated. Because of its simplicity, versatility and high accuracy, the method can be applied to study linear and nonlinear boundary value problems arising in other engineering applications.
This paper deals with the study of magneto-hydrodynamic (MHD) squeeze flow of an electrically conducting fluid between two infinite, parallel disks. The lower disk is stationary and permeable with (suction or injection). The upper disk is impermeable and moves toward the lower disk with a specified time dependent velocity. The use of a MHD fluid in lubrication prevents the adverse impact of temperature on the fluid viscosity when the system operates under extreme conditions. The problem considered is of general interest in the theory of magneto-hydrodynamic lubrication and other related applications. In particular, the results of the present investigation are directly applicable to the hydrodynamics of high temperature bearings lubricated with liquid metals. A number of theoretical and experimental investigations into magneto-hydrodynamic effects in lubrication have been reported. These include among other works of Hughes and Elco , Kuzma et al.  and Krieger et al. . These authors considered the electromagnetic force term in the Navier-Stokes equations but neglected some or all the inertia terms. When low viscosity lubricants are used to reduce energy losses in lubrication devices, inertial effects become important and must be included especially if the squeeze Reynolds number is not small.
In the present work we investigate the combined effect of inertia, electromagnetic forces, and surface suction or injection in a squeeze film between two parallel disks. This combination of effects in squeezing flow has not been studied previously. The plates are made of a nonconducting material. There is no externally applied electric field and the induced electric field is negligible. Since the problem defies an exact analytical solution, special techniques must be used to derive approximate analytical solution.
One such technique is the homotopy perturbation method (HPM) proposed and applied by He [4–6]. His pioneering work has prompted many authors [7–20] to apply HPM to solve a wide variety of homogeneous and nonhomogeneous linear and nonlinear problems. Over the past several years the accuracy of HPM has been repeatedly verified. In view of the success of HPM reported by many researchers, we were tempted to adopt this method to solve the present problem.
2. The Homotopy Perturbation Method
The basic idea embodied in the HPM and a brief summary of the method can be found in [4–6]. The convergence and stability of this method have been established in . The latest developments of the method can be found in a series of papers by pioneering researchers such as He [21–23] and Yıldırım [24–30]. Other papers that bear close relevance to the present work are those of Mahmood et al. , Mehmood and Ali , and Z. Z. Ganji and D. D. Ganji .
3. Solution Procedure
We consider axisymmetric incompressible flow between two parallel infinite disks, which at time , are spaced a distance apart and a magnetic field proportional to is applied perpendicular to the disks. The upper disk at is moving with velocity toward the stationary lower disk at . The axial coordinate is denoted by and the radial coordinate by . With the axial and radial velocities denoted by and , respectively, we introduce the following quantities:
into the continuity and momentum equations governing the flow. The continuity equation is automatically satisfied and the momentum equation is reduced to the following fourth-order, nonlinear, ordinary differential equation: where primes denote differentiation with respect to the similarity variable , , and . Here denotes the density, is the absolute viscosity, is kinematic viscosity, is the fluid electrical conductivity, is the magnetic field acting normal to the plates, is the distance separating the plates at time and a constant and has the units of .
The boundary conditions are given by where A is a constant such that corresponds to suction and to injection at the lower stationary disk.
According to the HPM, we can construct a homotopy of (2) as follows: where primes denote differentiation with respect to η and asterisks have been omitted for convenience.
We consider a three term-solution for as follows:
Assuming and substituting f from (5) into (4) and after some algebraic manipulation, we obtain the following set of equations: We utilized the symbolic algebra package Maple 11 to solve (6) with the appropriate boundary conditions and obtained the following analytical solutions
Once is known in an analytical form, the definition of in (1) can be used to calculate the axial component of the velocity, . By differentiating with respect to , which is straightforward and can be easily performed by Maple 11, we use the definition of in (1) to compute the radial component of the velocity.
The axial and radial velocities are each functions of parameters S, M and A. Since the axial and radial velocities are proportional to and , respectively, they suffices to present the results for . These results are displayed in Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, and 12 for twelve different combinations of parameters S, M, and A. Figures 1–8 show that for both the cases of no injection and injection , the axial component of the velocity increases monotonically as increases. The velocity profile is not significantly affected by the increase in the fluid electrical conductivity and/or the magnetic field, that is, the increase in M from 0 to 5. A careful examination of Figures 1–8 also reveals that the effect of squeeze Reynolds number S on the axial velocity profiles is minimal.
Figures 9–12 show the effect of suction on the axial velocity distributions. In each case, the axial velocity decreases monotonically as increases. Siddiqui et al.  have observed the same pattern for the axial (normal in their case) velocity in the study of two-dimensional MHD squeeze flow between parallel plates without injection or suction. As in the case for injection (Figures 1–8), the velocity patterns are only minimally influenced by the changes in parameters M and S.
The results for which represents the radial component of the velocity appear in Figures 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, and 24. Figures 13–20 show the results for the radial velocity with no injection and injection for different combination values of M and S. In each case, the radial velocity increases as increases, reaches a peak value in the neighborhood of 0.5 and then decays to zero around . For the case of suction , which is illustrated in Figures 21–24, the radial velocity becomes negative and decreases as increases. The trough in radial velocity is attained in the neighborhood of . After that the velocity increases and reaches the zero value around .
Figure 25 has been prepared to illustrate the effect of varying the imposed magnetic field (Hartmann number M ) on the radial distribution of velocity for suction and injection . As the strength of the magnetic field is increased, that is, as the Hartmann number increases, the impact is felt most on the maximum radial velocity, more so for injection than for suction. The effect of squeeze Reynolds number S on the radial velocity distribution is shown in Figure 26 for a fixed value of . It can be seen that even with a hundred-fold increase in S, the radial velocity profile is not affected significantly.
For every case investigated (Figures 1–26), the HPM predictions have been compared with the corresponding direct numerical solutions (NSs) obtained by using Maple 11 software. This software uses a Fehlberg fourth-fifth order Runge-Kutta finite-difference method for the numerical solution of the boundary value problem . In each case, the HPM and NS are found to be consistently close. This consistent closeness strongly vouches for the accuracy of HPM. Equation (9) therefore provides simple and highly accurate analytical solutions for the problem studied in this paper and is useful for rapid engineering calculations.
He’s homotopy perturbation method (HPM) has been utilized to derive approximate analytical solutions for the radial and axial velocity distributions in magneto-hydrodynamic (MHD) squeeze flow between two parallel infinite disks where one disk is impermeable and the other is porous with either suction or injection of the fluid. The approximate solutions have been compared with the direct numerical solutions generated by the symbolic algebra package Maple 11 which uses a Fehlberg fourth-fifth order Runge-Kutta finite-difference method for solving nonlinear boundary value problems. The comparison showed that the HPM solutions are highly accurate and provide a rapid means of computing the flow velocities between the plates.
For both the cases of no injection and injection, the axial component of the velocity increases monotonically as the similarity variable increases. The velocity profiles are not significantly affected by the increase in the fluid electrical conductivity and/or the magnetic field. Similarly, the effect of squeeze Reynolds number S on the axial velocity profiles is minimal. The effect of imposed magnetic field on the maximum radial velocity is more pronounced for injection than for suction.
The authors express their gratitude to the reviewers for their valuable suggestions which were incorporated to enhance the value of the paper.
R. J. Krieger, H. J. Day, and W. F. Hughes, “The MHD hydrostatics thrust bearings—theory and experiments,” ASME Journal of Lubrication Technology, vol. 89, pp. 307–313, 1967.View at: Google Scholar
D. D. Ganji and A. Sadighi, “Application of He's homotopy-perturbation method to nonlinear coupled systems of reaction-diffusion equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 4, pp. 411–418, 2006.View at: Google Scholar
P. D. Ariel, T. Hayat, and S. Asghar, “Homotopy perturbation method and axisymmetric flow over a stretching sheet,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 4, pp. 399–406, 2006.View at: Google Scholar
D. D. Ganji and A. Sadighi, “Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 24–34, 2007.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
A. Yıldırım, “An algorithm for solving the fractional nonlinear Schrödinger equation by means of the homotopy perturbation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, no. 4, pp. 445–451, 2009.View at: Google Scholar
A. Yıldırım, “The homotopy perturbation method for solving the modified Korteweg-de Vries equation,” Zeitschrift für Naturforschung A, vol. 63a, no. 10-11, pp. 621–626, 2008.View at: Google Scholar
A. Yıldırım, “Exact solutions of nonlinear differential-difference equations by He's homotopy perturbation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 9, no. 2, pp. 111–114, 2008.View at: Google Scholar
M. Mahmood, M. A. Hossain, S. Asghar, and T. Hayat, “Application of homotopy perturbation method to deformable channel with wall suction and injection in a porous medium,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 9, no. 2, pp. 195–206, 2008.View at: Google Scholar
A. A. Mehmood and A. Ali, “An application of He's homotopy perturbation method in fluid mechanics,” International Journal of Nonlinear Sciences and Numerical Simulations, vol. 10, no. 2, pp. 239–246, 2009.View at: Google Scholar
Z. Z. Ganji and D. D. Ganji, “Approximate solutions of thermal boundary-layer problems in a semi-infinite flat plate by using He's homotopy perturbation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 9, no. 4, pp. 415–422, 2008.View at: Google Scholar