- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Advances in Mechanical Engineering

Volume 2012 (2012), Article ID 267170, 13 pages

http://dx.doi.org/10.1155/2012/267170

## Hydromagnetic Natural Convection Flow in a Fluid-Saturated Porous Medium with Nonuniform Heat Generation

School of Mathematical Sciences, Universiti Kebangsaan Malaysia, Selangor, 43600 UKM Bangi, Malaysia

Received 17 February 2012; Accepted 6 May 2012

Academic Editor: Oronzio Manca

Copyright © 2012 Habibis Saleh and Ishak Hashim. 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

Convective flows can develop naturally within the porous materials if they are subject to external heating with/without internal heat source. In this work, the heat is assumed to be generated internally within the porous enclosure at a rate proportional to a power of the temperature difference. This relation is an approximation of the state of some exothermic chemical reaction. An inclined magnetic field is imposed to the porous enclosure. Darcy, model is used to formulate the porous layer and finite difference method is applied to solve the governing equations. The obtained results indicate that strong internal heating can generate significant maximum fluid temperatures above the heated wall temperature, and location of the maximum fluid temperature moves towards the center of the top wall by strengthening the magnetic field. Local heating exponent does not have a major effect on the flow and temperature distributions as well as the heat transfer performance within the porous medium. The large magnetic fields, regardless of direction, are effective to suppress the convective flows and reduce the rate of heat transfer.

#### 1. Introduction

A porous medium means a material consisting of a solid matrix with an interconnected void [1]. The interconnectedness of the void (the pores) allows the flow of one or more fluids through the material. Examples of natural porous media are beach sand, sandstone, limestone, rye bread, wood, rocks, and soil. Fluid which flows in porous media has occupied the central stage in many fundamental heat transfer analysis and has received considerable attention over the last few decades. This interest is due to its wide range of applications, for example, high-performance insulation for buildings, chemical catalytic reactors, packed sphere beds, grain storage, and geophysical problems as frost heave. Porous media are also of interest in relation to the underground spread of pollutants, solar power collectors, and to geothermal energy systems. The literature concerning convective flow in porous media is abundant and representative studies in this may be found in the recent books by Ingham and Pop [2], Pop and Ingham [3], Ingham et al. [4], Ingham and Pop [5], Vafai [6], de Lemos [7], and Vadasz [8].

Convective flows can develop naturally within the porous materials if they are subject to some form of external heating. In some situations this material provides its own source of heat such as spontaneous combustion in coal stockpiles, heat removal from nuclear fuel debris in nuclear reactors, underground disposal of radioactive waste materials, and exothermic chemical reactions in packed-bed reactors. Moreover, this phenomenon can be encountered during the storage of agricultural products where heat is generated as a result of metabolism of the products [9]. Uncontrolled the spontaneous combustion cause a destruction and financial losses [10] and [11]. Thus, understanding and controlling the convective flows are clearly important for prevention of fires in the coal stockpiles.

The study of the interaction of the geomagnetic field with the fluid, such as oil, gas, and water, in the geothermal regions arises in geophysics. We can determine the flow and temperature distribution by solving the Navier-Stokes and energy equations. Early studies on hydromagnetic natural convection flow through a porous medium are due to Raptis et al. [12] and Raptis and Vlahos [13]. The effect of the transverse magnetic field within an inclined porous tall enclosure were investigated numerically by Vasseur et al. [14] and Bian et al. [15]. They found that the temperature and the velocity fields are significantly modified with the application of the magnetic. Khanafer and Chamkha [16] studied a transverse magnetic field and fluid heat generation effects. Grosan et al. [17] and Saleh et al. [18] extended this to include an inclined magnetic field. They reported an optimum reducing of the heat transfer performance that was obtained for a large magnetic field in the horizontal direction. Recently, analytical solutions of the magnetohydrodynamic problem were obtained by Raftari and Yildirim [19] and Raftari et al. [20]. The present work investigates the effect of an inclined magnetic field on natural convection in an enclosure filled with a porous medium saturated with an electrically conducting fluid having nonuniform internal heat generation. The special purpose is to find the optimum parameters in controlling the convective flows within the enclosure.

#### 2. Mathematical Formulation

We consider the steady, two-dimensional natural convection flow in a square region filled with an electrically conducting fluid-saturated porous medium (see Figure 1). The coordinate system employed is also depicted in this figure. The top and bottom surfaces of the convective region are assumed to be thermally insulated, the face is held at the constant reference temperature and the face has a constant prescribed temperature above . Heat is assumed to be generated internally within the porous medium at a rate proportional to , where is the local temperature. This relation, as explained by Mealey and Merkin [21], is an approximation of the state of some exothermic chemical reaction. Further, the angle of inclination of the magnetic field from the horizontal plane and measured positively in the counterclockwise direction is denoted by .

The viscous, radiation and Joule heating effects are neglected. The resulting convective flow is governed by the combined mechanism of buoyancy force, internal heat generation, and magnetic field. The magnetic Reynolds number is assumed to be small so that the induced magnetic field can be neglected compared to the applied magnetic field. Under these assumptions, the conservation equations for mass, momentum, energy, and electric transfer are as follows: where is the fluid velocity, is the fluid temperature, is the pressure, is the external magnetic field, is the electric current, is the electric potential, is the gravitational acceleration, is the permeability of the porous medium, is the effective thermal diffusivity, is the density, is the dynamic viscosity, is the coefficient of thermal expansion, is the reference temperature, is the local heat generation, is the local-heating exponent, is the specific heat at constant pressure, is the electrical conductivity, is the reference density, and is the associated electric field. As discussed by Garandet et al. [4, 5, 22] reduce to . The unique solution is since there is always an electrically insulating boundary around the enclosure. Thus, it follows that the electric field vanishes everywhere [23].

The above equations correspond to having the porous medium-modeled according to Darcy's law. Because of the complexity of pore geometries in a porous medium, Darcy's law has to be used to obtain any meaningful insights into the physics of flow in porous media. Darcy’s law has been verified by the results of many experiments. Theoretical backing for it has been obtained in various ways, with the aid of either deterministic or statistical models [1]. Vafai [6] have described that Darcy’s approximation is reasonable for the geothermal flows, except perhaps near boreholes [6]. Darcy’s law is valid only when the pore Reynolds number is of the order of 1 and for many practical applications, Darcy’s law is not valid, and boundary and inertial effects need to be accounted for [6]. For example, modeling fluid flows in the fuel cells. This is because the fuel cell consists of a number of distinct layers of different porosities. Darcy’s law assumes no effect of boundaries and the fluid velocity in Darcy's equation is determined by the permeability of the matrix. Therefore at the interface between the regions of different porosity in the fuel cell, particularly between the free fluid flow region, such as a gas flow channel and a permeable medium, then a discontinuity in the fluid velocity and/or the shear stress could emerge [5].

Eliminating the pressure term in [2] in the usual way, the governing equations [1–3] can be written as where is the magnitude of and is the kinematic viscosity of the fluid. The boundary conditions are Further, we introduce the following nondimensional variables: where . Introducing the stream function defined as and , and using expressions (11) in (7) and (9), we obtain the following partial differential equations in nondimensional form: subject to the boundary condition where is the Rayleigh number, is the internal heat generation parameter, and is the Hartmann number for the porous medium. It should be mentioned that corresponds to a horizontal magnetic field and corresponds to a vertical magnetic field, respectively.

Once we know the temperature we can obtain the rate of heat transfer from each of the vertical walls, which are given in terms of the mean Nusselt number as at the hot wall and at the cold wall.

#### 3. Numerical Method and Validation

An iterative finite difference procedure will be applied to solve (12), (13). Central difference method is applied for discretization of equations. Forward and backward difference schemes are used for the Neumann boundary conditions. The solution of algebraic equations was performed using Gauss-Seidel iteration. The unknown and were calculated until the following criteria of convergence were fulfilled: where is either or , represents the iteration number, and is the convergence criterion.

In this study, the convergence criterion is set at . The effect of grid resolution was examined in order to select the appropriate grid density as tabulated in Table 1 for , , , , and . The results indicate that an mesh can be used in the final computations. As a validation, our results for the mean Nusselt number compare well with that obtained by Revnic et al. [24] for case , , , and as shown in Table 2.

#### 4. Results and Discussion

In this section, we present numerical results for the streamlines and isotherms for various values of the magnetic field parameter, , the inclination angle, , the internal heat generation parameter, , and the exponent in the local-heating term as well as the Rayleigh number, . Mealey and Merkin [21] have investigated that the steady solutions exist only for a finite range of and . In this study, we performed calculation for range from 0 until 5 and from 0 until 300. The low is also parallel with the fact of coefficient permeability values from until inside the formula. Therefore, for physical reality, should be small.

Figure 2 shows the streamlines and isotherms for the no internal heating case () and horizontal magnetic field () with and a range of values of . In the absence of a magnetic field, the fluid motion as shown in Figure 2(a) is described as follows. Since the temperature of the left wall is higher than that of the fluid inside the enclosure, the wall transmits heat to the fluid and raises the temperature of fluid particles adjoining the left wall. As the temperature rises, the fluid moves from the left wall (hot) to the right wall (cold), falling along the cold wall, then rising again at the hot wall. This movement creates a clockwise vortex cell inside the enclosure and the isotherms start either from the hot wall or from the bottom wall and end at the top wall or cold wall. When the magnetic field is relatively strengthened (Figure 2(b)), the central streamlines cell are becoming upright and the isotherms stratification in the core diminishes. Further for large (Figure 2(c)), the core vortex is elongated vertically, that is, the flow circulation is progressively inhibited by the retarding effect of the electromagnetic body force (Lorentz force). While the isotherms are almost parallel to the vertical wall, this implies that conduction is dominant.

Figure 3 shows the streamlines and isotherms for the strong heating case () and horizontal magnetic field () with , and a range of values of . The situation for the absence of a magnetic field (), is similar to that of the previous case, but with a small increase in intensity of convective motion (as the value of each and increases). It notes that symbol denoted a secondary counterclockwise flow, but relative to small to be displayed in the figure. The internal heat generation enhances the flow near the cold wall and forces the streamlines in a denser distribution. On the other hand, the negative buoyancy caused by produces a vertical downward flow motion in the vicinity of the top corner hot wall. The maximum temperature increase above that on the heated wall, that is, . When the magnetic field is relatively strengthened (Figure 3(b)), the central streamline cells are also becoming upright and the maximum temperature drifts towards the center of the top wall. This is due to the fact that the hot fluid reaching the top left corner of the enclosure is unable to reject energy since the velocities are small. Finally, for large (Figure 3(c)), the core vortex is also elongated vertically and the isotherms are almost parallel to the vertical wall, this implies that again conduction is dominant.

The streamlines and isotherms for a range of values of with , , , are shown in Figure 4. The values of decreases slightly as is increased. The trends for moderate and shown in Figure 4 suggest that does not have a major effect on the flow and heat transfer within the porous medium.

Figure 5 shows the influence of the inclination angle of the applied magnetic field on flow and temperature distributions for , , , . The influence of on flow and temperature distributions is apparent from these figures. The flow fields become distorted with stronger convective flows developing on the vertical walls. The maximum temperature is above and is located on the top wall, though the value of is reduced as is increased, from for to for .

The effect of on the stream function and temperature are shown in Figure 6 for , , , . As is increased, the streamline pattern becomes mildly distorted. This in turn distorts the isotherms, giving higher temperatures in a region towards the top of the heated wall.

Figure 7 presents relationship between the mean Nusselt number along the hot wall and along the cold wall against the Rayleigh number for the case , , and several values of and . The variations of and with for several values of are shown in Figure 8 for the case , , . Naturally, the heat transfer is maximum in the absence of a magnetic field, because the convection is maximum for this situation. In general, and initially decrease steeply with . As the value of is made larger, the strength of the convective motion is progressively suppressed and for . For a fixed moderate value of , increasing has the effect of decreasing and increasing .

Figure 9 shows the variations of and with when is varied for the case , , . Generally, for all values of , increasing has the effect of increasing and decreasing , contrasting to the effect of increasing above. The inclination angle of the magnetic field has only a moderate effect on and for a fixed small value of as shown in Figure 10.

#### 5. Conclusions

The present numerical study exhibits many interesting features concerning the effect of the inclined magnetic fields on natural convection in square enclosure filled with a porous medium having nonuniform internal heat source. Detailed numerical results for flow field, temperature distribution, and heat transfer performance have been presented in graphical form. The results of the numerical analysis lead to the following conclusions.

Strong internal heating can generate significant maximum fluid temperatures above the heated wall temperature and location of the maximum fluid temperature moves towards the center of the top wall by strengthening the magnetic field. Local heating exponent does not have a major effect on the flow and temperature distributions as well as the heat transfer performance within the porous medium. In general, effect of magnetic field is to retard the convective flows circulation within the porous enclosure. Combination of the large of applied magnetic field in any direction having a strong internal heat source with a small local heating exponent was found to be most effective in suppressing the heat transfer performance along the hot wall. However, combination of the large of applied magnetic field in any direction having a low internal heat source with a high local heating exponent was found to be most effective in suppressing the heat transfer performance along the cold wall.

The findings may be useful to find an effective way in controlling the convective flows due to external and internal heating. The results of the problem are also useful in the study of movement of oil or gas and water through the reservoir of an oil or gas field in the migration of underground water.

#### Nomenclature

: | Applied magnetic field |

: | Specific heat at constant pressure |

: | Gravitational acceleration vector |

: | Internal heat generation |

: | Hartmann number |

: | Electric current |

: | Width of the cavity |

: | Permeability of the porous medium |

: | Mean Nusselt numbers |

: | Local heating exponent |

: | Rayleigh number |

: | Fluid temperature |

: | Reference temperature |

: | Prescribed temperature above reference temperature |

: | Velocity component in the -direction and -direction |

: | Dimensionless velocity component in the -directionand -direction |

: | Space coordinates |

: | Dimensionless space coordinates. |

#### Greek Symbols

: | Effective thermal diffusivity |

: | Thermal expansion coefficient |

: | Angle of inclination magnetic field |

: | Stream function |

: | Dimensionless temperature |

: | Density of fluid |

: | Density of fluid |

: | Dynamic viscosity |

: | Kinematic viscosity |

: | Electrical conductivity |

: | Electric potential. |

#### References

- D. A. Nield and A. Bejan,
*Convection in Porous Media*, Springer, New York, NY, USA, 2nd edition, 1999. - D. B. Ingham and I. Pop,
*Transport Phenomena in Porous Media*, Elsevier Science, Oxford, UK, 1998. - I. Pop and D. B. Ingham,
*Convective Heat Transfer: Mathematical and Computational Modelling of Viscous Fluids and Porous Media*, Pergamon, Oxford, UK, 2001. - D. B. Ingham, A. Bejan, E. Mamut, and I. Pop,
*Emerging Technologies and Techniques in Porous Media*, Kluwer, Dordrecht, The Netherlands, 2004. - D. Ingham and I. Pop,
*Transport Phenomena in Porous Media III*, Elsevier Science, Oxford, UK, 2005. - K. Vafai,
*Handbook of Porous Media*, Taylor and Francis, New York, NY, USA, 2nd edition, 2005. - M. J. S. deLemos,
*Turbulence in Porous Media: Modeling and Applications*, Elsevier, Oxford, UK, 2006. - P. Vadasz,
*Emerging Topics in Heat and Mass Transfer in Porous Media*, Springer, New York, NY, USA, 2008. - M. Haajizadeh, A. F. Ozguc, and C. L. Tien, “Natural convection in a vertical porous enclosure with internal heat generation,”
*International Journal of Heat and Mass Transfer*, vol. 27, no. 10, pp. 1893–1902, 1984. View at Scopus - K. Brooks and D. Glasser, “A simplified model of spontaneous combustion in coal stockpiles,”
*Fuel*, vol. 65, no. 8, pp. 1035–1041, 1986. View at Scopus - K. Brooks, V. Balakotaiah, and D. Luss, “Effect of natural convection on spontaneous combustion of coal stockpiles,”
*AIChE Journal*, vol. 34, no. 3, pp. 353–365, 1988. View at Scopus - A. Raptis, C. Massalas, and G. Tzivanidis, “Hydromagnetic free convection flow through a porous medium between two parallel plates,”
*Physics Letters A*, vol. 90, no. 6, pp. 288–289, 1982. View at Scopus - A. Raptis and J. Vlahos, “Unsteady hydromagnetic free convective flow through a porous medium,”
*Letters in Heat and Mass Transfer*, vol. 9, no. 1, pp. 59–64, 1982. View at Scopus - P. Vasseur, M. Hasnaoui, E. Bilgen, and L. Robillard, “Natural convection in an inclined fluid layer with a transverse magnetic field. Analogy with a porous medium,”
*Journal of Heat Transfer*, vol. 117, no. 1, pp. 121–129, 1995. View at Scopus - W. Bian, P. Vasseur, E. Bilgen, and F. Meng, “Effect of an electromagnetic field on natural convection in an inclined porous layer,”
*International Journal of Heat and Fluid Flow*, vol. 17, no. 1, pp. 36–44, 1996. View at Publisher · View at Google Scholar · View at Scopus - K. M. Khanafer and A. J. Chamkha, “Hydromagnetic natural convection from an inclined porous square enclosure with heat generation,”
*Numerical Heat Transfer. Part A*, vol. 33, no. 8, pp. 891–910, 1998. View at Scopus - T. Grosan, C. Revnic, I. Pop, and D. B. Ingham, “Magnetic field and internal heat generation effects on the free convection in a rectangular cavity filled with a porous medium,”
*International Journal of Heat and Mass Transfer*, vol. 52, no. 5-6, pp. 1525–1533, 2009. View at Publisher · View at Google Scholar · View at Scopus - H. Saleh, R. Roslan, and I. Hashim, “Natural convection in a porous trapezoidal enclosure with an inclined magnetic field,”
*Computers and Fluids*, vol. 47, no. 1, pp. 155–164, 2011. View at Publisher · View at Google Scholar · View at Scopus - B. Raftari and A. Yildirim, “The application of homotopy perturbation method for MHD flows of UCM fluids above porous stretching sheets,”
*Computers and Mathematics with Applications*, vol. 59, no. 10, pp. 3328–3337, 2010. View at Publisher · View at Google Scholar · View at Scopus - B. Raftari, S. T. Mohyud-Din, and A. Yildirim, “Solution to the MHD flow over a non-linear stretching sheet by homotopy perturbation method,”
*Science China: Physics, Mechanics and Astronomy*, vol. 54, no. 2, pp. 342–345, 2011. View at Publisher · View at Google Scholar · View at Scopus - L. R. Mealey and J. H. Merkin, “Steady finite Rayleigh number convective flows in a porous medium with internal heat generation,”
*International Journal of Thermal Sciences*, vol. 48, no. 6, pp. 1068–1080, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. P. Garandet, T. Alboussiere, and R. Moreau, “Buoyancy driven convection in a rectangular enclosure with a transverse magnetic field,”
*International Journal of Heat and Mass Transfer*, vol. 35, no. 4, pp. 741–748, 1992. View at Scopus - S. Alchaar, P. Vasseur, and E. Bilgen, “Natural convection heat transfer in a rectangular enclosure with a transverse magnetic field,”
*Journal of Heat Transfer*, vol. 117, no. 3, pp. 668–673, 1995. View at Scopus - C. Revnic, T. Grosan, I. Pop, and D. B. Ingham, “Magnetic field effect on the unsteady free convection flow in a square cavity filled with a porous medium with a constant heat generation,”
*International Journal of Heat and Mass Transfer*, vol. 54, no. 9-10, pp. 1734–1742, 2011. View at Publisher · View at Google Scholar · View at Scopus