Volume 2013 (2013), Article ID 183079, 13 pages
http://dx.doi.org/10.1155/2013/183079
Research Article

## Magnetic Field Effect on Natural Convection in a Porous Cavity Heating from below and Salting from Side

1School of Mathematical Sciences, National University of Malaysia, 43600 Bangi, Selangor, Malaysia
2Department of Mechanical, Materials and Manufacturing Engineering, Faculty of Engineering, The University of Nottingham, 43500 Semenyih, Selangor, Malaysia

Received 10 July 2012; Accepted 17 December 2012

Copyright © 2013 A. A. Altawallbeh 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.

#### Abstract

The effect of magnetic field on double-diffusive natural convection in a square cavity filled with a fluid-saturated porous medium is studied numerically. The bottom wall is fully heated at a constant temperature, and the top wall is maintained at a constant cold temperature. The right wall is fully salted to a high concentration, while the left wall is fully salted at a lower concentration than the right one. A magnetic force is applied on the cavity along the gravity force direction. The Darcy model is used for the mathematical formulation of the fluid flow through porous media. The governing equations for heat and mass transfer are solved using the finite volume method. The governing parameters of the present study are Rayleigh number (Ra), Lewis number (Le), buoyancy ratio (), and Hartmann number (Ha). The numerical solutions were studied in the range of −, , , and . The results were discussed considering the effect of these parameters on the heat and mass transfer processes. The results were presented in terms of streamlines, isotherms, isoconcentration, average Nusselt number, and average Sherwood number for different values of the governing parameters. In general, it has been found that the increase of magnetic force has an effect to retard the strength of the flow inside the cavity and reduce the heat and mass transfer processes. For high Hartmann number, the flow is almost suppressed.

#### 1. Introduction

Double-diffusion natural convection inside cavities filled with a fluid-saturated porous medium has received considerable attention in recent years, on account of its wide applications in natural systems and industries, for example, the effect of contaminants on lakes and underground water, atmospheric pollution, chemical processes, food processing, energy storage, and materials processing, to name but a few. The study of double-diffusive convection in a fluid-saturated porous medium under the influence of magnetic field has importance in many fields of science, engineering, and technology. It is useful in commercial production of the magnetic fluids, in chemical engineering and in the performance of petroleum reservoirs, where the geothermal areas are influenced by the magnetic field of the earth. The investigation is also useful in geophysics, to study the earth’s core, where the molting fluid is conducting, which can become convectively unstable as a result of differential diffusion. The process of manufacturing materials in industrial problems and microelectronic heat transfer devices involves an electrically conducting fluid subjected to a magnetic field. In that case the fluid experiences a Lorentz force, and its effect is to reduce the flow velocities [1]. This in turn affects the rate of heat and mass transfer. It is, therefore, important to study the detailed characteristics of transport phenomena in such a process so that good quality product can be developed with improved design in the manufacturing processes.

In double diffusion, the buoyancy force is affected not only by the difference of temperature, but also by the difference of concentration of the fluid. A detailed discussion on double diffusion natural convection can be found in the books by Nield and Bejan [2], Ingham and Pop [3], and Vafai [4]. Magneto-hydrodynamics (MHD) is the science of the motion of electrically conducting fluids under the influence of applied magnetic forces. The symbiotic interaction between the fluid velocity field and the electromagnetic forces gives rise to flow scenarios; the magnetic field affects the motion. Electromagnetic field has an important influence on the hydrodynamics. One of the main purposes of the electromagnetic control is to stabilize the flow and suppress oscillatory instabilities [5].

The phenomenon of double-diffusive natural convection in pure fluids and fluid-saturated porous media has received more attention in the literature due its important applications as mentioned before. For some related problems with just a heat effect (pure convection) from below, see, for example, Elder [6], Robillard et al. [7], Sezai and Mohamad [8], and Calcagni et al. [9], and for heat and mass transfer effect see, for example, Bejan and Khair [10]. The review of the literature on double-diffusive natural convection in fluids filled with a porous medium shows that most of studies focus on the thermosolutal effect inside cavities or enclosures where different boundary conditions apply. Studies by Bahloul et al. [11], Zhao et al. [12], Liu et al. [13], Zhao et al. [14], and Mamou et al. [15], for example, used boundary conditions of heating and salting from side. Other studies by Bourich et al. [16] and Mansour et al. [17] investigated heating from below and salting from sides. Some studies, for example, Zhao et al. [18] and Alloui et al. [19], looked at double-diffusive natural convection in boundary conditions of heating and salting from below. Many studies used a Darcy model to formulate the momentum equation; see for example, Liu et al. [13] Zhao et al. [14], Mamou et al. [15], Alloui et al. [19], Trevisan and Bejan [20], Mamou et al. [21], Bourich et al. [22], Saeid and Pop [23], Saleh and Hashim [24], and Akbal and Baytaş [25]. For non-Darcy, see for instance, Zhao et al. [12], Zhao et al. [18], and Goyeau et al. [26]. For an analytical study, see, for example, Bahloul et al. [11], Trevisan and Bejan [27], Amahmid et al. [28], and Bhadauria and Srivastava [29].

In the present work, the effect of magnetic field on the heat and mass transfer in a square porous cavity heated from below and salted from side is considered. Related studies can be found in, for example, Mahmud and Fraser [1], Revnic et al. [5], Garandet et al. [30], Alchaar et al. [31], Chamkha and Al-Naser [32], Ece and Buyuk [33], Saleh and Hashim [24], Teamah et al. [34], and Grosan et al. [35]. Mahmud and Fraser [1] studied magnetohydrodynamic-free convection and entropy generation in a square porous cavity with heating from side and adiabatic horizontal walls. They found that increases in the value of magnetic force have a tendency to retard the fluid motion inside the cavity, and the entropy generation rate is decreased in magnitude as the magnetic force is introduced and strengthened. Teamah et al. [34] discussed double-diffusive natural convection in an inclined rectangular enclosure in the presence of magnetic field with heating and salting from side. Grosan et al. [35] discussed the steady magnetohydrodynamics-free convection in a fluid-saturated porous medium with internal heat generation with an inclined magnetic field. They found that the convection modes within the cavity were found to depend upon both the strength and the inclination of the magnetic field. In addition, the applied magnetic field in the horizontal direction was found to be most effective in suppressing the convection flow for a stronger magnetic field in comparison with the vertical direction. Revnic et al. [5] studied the effect of an inclined magnetic field and heat generation on unsteady natural convection in a square cavity filled with a fluid-saturated porous medium with heating from side. They found that, for a fixed value of the Rayleigh number, the average Nusselt number decreases to the steady case when the Hartmannn number is increased. In addition, for all values of the governing parameters considered, the average Nusselt number reaches the steady-case faster when the magnetic field is parallel to the vertical walls of the cavity.

Most of those related papers discussed natural convection with heating the vertical walls and without salting source. The aim of this work is to study the effect of a magnetic field on double-diffusive natural convection in a square cavity filled with a fluid-saturated porous medium with heating from below and salting from side. Parametric study is carried out to study the effect of the following parameters: Hartmann number (), buoyancy ratio (), Lewis number (), and Rayleigh number () on the heat and mass transfer processes. The results were presented in terms of streamlines, isotherms, iso-concentration, average Nusselt number, and average Sherwood number for different cases in the parametric study.

#### 2. Mathematical Formulation

The physical model under consideration is a square cavity of length filled with an electrically conducting fluid-saturated porous medium as shown in Figure 1. The bottom and the right walls are fully heated and fully salted, respectively, while the top and the left walls are cooled and salted to a constant temperature and concentration, respectively. The two vertical walls are adiabatic and impermeable. A magnetic force affects the fluid in the cavity along the gravity force direction. The third dimension is assumed to be large enough to consider a two-dimensional simplification. The full elliptic equations are considered, while, the simplified boundary layer (parabolic) equations are not considered. The porous medium is assumed to be in local thermal equilibrium. In addition, the properties of the fluid and the porous medium are assumed to be constants. To simplify the problem, the Soret and Duffour effects are neglected. The Darcy model is employed, which is valid for the low values of Rayleigh number, and the Boussinesq approximation for the relationship of the density and the temperature and concentration is adopted: There are many industrial applications of the electrically conducting fluids in double diffusion of cross-thermosolute gradients in the presence of a magnetic forces, for example, crystal growth, electronic packages, metallurgical applications involving continuous casting, and solidification of metal alloys and others.

Figure 1: Physical model with coordinate system.

The mathematical model can be derived based on the following dimensionless variables: The dimensionless conservation equations for mass, momentum, energy, and concentration for the two-dimensional flow are where , , , . The velocity components and defined in terms of the stream function as and . The dimensionless boundary conditions are as follows: The average Nusselt number on the bottom wall and the average Sherwood number on the right wall are calculated using the following formulas:

#### 3. Numerical Method and Validation

The finite volume method (FVM) was used to solve the governing equations (3) subjected to the boundary conditions (4). The power law scheme is applied for convection-diffusion formulation. The details of this method can be found in the books by Patankar [36] and Versteeg and Malalasekera [37]. The solution is initialized by assuming stagnant fluid in the cavity with . The line-by-line iteration using the tri diagonal matrix algorithm (TDMA) was used to solve the discretized equations. The iterations are continued until the following condition is satisfied: where refers to , or and denotes the iteration step. The effect of the grid size was tested in terms of the average Nusselt number and average Sherwood number with and Lewis number , with two values of buoyancy ratio (, and ) as shown in Table 1. It is found that is an appropriate size for most calculations in the present problem. A good agreement was found between the present results and the results presented in the literature as shown in Table 2. The average Nusselt number () is compared with a literature using two different values of with heating from below and without the effect of concentration buoyancy force (). The results presented in Table 3 show a good agreement between the present results and those from the literature. The calculated average Nusselt number with various values of Hartmann number is presented in Table 4. Again a good agreement was found between the present results and the results presented by Mahmud and Fraser [1]. Based on the above comparisons, the numerical method and the developed code are accurate and valid for the parametric study.

Table 1: Grid size test for , with two values of .
Table 2: Comparison of for side heating with some previous results.
Table 3: Code validation in the case of pure convection () and with the heating from below.
Table 4: Comparison of for side heating with Mahmud and Fraser [1] results with .

#### 4. Results and Discussion

The mathematical formulation shows that there are four parameters affecting the heat and mass transfer process when a magnetic field is applied on the fluid in a square cavity filled with a porous medium. These parameters are , and . The numerical simulation was conducted and limited under different ranges of these parameters: , , , and . The results were presented in terms of streamlines, isotherms, iso-concentration, average Nusselt number, and average Sherwood number.

##### 4.1. Effect of the Buoyancy Ratio Parameter ()

In this section the effect of buoyancy ratio parameter () in double-diffusive natural convection is considered. The numerical simulation was conducted using the following buoyancy ratio , and , with different values of Hartmann number , and . The values of Rayleigh number and Lewis number were fixed with and .

Streamlines, isotherms, and iso-concentration inside the cavity shown in Figure 2 are for different values of Hartmann number, and with . Figure 2(a) shows the case without any effect of magnetic field . Isotherms spread wide inside the cavity. The buoyancy force enhances the fluid to flow rising up along the right wall until the top wall, then circulating and going down along the left wall making a circulation with counterclockwise direction as shown in the streamlines with . Since , iso-concentration spreads under the effect of thermal buoyancy force and follow the same direction. At the flow became slower and as shown in Figure 2(b). The spread of isotherms is not wide and the high temperatures are still in the bottom half of the cavity. The shape of isotherms is different from the previous case, which change the type of solution from one-circulation cell in the previous case to a triple-circulation cell in this case, with two circulations in counterclockwise direction and one circulation in the middle in clockwise direction. It is noticeable that a big difference occurred between the previous values of and when and the present when , where the previous and are , respectively, while the present and are and respectively.

Figure 2: (left) Streamlines, (middle) isotherms, and (right) iso-concentration in the case of , , . (a) , , , , (b) , , , .

Figure 3 shows the streamlines, isotherms, and iso-concentration for with different values of Hartmann number (, and ). When (see Figure 3(a)), both isotherms and iso-concentration spread inside the cavity and cooperating together to make a circulation with counterclockwise direction with , , and . At the introduction of magnetic force (), isotherms and iso-concentration still spread widely inside the cavity and working together to enhance the fluid to flow rising up along the right wall, then circulating and going down along the left wall making a circulation with counterclockwise direction. It is noticeable that the magnetic force tends to retard the fluid flow, where , , and as shown in Figure 3(b). As the magnetic force increased to , as shown in Figure 3(c), the effect of magnetic force is relatively big. The variation of isotherms and iso-concentration becomes less intensive and the strength of flow becomes very small with . Hence, the processes of heat and mass transfer decreased with and .

Figure 3: (left) Streamlines, (middle) isotherms, and (right) iso-concentration in the case of , , . (a) , , , , (b) , , , , and (c) , , , .

Streamlines, isotherms, and iso-concentration shown in Figure 4 are for and for different values of Hartmann number , and . The buoyancy ratio has the same value as the previous case but in the opposite direction. It is noticeable from Figures 4(a)4(c) that the flow has the same behavior as the previous case but in the opposite direction, and this is due to the opposite direction of buoyancy ratio. The strength of the flow, heat transfer, and mass transfer have almost the same values with the previous case as shown in Figures 4(a)4(c). The circulations in all cases reverse their directions to clockwise direction, where the thermal and solutal buoyancy forces cooperate together to enhance the fluid to rise up along the left wall and then circulating and going down along the right wall making a circulation with clockwise direction with small skewness to the top right corner and the bottom left corner.

Figure 4: (left) Streamlines, (middle) isotherms, and (right) iso-concentration in the case of , , . (a) , , , , (b) , , , , and (c) , , , .

The streamlines, isotherms, and iso-concentration are shown in Figures 5(a), 5(b), and 5(c) for and , and , respectively. When , both isotherms and iso-concentration spread inside the cavity as shown in Figure 5(a). Both buoyancy forces are working together to rise the fluid up along the right wall and then circulating and going down along the left wall making a circulation with counterclockwise direction with , and . The flow is progressively inhibited due to the retarding effect of magnetic force as shown in Figures 5(b) and 5(c). A sharp decreasing in the flow circulation strength is noticeable when the magnetic field is strengthened to as shown in Figure 5(c), with and , and . The spread and distortion of isotherms and iso-concentration become relatively less, but still working together to make a circulation in a counterclockwise direction with two egg-shaped cores appeared for as shown in Figure 5(c).

Figure 5: (left) Streamlines, (middle) isotherms, and (right) iso-concentration in the case of , , . (a) , , , , (b) , , , , and (c) , , , .
##### 4.2. Effect of Lewis Number (Le)

In this section the effect of Lewis number on double-diffusive natural convection is considered. The numerical simulation was conducted under different values of Lewis number ( and ), with different values of Hartmann number (, 1 and ). The values of Rayleigh number and buoyancy ratio were fixed with and .

Figure 6 presents the streamlines, isotherms, and iso-concentration for and different values of Hartmann number and . For the case without the effect of magnetic force (see Figure 6(a)), isotherms and iso-concentration are cooperating together to rise the fluid up along the right wall to the top wall an then circulating and going down along the left wall making a circulation with counterclockwise direction. The strength of the circulation is relatively high with . Distortion and spread of isotherms are clear and wide in comparison with iso-concentration where and , and this is due to the small value of Lewis number. As the introduction of magnetic force and (see Figures 6(b) and 6(c)), the flow is progressively inhibited due to retarding effect of magnetic force. When , the circulation strength is decreased and . While , a sharp decreasing occurred in the circulation strength with . In this case the flow is almost driven by the thermal buoyancy force making a circulation with counterclockwise direction. A one egg-shaped core appeared. Distortion and spread of isotherms and iso-concentration are decreased, as shown in Figure 6(c). Due to a small value of Lewis number, iso-concentration lines became almost parallel and close to the pure diffusion case as shown in Figure 6(c). Heat and mass transfer also decreased and reached to small values when with and .

Figure 6: (left) Streamlines, (middle) isotherms, and (right) iso-concentration in the case of , , . (a) , , , , (b) , , , , (c) and , , , .

Streamlines, isotherms, and iso-concentration are presented in Figure 7 for with different values of Hartmann number ( and ). For , both isotherms and iso-concentration are spread widely inside the cavity, and cooperating together to enhance the fluid to rise up along the right wall to the top wall and then circulating and going down along the left wall making a circulation with counterclockwise direction. The strength of the circulation is relatively high with . Figures 7(b) and 7(c) illustrate the flow under the effect of magnetic force ( and ). The flow is progressively inhibited due to retarding effect of magnetic force on the fluid flow. Distortion of isotherms and iso-concentration is a start to disappear with increasing Hartmann number. As Hartmann number increases, isotherms and iso-concentration lines inside the cavity approache more and more towards the conduction-like distribution pattern as shown in Figure 7(c). A sharp decreasing of circulation strength where the becomes for . In addition, a two egg-shaped cores appeared as shown in Figure 7(c).

Figure 7: (left) Streamlines, (middle) isotherms, and (right) iso-concentration in the case of , , . (a) , , , , (b) , , , , and (c) , , , .
##### 4.3. Heat and Mass Transfer

In this section the heat and mass transfer processes are discussed and presented in terms of average Nusselt number () and average Sherwood number (). The numerical simulation was conducted with different values of the governing parameters. In each case some of the parameters are varied and the others are fixed.

Figure 8 illustrates the relation between buoyancy ratio () and average Nusselt number () and average Sherwood number (). The values of are varied with and Hartmann number has different values where , and , while the other parameters were fixed with and . A mathematical explanation of the selected range of parameters was studied by Mahmud and Fraser [1]. The relation between and is shown in Figure 8(a). It is noticeable that the highest values of occurred for where there is no effect of the magnetic field on the fluid flow. As increases the values of decreased until reaching the smallest values when and , where the value of is almost , which means that the thermal pure diffusion case occurred. It is observed that, for and , the values of are close to its asymptotic value and the flow is almost suppressed. Figure 8(b) shows the relation between and . Similarly, as increases, the value of increases for any value of . On the other hand, as is increased, the values of decrease for any values of . It is observed that when and , the mass transfer process is suppressed and average Sherwood number is close to 1.

Figure 8: The effect of buoyancy ratio () on (a) the average Nusselt number , and (b) the average Sherwood number with different values of .

Figure 9 illustrates the relation between average Nusselt number () and average Sherwood number () with Lewis number . The values of and are changeable, where and , and 10 while and were fixed at 1 and 100, respectively.

Figure 9: The effect of Lewis number () on (a) the average Nusselt number and (b) the average Sherwood number with different values of .

The relation between and is shown in Figure 9(a). As Lewis number increases, the values of average Nusselt number decreased for all values of Hartmann number. For high values of Hartmann number and , the process of heat transfer is almost suppressed and the average Nusselt number is close to its asymptotic value .

The relation between and is shown in Figure 9(b). A direct relation is observed between and the values of for all values of . For , the mass transfer process is very low and the values of are close to 1. Also, the same behavior is observed for and for small values of Lewis number .

Figure 10 illustrates the relation between average Nusselt number and average Sherwood number with Rayleigh number, under the variation of Hartmann number. The Rayleigh number and Hartmann number vary: and , and . The buoyancy ratio and Lewis number were fixed at and in the results presented in Figure 10. A direct relation is observed from Figure 10 between and with for all values of Hartmann number. From Figure 10(a), very small values of are found when , and the values are close to its asymptotic value . Otherwise, the heat transfer process is active and is sometimes relatively high. The average Sherwood number has also a direct relation with the Rayleigh number as shown in Figure 10(b) for each value of the Hartmann number. For , the mass transfer process is very low, and the values of are close to 1 for . Relatively high values of occurred for , and .

Figure 10: The effect of Rayleigh number () on (a) the average Nusselt number and (b) the average Sherwood number with different values of .

As a result, as the Hartmann number increases, the values of the average Nusselt number and average Sherwood number decrease for each value of the buoyancy ratio, Lewis number, and Rayleigh number. In some cases as discussed before with high Hartmann number , they become close to 1 and approach to the pure diffusion case.

#### 5. Conclusion

The double-diffusive natural convection in a square cavity filled with a porous medium is studied numerically. The bottom and right walls are fully heated and fully salted, respectively. The top and left walls are heated and salted, respectively, with high temperature and high concentration on the bottom and right walls, respectively. The magnetic field affects the fluid flow in the cavity on the direction of gravity acceleration. The results are generated and presented in terms of the average Nusselt number and average Sherwood number for different values of buoyancy ratio, Lewis number, and Rayleigh number with a combined effect of Hartmann number. Streamlines, isotherms, and iso-concentrations are also presented for some cases. In general, it is found that the increase of magnetic force has an effect to retard the strength of the flow inside the cavity. In addition, the distribution of isotherms and iso-concentration becomes less intensive and for some cases close to parallel lines. As the Hartmann number increases the heat, mass transfer processes are decreased, where the values of the average Nusselt number and average Sherwood number decreased to small values. For high values of Hartmann number, the heat and mass transfer are suppressed, and the values of average Nusselt number and average Sherwood number approach their asymptotic values. For any fixed value of the Hartmann number, an increasing of buoyancy ratio and Rayleigh number will enhance the heat and mass transfer processes, while increasing the Lewis number will decrease heat transfer and enhance mass transfer for any fixed value of the Hartmann number.

#### Nomenclature

 : The magnitude of the external magnetic field : Dimensional concentration : Mass diffusivity : Gravitational acceleration : Hartmann number : Permeability of porous medium : Lewis number : Buoyancy ratio parameter : Length of the cavity : Nusselt number : Rayleigh number : Dimensionless concentration : Sherwood number : Temperature : Dimensional velocity components in and directions, respectively : Dimensionless velocity components in and directions, respectively : Dimensional space coordinates : Dimensionless space coordinates.
Greek Symbols
 : Thermal diffusivity of porous medium : Concentration expansion coefficient : Thermal expansion coefficient : Dynamic viscosity : Dimensionless temperature : Kinematic viscosity : Dimensionless stream function : Density : The electrical conductivity.

#### References

1. S. Mahmud and R. A. Fraser, “Magnetohydrodynamic free convection and entropy generation in a square porous cavity,” International Journal of Heat and Mass Transfer, vol. 47, no. 14–16, pp. 3245–3256, 2004.
2. D. A. Nield and A. Bejan, Convection in Porous Media, Springer, New York, NY, USA, 3rd edition, 2006.
3. D. B. Ingham and I. Pop, Transport Phenomena in Porous Media, vol. 3, Pergamon, Oxford, UK, 2005.
4. K. Vafai, Handbook of Porous Media, Taylor & Francis, LLC, U.S., Boca Raton, Fla, USA, 2nd edition, 2005.
5. 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.
6. J. W. Elder, “Transient convection in a porous medium,” The Journal of Fluid Mechanics, vol. 47, pp. 609–623, 1967.
7. L. Robillard, C. H. Wang, and P. Vasseur, “Multiple steady states in a confined porous medium with localaized heating from below,” Numerical Heat Transfer, vol. 13, no. 1, pp. 91–110, 1988.
8. I. Sezai and A. A. Mohamad, “Natural convection from a discrete heat source on the bottom of a horizontal enclosure,” International Journal of Heat and Mass Transfer, vol. 43, no. 13, pp. 2257–2266, 2000.
9. B. Calcagni, F. Marsili, and M. Paroncini, “Natural convective heat transfer in square enclosures heated from below,” Applied Thermal Engineering, vol. 25, no. 16, pp. 2522–2531, 2005.
10. A. Bejan and K. R. Khair, “Heat and mass transfer by natural convection in a porous medium,” International Journal of Heat and Mass Transfer, vol. 28, no. 5, pp. 909–918, 1985.
11. A. Bahloul, N. Boutana, and P. Vasseur, “Double-diffusive and soret-induced convection in a shallow horizontal porous layer,” The Journal of Fluid Mechanics, vol. 491, pp. 325–352, 2003.
12. F. Y. Zhao, D. Liu, and G. F. Tang, “Natural convection in a porous enclosure with a partial heating and salting element,” International Journal of Thermal Sciences, vol. 47, no. 5, pp. 569–583, 2008.
13. D. Liu, F. Y. Zhao, and G. F. Tang, “Thermosolutal convection in saturated porous enclosure with concentrated energy and solute sources,” Energy Conversion and Management, vol. 49, no. 1, pp. 16–31, 2008.
14. F. Y. Zhao, D. Liu, and G. F. Tang, “Free convection from one thermal and solute source in a confined porous medium,” Transport in Porous Media, vol. 70, no. 3, pp. 407–425, 2007.
15. M. Mamou, P. Vasseur, and E. Bilgen, “Multiple solutions for double-diffusive convection in a vertical porous enclosure,” International Journal of Heat and Mass Transfer, vol. 38, no. 10, pp. 1787–1798, 1995.
16. M. Bourich, A. Amahmid, and M. Hasnaoui, “Double diffusive convection in a porous enclosure submitted to cross gradients of temperature and concentration,” Energy Conversion and Management, vol. 45, no. 11-12, pp. 1655–1670, 2004.
17. A. Mansour, A. Amahmid, M. Hasnaoui, and M. Bourich, “Multiplicity of solutions induced by thermosolutal convection in a square porous cavity heated from below and submitted to horizontal concentration gradient in the presence of soret effect,” Numerical Heat Transfer A, vol. 49, no. 1, pp. 69–94, 2006.
18. F. Y. Zhao, D. Liu, and G. F. Tang, “Natural convection in an enclosure with localized heating and salting from below,” International Journal of Heat and Mass Transfer, vol. 51, no. 11-12, pp. 2889–2904, 2008.
19. Z. Alloui, L. Dufau, H. Beji, and P. Vasseur, “Multiple steady states in a porous enclosure partially heated and fully salted from below,” International Journal of Thermal Sciences, vol. 48, no. 3, pp. 521–534, 2009.
20. O. V. Trevisan and A. Bejan, “Natural convection with combined heat and mass transfer buoyancy effects in a porous medium,” International Journal of Heat and Mass Transfer, vol. 28, no. 8, pp. 1597–1611, 1985.
21. M. Mamou, P. Vasseur, and E. Bilgen, “A Galerkin finite-element study of the onset of double-diffusive convection in an inclined porous enclosure,” International Journal of Heat and Mass Transfer, vol. 41, no. 11, pp. 1513–1529, 1998.
22. M. Bourich, M. Hasnaoui, and A. Amahmid, “Double-diffusive natural convection in a porous enclosure partially heated from below and differentially salted,” International Journal of Heat and Fluid Flow, vol. 25, no. 6, pp. 1034–1046, 2004.
23. N. H. Saeid and I. Pop, “Natural convection from a discrete heater in a square cavity filled with a porous medium,” Journal of Porous Media, vol. 8, no. 1, pp. 55–63, 2005.
24. H. Saleh and I. Hashim, “Hydromagnetic natural convection flow in a fluid-saturated porous medium with nonuniform heat generation,” Advances in Mechanical Engineering, vol. 2012, Article ID 267170, 13 pages, 2012.
25. S. Akbal and F. Baytaş, “Effects of non-uniform porosity on double diffusive natural convection in a porous cavity with partially permeable wall,” International Journal of Thermal Sciences, vol. 47, no. 7, pp. 875–885, 2008.
26. B. Goyeau, J. P. Songbe, and D. Gobin, “Numerical study of double-diffusive natural convection in a porous cavity using the Darcy-Brinkman formulation,” International Journal of Heat and Mass Transfer, vol. 39, no. 7, pp. 1363–1378, 1996.
27. O. V. Trevisan and A. Bejan, “Mass and heat transfer by natural convection in a vertical slot filled with porous medium,” International Journal of Heat and Mass Transfer, vol. 28, no. 3, pp. 403–415, 1986.
28. A. Amahmid, M. Hasnaoui, M. Mamou, and P. Vasseur, “Boundary layer flows in a vertical porous enclosure induced by opposing buoyancy forces,” International Journal of Heat and Mass Transfer, vol. 42, no. 19, pp. 3599–3608, 1999.
29. B. S. Bhadauria and A. K. Srivastava, “Magneto-double diffusive convection in an electrically conducting-fluid-saturated porous medium with temperature modulation of the boundaries,” International Journal of Heat and Mass Transfer, vol. 53, no. 11-12, pp. 2530–2538, 2010.
30. 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.
31. S. Alchaar, P. Vasseur, and E. Bilgen, “Effect of a magnetic field on natural convection in a shallow cavity heated from below,” Chemical Engineering Communications, vol. 134, pp. 195–209, 1995.
32. A. J. Chamkha and H. Al-Naser, “Hydromagnetic double-diffusive convection in a rectangular enclosure with uniform side heart and mass fluxes and opposing temperature and concentration gradients,” International Journal of Thermal Sciences, vol. 41, no. 10, pp. 936–948, 2002.
33. M. C. Ece and E. Buyuk, “Natural-convection flow under a magnetic field in an inclined rectangular enclosure heated and cooled on adjacent walls,” Fluid Dynamics, vol. 38, Article ID 564590, 2006.
34. M. A. Teamah, A. F. Elsafty, and E. Z. Massoud, “Numerical simulation of doublediffusive natural convective flow in an inclined rectangular enclosure in the presence of magnetic field and heat source,” International Journal of Thermal Sciences, vol. 52, pp. 161–175, 2012.
35. 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.
36. S. V. Patankar, Numerical Heat Transfer and Fluid Flow, Hemisphere, Washington, DC, USA, 1st edition, 2007.
37. H. K. Versteeg and W. Malalasekera, An Introduction to Computational Fluid Dynamics, Pearson Education, Essex, UK, 2nd edition, 2007.
38. A. C. Baytas and I. Pop, “Free convection in oblique enclosures filled with a porous medium,” International Journal of Heat and Mass Transfer, vol. 42, no. 6, pp. 1047–1057, 1999.
39. N. H. Saeid and A. A. Mohamad, “Natural convection in a porous cavity with spatial sidewall temperature variation,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 15, no. 6, pp. 555–566, 2005.
40. K. L. Walker and G. M. Homsy, “Convection in a porous cavity,” The Journal of Fluid Mechanics, vol. 87, no. 3, pp. 449–474, 1978.