Research Article  Open Access
Reza Hosseini, Sadegh Poozesh, Saeed Dinarvand, "MHD Flow of an Incompressible Viscous Fluid through Convergent or Divergent Channels in Presence of a High Magnetic Field", Journal of Applied Mathematics, vol. 2012, Article ID 157067, 12 pages, 2012. https://doi.org/10.1155/2012/157067
MHD Flow of an Incompressible Viscous Fluid through Convergent or Divergent Channels in Presence of a High Magnetic Field
Abstract
The flow of an incompressible electrically conducting viscous fluid in convergent or divergent channels under the influence of an externally applied homogeneous magnetic field is studied both analytically and numerically. NavierStokes equations of fluid mechanics and Maxwell’s electromagnetism equations are reduced into highly nonlinear ordinary differential equation. The resulting nonlinear equation has been solved analytically using a very efficient technique, namely, differential transform method (DTM). The DTM solution is compared with the results obtained by a numerical method (shooting method, coupled with fourthorder RungeKutta scheme). The plots have revealed the physical characteristics of flow by changing angles of the channel, Hartmann and Reynolds numbers.
1. Introduction
The incompressible viscous fluid flow through convergent or divergent channels is one of the most applicable cases in many applications such as aerospace, chemical, civil, environmental, mechanical, and biomechanical engineering as well as in understanding rivers and canals. Jeffery [1] and Hamel [2] have carried out the mathematical formulations of this problem in 1915 and 1916, respectively. If we simplify NavierStokes equations in the particular case of twodimensional flow through a channel with inclined walls, finally we can reach JefferyHamel problem [3–6]. JefferyHamel flows have been extensively studied by several authors and discussed in many textbooks, for example, [7–11], and so forth.
The study of electrically conducting viscous fluid that flows through convergent or divergent channels under the influence of an external magnetic field not only is fascinating theoretically but also finds applications in mathematical modeling of several industrial and biological systems. A possible practical application of the theory we envisage is in the field of industrial metal casting, the control of molten metal flows. Another area in which the theoretical study may be of interest is in the motion of liquid metals or alloys in the cooling systems of advanced nuclear reactors [12]. Clearly, the motion in the region with intersecting walls may represent a local transition between two parallel channels with different crosssections, a widening or a contraction of the flow. The first recorded use of the word magnetohydrodynamics (MHD) is by Bansal [13]. The theory of MHD is inducing current in a moving conductive fluid in the presence of magnetic field which creates force on electrons of the conductive fluid and also changes the magnetic field itself. A survey of magnetohydrodynamics studies in the mentioned technological field can be found in [14]. The problem is basically an extension of classical JefferyHamel flows of ordinary fluid mechanics to MHD. In the MHD solution an external magnetic field acts as a control parameter for both convergent and divergent channel flows. Here, beside the flow Reynolds number and the channel angular widths, at least an additional dimensionless parameter appears, namely, the Hartman number. Hence, a much larger variety of solutions than in the classical problem are expected.
Zhou was the first one who introduced DTM as an efficient method to apply for electrical circuits in his paper entitled “Differential transformation and its application for electrical circuits,” [15]. It was used to solve both linear and nonlinear initial value problems in electric circuit analysis. The differential transform method (DTM) is an analytical method for solving ordinary differential equations, partial differential and integral equations. The method provides us with easily computable components and the solution is obtained in terms of convergent series. The main advantages of this method, compared to other analytic methods are controllable accuracy, and high efficiency, which is exhibited by the rapid convergence of the solution. The DTM gives exact values of the th derivative of an analytic function at a point in terms of known and unknown boundary conditions. This method constructs, for differential equations, an analytical solution in the form of a polynomial. It is different from the traditional highorder Taylor series method, which requires symbolic computations of the necessary derivatives of the data functions. The disadvantage of Taylor series method is that this method computationally takes long time for large orders. The DTM is an iterative procedure for obtaining analytic Taylor series solutions of differential equations in a fast manner. This methodology introduces a promising approach for many applications in various domains of nonlinear problems. Various applications of DTM can be found in [15–22].
The aim of the present work is to investigate the MHD flow through convergent or divergent channels in presence of a high magnetic field, by the differential transform method (DTM). The governing highly nonlinear equation of this problem is also solved numerically by shooting method, coupled with fourthorder RungeKutta scheme.
2. Problem Statement and Mathematical Formulation
Consider a system of cylindrical polar coordinates, where the steady twodimensional flow of an incompressible conducting viscous fluid from a source or sink at channel walls lie in planes and intersect in axis. The schematic diagram of problem is illustrated in Figure 1. We assumed that; it means that there are no changes with respect to direction; thus the motion is purely in radial direction and merely depends on and [3, 23], and there is no magnetic field along axis. The reduced form of continuity, NavierStokes and Maxwell’s equations in polar system are where is the electromagnetic induction strength, the conductivity of the fluid, the velocity along radial direction, the fluid pressure, the coefficient of kinematic viscosity, and the fluid density. From (2.1), Using dimensionless parameters, where is the semiangle between the two inclined walls. Substituting (2.5) into (2.2) and (2.3) and eliminating , one can obtain an ordinary differential equation for the normalized function profile [3]: where the Reynolds number is Then we have boundary condition as follows: The Hartmann number is
3. Analytical Approximations by Means of DTM
The differential transformation method is an analytical method for a vast variety of differential equations including ODEs and PDEs [15]. This method uses polynomials form to approximate the exact solutions. We now take a brief review to the DTM. The differential transform of the th derivative of function is defined as follows [24]: where is the base function and is the transformed function. The differential inverse transform of is defined as Equations (3.1) and (3.2) give the following: This shows that differential transform is derived from Taylor series expansion, but the method does not evaluate the derivatives symbolically. However, relative derivatives are calculated by an iterative way which is described by the transformed equations of the base function. We approximate by a finite series and (3.2) can be written as, The main steps of the DTM are the following. First, we apply the differential transform (3.1) to the given differential equation or a system of differential equations to obtain a recursive relation. Second, solving the recursive relation and then using the differential inverse transform (3.2) we obtain the solution of the problem.
Using (3.1) and (3.2) the following theorems can be deduced.
Theorem 3.1. If , then .
Theorem 3.2. If , then, where is a constant.
Theorem 3.3. If , then .
Theorem 3.4. If , then
Theorem 3.5. If , then,
Taking differential transform from (2.6), one can obtain: The boundary conditions, (2.8), are transformed into:
4. Numerical Simulation
Many categories of numerical methods have been developed for the numerical solution of differential equations; among them are the wellknown fitted RungeKutta (RK) and RungeKuttaNyström (RKN) methods [24]. In the last decade exponentially fitted RK and RKN methods have been constructed by many authors. Simos [25, 26] and Vanden Berghe et al. [27, 28] first constructed exponentially fitted RK methods. Also exponentially fitted RKN methods have been studied by Simos [29], Van de Vyver [30], and Franco [31].
In recent years, modern numerical techniques have been successfully employed to solve various types of mathematics and engineering problems. For example, Anastassi and Simos have applied Phasefitted RungeKuttaNystrom method for the numerical solution of initial value problems with oscillating solutions [32]. Simos has analyzed the Schrodinger equation by a new Numerovtype method [33]. The solution of the Schrodinger equation and related problems by means of the optimized explicit RungeKutta method with increased phaselag order has been presented by Kosti et al. [34]. Stavroyiannis investigated the linear periodic initial value problems by nonlinear explicit twostep fourthalgebraicorder method [35]. The symmetric eightstep predictorcorrector method was applied for the radial Schrodinger equation and related IVPs with oscillating solutions by Panopoulos et al. [36]. Briefly speaking, this shows the great potential of modern numerical techniques for solving different problems in science and engineering.
The shooting method works by considering the boundary conditions as a multivariable function of initial conditions at some point, reducing the boundary value problem to finding the initial conditions that give a root. The advantage of the shooting method is that it takes advantage of the speed and adaptivity of methods for initial value problems. The basic concept of the shooting method can be obtained from [37]. In this paper the shooting method, coupled with fourthorder RungeKutta scheme, is used for solving the problem of MHD flow through convergent or divergent channel in presence of a high magnetic field. An important point to note is that the numerical results have been obtained using Mathematica software. It is worth to mentioning that authors intend to apply some modern numerical techniques for the present problem and to compare the efficiency of different methods in the future.
5. Results and Discussions
The objective of the present study was to apply the differential transform method to obtain an explicit analytic solution of MHD flow through convergent or divergent channels in presence of a high magnetic field. The magnetic field plays its role in nondimensional parameter, namely, the Hartmann number. Figure 2 that shows velocity profile variation for a convergent channel with fixed Re number reveals the fact that by increasing magnetic field the velocity profile becomes flat and thickness of boundary layer decreases. In fact magnetic field induces a force in opposite of the momentum’s direction that stabilizes the velocity profile. There is a reverse behavior in velocity profile as shown in Figure 3 for divergent channel, where we see that, by decreasing Hartman number, the velocity profile becomes flat and thickness of boundary layer decreases. As demonstrated in Figure 4 for pure inflow regime (negative Re number) in divergent channel, when Hartmann number increases we have a little change in velocity profile. Comparison of Figures 3 and 5 shows moderate increases in the velocity with increasing Hartmann numbers at small angle and differences between velocity profiles are more noticeable at greater angle. Figures 6 and 7 show the velocity profile in divergent and convergent channels, respectively. In Figure 6, we can see that as increases the effect of walls on fluid flow decreases when we move away from them which lead to an increase of velocity, while there is a reverse behavior in velocity profile as shown in Figure 7, where the values of decrease. Figures 8 and 9 illustrate the effect of Reynolds numbers on the fluid velocity for fixed Hartmann numbers. According to Figure 8, for inflow regime as Re number increases, in a particular distance from the wall, velocity decreases and effects of backflow phenomenon are revealed. There is a reverse behavior in Figure 9 where as Re number increases, value of local velocity becomes more close to value of centerline velocity. We can infer from Figures 8 and 9 for inflow regime, back flow is prevented in the case of convergent channels but is possible for large Reynolds numbers in the case of divergent channels. Figures 8 and 9 show that there is a reverse condition for outflow regime.
For comparison, a few limited cases of the DTM solutions are compared with the numerical results. The comparison between the numerical results and DTM solution for velocity when and is shown in Table 1. The error bar shows an acceptable agreement between the results observed, which confirms the validity of the DTM. In this table error is introduced as .

6. Final Remarks
In this investigation, the flow of an incompressible electricallyconducting viscous fluid in convergent or divergent channels under the influence of an externally applied homogeneous magnetic field is studied both analytically and numerically. The governing highly nonlinear equation has been solved analytically using the differential transform method (DTM) and numerically by shooting method, coupled with fourthorder RungeKutta scheme. The behavior of the DTM solution is in good agreement with the numerical simulation. Graphical results are presented to investigate the influence of the angles of the channel, Hartmann number and Reynolds number on the velocity profiles. From the presented analysis, the following observations are noted.(i)For a convergent channel with fixed number, with increasing the magnetic field, the velocity profile becomes flat and thickness of boundary layer decreases. While there is a reverse behavior in velocity profile for divergent channel.(ii)For a divergent channel, increasing in the Reynolds numbers leads to adverse pressure gradient causing velocity reduction near the walls.(iii)For pure outflow regime (positive Re number), there are moderate increases in the velocity with increasing Hartmann number.(iv)For divergent channels, as increases the effect of walls on fluid flow decreases when we move away from them which lead to an increase of velocity, while there is a reverse behavior in velocity profile for convergent channels.(v)For inflow regime (negative Re number), as Re number increases, in a particular distance from the wall, velocity decreases and effects of backflow phenomenon are revealed.
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,” Jahresber Deutsch MathVerein, vol. 25, pp. 34–60, 1916. View at: Google Scholar
 H. Schlichting, Boundary Layer Theory, McGrawHill, New York, NY, USA, 1955.
 L. Rosenhead, “The steady twodimensional radial flow of viscous fluid between two inclined plane walls,” Proceedings of the Royal Society A, vol. 175, pp. 436–467, 1940. View at: Google Scholar
 K. Batchelor, An Introduction to Fluid Dynamics, Cambridge University Press, 1967.
 I. J. Sobey and P. G. Drazin, “Bifurcations of twodimensional channel flows,” Journal of Fluid Mechanics, vol. 171, pp. 263–287, 1986. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. E. Fraenkel, “Laminar flow in symmetrical channels with slightly curved walls. I. On the JefferyHamel solutions for flow between plane walls,” Proceedings of the Royal Society A, vol. 267, pp. 119–138, 1962. View at: Google Scholar  Zentralblatt MATH
 O. D. Makinde and P. Y. Mhone, “HermitePadé approximation approach to MHD JefferyHamel flows,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 966–972, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Hamadiche, J. Scott, and D. Jeandel, “Temporal stability of JefferyHamel flow,” Journal of Fluid Mechanics, vol. 268, pp. 71–88, 1994. View at: Google Scholar
 A. McAlpine and P. G. Drazin, “On the spatiotemporal development of small perturbations of JefferyHamel flows,” Fluid Dynamics Research, vol. 22, no. 3, pp. 123–138, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. M. Moghimi, G. Domairry, S. Soleimani, E. Ghasemi, and H. Bararnia, “Application of homotopy analysis method to solve MHD JefferyHamel flows in nonparallel walls,” Advances in Engineering Software, vol. 42, no. 3, pp. 108–113, 2011. View at: Publisher Site  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
 L. Bansal, Magnetofluiddynamics of Viscous Fluids, Jaipur Publishing House, Jaipur, India, 1994.
 R. Moreau, Magnetohydrodynamics, vol. 3 of Fluid Mechanics and its Applications, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1990.
 J. K. Zhou, Differential Transformation and Its Applications for Electrical Circuits, Huazhong University Press, Wuhan, China, 1986.
 J. Biazar and M. Eslami, “Analytic solution for Telegraph equation by differential transform method,” Physics Letters. A, vol. 374, no. 29, pp. 2904–2906, 2010. View at: Publisher Site  Google Scholar
 M. S. Abel, N. Mahesha, and J. Tawade, “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, no. 8, pp. 3430–3441, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 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. Arikoglu and I. Ozkol, “Solution of boundary value problems for integrodifferential equations by using differential transform method,” Applied Mathematics and Computation, vol. 168, no. 2, pp. 1145–1158, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  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
 A. S. V. Ravi Kanth and K. Aruna, “Differential transform method for solving the linear and nonlinear KleinGordon equation,” Computer Physics Communications, vol. 180, no. 5, pp. 708–711, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. M. Odibat, C. Bertelle, M. A. AzizAlaoui, and G. H. E. Duchamp, “A multistep differential transform method and application to nonchaotic or chaotic systems,” Computers & Mathematics with Applications, vol. 59, no. 4, pp. 1462–1472, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. K. Rathy, An Introduction to Fluid Dynamics, Oxford & IBH Publishing Co., New Delhi, India, 1976.
 J. R. Cash and S. Girdlestone, “Variable step RungeKuttaNyström methods for the numerical solution of reversible systems,” Journal of Numerical Analysis, Industrial and Applied Mathematics, vol. 1, no. 1, pp. 59–80, 2006. View at: Google Scholar  Zentralblatt MATH
 T. E. Simos, “An exponentiallyfitted RungeKutta method for the numerical integration of initialvalue problems with periodic or oscillating solutions,” Computer Physics Communications., vol. 115, no. 1, pp. 1–8, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. E. Simos, “Exponentially fitted RungeKutta methods for the numerical solution of the Schrödinger equation and related problems,” Computational Materials Science, vol. 18, no. 34, pp. 315–332, 2000. View at: Google Scholar
 L. Gr. Ixaru and G. Vanden Berghe, Exponential Fitting, vol. 568 of Mathematics and Is Applications, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2004.
 G. Vanden Berghe, H. De Meyer, M. Van Daele, and T. Van Hecke, “Exponentiallyfitted explicit RungeKutta methods,” Computer Physics Communications, vol. 123, no. 1–3, pp. 7–15, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. E. Simos, “Exponentiallyfitted RungeKuttaNyström method for the numerical solution of initialvalue problems with oscillating solutions,” Applied Mathematics Letters, vol. 15, no. 2, pp. 217–225, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Van De Vyver, “A symplectic exponentially fitted modified RungeKuttaNyström method for the numerical integration of orbital problems,” New Astronomy, vol. 10, no. 4, pp. 261–269, 2005. View at: Publisher Site  Google Scholar
 J. M. Franco, “Exponentially fitted explicit RungeKuttaNyström methods,” Journal of Computational and Applied Mathematics, vol. 167, no. 1, pp. 1–19, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. F. Papadopoulos, Z. A. Anastassi, and T. E. Simos, “A phasefitted RungeKuttaNyström method for the numerical solution of initial value problems with oscillating solutions,” Computer Physics Communications, vol. 180, no. 10, pp. 1839–1846, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 T. E. Simos, “A new Numerovtype method for the numerical solution of the Schrödinger equation,” Journal of Mathematical Chemistry, vol. 46, no. 3, pp. 981–1007, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. A. Kosti, Z. A. Anastassi, and T. E. Simos, “An optimized explicit RungeKutta method with increased phaselag order for the numerical solution of the Schrödinger equation and related problems,” Journal of Mathematical Chemistry, vol. 47, no. 1, pp. 315–330, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Stavroyiannis and T. E. Simos, “A nonlinear explicit twostep fourth algebraic order method of order infinity for linear periodic initial value problems,” Computer Physics Communications, vol. 181, no. 8, pp. 1362–1368, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 G. A. Panopoulos, Z. A. Anastassi, and T. E. Simos, “A symmetric eightstep predictorcorrector method for the numerical solution of the radial Schrödinger equation and related IVPs with oscillating solutions,” Computer Physics Communications, vol. 182, no. 8, pp. 1626–1637, 2011. View at: Publisher Site  Google Scholar
 S. M. Roberts and J. S. Shipman, TwoPoint Boundary Value Problems: Shooting Methods, Modern Analytic and Computational Methods in Science and Mathematics, No. 3, American Elsevier, New York, NY, USA, 1972.
Copyright
Copyright © 2012 Reza Hosseini 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.