Research Article  Open Access
Analytical Approximate Solution of Nonlinear Differential Equation Governing JefferyHamel Flow with High Magnetic Field by Adomian Decomposition Method
Abstract
The magnetohydrodynamic JefferyHamel flow is studied analytically. The traditional NavierStokes equation of fluid mechanics and Maxwell's electromagnetism governing equations reduce to nonlinear ordinary differential equations to model this problem. The analytical tool of Adomian decomposition method is used to solve this nonlinear problem. The velocity profile of the conductive fluid inside the divergent channel is studied for various values of Hartmann number. Results agree well with the numerical (RungeKutta method) results, tabulated in a table. The plots confirm that the method used is of high accuracy for different α, Ha, and Re numbers.
1. Introduction
The flow of fluid through a divergent channel is called JefferyHamel flow since introducing this problem by Jeffery [1] and Hamel [2] in 1915 and 1916, respectively. On the other hand, the term of magnetohydrodynamic (MHD) was first introduced by Bansal [3] in 1994. The theory of magnetohydrodynamics is inducing current in a moving conductive fluid in presence of magnetic field; such induced current results in force on ions of the conductive fluid. The theoretical study of magnetohydrodynamic (MHD) channel has been a subject of great interest due to its extensive applications in designing cooling systems with liquid metals, MHD generators, accelerators, pumps, and flow meters [4–7].
In fluid mechanics most of the problems are nonlinear. It is very important to develop efficient methods to solve them. Up to now, it is very difficult to obtain analytical approximations of nonlinear partial differential equations even though there are highperformance computers and computation software. The small disturbance stability of Magnetohydrodynamic stability of plane Poiseuille flow has been investigated by Makinde and Motsa [8] and Makinde [9] for generalized plane Couette flow. Their results show that magnetic field has stabilizing effects on the flow. Considerable efforts have been done to study the MHD theory for technological application of fluid pumping system in which electrical energy forces the working conductive fluid. Damping and controlling of electrically conducting fluid can be achieved by means of an electromagnetic body force (Lorentz force) produced by the interaction of an applied magnetic field and an electric current that usually is externally supplied. Harada et al. [10] studied the fundamental characteristics of linear Faraday MHD theoretically and numerically. In 2005, Anwari et al. [11] continued the Harada et al. [10] work numerically and theoretically, for various loading configurations. Jang and Lee [12] emphasized on the idea that, in such problems, the moving ions drag the bulk fluid with themselves and such MHD system induces continues pumping of conductive fluid without any moving part. Homsy et al. [13] worked and developed the same idea mentioned above. The purpose of the current work is to study the mechanics of the fluid through a divergent channel in presence of electromagnetic field (Figure 1).
The Adomian decomposition method (ADM) is used to solve a wide range of physical problems. One of the semiexact methods which does not need linearization or discretization is Adomian decomposition method, and several modifications have improved its ability [14–16].
An advantage of this method is that it can provide analytical approximation or an approximated solution to a rather wide class of nonlinear (and stochastic) equations without linearization, perturbation, closure approximation, or discretization methods. Unlike the common methods, that is, weak nonlinearity and small perturbation which change the physics of the problem due to simplification, ADM gives the approximated solution of the problem without any simplification. Thus, its results are more realistic [17]. ADM abilities have attracted many authors to use this method for solving fluid dynamic problems.
Jafari and DaftardarGejji [18] presented a modified ADM to solve a system of nonlinear equations, which yielded a series of solutions with faster accelerated convergence than the series obtained by the standard ADM. Bulut et al. [19] studied viscous incompressible flow through orifice. Allen and Syam [20] and Wang [21] investigated nonhomogeneous and classical Blasius equation by ADM. Soh [22] applied ADM to solve thin film equation. Hashim [23] presented the Adomian decomposition method for solving BVPs for fourthorder integrodifferential equations and the Blasius equation [24]. Kechil and Hashim [25] presented a nonperturbative solution of freeconvective boundarylayer equation by ADM. Chang [26] presented a decomposition solution for fins with temperaturedependent surface heat flux. Arslanturk [27] inspected on the fins efficiency of convective straight fins with temperaturedependent thermal conductivity using the decomposition method. ADM also has been used by several researchers to solve a wide range of physical problems in various engineering fields such as fluid flow and porous media simulation [28–31] and other nonlinear systems [32–38]. In recent years some researchers used new methods to solve these kinds of problems [39–41].
In this paper, we have applied ADM to find the approximate solutions of nonlinear differential equations governing the MHD JefferyHamel flow, and a comparison between the results and the numerical solution has been provided. The numerical results of this problem are done using Maple 12.
2. Governing Equations
Consider a system of cylindrical polar coordinates which steady twodimensional flow of an incompressible conducting viscous fluid from a source or sink at channel walls, lie in planes, and intersect in the axis of , assuming purely radial motion which means that there is no change in the flow parameter along the direction. The flow depends on, , and further assume that there is no magnetic field in the direction. In the reduced form of continuity, NavierStokes and Maxwell’s equations are [40] where is the electromagnetic induction, the conductivity of the fluid, the velocity along radial direction, the fluid pressure, the coefficient of kinematic viscosity, and the fluid density. Considering for purely radial flow, one can define the velocity parameter as Introducing the as the dimensionless degree, the dimensionless form of the velocity parameter can be obtained by dividing that to its maximum value as Substituting (2.5) into (2.2) and (2.3) and eliminating , one can obtain the ordinary differential equation for the normalized function profile as [10] with the following reduced form of boundary condition We introduce the Reynolds number and the Hartmann number based on the electromagnetic parameter as follows, respectively:
3. Fundamentals of Adomian Decomposition Method (ADM)
Consider equation , where represents a general nonlinear ordinary or partial differential operator including both linear and nonlinear terms. The linear terms are decomposed into , where is easily invertible (usually the highest order derivative) and is the remainder of the linear operator. Thus, the equation can be written as [14]: where indicates the nonlinear terms. By solving this equation for , since is invertible, we can write If is a secondorder operator, is twofold in definite integral. By solving (3.2), we have [14] where and are constants of integration and can be found from the boundary or initial conditions. Adomian method assumes that the solution can be expanded into infinite series as Also, the nonlinear term will be written as where are the special Adomian polynomials. By specified , the next component of can be determined: Finally, after some iteration and getting sufficient accuracy, the solution of the equation can be expressed by (3.5).
In (3.6), the Adomian polynomials can be generated by several means. Here we used the following recursive formulation [41]: Since the method does not resort to linearization or assumption of weak nonlinearity, the solution generated is in general more realistic than those achieved by simplifying the model of the physical problem.
4. Application
According to (3.1), (2.6) must be written as follows: where the differential operator is given by . Assume that the inverse of the operator exists and it can be integrated from 0 to , that is, .
Operating with on (4.1) and after exerting boundary condition on it, we have where ADM introduced the following expression: To determine the components of the is defined by applying the boundary condition of (2.7), and by assuming ,
and can be determined in a similar way from (4.5).
Using , thus According to (4.8), the accuracy of ADM solution increases by increasing the number of solution terms (). For the complete solution of (4.8), should be determined, with boundary condition of .
5. Results and Discussion
In this study the objective was to apply Adomian decomposition method to obtain an explicit analytic solution of the MHD JefferyHamel problem. The magnetic field acts as a control parameter such as the flow Reynolds number and the angle of the walls, in MHD JefferyHamel problems. There is an additional nondimensional parameter that determines the solutions, namely, the Hartmann number. Table 1 shows the value of constant for different , Ha, and Re numbers at the divergent channel.

For comparison, a few limited cases of the ADM solutions are compared with the numerical results. The comparison between the numerical results and ADM solution for velocity when and is shown in Table 2. The error bar shows an acceptable agreement between the results observed, which confirms the validity of the ADM. In these tables the error is introduced as follows: %Error of ADM for at different steps when , , and can be seen in Figure 2; it shows that at first we have large error then it decreases and after 12 steps the error becomes minimized.

Figures 3, 4, and 5 show the magnetic field effect on the velocity profiles for divergent channels. There are good agreements between the numerical solution obtained by the fourthorder RungeKutta method and the differential transformation method.
Under magnetic field the Lorentz force effect is in opposite of the momentum’s direction that stabilizes the velocity profile.
The results show moderate increases in the velocity with increasing Hartmann numbers at small angle () and differences between velocity profiles are more noticeable at greater angles. Backflow is excluded in converging channels [42] but it may occur for large Reynolds numbers in diverging channels. For specified opening angle, after a critical Reynolds number, we observe that separation and backflow are started.
Figures 6–8 show the magnetic field effects at constant and different Reynolds numbers. At , with increasing Hartmann number the velocity profile becomes flat and thickness of boundary layer decreases, but at this Reynolds number no backflow is observed as shown in Figure 6. It can be seen in Figure 9 that without magnetic field at , the backflow starts and that with increasing Hartmann number this phenomenon eliminates.
By increasing Reynolds number the backflow expands and so greater magnetic field is needed in order to eliminate it. As shown in Figures 7 and 8 at , the back flow eliminates at while at , this occurs at .
6. Conclusion
In this paper, magnetohydrodynamic JefferyHamel flow has been solved via a sort of analytical method, Adomian decomposition method (ADM). Also this problem is solved by a numerical method (the RungeKutta method of order 4), and some conclusions are summarized as follows.(a)Adomian decomposition method is a powerful approach for solving MHD JefferyHamel flow in high magnetic field, and it can be observed that there is a good agreement between the present and numerical results. (b)Increasing Reynolds numbers leads to adverse pressure gradient which causes velocity reduction near the walls.(c)Increasing Hartmann number will lead to backflow reduction. In greater angles or Reynolds numbers high Hartmann number is needed for the reduction of backflow.
Nomenclature
:  Magnetic field () 
:  Constant 
:  Linear term 
:  Nonlinear term 
:  The remainder of linear operator 
:  General nonlinear operator 
:  Adomian polynominal 
:  Dimensionless velocity 
:  Hartmann number 
:  Density 
:  Pressure term 
:  Reynolds number 
:  Cylindrical coordinates 
:  Maximum value of velocity 
:  Velocity components along , axes, respectively. 
Greek Symbols
:  Kinematic viscosity 
:  Viscous stresses 
:  Angle of the channel 
:  Any angle 
:  Dimensionless angle. 
References
 G. B. Jeffery, “The twodimensional steady motion of a viscous fluid,” Philosophical Magazine, vol. 6, pp. 455–465, 1915. View at: Google Scholar
 G. Hamel, “Spiralförmige Bewgungen Zäher Flüssigkeiten,” Jahresbericht der Deutschen MathematikerVereinigung, vol. 25, pp. 34–60, 1916. View at: Google Scholar
 L. Bansal, Magnetofluiddynamics of Viscous Fluids, Jaipur Publishing House, Jaipur, India, 1994.
 J. E. Cha, Y. C. Ahn, and M. H. Kim, “Flow measurement with an electromagnetic flowmeter in twophase bubbly and slug flow regimes,” Flow Measurement and Instrumentation, vol. 12, no. 56, pp. 329–339, 2002. View at: Publisher Site  Google Scholar
 M. Tendler, “Confinement and related transport in extrap geometry,” Nuclear Instruments and Methods In Physics Research, vol. 207, no. 12, pp. 233–240, 1983. View at: Google Scholar
 J. Mossino, “Some nonlinear problems involving a free boundary in plasma physics,” Journal of Differential Equations, vol. 34, no. 1, pp. 114–138, 1979. View at: Google Scholar
 R. Nijsing and W. Eifler, “A computational analysis of transient heat transfer in fuel rod bundles with single phase liquid metal cooling,” Nuclear Engineering and Design, vol. 62, no. 1–3, pp. 39–68, 1980. View at: Google Scholar
 O. D. Makinde and S. S. Motsa, “Hydromagnetic stability of plane poiseuille flow using chebyshev spectral collocation method,” Journal of Institute of Mathematics and Computer Sciences, vol. 12, no. 2, pp. 175–183, 2001. View at: Google Scholar
 O. D. Makinde, “Magnetohydrodynamic stability of planePoiseuille flow using multideck asymptotic technique,” Mathematical and Computer Modelling, vol. 37, no. 34, pp. 251–259, 2003. View at: Google Scholar
 N. Harada, J. Ikewada, and Y. Terasaki, “Basic studies on an MHD accelerator,” The American Institute of Aeronautics and Astronautics, vol. 2175, 2002. View at: Google Scholar
 M. Anwari, N. Harada, and S. Takahashi, “Performance study of a magnetohydrodynamic accelerator using airplasma as working gas,” Energy Conversion and Management, vol. 46, no. 1516, pp. 2605–2613, 2005. View at: Publisher Site  Google Scholar
 J. Jang and S. S. Lee, “Theoretical and experimental study of MHD (magnetohydrodynamic) micropump,” Sensors and Actuators, vol. 80, no. 1, pp. 84–89, 2000. View at: Publisher Site  Google Scholar
 A. Homsy, S. Koster, J. C. T. Eijkel et al., “A high current density DC magnetohydrodynamic (MHD) micropump,” The Royal Society of Chemistry, Lab Chip, vol. 5, no. 4, pp. 466–471, 2005. View at: Publisher Site  Google Scholar
 G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988. View at: Google Scholar
 H. Jafari and V. DaftardarGejji, “Revised Adomian decomposition method for solving systems of ordinary and fractional differential equations,” Applied Mathematics and Computation, vol. 181, no. 1, pp. 598–608, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 S. Ghosh, A. Roy, and D. Roy, “An adaptation of Adomian decomposition for numericanalytic integration of strongly nonlinear and chaotic oscillators,” Computer Methods in Applied Mechanics and Engineering, vol. 196, no. 4–6, pp. 1133–1153, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988. View at: Google Scholar
 H. Jafari and V. DaftardarGejji, “Revised Adomian decomposition method for solving a system of nonlinear equations,” Applied Mathematics and Computation, vol. 175, no. 1, pp. 1–7, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 H. Bulut, M. Ergut, V. Asil, and R. H. Bokor, “Numerical solution of a viscous incompressible flow problem through an orifice by Adomian decomposition method,” Applied Mathematics and Computation, vol. 153, no. 3, pp. 733–741, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 F. M. Allan and M. I. Syam, “On the analytic solutions of the nonhomogeneous Blasius problem,” Journal of Computational and Applied Mathematics, vol. 182, no. 2, pp. 362–371, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Wang, “A new algorithm for solving classical Blasius equation,” Applied Mathematics and Computation, vol. 157, no. 1, pp. 1–9, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 C. W. Soh, “Nonperturbative semianalytical sourcetype solutions of thinfilm equation,” Applied Mathematics and Computation, vol. 174, no. 2, pp. 1576–1585, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 I. Hashim, “Adomian decomposition method for solving BVPs for fourthorder integrodifferential equations,” Journal of Computational and Applied Mathematics, vol. 193, no. 2, pp. 658–664, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 I. Hashim, “Comments on "a new algorithm for solving classical Blasius equation,” The Journal of Computational and Applied Mathematics, vol. 182, pp. 362–371, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 S. A. Kechil and I. Hashim, “Nonperturbative solution of freeconvective boundarylayer equation by Adomian decomposition method,” Physics Letters A, vol. 363, no. 12, pp. 110–114, 2007. View at: Publisher Site  Google Scholar
 M. H. Chang, “A decomposition solution for fins with temperature dependent surface heat flux,” International Journal of Heat and Mass Transfer, vol. 48, no. 9, pp. 1819–1824, 2005. View at: Publisher Site  Google Scholar
 C. Arslanturk, “A decomposition method for fin efficiency of convective straight fins with temperaturedependent thermal conductivity,” International Communications in Heat and Mass Transfer, vol. 32, no. 6, pp. 831–841, 2005. View at: Publisher Site  Google Scholar
 F. M. Allan and M. I. Syam, “On the analytic solutions of the nonhomogeneous Blasius problem,” Journal of Computational and Applied Mathematics, vol. 182, no. 2, pp. 362–371, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. Kaya and A. Yokus, “A decomposition method for finding solitary and periodic solutions for a coupled higherdimensional Burgers equations,” Applied Mathematics and Computation, vol. 164, no. 3, pp. 857–864, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Pamuk, “Solution of the porous media equation by Adomian's decomposition method,” Physics Letters A, vol. 344, no. 2–4, pp. 184–188, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 F. M Allan and K. AlKhaled, “An approximation of the analytic solution of the shock wave equation,” Journal of Computational and Applied Mathematics, vol. 192, no. 2, pp. 301–309, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 V. DaftardarGejji and H. Jafari, “An iterative method for solving nonlinear functional equations,” Journal of Mathematical Analysis and Applications, vol. 316, no. 2, pp. 753–763, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. Lesnic, “Decomposition methods for nonlinear, noncharacteristic Cauchy heat problems,” Communications in Nonlinear Science and Numerical Simulation, vol. 10, no. 6, pp. 581–596, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Zhu, Q. Chang, and S. Wu, “A new algorithm for calculating Adomian polynomials,” Applied Mathematics and Computation, vol. 169, no. 1, pp. 402–416, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 X. G. Luo, “A twostep Adomian decomposition method,” Applied Mathematics and Computation, vol. 170, no. 1, pp. 570–583, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 X. Zhang, “A modification of the Adomian decomposition method for a class of nonlinear singular boundary value problems,” Journal of Computational and Applied Mathematics, vol. 180, no. 2, pp. 377–389, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 D. Kaya and A. Yokus, “A comparison of partial solutions in the decomposition method for linear and nonlinear partial differential equations,” Mathematics and Computers in Simulation, vol. 60, no. 6, pp. 507–512, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. I. Axford, “The mahnetohydrodynamic JefferyHamel problem for a weakly conducting fluid,” The Quarterly Journal of Mechanics and Applied Mathematics, vol. 14, pp. 335–351, 1961. View at: Google Scholar
 M. Esmaeilpour and D. D. Ganji, “Solution of the JefferyHamel flow problem by optimal homotopy asymptotic method,” Computers and Mathematics with Applications, vol. 59, no. 11, pp. 3405–3411, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. Z. Ganji, D. D. Ganji, and M. Esmaeilpour, “Study on nonlinear JefferyHamel flow by He's semianalytical methods and comparison with numerical results,” Computers and Mathematics with Applications, vol. 58, no. 1112, pp. 2107–2116, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 O. D. Makinde and P. Y. Mhone, “Temporal stability of small disturbances in MHD JefferyHamel flows,” Computers and Mathematics with Applications, vol. 53, no. 1, pp. 128–136, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. Z. Ganji, D. D. Ganji, and M. Esmaeilpour, “Study on nonlinear JefferyHamel flow by He's semianalytical methods and comparison with numerical results,” Computers and Mathematics with Applications, vol. 58, no. 1112, pp. 2107–2116, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2011 D. D. Ganji 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.