Research Article  Open Access
Using Differential Transform Method and Padé Approximant for Solving MHD Flow in a Laminar Liquid Film from a Horizontal Stretching Surface
Abstract
The purpose of this study is to approximate the stream function and temperature distribution of the MHD flow in a laminar liquid film from a horizontal stretching surface. In this paper DTMPadé method was used which is a combination of differential transform method (DTM) and Padé approximant. The DTM solutions are only valid for small values of independent variables. Comparison between the solutions obtained by the DTM and the DTMPadé with numerical solution (fourthorder Runge–Kutta) revealed that the DTMPadé method is an excellent method for solving MHD boundarylayer equations.
1. Introduction
As the researches indicated, the nonlinear equations are one of the most important phenomena across the world. Nonlinear phenomena have important effects on applied mathematics, physics, and issues related to engineering. Then the variation of each parameter depends on different factors. The importance of obtaining the exact or approximate solutions of nonlinear partial differential equations (NLPDEs) in physics and mathematics is still a big problem that needs new methods to discover new exact or approximate solutions. Most of nonlinear equations do not have a precise analytic solution; so numerical methods have largely been used to handle these equations. There are also some analytic techniques for nonlinear equations. Some of the classic analytic methods are Lyapunov’s artificial small parameter method [1], perturbation techniques [2–11], and δexpansion method [12]. In the recent years, many authors mainly had paid attention to study solutions of nonlinear partial differential equations by using various methods. Among these are the Adomian decomposition method (ADM) [13, 14], tanh method, homotopy perturbation method (HPM), sinhcosh method, HAM, the DTM, and variational iteration method (VIM) [15, 16].
Magnetohydrodynamics (MHD) is the study of the interaction of conducting fluids with electromagnetic phenomena. The flow of an electrically conducting fluid in the presence of a magnetic field is of importance in various areas of technology and engineering such as MHD power generation, MHD flow meters, and MHD pumps [17–22]. The viscous flow due to stretching boundary is important in extrusion processes where sheet material is pulled out of an orifice with increasing velocity. If the boundary velocity is linear with respect to a fixed point, exact solutions of the NavierStokes equations may be obtained [23, 24].
In recent years the analysis of fluid flow across a thin liquid film has attracted the attention of a number of researchers because of its possible applications in many branches of science and technology. The knowledge of flow and heat transfer within a thin liquid film is crucial in understanding the coating process and design of various heat exchangers and chemical processing equipments. Other applications include wire and fiber coating, food stuff processing, reactor fluidization, and transpiration cooling. The prime aim in almost every extrusion applications is to maintain the surface quality of the extrudate. All coating processes demand a smooth glossy surface to meet the requirements for best appearance and optimum service properties such as low friction, transparency, and strength. The problem of extrusion of thin surface layers needs special attention to gain some knowledge for controlling the coating product efficiently.
The motivation of this letter is the use of the differential transform method and Padé approximant to construct analytical approximate solutions of the Heat transfer in a liquid film over an unsteady stretching surface with viscous dissipation in presence of external magnetic field. The concept of differential transform method was first introduced by Zhou [25] in 1986 and it was used to solve both linear and nonlinear initial value problems in electric circuit analysis. The main advantage of this method is that it can be applied directly to NLPDEs without requiring linearization, discretization, or perturbation. It is a semianalyticalnumerical technique that formulizes Taylor series in a very different manner. This method constructs, for differential equations, an analytical solution in the form of a polynomial. Not like the traditional highorder Taylor series method that requires symbolic computation, the DTM is an iterative procedure for obtaining Taylor series solutions. Another important advantage is that this method reduces the size of computational work while the Taylor series method is computationally taken long time for large orders. This method is well addressed in [26–33].
2. Mathematical Formulation
2.1. Governing Equations and Boundary Conditions
Consider a thin elastic liquid film of uniform thickness lying on the horizontal stretching sheet. The axis is chosen in the direction along which the sheet is set to motion and the axis is taken perpendicular to it. The fluid motion within the film is primarily caused solely by stretching of the sheet. The sheet is stretched by the action of two equal and opposite forces along the axis. The sheet is assumed to have velocity and the flow field is exposed to the influence of an external transverse magnetic field of strength B. We have neglected the effect of latent heat due to evaporation by assuming the liquid to be nonvolatile. Further the buoyancy is neglected due to the relatively thin liquid film, but it is not so thin that intermolecular forces come into play. The velocity and temperature fields of the liquid film obey the following boundarylayer equations: The pressure in the surrounding gas phase is assumed to be uniform and the gravity force gives rise to a hydrostatic pressure variation in the liquid film. In order to justify the boundarylayer approximation, the length scale in the primary flow direction must be significantly larger than the length scale in the cross stream direction. We choose the representative measure of the film thickness to be so that the scale ratio is large enough. This choice of length scale enables us to employ the boundarylayer approximations. Further it is assumed that the induced magnetic field is negligibly small. The associated boundary conditions are given by At this juncture we make a note that the mathematical problem is implicitly formulated only for . Further it is assumed that the surface of the planar liquid film is smooth so as to avoid the complications due to surface waves. The influence of interfacial shear due to the quiescent atmosphere, in other words the effect of surface tension, is assumed to be negligible. The viscous shear stress and the heat flux vanish at the adiabatic free surface (at ).
Let us consider a thin elastic sheet which emerges from a narrow slit at the origin of a Cartesian coordinate system. The continuous sheet at is parallel with the axis and moves in its own plane with the velocity where and are both positive constants with dimension per time. The surface temperature of the stretching sheet is assumed to vary with the distance from the slit as where is the temperature at the slit and can be taken as a constant reference temperature such that . The term can be recognized as the local Reynolds number based on the surface velocity . The expression (2.7) for the velocity of the sheet reflects that the elastic sheet which is fixed at the origin is stretched by applying a force in the positive direction and the effective stretching rate increases with time as With the same analogy the expression for the surface temperature given by (2.8) represents a situation in which the sheet temperature decreases from at the slit in proportion to and such that the amount of temperature reduction along the sheet increases with time. The applied transverse magnetic field is assumed to be of variable kind and is chosen in its special form as The particular form of the expressions for , and are chosen so as to facilitate the construction of a new similarity transformation which enables in transforming the governing partial differential equations of momentum and heat transfer into a set of nonlinear ordinary differential equations.
2.2. Similarity Transformations
We now introduce dimensionless variables and and the similarity variable as The physical stream function automatically assures mass conversion given in (2.1). The velocity components are readily obtained as The system of partial differential equations (2.1)–(2.3) with boundary conditions (2.4)–(2.6) transforms exactly into a set of ordinary differential equations and their associated boundary conditions [34]: Here is the dimensionless measure of the unsteadiness and the prime indicates differentiation with respect to Further, denotes the value of the similarity variable at the free surface so that (2.12) gives Yet is an unknown constant, which should be determined as an integral part of the boundary value problem. The rate at which film thickness varies can be obtained differentiating equation (2.20) with respect to in the form Thus the kinematics constraint at given by (2.6) transforms into the free surface condition (2.21). It is noteworthy that the momentum boundarylayer equation defined by (2.15) subject to the relevant boundary conditions (2.17)–(2.19) is decoupled from the thermal field; on the other hand the temperature field is coupled with the velocity field Since the sheet is stretched horizontally, the convection least affects the flow (i.e., buoyancy effect is negligibly small) and hence there is a oneway coupling of velocity and thermal fields.
The local skin friction coefficient, which of practical importance, is given by and the heat transfer between the surface and the fluid conventionally expressed in dimensionless form as a local Nusselt number is given by where denotes the local Reynolds number and denotes the same reference temperature (temperature difference) as in (2.8). We now march on to find the solution of the boundary value problem (2.15)–(2.19).
3. The Differential Transform Method
Transformation of the derivative of a function in one variable is as follows [35]: and the inverse transformation is defined by [36] where is the differential transform of
Taking differential transform of (2.15)(2.16) by using the related definitions in Table 1, we obtain and by using the DTM, the boundary conditions ((2.17) and (2.18)) are transformed into a recurrence equation that finally leads to the solution of a system of algebraic equations. We can consider the boundary conditions (3.4) as follows:

The differential transform of the boundary conditions is as follows:
4. Results and Discussion
System of nonlinear ordinary differential equations (3.4) with transformed boundary conditions was solved analytically using the DTM and the DTMPadé. As there are not exact solutions for nonlinear equations, a comparison of analytical and numerical solutions was conducted. In numerical solutions, the fourthorder RungeKutta method was applied. It should be noted that the solution exists only for small value of unsteadiness parameter
After finding the DTM solutions for (3.4), the Padé approximant must be applied. Ismail and Abd Rabboh [37] presented a restrictive Padé approximation for the generalized Fisher and BurgerFisher equations. The Padé approximants [38] that often show superior performance over series approximations provide a successful tool and promising scheme for identical applications. For analytical solution, the convergence analysis was performed, and in (3.3), the value is selected equal to 20. The order of Padé approximation [L,M], has sufficient accuracy; on the other hand, if the order of Padé approximation increases, the accuracy of the solution increases but sometimes increasing the order of the DTMPadé significantly increases the volume of computations. The detail of this method is well addressed in [39].
Figure 1 presents the analytical and numerical solutions of dimensionless stream function for , , and two different values of unsteadiness parameter , respectively. It observed that increases by increasing similarity variable Figure 2 indicates that increasing values of decreases the velocity profiles . Figures 3, 4, 5, and 6 display the variations of dimensionless temperature for different values of prandtl number and two different values of Comparison of the solutions obtained by the DTM and the DTMPadé revealed that DTMPadé results are in excellent agreement with that of numerical solutions [36]. In fact, using the DTMPadé aids to convergence of the DTM impressively.
5. Conclusion
In the present study, the DTM and DTMPadé were used to find analytical solutions of magnetohydrodynamics boundarylayer equations. It was found that DTMPadé is a powerful method for solving problems consisting of systems of nonlinear differential equations. In this paper, a reliable algorithm is presented based on the DTM and DTMPadé to solve some nonlinear equations. Comparison between the solutions obtained by the DTM and DTMPadé with numerical solution (fourthorder RungeKutta) remarked that the accuracy of DTMPadé is very good. The method has been applied directly without requiring linearization, discretization, or perturbation. The obtained results certify the reliability of the algorithm and give it a wider applicability to nonlinear differential equations.
Nomenclature
:  Stretching rate 
:  Sheet velocity 
:  Horizontal coordinate 
:  Vertical coordinate 
:  Horizontal velocity component 
:  Vertical velocity component 
:  Temperature 
:  Time 
:  Film thickness 
:  Unsteadiness parameter, 
:  Specific heat 
:  Dimensionless stream function, equation(2.10) 
:  Prandtl number, 
:  Eckert number, 
:  Magnetic parameter, 
:  Heat flux, 
:  Local Reynolds number, 
:  Local Nusselt number, equation (2.23). 
Greek symbols
:  Constant 
:  Dimensionless film thickness 
:  Similarity variable, equation(2.12) 
:  Dimensionless temperature, equation(2.11) 
:  Thermal diffusivity 
:  Dynamic viscosity 
:  Kinematic viscosity 
:  Density 
:  Shear stress, 
:  Stream function . 
Subscripts
:  Origin 
:  Reference value 
:  Sheet 
:  Local value. 
Superscripts
':  First derivative 
”:  Second derivative 
”':  Third derivative. 
Acknowledgment
The authors thank anonymous reviewers for their valuable comments.
References
 A. M. Lyapunov, The General Problem of the Stability of Motion, Taylor & Francis, London, UK, 1992, English translation. View at: MathSciNet
 J.H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 34, pp. 257–262, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Dehghan and F. Shakeri, “Solution of an integrodifferential equation arising in oscillating magnetic fields using He's homotopy perturbation method,” Progress in Electromagnetic Research, vol. 78, pp. 361–376, 2008. View at: Google Scholar
 M. Dehghan and F. Shakeri, “Use of He's homotpy perturbation method for solving a partial differential equation arising in modeling of flow in porous media,” Journal of Porous Media, vol. 11, pp. 765–778, 2008. View at: Google Scholar
 M. Dehghan and F. Shakeri, “The numerical solution of the second Painlevé equation,” Numerical Methods for Partial Differential Equations, vol. 25, no. 5, pp. 1238–1259, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Dehghan and J. Manafian, “The solution of the variable coefficients fourthorder parabolic partial differential equations by the homotopy perturbation method,” Zeitschrift für Naturforschung A, vol. 64a, pp. 411–419, 2009. View at: Google Scholar
 F. Soltanian, M. Dehghan, and S. M. Karbassi, “Solution of the differentialalgebraic equations via homotopy perturbation method and their engineering applications,” International Journal of Computer Mathematics. In press (2009) DOI 10.1080. View at: Google Scholar
 A. Saadatmandi, M. Dehghan, and A. Eftekhari, “Application of He's homotopy perturbation method for nonlinear system of secondorder boundary value problems,” Nonlinear Analysis: Real World Applications, vol. 10, no. 3, pp. 1912–1922, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Shakeri and M. Dehghan, “Inverse problem of diffusion equation by He's homotopy perturbation method,” Physica Scripta, vol. 75, pp. 551–556, 2007. View at: Google Scholar
 M. Dehghan and R. Salehi, “A seminumeric approach for solution of the Eikonal partial differential equation and its applications,” Numerical Methods for Partial Differential Equations, vol. 26, pp. 702–722, 2009. View at: Publisher Site  Google Scholar
 M. Dehghan and F. Shakeri, “Solution of a partial differential equation subject to temperature overspecification by He's homotopy perturbation method,” Physica Scripta, vol. 75, no. 6, pp. 778–787, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. V. Karmishin, A. I. Zhukov, and V. G. Kolosov, Methods of Dynamics Calculation and Testing for ThinWalled Structures, Mashinostroyenie, Moscow, Russia, 1990.
 M. Dehghan, M. Shakourifar, and A. Hamidi, “The solution of linear and nonlinear systems of Volterra functional equations using AdomianPade technique,” Chaos, Solitons and Fractals, vol. 39, no. 5, pp. 2509–2521, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Dehghan and F. Shakeri, “The use of the decomposition procedure of Adomian for solving a delay differential equation arising in electrodynamics,” Physica Scripta, vol. 78, Article ID 065004, 11 pages, 2008. View at: Publisher Site  Google Scholar
 J. H. He, “A new approach to nonlinear partial differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 2, no. 4, 1997. View at: Google Scholar
 M. M. Rashidi and H. Shahmohamadi, “Analytical solution of threedimensional Navier–Stokes equations for the flow near an infinite rotating disk,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 7, pp. 2999–3006, 2009. View at: Google Scholar
 T. Hayat, C. Fetecau, and M. Sajid, “Analytic solution for MHD transient rotating flow of a second grade fluid in a porous space,” Nonlinear Analysis: Real World Applications, vol. 9, no. 4, pp. 1619–1627, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Hayat, T. Javed, and M. Sajid, “Analytic solution for MHD rotating flow of a second grade fluid over a shrinking surface,” Physics Letters A, vol. 372, no. 18, pp. 3264–3273, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 M. M. Abdelkhalek, “Heat and mass transfer in MHD flow by perturbation technique,” Computational Materials Science, vol. 43, pp. 384–391, 2008. View at: Google Scholar
 M. S. Abel and M. M. Nandeppanavar, “Heat transfer in MHD viscoelastic boundary layer flow over a stretching sheet with nonuniform heat source/sink,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 2120–2131, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 A. Ishak, R. Nazar, and I. Pop, “MHD boundarylayer flow of a micropolar fluid past a wedge with constant wall heat flux,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 1, pp. 109–118, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 K. V. Prasad, D. Pal, and P. S. Datti, “MHD powerlaw fluid flow and heat transfer over a nonisothermal stretching sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, pp. 2178–2189, 2009. View at: Google Scholar
 C. Y. Wang, “Analysis of viscous flow due to a stretching sheet with surface slip and suction,” Nonlinear Analysis: Real World Applications, vol. 10, no. 1, pp. 375–380, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Y. Wang, “Exact solutions of the steadystate NavierStokes equations,” Annual Review of Fluid Mechanics, vol. 23, pp. 159–177, 1991. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 J. K. Zhou, Differential Transformation and Its Applications for Electrical Circuits, Huazhong University Press, Wuhan, China, 1986.
 C. K. Chen and S. H. Ho, “Solving partial differential equations by twodimensional differential transform method,” Applied Mathematics and Computation, vol. 106, no. 23, pp. 171–179, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M.J. Jang, C.L. Chen, and Y.C. Liu, “Twodimensional differential transform for partial differential equations,” Applied Mathematics and Computation, vol. 121, no. 23, pp. 261–270, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 I. H. A.H. Hassan, “Different applications for the differential transformation in the differential equations,” Applied Mathematics and Computation, vol. 129, no. 23, pp. 183–201, 2002. 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
 F. Ayaz, “Solutions of the system of differential equations by differential transform method,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 547–567, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Kurnaz, G. Oturanç, and M. E. Kiris, “$n$dimensional differential transformation method for solving PDEs,” International Journal of Computer Mathematics, vol. 82, no. 3, pp. 369–380, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 I. H. A.H. Hassan, “Comparison differential transformation technique with Adomian decomposition method for linear and nonlinear initial value problems,” Chaos, Solitons and Fractals, vol. 36, no. 1, pp. 53–65, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Kangalgil and F. Ayaz, “Solitary wave solutions for the KdV and mKdV equations by differential transform method,” Chaos, Solitons and Fractals, vol. 41, no. 1, pp. 464–472, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Subhas Abel and N. Mahesha, “Heat transfer in a liquid film over an unsteady stretching surface with viscous dissipation in presence of external magnetic field,” Applied Mathematical Modelling, vol. 33, pp. 3430–3441, 2009. View at: Google Scholar
 I. H. AbdelHalim Hassan, “Comparison differential transformation technique with Adomian decomposition method for linear and nonlinear initial value problems,” Chaos, Solitons and Fractals, vol. 36, no. 1, pp. 53–65, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. M. Rashidi and E. Erfani, “New analytical method for solving Burgers' and nonlinear heat transfer equations and comparison with HAM,” Computer Physics Communications, vol. 180, pp. 1539–1544, 2009. View at: Google Scholar
 H. N. A. Ismail and A. A. Abd Rabboh, “A restrictive Padé approximation for the solution of the generalized Fisher and BurgerFisher equations,” Applied Mathematics and Computation, vol. 154, no. 1, pp. 203–210, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A.M. Wazwaz, “Analytical approximations and Padé approximants for Volterra's population model,” Applied Mathematics and Computation, vol. 100, no. 1, pp. 13–25, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. M. Rashidi, “The modified differential transform method for solving MHD boundarylayer equations,” Computer Physics Communications, vol. 180, pp. 2210–2217, 2009. View at: Google Scholar
Copyright
Copyright © 2010 Mohammad Mehdi Rashidi and Mohammad Keimanesh. 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.