Research Article  Open Access
Combined AnalyticalNumerical Solution for MHD Viscous Flow over a Stretching Sheet
Abstract
We studied the twodimensional flow of viscous and electrically conducting fluid over a stretching sheet under the influence of constant magnetic field. Approximate analytical solution of governing nonlinear boundary layer equation via optimal homotopy asymptotic method (OHAM) is obtained. For numerical comparison we used RungeKuttaFehlberg fourthfifthorder method. The effect of different parameters on fluid flow is analyzed. It is found that the OHAM solution is very close to the numerical solution for different assigned values of parameters; this thus indicates the feasibility of the proposed method (OHAM).
1. Introduction
Magnetohydrodynamics (MHD) is the study of interaction of conducting fluids with electromagnetic phenomena. The flow of an electrically conducting fluid in the presence of magnetic field has many applications in engineering. Also the flow influence by a moving boundary is of crucial importance in the extrusion processes in chemical industries [1, 2].
Sakiadis [3, 4] is a pioneer in this area who has investigated the boundary layer flow with uniform speed over continuously stretching surface. Later on, Tsou et al. [5] experimentally verified the work of Sakiadis. Crane [6] studied the steadystate twodimensional boundary layer flow caused by a stretching sheet whose velocity varies linearly with the distance from a fixed point on the sheet. Chiam [7] and Dandapat and Gupta [8] considered the motion of micropolar and powerlaw fluids, respectively. Most attention has so far been devoted to the analysis of flow of viscoelastic fluids [9–13] and the joint effect of viscoelasticity and magnetic field has been worked out by Ariel [14]. Beside this, Liao [15, 16] has provided solution for impermeable and permeable stretching sheets, which shows that multiple solutions for the stretching surfaces are possible under definite conditions.
The optimal homotopy asymptotic method is a powerful approximate analytical technique that is straightforward to use and does not require the existence of any small or large parameter. Optimal homotopy asymptotic method (OHAM) is employed to construct the series solution of the problem. This method is a consistent analytical tool and it has already been applied to a number of nonlinear differential equations arising in science and engineering [17–19]. So far, as we are aware, there have been no solutions for MHD viscous flow over stretching sheet via OHAM.
This paper is organized as follows. First, in Section 2, we formulate the problem. In Section 3 we present basic principles of OHAM. The OHAM solution for MHD viscous flow is given in Section 4. In Section 5, we analyze the comparison of the solutions using OHAM with numerical method (NM). Section 6 is devoted to the concluding remarks.
2. Governing Equation
We consider the steady twodimensional flow of an incompressible viscous and electrically conducting fluid over a stretching sheet. A uniform magnetic field is applied and further assumes that the induced magnetic field is negligible due to small magnetic Reynolds number fluid. Under such assumption, the MHD boundary layer flow over a flat plate is governed by where and are the velocity components in the and directions, respectively, is the kinematic viscosity, is fluid density, is the electrical conductivity, and is the magnetic field strength, where .
The boundary conditions are given below: We introduce dimensionless variables for nondimensionalized form of momentum and energy equations [13]: Using (3), the governing equations can be reduced to nonlinear differential equation, where is a function of the similarity variable : subject to the boundary conditions as where is the nondimensional parameter and is the magnetic parameter.
3. Basic Principles of OHAM
We review the basic principles of OHAM as developed in [17–20] in the following five steps.
(i) Let us consider the following differential equation: where is problem domain, , where , are linear and nonlinear operator, is an unknown function, and is a known function.
(ii) Construct an optimal homotopy equation as where is an embedding parameter, is auxiliary function on which the convergence of the solution is greatly dependent. The auxiliary function also adjusts the convergence domain and controls the convergence region.
(iii) Expand in Taylor’s series about ; one has an approximate solution: Many researchers [19, 20] have observed that the convergence of the series (9) depends upon , ; if it is convergent, then we obtain
(iv) Substituting (10) into (8), we have the following residual: If , then will be the exact solution. For nonlinear problems, generally this will not be the case. For determining , , Galerkin’s Method, Ritz Method, or the method of least squares can be used.
(v) Finally, substitute these constants in (11) and one can get the approximate solution.
4. RungeKuttaFehlberg Method
The RungeKuttaFehlberg method (RKF45) is a technique to resolve the below problem: It is called RKF45 because the fourthorder method with five stages is used together with a fifthorder method with six stages that uses all of the points of the first one. It has a procedure to determine if the proper step size is being used. At each step, two different approximations for the solution are made and compared. If the two answers are in close agreement, the approximation is accepted. If the two answers do not agree to a specified accuracy, the step size is reduced. If the answers agree to more significant digits than required, the step size is increased. Each step requires the use of the following six values.
We need the definitions of the following: Then an approximation to the solution of the initial valued problem (I.V.P) is made using a RungeKutta method of order 4: where the four function values , , and are used.
5. Solution of the Problem via OHAM
According to the OHAM, applying (8) to (4), where primes denote differentiation with respect to .
We consider and as the following: Using (16) in (15) and some simplifying and rearranging of the terms based on the powers of , the zeroth, first, and secondorder problems are as follows.
5.1. ZerothOrder Problem
We have with boundary conditions Its solution is
5.2. FirstOrder Problem
We have with boundary conditions Its solution is
5.3. SecondOrder Problem
We have with boundary conditions Its solution becomes
5.4. ThirdOrder Problem
We have with boundary conditions Its solution is a long expression; the few terms are We obtain the fourterm solution using OHAM for : We use the method of least squares to obtain the unknown convergent constants , , and in (29); for particular case, if and , then the values of , , and are given as ; ; .
6. Results and Discussion
Tables 1 and 2 show the comparison of OHAM results with numerical RKF45 and RK4 for different values of parameters. It is noteworthy to mention here that the low error of OHAM is remarkable, while the effectiveness of the proposed method (OHAM) can be concluded from Figures 1(a), 1(b), 1(c), and 1(d), Figures 2(a), 2(b), 2(c), and 2(d), and Figures 3(a) and 3(b), where we have presented the effect of parameters and for several assigned values on . The effect of parameters on dimensionless velocity profile is observed in Figure 4, where the boundary layer thickness decreases with the increase in magnetic parameter keeping fixed the other parameter .


(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
7. Conclusion
In this study, we have successfully applied the optimal homotopy asymptotic method for MHD viscous flow over a stretching sheet. Both numerical and approximate analytical results are obtained for the problem. The results obtained via OHAM are presented in tabular and graphical forms. An excellent agreement is observed in OHAM and numerical solution. It was found that OHAM results are matched well with numerical solution. The approximate solution obtained using OHAM is also consistent with the numerical solution for various values of and . The objective of this work was to show the usefulness of OHAM.
The proposed method is simple in applicability, as it does not require discretization like numerical methods. Furthermore, this method provides a convenient way to control the convergence by optimally determining the auxiliary constants. Additionally, this method converges rapidly at lower order of approximations. Therefore, OHAM shows its latent strength and potential for the solution of nonlinear physical problems in reallife applications.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 T. Altan, S. Oh, and H. Gegel, Metal Forming Fundamentals and Applications, American Society of Metals, Metals Park, Ohio, USA, 1979.
 Z. Tadmor and I. Klein, Engineering Principles of Plasticating Extrusion, Polymer Science and Engineering Series, Van Nostrand Reinhold, New York, NY, USA, 1970.
 B. C. Sakiadis, “Boundarylayer behavior on continuous solid surface: I. Boundarylayer equations for twodimensional and axisymmetric flow,” AIChE Journal, vol. 7, pp. 26–28, 1961. View at: Google Scholar
 B. C. Sakiadis, “Boundarylayer behavior on continuous solid surface: II. Boundarylayer equations for twodimensional and axisymmetric flow,” AIChE Journal, vol. 7, pp. 221–225, 1961. View at: Google Scholar
 F. K. Tsou, E. M. Sparrow, and R. J. Goldstein, “Flow and heat transfer in the boundary layer on a continuous moving surface,” International Journal of Heat and Mass Transfer, vol. 10, no. 2, pp. 219–235, 1967. View at: Google Scholar
 L. J. Crane, “Flow past a stretching plate,” Zeitschrift für angewandte Mathematik und Physik, vol. 21, no. 4, pp. 645–647, 1970. View at: Publisher Site  Google Scholar
 T. C. Chiam, “Micropolar fluid flow over a stretching sheet,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 62, pp. 565–568, 1982. View at: Google Scholar
 B. S. Dandapat and A. S. Gupta, “Flow and heat transfer in a viscoelastic fluid over a stretching sheet,” International Journal of NonLinear Mechanics, vol. 24, no. 3, pp. 215–219, 1989. View at: Google Scholar
 P. D. Ariel, “Computation of flow of viscoelastic fluids by parameter differentiation,” International Journal for Numerical Methods in Fluids, vol. 15, no. 11, pp. 1295–1312, 1992. View at: Google Scholar
 P. D. Ariel, “On the second solution of flow of viscoelastic fluid over a stretching sheet,” Quarterly of Applied Mathematics, vol. 53, no. 4, pp. 629–632, 1995. View at: Google Scholar
 M. Sajid and T. Hayat, “The application of homotopy analysis method for MHD viscous flow due to a shrinking sheet,” Chaos, Solitons and Fractals, vol. 39, no. 3, pp. 1317–1323, 2009. View at: Publisher Site  Google Scholar
 M. Sajid, T. Javed, and T. Hayat, “MHD rotating flow of a viscous fluid over a shrinking surface,” Nonlinear Dynamics, vol. 51, no. 12, pp. 259–265, 2008. View at: Publisher Site  Google Scholar
 M. Fathizadeh, M. Madani, Y. Khan, N. Faraz, A. Yildirim, and S. Tutkun, “An effective modification of the homotopy perturbation method for MHD viscous flow over a stretching sheet,” Journal of King Saud UniversityScience, vol. 25, no. 2, pp. 107–113, 2011. View at: Google Scholar
 P. D. Ariel, “MHD flow of a viscoelastic fluid past a stretching sheet with suction,” Acta Mechanica, vol. 105, pp. 49–56, 1994. View at: Publisher Site  Google Scholar
 S. Liao, “A new branch of solutions of boundarylayer flows over an impermeable stretched plate,” International Journal of Heat and Mass Transfer, vol. 49, no. 12, pp. 2529–2539, 2005. View at: Publisher Site  Google Scholar
 S. Liao, “A new branch of solutions of boundarylayer flows over a permeable stretching plate,” International Journal of NonLinear Mechanics, vol. 42, no. 6, pp. 819–830, 2007. View at: Publisher Site  Google Scholar
 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
 V. Marinca and N. Herişanu, “Application of optimal 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
 V. Marinca and N. Herişanu, “Determination of periodic solutions for the motion of a particle on a rotating parabola by means of the optimal homotopy asymptotic method,” Journal of Sound and Vibration, vol. 329, no. 9, pp. 1450–1459, 2010. View at: Publisher Site  Google Scholar
 F. Mabood, W. A. Khan, and A. I. M. Ismail, “Optimal homotopy asymptotic method for flow and heat transfer of a viscoelastic fluid in an axisymmetric channel with a porous wall,” PLoS ONE, vol. 8, no. 12, Article ID e83581. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Fazle Mabood and Waqar A. Khan. 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.