`ISRN Mechanical EngineeringVolumeΒ 2011, Article IDΒ 505673, 10 pageshttp://dx.doi.org/10.5402/2011/505673`
Research Article

## Effect of Temperature-Dependent Variable Viscosity on Magnetohydrodynamic Natural Convection Flow along a Vertical Wavy Surface

Department of Mathematics, Bangladesh University of Engineering and Technology, Dhaka 1000, Bangladesh

Received 17 January 2011; Accepted 21 February 2011

Academic Editors: C.-Y.Β Huang, G.-J.Β Wang, and M.Β Al-Nimr

Copyright Β© 2011 Nazma Parveen and Md. Abdul Alim. 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.

#### Abstract

The effect of temperature dependent variable viscosity on magnetohydrodynamic (MHD) natural convection flow of viscous incompressible fluid along a uniformly heated vertical wavy surface has been investigated. The governing boundary layer equations are first transformed into a nondimensional form using suitable set of dimensionless variables. The resulting nonlinear system of partial differential equations are mapped into the domain of a vertical flat plate and then solved numerically employing the implicit finite difference method, known as Keller-box scheme. The numerical results of the surface shear stress in terms of skin friction coefficient and the rate of heat transfer in terms of local Nusselt number, the stream lines and the isotherms are shown graphically for a selection of parameters set consisting of viscosity parameter (), magnetic parameter (), and Prandtl number (Pr). Numerical results of the local skin friction coefficient and the rate of heat transfer for different values are also presented in tabular form.

#### 1. Introduction

Laminar natural convection boundary layer flow and heat transfer problem from a vertical wavy surface gets a great deal of attention in various branches of engineering. If the surface is roughened, the flow is disturbed by the surface and this alters the rate of heat transfer. These types of roughened surface are taken into account in several heat transfer collectors, flat plate condensers in refrigerators and heat exchanger. One common example of a heat exchanger is the radiator used in cars, in which the heat generated from engine transferred to air flowing through the radiator. Machine-roughened surface enhanced heat transfer and the interface between concurrent or countercurrent two-phase flow is another example remotely related to this problem. Such an interface is always wavy and momentum transfer across it is by no means similar to that across a smooth, flat surface and neither is the heat transfer.

The study of the flow of electrically conducting fluid in the presence of magnetic field is also important from the technical point of view and such types of problems have received much attention by many researchers. The specific problem selected for study is the flow and heat transfer in an electrically conducting fluid adjacent to the surface. The surface is maintained at a uniform temperature , which may either exceed the ambient temperature or may be less than . When , an upward flow is established along the surface due to free convection, where as for , there is a down flow. The interaction of the magnetic field and the moving electric charge carried by the flowing fluid induces a force, which is proportional to the magnitude of the longitudinal velocity and acts in the opposite direction is also very small. Additionally, a magnetic field of strength acts normal to the surface. Consequently, the influence of the magnetic field on the boundary layer is exerted only through induced forces within the boundary layer itself, with no additional effects arising from the free stream pressure gradient. The problem of magnetohydrodynamic free convection in a strong cross field was investigated by Kuiken [1]. Also the effect of magnetic field on free convection heat transfer has been studied by Sparrow and Cess [2].

The viscosity of the fluid is proportional to a linear function of temperature which was proposed by Charraudeau [3]. The effects of surface waviness on the natural convection boundary layer flow of a Newtonian fluid have studied by Yao [4] and Moulic and Yao [5, 6]. In these papers they used an extended Prantdlβs transposition theorem and a finite-difference scheme. Yao proposed a simple transformation to study the natural convection heat transfer for an isothermal vertical sinusoidal surface. This simple coordinate transformation method to change the wavy surface into a flat plate. Transient-free convection flow with temperature dependent viscosity in a fluid saturated porous media has shown by Mehta and Sood [7]. In this paper they found that the flow characteristics substantially change when the effect of temperature dependent viscosity is considered. The effects of temperature dependent viscosity on the free and mixed convection flow from a vertical flat plate in the region near the leading edge have studied by Kafoussius and Williams [8] and Kafoussius and Rees [9]. Rees and Pop [10β12] investigated the natural convection boundary layer induced by vertical and horizontal wavy surface exhibiting small amplitude waves embedded in a porous medium. Hady et al. [13] investigated the mixed convection boundary layer flow on a continuous flat plate with variable viscosity. Alam et al. [14] have also studied the problem of free convection from a wavy vertical surface in presence of a transverse magnetic field. On the other hand, the combined effects of thermal and mass diffusion on the natural convection flow of a viscous incompressible fluid along a vertical wavy surface have been investigated by Hossain and Rees [15]. In this paper the effect of waviness of the surface on the heat and mass flux is investigated in combination with the species concentration for a fluid having Prandtl number equal to 0.7. The natural convection heat and mass transfer near a vertical wavy surface with constant wall temperature and concentration in a porous medium studied by Cheng [16]. Hossain et al. [17] investigated the natural convection flow past a permeable wedge for the fluid having temperature dependent viscosity and thermal conductivity. Free convection flow with variable viscosity and thermal diffusivity along a vertical plate in the presence of magnetic field investigated by Elbashbeshy [18]. Munir et al. [19, 20] investigated natural convection with variable viscosity and thermal conductivity along a vertical wavy cone. The problems of natural convection of fluid with temperature dependent viscosity along a heated vertical wavy surface have studied by Kabir et al. [21]. Molla et al. [22] have studied natural convection flow along a vertical wavy surface with uniform surface temperature in presence of heat generation/absorption. Nasrin and Alim [23] investigated magnetohydrodynamic free convection flow along a vertical flat plate with thermal conductivity and viscosity depending on temperature. Numerical study on a vertical plate with variable viscosity and thermal conductivity has investigated by Palani and Kim [24]. Aldoss et al. [25] investigated MHD mixed convection from a horizontal circular cylinder. Al-Nimr and Alkam [26] have studied magneto-hydrodynamics transient free convection in open-ended vertical annuli. Al-Nimr and Hader [27] also investigated MHD free convection flow in open-ended vertical porous channels. Damseh et al. [28] investigated entropy generation during fluid flow in a channel under the effect of transverse magnetic field. From the above investigations it is found that variation of viscosity with temperature of magnetic field is an interesting macroscopic physical phenomenon in fluid mechanics.

The present study is to incorporate the idea that the effects of temperature dependent viscosity in presence of magnetic field of electrically conducting fluid with free convection boundary layer flow along a vertical wavy surface. However, it is known that this physical property (viscosity) may be change significantly with temperature. To predict accurately the flow behavior, it is necessary to take into account of viscosity. The governing partial differential equations are reduced to locally non-similar partial differential forms by adopting some appropriate transformations. The transformed boundary layer equations are solved numerically using implicit finite difference scheme together with the Keller-box technique [29]. We have focused our attention on the evolution of the surface shear stress in terms of local skin friction coefficient and the rate of heat transfer in terms of local Nusselt number, the stream lines and the isotherms for selected values of parameters consisting of the magnetic parameter , Prandtl number Pr and the viscosity variation parameter .

#### 2. Formulation of the Problem

The steady two-dimensional laminar free convection boundary layer flow of a viscous incompressible and electrically conducting fluid along a vertical wavy surface in presence of uniform transverse magnetic field of strength with temperature dependent variable viscosity is considered. It is assumed that the wavy surface is electrically insulated and is maintained at a uniform temperature . Far above the wavy plate, the fluid is stationary and is kept at a temperature , where . The boundary layer analysis outlined below allows being arbitrary, but our detailed numerical work assumed that the surface exhibits sinusoidal deformations. The wavy surface may be described by where is the characteristic length associated with the wavy surface.

The geometry of the wavy surface and the two-dimensional Cartesian coordinate system are shown in Figure 1.

Figure 1: The coordinate system and the physical model.

The governing equations of such flow under the usual boundary layer and Boussinesq approximation can be written as where are the dimensional coordinates along and normal to the tangent of the surface and are the velocity components parallel to , is the Laplacian operator, is the acceleration due to gravity, is the dimensional pressure of the fluid, is the density, is the strength of magnetic field, is the electrical conduction, is the thermal conductivity, is the coefficient of thermal expansion, is the viscosity of the fluid depending on temperature of the fluid in the boundary layer region and is the specific heat due to constant pressure.

The boundary conditions relevant to the above problem are where is the surface temperature, is the ambient temperature of the fluid and is the pressure of fluid outside the boundary layer.

There are very few forms of viscosity variation available in the literature. Among them we have considered that one which is appropriate for liquid introduced by Hossain et al. [17] as follows: where is the viscosity of the ambient fluid and is a constant evaluated at the film temperature of the flow .

Following Yao [4], we now introduce the following nondimensional variables: where is the dimensionless temperature function and are the dimensionless velocity components parallel to and is the kinematic viscosity.

Introducing the above dimensionless dependent and independent variables into (2), the following dimensionless form of the governing equations are obtained after ignoring terms of smaller orders of magnitude in , the Grashof number defined in(5) In the above equations Pr, , and are, respectively, known as the Prandtl number, dimensionless viscosity variation parameter and dimensionless magnetic parameter, which are defined as It can easily be seen that the convection induced by the wavy surface is described by (6)β(9). We further notice that, (8) indicates that the pressure gradient along the -direction is , which implies that lowest-order pressure gradient along -direction can be determined from the inviscid flow solution. For the present problem this pressure gradient is zero. Equation (8) further shows that is and is determined by the left-hand side of this equation. Thus, the elimination of from (7) and (8) leads to The corresponding boundary conditions for the present problem then turn into Now we introduce the following transformations to reduce the governing equations to a convenient form: where is the pseudo similarity variable and is the stream function that satisfies (6) and is defined by Introducing the transformations given in (13) and into (11) and (9), we have The boundary conditions (12) now take the following form: In the above equations prime denotes the differentiation with respect to .

#### 3. Local Skin-Friction Coefficient and the Local Rate of Heat Transfer

In practical applications, the physical quantities of principle interest are the shearing stress and the rate of heat transfer in terms of the skin friction coefficient and Nusselt number , respectively, which can be written as where Using the transformations (13) and (19) into (18), the local skin friction coefficient and the rate of heat transfer in terms of the local Nusselt number take the following form:

#### 4. Method of Solution

Along with the boundary condition (17), the solution of the parabolic differential equations (15) and (16) will be found by using the implicit finite difference method together with Keller- box scheme [29], which is well documented by Cebeci and Bradshaw [30]. This method has been extensively used recently by Hossain et al. [14, 15, 17, 19β22]. To apply the aforementioned method, (15) and (16) their boundary condition (17) are first converted into the following system of first-order equations. For this purpose we introduce new dependent variables , , , and , , so that the transformed momentum and energy equations can be written as where , and and the boundary conditions (17) are The finite difference approximation according to box method then reduced the system of equation into nonlinear system of difference equation. These are as follows: The above equations are to be linearized by using Newtonβs quasilinearization method. Then linear algebraic equations can be written in a block matrix, which forms a coefficient matrix. The whole procedure, namely reduction to first-order followed by central difference approximations, Newtonβs quasilinearization method and the block Thomas algorithm, is well known as Keller-box method.

#### 5. Results and Discussion

In this paper the effect of temperature dependent variable viscosity on magnetohydrodynamic natural convection flow of viscous incompressible fluid along a uniformly heated vertical wavy surface has been investigated. Although there are four parameters of interest in the present problem, our main aim is to determine the effects of varying , the temperature dependent viscosity, the strength of magnetohydrodynamic field and Prandtl number Pr.

Numerical values of local shearing stress and the rate of heat transfer are calculated in terms of the skin friction coefficient and Nusselt number from (20) for natural convection flow of a variable viscosity electrically conducting fluid in the presence of magnetic field for a wide range of the axial distance variable starting from the leading edge. These are shown in tabular form in Table 1 and graphically in Figures 4β7 for different values of the aforementioned parameters , and Pr. From this table it can be concluded that if the value of Pr increases the values of the skin-friction coefficient decreases and the rate of heat transfer in terms of Nusselt number increases for any value of and . It is also shown that if the value of increases the values of decreases while it increases of and the rate of heat transfer in terms of Nusselt number decreases for any value of and .

Table 1: Comparison of Skin friction coefficient and the local rate of heat transfer against for the variation of Prandtl number Pr with and without effects of magnetic parameter and viscosity parameter with .

Figure 2 illustrate the effects of temperature dependent variable viscosity parameter and magnetic parameter on the development of streamlines which are plotted for the amplitude of the wavy surface and Prandtl number . When and , the problem discussed by Yao [4], where the viscosity is independent of temperature and in absence of magnetic parameter as shown in Figure 2(a). In this case . In Figure 2(b) it is found that an increase in the value of causes the effects of the wavy surface to be attenuated and the boundary layer becomes thinner where . For the case of fluid with constant viscosity () and in presence of magnetic field (i.e., ) which was studied by Alam et al. [14]. The magnetic field acting along the horizontal direction retards the fluid velocity. For this there creates a lorentz force by the interaction between the applied magnetic field and flow field. This force acts against the fluid flow and reduce the velocity distribution. On the other hand, when the value of increases similar thing happens and the maximum value of , that is, is 8.08 which is shown in Figure 2(c). The combined effects of and are shown in Figure 2(d). In this case the maximum value of is 5.23 decreases comparatively in presence of magnetic parameter with variable viscosity parameter .

Figure 2: Streamlines for (a) , (b) , (c) , (d) , while , and .

The influence of the variable viscosity parameter and the magnetic parameter on the isotherms profile for and are shown in Figure 3. As mentioned before, owing to the presence of the variable viscosity parameter and the magnetic field effects ( > 0 and ), the thermal state of the fluid increases, causing the thermal boundary layer to increase. In the down stream region the temperature distribution is negligible in this case. For absence of the variable viscosity parameter and the magnetic field , it can be observed that the opposite phenomenon happens.

Figure 3: Isotherms for (a) , (b) , (c) , (d) , while , and .
Figure 4: Variation of (a) skin friction coefficient and (b) rate of heat transfer against for different values of Prandtl number Pr while , , and .
Figure 5: Variation of (a) skin friction coefficient and (b) rate of heat transfer against for different values of Prandtl number Pr while and , and .
Figure 6: Variation of (a) skin friction coefficient and (b) rate of heat transfer against for different values of while , , and .
Figure 7: Variation of (a) skin friction coefficient and (b) rate of heat transfer against for different values of magnetic parameter while , , and .

The variation of the local skin friction coefficient and local rate of heat transfer for different values of Prandtl number Pr for and , are shown in Figures 4 and 5 while . The skin friction coefficient and local rate of heat transfer varies according to the slope of the wavy surface. This is due to the alignment of the buoyancy force , as shown in (15), which drives the flow tangentially to the wavy surface. It can be observed from Figure 4 that without effects of variable viscosity parameter and the magnetic parameter the skin friction coefficient and their amplitude reduce at a great extent for increasing values of the Prandtl number Pr and the rate of heat transfer increases. The local skin friction coefficient decreases by 32.94% and the rate of heat transfer increases by 52.65% as Pr increases from 0.7 to 7.0. Figure 5 it is noted that for the influence of the parameters and , the decreasing skin friction coefficient becomes slower in the downstream region. On the other hand, increasing the values of Prandtl number Pr speed up the decay of the temperature field away from the heated surface with a consequent increase in the rate of heat transfer and a reduction in the thermal boundary layer thickness. In this case the local skin friction coefficient decreases by 36.31% and the rate of heat transfer are increases by 49.31% as Pr increases from Pr = 0.7 to 7.0 which occur at point .

The effect of on the surface shear stress in terms of the local skin friction coefficient and the rate of heat transfer in terms of the local Nusselt number are depicted graphically in Figure 6 when , and . From this figure it can be noted that an increase in the variable viscosity variation parameter , the skin friction coefficient increases monotonically along the upward direction of the plate and to decrease of the heat transfer rates. Here it is concluded that for high viscous fluid the skin friction is large and the corresponding rate of heat transfer is slow. In Figure 6(a), the maximum values of local skin friction coefficient are 0.86640 and 1.54688 for and 5.0, respectively, which occur at the same point and it is seen that the local skin friction coefficient increases by 43.99% as increases from 0.0 to 5.0. Increasing values of lead to increase the amplitude of the . Again Figure 6(b) it is found that the rate of heat transfer decreases by 21.22% due to the increased value of .

The influence of the parameter , on the reduced local skin friction coefficient and local rate of heat transfer are illustrated in Figure 7 for Prandtl number , amplitude of wavy surface and viscosity parameter . From Figure 7 it observed that an increase in the magnetic parameter leads to decrease the local skin friction coefficient and local rate of heat transfer at different position of . The skin friction coefficient and the rate of heat transfer coefficient decreases by 23.95% and 12.01%, respectively, as increases from 0.0 to 2.0. The magnetic field acts against the flow and reduces the skin friction and the rate of heat transfer.

#### 6. Conclusion

For different values of relevant physical parameters including the magnetic parameter , the effect of natural convection flow of viscous incompressible fluid with temperature dependent variable viscosity along a uniformly heated vertical wavy surface has been studied. New variables to transform the complex geometry into a simple shape and were used a very efficient implicit finite difference method known as Keller-box scheme was employed to solve the boundary later equations. From the present investigation the following conclusions may be drawn: (i)The skin friction coefficient decreases for an increase of the Prandtl number Pr, over the whole boundary layer but significantly increase the rate of heat transfer. When with and without the effects of magnetic parameter and viscosity parameter are included. (ii)The effect of increasing viscosity parameter results in decreasing the local rate of heat transfer and increasing the local skin friction coefficient . (iii)An increase in the values of leads to decrease the local skin friction coefficient and the local rate of heat transfer .

#### Nomenclature

 πΆππ₯: Local skin friction coefficient πΆπ: Specific heat at constant pressure (JΒ·kgβ1Β·Kβ1) π: Dimensionless stream function π: Acceleration due to gravity (mΒ·s-2) Gr: Grashof number π: Thermal conductivity (WΒ·mβ1Β·Kβ1) πΏ: Characteristic length associated with the wavy surface (m) π: Magnetic parameter Nuπ₯: Local Nusselt number π: Pressure of the fluid (NΒ·mβ2) Pr: Prandtl number π: Temperature of the fluid in the boundary layer (K) ππ€: Temperature at the surface (K) πβ: Temperature of the ambient fluid (K) π’,π£: Dimensionless velocity components along the (π₯,π¦) axes (mΒ·sβ1) π₯,π¦: Axis in the direction along and normal to the tangent of the surface.

#### Greek Symbols

 πΌ: Amplitude of the surface waves π½: Volumetric coefficient of thermal expansion (Kβ1) π½0: Applied magnetic field strength π: Viscosity variation parameter π: Dimensionless similarity variable π: Dimensionless temperature function π: Stream function (m2Β·sβ1) π: Viscosity of the fluid (kgΒ·mβ1Β·sβ1) πβ: Viscosity of the ambient fluid π: Kinematic viscosity (m2Β·sβ1) π: Density of the fluid (kgΒ·mβ3) π0: Electrical conductivity ππ€: Shearing stress π(π₯): Surface profile function defined in (1).

#### Subscripts

 π€: Wall conditions β: Ambient conditions π₯: Differentiation with respect to π₯.

#### References

1. H. K. Kuiken, βMagneto-hydrodynamic free convection in a strong cross field,β Journal of Fluid Mechanics, vol. 40, no. 1, pp. 21β38, 1970.
2. E. M. Sparrow and R. D. Cess, βThe effect of a magnetic field on free convection heat transfer,β International Journal of Heat and Mass Transfer, vol. 3, no. 4, pp. 267β274, 1961.
3. J. Charraudeau, βInfluence de gradients de properties physiques en convection force application au cas du tube,β International Journal of Heat and Mass Transfer, vol. 18, no. 1, pp. 87β95, 1975.
4. L. S. Yao, βNatural convection along a vertical wavy surface,β Journal of Heat Transfer, vol. 105, no. 3, pp. 465β468, 1983.
5. S. G. Moulic and L. S. Yao, βNatural convection along a wavy surface with uniform heat flux,β Journal of Heat Transfer, vol. 111, pp. 1106β1108, 1989.
6. S. G. Moulic and L. S. Yao, βMixed convection along wavy surface,β Journal of Heat Transfer, vol. 111, pp. 974β979, 1989.
7. K. N. Mehta and S. Sood, βTransient free convection flow with temperature dependent viscosity in a fluid saturated porous medium,β International Journal of Engineering Science, vol. 30, no. 8, pp. 1083β1087, 1992.
8. N. G. Kafoussias and E. W. Williams, βThe effect of temperature-dependent viscosity on free-forced convective laminar boundary layer flow past a vertical isothermal flat plate,β Acta Mechanica, vol. 110, no. 1–4, pp. 123β137, 1995.
9. N. G. Kafoussias and D. A. S. Rees, βNumerical study of the combined free-forced convective laminar boundary layer flow past a vertical isothermal flat plate with temperature-dependent viscosity,β Acta Mechanica, vol. 127, no. 1–4, pp. 39β50, 1995.
10. D. A. S. Rees and I. Pop, βNote on free convection along a vertical wavy surface in a porous medium,β Journal of Heat Transfer, vol. 116, no. 2, pp. 505β508, 1994.
11. D. A. S. Rees and I. Pop, βFree convection induced by a horizontal wavy surface in a porous medium,β Fluid Dynamics Research, vol. 14, no. 4, pp. 151β166, 1994.
12. D. A. S. Rees and I. Pop, βFree convection induced by a vertical wavy surface with uniform heat flux in a porous medium,β Journal of Heat Transfer, vol. 117, no. 2, pp. 547β550, 1995.
13. F. M. Hady, A. Y. Bakier, and R. S. R. Gorla, βMixed convection boundary layer flow on a continuous flat plate with variable viscosity,β International Journal of Heat and Mass Transfer, vol. 31, no. 3, pp. 169β172, 1996.
14. K. C. A. Alam, M. A. Hossain, and D. A. S. Rees, βMagnetohydrodynamic free convection along a vertical wavy surface,β International Journal of Applied Mechanics and Engineering, vol. 1, pp. 555β566, 1997.
15. M. A. Hossain and D. A. S. Rees, βCombined heat and mass transfer in natural convection flow from a vertical wavy surface,β Acta Mechanica, vol. 136, no. 3, pp. 133β141, 1999.
16. C. Y. Cheng, βNatural convection heat and mass transfer near a vertical wavy surface with constant wall temperature and concentration in a porous medium,β International Communications in Heat and Mass Transfer, vol. 27, no. 8, pp. 1143β1154, 2000.
17. M. A. Hossain, M. S. Munir, and D. A. S. Rees, βFlow of viscous incompressible fluid with temperature dependent viscosity and thermal conductivity past a permeable wedge with uniform surface heat flux,β International Journal of Thermal Sciences, vol. 39, no. 6, pp. 635β644, 2000.
18. E. M. A. Elbashbeshy, βFree convection flow with variable viscosity and thermal diffusivity along a vertical plate in the presence of the magnetic field,β International Journal of Engineering Science, vol. 38, no. 2, pp. 207β213, 2000.
19. M. S. Munir, M. A. Hossain, and I. Pop, βNatural convection flow of a viscous fluid with viscosity inversely proportional to linear function of temperature from a vertical wavy cone,β International Journal of Thermal Sciences, vol. 40, no. 4, pp. 366β371, 2001.
20. M. S. Munir, M. A. Hossain, and I. Pop, βNatural convection with variable viscosity and thermal conductivity from a vertical wavy cone,β International Journal of Thermal Sciences, vol. 40, no. 5, pp. 437β443, 2001.
21. S. Kabir, M. A. Hossain, and D. A. S. Rees, βNatural convection of fluid with variable viscosity from a heated vertical wavy surface,β Zeitschrift fur Angewandte Mathematik und Physik, vol. 53, no. 1, pp. 48β57, 2002.
22. M. M. Molla, M. A. Hossain, and L. S. Yao, βNatural convection flow along a vertical wavy surface with uniform surface temperature in presence of heat generation/absorption,β International Journal of Thermal Sciences, vol. 43, no. 2, pp. 157β163, 2004.
23. R. Nasrin and M. A. Alim, βMHD free convection flow along a vertical flat plate with thermal conductivity and viscosity depending on temperature,β Journal of Naval Architecture and Marine Engineering, vol. 6, pp. 72β83, 2009.
24. G. Palani and K. Y. Kim, βNumerical study on a vertical plate with variable viscosity and thermal conductivity,β Archive of Applied Mechanics, vol. 80, pp. 711β725, 2009.
25. T. K. Aldoss, Y. D. Ali, and M. A. Al-Nimr, βMHD mixed convection from a horizontal circular cylinder,β Numerical Heat Transfer A, vol. 30, no. 4, pp. 379β396, 1996.
26. M. A. Al-Nimr and M. K. Alkam, βMagnetohydrodynamics transient free convection in open-ended vertical annuli,β Journal of Thermophysics and Heat Transfer, vol. 13, no. 2, pp. 256β265, 1999.
27. M. A. Al-Nimr and M. A. Hader, βMHD free convection flow in open-ended vertical porous channels,β Chemical Engineering Science, vol. 54, no. 12, pp. 1883β1889, 1999.
28. R. A. Damseh, M. Q. Al-Odat, and M. A. Al-Nimr, βEntropy generation during fluid flow in a channel under the effect of transverse magnetic field,β Heat and Mass Transfer, vol. 44, no. 8, pp. 897β904, 2008.
29. H. B. Keller, βNumerical methods in boundary layer theory,β Annual Review of Fluid Mechanics, vol. 10, pp. 417β433, 1978.
30. T. Cebeci and P. Bradshaw, Physical and Computational Aspects of Convective Heat Transfer, Springer, New York, NY, USA, 1984.