Numerical and Soft Computing Methods for Characteristic Value Problems of ODE and ODEs Systems
View this Special IssueResearch Article  Open Access
Inayat Ullah, Hamid Khan, M. T. Rahim, "Approximation of First Grade MHD Squeezing Fluid Flow with Slip Boundary Condition Using DTM and OHAM", Mathematical Problems in Engineering, vol. 2013, Article ID 816262, 9 pages, 2013. https://doi.org/10.1155/2013/816262
Approximation of First Grade MHD Squeezing Fluid Flow with Slip Boundary Condition Using DTM and OHAM
Abstract
A steady twodimensional axisymmetric flow of an incompressible viscous fluid under the influence of a uniform transverse magnetic field with slip boundary condition is studied. An ordinary nonlinear differential equation is formed by transforming the NavierStokes equations using the transformation . Differential transform and optimal homotopy analysis methods have been used to obtain the solutions by varying pertinent flow parameters. By using residuals in each case, the validity of solutions is established. Excellent results are obtained by using the proposed schemes. The influence of different parameters on the flow is shown through graphs.
1. Introduction
Squeezing flows induce when normal stresses or vertical velocities are applied externally by means of a mobile boundary [1]. Stefan [2] carried out the fundamental research in this field. Engineers studied the analysis of Newtonian fluid squeezed between two infinite planar plates [3, 4]. There are many applications of squeezing flow in food industry, especially in chemical engineering [4]. Polymer processing, compression, and injection modeling are the practical examples of squeezing flow. Thin Newtonian liquid films squeezing between two plates were studied by Grimm [5]. Squeezing flow under the influence of magnetic field is widely applied to bearing with liquidmetal lubrication [6–9]. Islam et al. [10] studied squeezing fluid flow between the two infinite parallel plates in a porous medium channel.
In case of many polymeric liquids when the weight of molecule is high, then they show slip at the boundary. The noslip boundary condition is not applicable in this case. In many cases such as thin film problems, rarefied fluid problems, fluids containing concentrated suspensions, and flow on multiple interfaces, the noslip boundary condition fails to work. Navier [11], for the first time, proposed the general boundary condition which demonstrates the fluid slip at the surface. The difference of fluid velocity and velocity of the boundary is proportional to the shear stress at that boundary. The proportionality constant is named the slip parameter having length as its dimension.
The slip condition is of great importance especially when fluids with elastic character are under consideration [12]. Newtonian fluid was considered by Ebaid [13] to study the effects of magnetic field and wall slip conditions on the peristaltic transport in an asymmetric channel. It has great importance in medical sciences, particularly in polishing artificial heart valves and internal cavities in many manufactured parts achieved by embedding such fluids with abrasives [14]. The influence of slip on the peristaltic motion of thirdorder fluid in asymmetric channel is studied by Hayat et al. [15]. The effects of slip condition on the rotating flow of a third grade fluid in a nonporous medium are investigated by Hayat and Abelman [16]. The extension of the work in [16] to a porous medium and obtaining the numerical solutions for the steady magnetohydrodynamics flow of a third grade fluid in a rotating frame is presented by Abelman et al. [17].
Various perturbation techniques can be applied to solve nonlinear differential equations [18]. The limitations of these methods are based on the assumption of small parameter and there is no proper way to select this parameter. The idea of homotopy was combined with perturbation in the last decade. Liao and He did the fundamental work on it. Homotopy analysis method (HAM) was proposed by Liao [19–21]. Homotopy perturbation method was introduced by He [22–24]. Marinca et al. [25, 26] introduced OHAM for approximate solution of nonlinear problems of thin film flow of a fourthgrade fluid down a vertical cylinder and for the study of the behavior of nonlinear mechanical vibration of electrical machines. It is observed that HPM and HAM are the special cases of OHAM [27, 28].
Differential transform method (DTM) was initially introduced by Zhou in 1986 [29]. Çatal [30–33] studied solution of free vibration equations of beam on elastic soil, buckling analysis of partially embedded pile in elastic soil, analysis of free vibration of beam on elastic soil, and response of forced EulerBernoulli beams by using DTM. Ayaz [34] applied this method to differential algebraic equations (DAEs) of index1. Liu and Song [35] investigated that DTM is effective in case of index2 DAEs but not in index3. Ayaz [36] studied the applications of twodimensional DTM in case of partial differential equations.
The objective of this research paper is to use OHAM and DTM for studying MHD squeezing flow with slip boundary condition between two infinite plates approaching each other slowly. In Section 2, the basic equations are derived and the model is reduced to a nonlinear boundary value problem. In Sections 3 and 4, the basic ideas of DTM and its use for our problem are presented while the idea of OHAM and its application to our problem are in Sections 5 and 6, respectively.
2. Problem Formulation
Let us consider, in the presence of a magnetic field, a squeezing flow of an incompressible Newtonian fluid with constant density and viscosity , squeezed between two large planar parallel plates separated by a small distance approaching each other with a low constant velocity (Figure 1). Assume that the flow is quasisteady [1, 3], and the NavierStokes equations governing such flow when inertial terms are retained are where denotes the material time derivative, is the Cauchy stress tensor given by with . is the total magnetic field given by . and represent the imposed and induced magnetic fields, respectively. The modified Ohm’s law and Maxwell’s equations (see [37] and the references therein), in the absence of displacement currents, are Here is the electric current density, represents the electrical conductivity, the electric field, and the magnetic permeability. If , , and are constant, is negligible as compared to , is perpendicular to so that the Reynold number is small with no electric field in the fluid flow region and then the magneto hydrodynamic force involved can be written as Assuming that the plates are nonconducting and the magnetic field is applied along the axis. The gap distance between the plates changes slowly with time for small values of the velocity so that it can be taken as constant. An axisymmetric flow in cylindrical coordinates with axis perpendicular to plates and at the plates. For axial symmetry, is represented by . In view of negligible body forces with no tangential velocity, NavierStokes equations [1, 6, 10] in cylindrical coordinates are where Introducing the stream function , we have, Eliminating from (4), we have where Using the transformation , (7) can be written as subject to the slip boundary conditions The nondimensional parameters are , , and . Omitting the *, (9) and (10) become with and , are the Reynolds and Hartmann numbers, respectively.
3. Differential Transform Method
Onedimensional differential transform of a function is defined as follows [29–35]: is also called the function. The inverse transform of is is defined as Combining (12) and (13), we can write Usually we express in a finite series, that is, which implies that is negligibly small.
Some fundamental theorems on onedimensional differential transform are as follows.
Theorem 1. If , then .
Theorem 2. If , then .
Theorem 3. If , then .
Theorem 4. If , then
4. Analysis of Differential Transform Method
For fourthorder boundary value problem [10], with boundary conditions The differential transform is given by with transformed boundary conditions Using (20) and (21), for is evaluated to obtain the solution up to as follows:
5. Application of DTM to Our Problem
Here, we consider an incompressible Newtonian fluid, squeezed between two large planar, parallel plates which are separated by a small distance and moving towards each other with velocity . For small values of , the gap distance between the plates changes slowly with time , so that it can be taken as constant, and the flow is steady as with slip boundary conditions The differential transform of (23) is given by with transformed boundary conditions where and are unknowns to be determined later [38]. Using (25) and (26), the values of , are As , have been nondimensionalized, so let us fix them one to find the values of and from the following system: which leads us to the following values: Using (22), the approximate solution of the BVP is as follows: This is a series solution and we note that the terms approach to zero if we go forward which is the beauty of DTM and other series solutions.
6. Basic Idea of OHAM
Let us consider a boundary value problem [25–27] along with boundary conditions where denotes the independent variable, , are unknown and known functions and , , are linear, nonlinear, and boundary operators, respectively.
According to OHAM a homotopy, satisfying for and , is constructed with nonzero auxiliary function for with and for an unknown function . For the equalities and hold. So it is concluded that the solution approaches from to as varies from to . For , We choose the auxiliary function such that where are the convergence controlling constants to be determined. Expanding in Taylor’s series about to get the approximation substituting (36) into (33), and equating the coefficients of like powers of , we get the following.
Zerothorder problem is Firstorder problem is Secondorder problem is Generally, with boundary conditions where is the coefficient of in the expansion of about as The convergence of the series (36) depends upon . For convergence at , the thorder approximation is Substituting (43) into (31), the expression for the residual is will be the exact solution if , but in case of nonlinearity, it does not happen generally. To find the optimal values of the constants , there are different methods to apply. We apply the method of least square as follows: Minimizing this function, we have where and are chosen in the domain of the problem for locating the suitable . The approximate solution is well determined for these values of .
7. Application of OHAM to Our Problem
Using , , and in (11), we have the following.
Zerothorder problem is Firstorder problem is Secondorder problem is Thirdorder problem is Fourthorder problem is Fifthorder problem is By considering the fifthorder solution, we have The residual of the problem is In order to find , , we apply the method of least square as follows: Solving (56) for , , we get Using these values of , , the approximate solution is given by: As , the components of vanish when going onward. In Table 1, the solutions obtained by OHAM and DTM along with their residual are presented. Table 2 is constructed to discuss residuals for different values of .


8. Numerical Results and Discussion
Here we analyzed, under a transverse magnetic field, the squeezing fluid flow between two parallel infinite plates. DTM and OHAM are applied to obtain the solution of reduced nonlinear boundary value problem. The convergence region can easily be adjusted using OHAM as compared to DTM in this case. Differential transform method which is the straight forward application of Taylor’s series suffers in general and particularly in initial value problems. Beside this fact, DTM leads easily to closed form solution. OHAM does not require discretization like other numerical and approximate methods. This technique requires less computational effort and the convergence of this method to the exact is fast enough. For computational purposes, we use Mathematica software.
In order to get an insight into the problem, it is important to discuss the influence of different flow parameters on the velocity profile of the fluid. Figure 2 represents the influence of Reynold’s number on the velocity of the fluid. It is observed from Figure 2 that by increasing the velocity of the fluid decreases. The influence of Hartmann number on the velocity of the fluid is shown in Figure 3. It is concluded that the velocity increases by increasing in case of . However, from Figure 4, the velocity decreases due to an increase in for . This observation is in agreement with the literature (see [39, Figure 5]). Figure 5 shows the influence of the slip parameter on the fluid velocity. By increasing , it is observed that the velocity of the fluid increases.
References
 T. C. Papanastasiou, G. C. Georgiou, and A. N. Alexandrou, Viscous Fluid Flow, CRC Press, 1994.
 M. J. Stefan, “Versuch Uber die scheinbare adhesion,” Akademie Der Wissenschaften in Wien, vol. 69, article 713, 1874. View at: Google Scholar
 Q. K. Ghori, M. Ahmed, and A. M. Siddiqui, “Application of homotopy perturbation method to squeezing flow of a newtonian fluid,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 179–184, 2007. View at: Google Scholar
 X. J. Ran, Q. Y. Zhu, and Y. Li, “An explicit series solution of the squeezing flow between two infinite parallel plates,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 1, pp. 119–132, 2007. View at: Publisher Site  Google Scholar
 R. J. Grimm, “Squeezing flows of Newtonian liquid films an analysis including fluid inertia,” Applied Scientific Research, vol. 32, no. 2, pp. 149–166, 1976. View at: Publisher Site  Google Scholar
 W. F. Hughes and R. A. Elco, “Magnetohydrodynamic lubrication flow between parallel rotating disks,” Journal of Fluid Mechanics, vol. 13, pp. 21–32, 1962. View at: Publisher Site  Google Scholar  MathSciNet
 S. Kamiyama, “Inertia Effects in MHD hydrostatic thrust bearing,” Transactions ASME, vol. 91, pp. 589–596, 1969. View at: Google Scholar
 E. A. Hamza, “Magnetohydrodynamic squeeze film,” Journal of Tribology, vol. 110, no. 2, pp. 375–377, 1988. View at: Google Scholar
 S. Bhattacharyya and A. Pal, “Unsteady MHD squeezing flow between two parallel rotating discs,” Mechanics Research Communications, vol. 24, no. 6, pp. 615–623, 1997. View at: Google Scholar
 S. Islam, H. Khan, I. A. Shah, and G. Zaman, “An axisymmetric squeezing fluid flow between the two infinite parallel plates in a porous medium channel,” Mathematical Problems in Engineering, vol. 2011, Article ID 349803, 10 pages, 2011. View at: Publisher Site  Google Scholar
 C.L.M.H. Navier, “Sur les lois de l’équilibre et du mouvement des corps solides élastiques,” Bulletin des Sciences par la Société Philomatique de Paris, pp. 177–181, 1823. View at: Google Scholar
 C. le Roux, “Existence and uniqueness of the flow of secondgrade fluids with slip boundary conditions,” Archive for Rational Mechanics and Analysis, vol. 148, no. 4, pp. 309–356, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Ebaid, “Effects of magnetic field and wall slip conditions on the peristaltic transport of a Newtonian fluid in an asymmetric channel,” Physics Letters A, vol. 372, no. 24, pp. 4493–4499, 2008. View at: Publisher Site  Google Scholar
 L. J. Rhooades, R. Resnic, T. O' Bradovich, and S. Stegman, “Abrasive flow machining of cylinder heads and its positive effects on performance and cost characteristics,” Tech. Rep., Dearborn, Mich, USA, 1996. View at: Google Scholar
 T. Hayat, M. U. Qureshi, and N. Ali, “The inuence of slip on the peristaltic motion of third order fluid in an asymmetric channel,” Physics Letters A, vol. 372, pp. 2653–2664, 2008. View at: Google Scholar
 T. Hayat and S. Abelman, “A numerical study of the influence of slip boundary condition on rotating flow,” International Journal of Computational Fluid Dynamics, vol. 21, no. 1, pp. 21–27, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Abelman, E. Momoniat, and T. Hayat, “Steady MHD flow of a third grade fluid in a rotating frame and porous space,” Nonlinear Analysis: Real World Applications, vol. 10, no. 6, pp. 3322–3328, 2009. View at: Publisher Site  Google Scholar
 A. H. Nayfeh, Problems in Perturbation, A WileyInterscience Publication, John Wiley & Sons, New York, NY, USA, 1985. View at: MathSciNet
 S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, 1992.
 S. Liao, Beyond Perturbation. Introduction to the Homotopy Analysis Method, vol. 2 of CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2004. View at: MathSciNet
 S. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 12, pp. 87–88, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695–700, 2005. View at: Publisher Site  Google Scholar
 J.H. He, “Comparison of homotopy perturbation method and homotopy analysis method,” Applied Mathematics and Computation, vol. 156, no. 2, pp. 527–539, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. Marinca, N. Herişanu, and I. Nemeş, “Optimal homotopy asymptotic method with application to thin film flow,” Central European Journal of Physics, vol. 6, no. 3, pp. 648–653, 2008. View at: Publisher Site  Google Scholar
 V. Marinca and N. Herişanu, “Application of optimal homotopy asymptotic Method for solving nonlinear equations arising in heat transfer,” International Communications in Heat and Mass Transfer, vol. 35, no. 6, pp. 710–715, 2008. View at: Publisher Site  Google Scholar
 J. Ali, S. Islam, S. Islam, and G. Zaman, “The solution of multipoint boundary value problems by the optimal homotopy asymptotic method,” Computers & Mathematics with Applications, vol. 59, no. 6, pp. 2000–2006, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. Marinca and N. Herisanu, “Optimal homotopy perturbation method for strongly nonlinear differential equations,” Nonlinear Science Letters A, vol. 1, no. 3, pp. 273–280, 2010. View at: Google Scholar
 J. K. Zhou, Deferential Transformation and Its Application for Electrical Circuits, Huazhong University Press, Wuhan, China, 1986.
 S. Çatal, “Solution of free vibration equations of beam on elastic soil by using differential transform method,” Applied Mathematical Modelling, vol. 32, no. 9, pp. 1744–1757, 2008. View at: Publisher Site  Google Scholar
 S. Çatal and H. H. Çatal, “Buckling analysis of partially embedded pile in elastic soil using differential transform method,” Structural Engineering and Mechanics, vol. 24, no. 2, pp. 246–269, 2006. View at: Google Scholar
 S. Çatal, “Analysis of free vibration of beam on elastic soil using differential transform method,” Structural Engineering and Mechanics, vol. 24, no. 1, pp. 51–63, 2006. View at: Google Scholar
 S. Catal, “Response of forced EulerBernoulli beams using differential transform method,” Structural Engineering and Mechanics, vol. 42, no. 1, pp. 95–119, 2012. View at: Google Scholar
 F. Ayaz, “Applications of differential transform method to differentialalgebraic equations,” Applied Mathematics and Computation, vol. 152, no. 3, pp. 649–657, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 H. Liu and Y. Song, “Differential transform method applied to high index differentialalgebraic equations,” Applied Mathematics and Computation, vol. 184, no. 2, pp. 748–753, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Ayaz, “On the twodimensional differential transform method,” Applied Mathematics and Computation, vol. 143, no. 23, pp. 361–374, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. R. Mohyuddin and T. Götz, “Resonance behaviour of viscoelastic fluid in Poiseuille flow in the presence of a transversal magnetic field,” International Journal for Numerical Methods in Fluids, vol. 49, no. 8, pp. 837–847, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. H. Chang and I. L. Chang, “A new algorithm for calculating onedimensional differ ential transform of nonlinear functions,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 799–808, 2008. View at: Publisher Site  Google Scholar
 S. Islam, M. Ullah, G. Zaman, and M. Idrees, “Approximate solutions to MHD squeezing fluid flow,” Journal of Applied Mathematics & Informatics, vol. 29, no. 56, pp. 1081–1096, 2011. View at: Google Scholar  MathSciNet
Copyright
Copyright © 2013 Inayat Ullah et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.