Journal of Applied Mathematics

VolumeΒ 2012Β (2012), Article IDΒ 594701, 29 pages

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

## A Finite Element Study of Double Diffusive Mixed Convection in a Concentration Stratified Darcian Fluid Saturated Porous Enclosure under Injection/Suction Effect

^{1}Department of Mathematics and Statistics, Indian Institute of Technology Kanpur, Kanpur 208016, India^{2}Department of Applied Mathematics, Defence Institute of Advanced Technology, Pune 411025, India

Received 9 April 2012; Accepted 23 April 2012

Academic Editor: MehmetΒ Sezer

Copyright Β© 2012 B. V. Rathish Kumar and S. V. S. S. N. V. G. Krishna Murthy. 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

Numerical investigation of mixed convection flow in a concentration-stratified fluid-saturated vertical square porous enclosure is investigated by Galerkin finite element method. The forced flow conditions are imposed by providing an inlet at the bottom wall and an outlet with a suction on the top wall. The free convection is induced by introducing a hot but isothermal temperature on the left vertical wall together with Boussinesq approximation on density. Numerical results are presented by tracing the cumulative heat and mass fluxes, streamline and isotherms of the fluid for a wide range of governing parameters such as suction/injection velocity ββ, suction/injection width (), mass stratification parameters (), Rayleigh number βRaβ, buoyancy ratio ββ, and Lewis number βLeβ.

#### 1. Introduction

A wide range of scientific and engineering applications such as packed bed reactors, porous insulation, beds of fossil fuels, nuclear waste disposal, usage of porous conical bearings in lubrication technology, fibrous insulation systems, grain storage, chemical catalytic reactors, solar collectors, heat exchangers, food processing, energy efficient drying process, geophysics and energy related problems, requires a detailed understanding of convective heat transfer process in a fluid-saturated porous media. These applications are reviewed in recent books by [1β6], provided a detailed coverage of research works related to convective heat transfer in Porous media. In particular applications such as cooling of radioactive waste containers, solar power collectors, electronic equipment cooling, thermal insulation of buildings, have drawn considerable attention in last several years onto flow through porous enclosures. Some of the industrial technical processes like solar central receivers exposed to wind currents, electronic devices cooled by fans, heat exchangers, microelectronic devices packed with micro-pin-fin, vented enclosures filled with microspheres require a detailed study of mixed convection process in porous enclosures under different conditions.

Mixed convection is a fundamentally significant heat transfer mechanism that occurs in a selection of industrial and technological applications. For example, it is of great interest to thermal designers in the field of electronic cooling. Given that its occurrence is frequent, there is a necessity to understand the physics of this fundamental mechanism at board level. It is motivated both by a desire to gain an understanding the fundamental physics of mixed convection and to solve design problems in scientific and electronic applications where mixed convection often occurs. Apart from this situation there are many more industrial and natural problems like microwave heating of dielectric materials using a rectangular wave guide, sintered porous heat sink for cooling of high-powered microprocessors for its applications, biomass growth on the hydraulic properties of saturated porous media, that is, bioclogging and so forth, where a thorough understanding of mixed convection process is essential. Experimental studies need to be supplemented by theoretical support. This certainly calls for a complete theoretical investigation of mixed convection process in porous media.

A comprehensive review of several analytical, numerical, and experimental results has been published on mixed convection process in porous enclosures that are reported in the literature. Reference [7] has numerically studied mixed convection in a vented enclosure with continuum fluid and isothermal vertical surface. Bansod and Jadhav [8] presented the mixed convection heat and mass transfer near a vertical surface in a stratified porous medium using an integral method. Chamkha [9, 10] in his mathematical model had additionally included transverse magnetic effects too. W. J. Chang and W. L. Chang [11] looked at mixed convection in a vertical tube partially filled with porous medium. Chen [12], Chou and Chung [13], and Prasad et al. [14] have analytically analyzed mixed convection with the free convection induced by either a horizontal or a vertical hot plate enclosed in a semi-infinite porous medium under boundary layer assumptions. C. K. Chen and C. H. Chen [15] have analyzed the non-Darcian and nonuniform porosity effects on conjugate mixed convection heat transfer from a plate fin in porous media under boundary layer assumptions, employing one-dimensional heat conduction equation for closing the conjugate heat transfer model. El-Amin et al. [16] numerically investigated the interaction of mixed convection with thermal radiation of an optical dense viscous fluid adjacent to an isothermal cone imbedded in a porous medium with Rossel and diffusion approximation incorporating the variation of permeability and thermal conductivity. Elbashbeshy [17] solved the problem of mixed convection along a vertical plate embedded in a non-Darcian porous medium with suction and injection under boundary layer assumptions using similar variables defined based on scale analysis.

Hassanien and Al-arabi [18] have simplified the boundary layer equations using scale analysis and numerically solved the problem of non-Darcian unsteady mixed convection flow near the stagnation point on heated vertical surface embedded in a porous medium with thermal radiation, and variable viscosity has been considered. Jayanthi and Kumari [19] have analytically studied the effect of variable viscosity on non-Darcian free or mixed convection flow on a vertical surface in a non-Newtonian fluid saturated porous medium. Khanafer and Vafai [20] investigated the double diffusive mixed convection in a square enclosure filled with a porous medium in the presence of moving boundary which was formulated and solved numerically by finite-volume approach along with the alternating direction implicit scheme. Recently, Wong and Saied [21] have studied numerically the non-Darcian effects, based on Brinkman-Forchheimer-extended Darcy model, on jet impingement cooling in a horizontal porous layer in a mixed convection regime. They have employed finite volume method in conjunction with power-law scheme for convection-diffusion model. Lai [22] explained in detail about the mixed convection in saturated porous media. Kumar et al. [23] have studied mixed convection flow in a vertical square enclosure filled with a non-Darcian fluid saturated homogeneous porous medium by providing an inlet at the bottom wall and an outlet with a suction on the top wall numerically by using the Galerkin finite element method. Nazar et al. [24] analyzed mixed convection boundary layer flow past a horizontal circular cylinder in a porous medium based on Brinkman model.

Singh et al. [25] reported on mixed convection process in a flow past a porous vertical plate bounded by a porous medium in a rotating system in the presence of magnetic filed. Shohel Mahmud and Pop [26] have numerically extended the study of [7] to Darcian porous square vented enclosure using finite volume approach. Srinivasacharya and RamReddy [27] have analysed the effects of thermal and solutal stratification on mixed convection along a vertical plate embedded in a micropolar fluid saturated non-Darcy porous medium. kurdyumov and Linan [28] presented an analysis for the steady, two-dimensional, free convection around line sources of heat and heated cylinders in unbounded saturated porous media. Vafai and Hadim [29] have presented an overview on the computational studies of heat transfer in porous media both under natural and mixed convection processes. Waheed et al. [30] have studied numerically the buoyancy- and the shear-driven flow induced by a hot plate moving through the horizontal midplane of a rectangular enclosure filled with fluid-saturated porous medium.

So far not much work has been reported on Darcian mixed convection in vertical square porous enclosure. In particular mixed convection in a vertical porous enclosure under simultaneous suction/injection effects on opposite walls, which becomes very relevant in the context of electronic devices, environmental chamber for bacterial culture preservation, and so forth, has not been considered so far. In the present study, attention will be focused on a new problem of steady mixed convection heat and mass transfer with concentration stratification inside a vertical square fluid saturated porous enclosure with fluid injection at the bottom wall and a suction at the top wall. Thus forced convection is imposed by this combination of suction/injection flow conditions. The free convection is induced by the hot and isothermal left vertical wall. Detailed numerical simulations are carried out by Galerkin finite element method for a wide range of parameters such as Rayleigh number ββ, suction/injection flow speed ββ, suction/injction width (), mass stratification parameters (), buoyancy ratio ββ, Lewis number ββ. Flow and temperature distribution is analyzed by tracing streamlines, isotherms, local/cumulative heat fluxes.

#### 2. Mathematical Formulation

We consider a two-dimensional square enclosure filled with a fluid saturated porous medium with the left vertical wall at the uniform temperature and the surface concentration as which are assumed to be higher than the ambient temperature and concentration, respectively, and the right vertical wall is kept at uniform ambient temperature and with the ambient species concentration as is defined as and shown in Figure 1: other two walls are adiabatic, excluding the inlet/outlet portion at the bottom/top wall. The other two walls are excluding the suction/injection windows, which are permeable with zero normal thermal gradients. Here the forced flow conditions are brought in by the combination of suction and injection effects on bottom and top walls, respectively. The inflow is considered to be at ambient temperature. is the width of the inlet and outlet and is the height of the cavity. The porous medium is assumed to be isotropic with permeability , and Boussinesq approximation is valid in the momentum equation for the density term. Dispersion effect has been neglected. and are the horizontal and vertical components, respectively. The physical configuration and the coordinate system used in the current study are shown in Figure 2. The equations governing the steady-state double diffusive Darcy mixed convection flow along a heated vertical wall under concentration stratification with the above mentioned assumptions, nondimensional equations [31β35], governing the flow, heat, and mass transport in a fluid saturated porous media, can be written in terms of the following nondimensional governing equations: with the following boundary conditions:

The nondimensional variables are defined as follows:

The nondimensional stream function is defined as The local and cumulative heat flux () along the left vertical wall is computed by the relation for Darcy model: Similarly, the local and cumulative mass flux () along the vertical wall is given by the relation for Darcy model: where is the running coordinate inββββdirection and the upper limits of the integration, that is, in the above two expressions, give the global heat flux or Nusselt number () and global mass flux or Sherwood number (), respectively.

#### 3. Solution Methodology: Finite Element Formulation

Governing equations (2.2) together with boundary conditions (2.3) has been solved using Galerkin finite element method. Following is the outline of finite element formulation used in the current study. Let denote the domain of interest and let be the boundary of the domain. The discretized representation of is given by , where denotes a typical bilinear element of the discretized domain and is the total number of such elements. The discretized elements are fully disjoint, that is, . The discretized representation of the field variables on a typical bilinear element is where denotes the standard bilinear interpolation function on a typical element .

Consider the Galerkin weighted residual forms of governing equation (2.2) on : Here, is the weight function that is chosen as equal to the interpolation function. Rewriting (3.2) in the weak form and on introducing the element-level discretized representation for the field variables, that is, (3.1), one would arrive at the following element level matrix equation: where Here, Here, denotes the normal to the boundary of the porous cavity. The nonlinear global system obtained by assembling the local element matrix (3.3) is solved iteratively using frontal method for nonlinear systems to an accuracy of on the relative error of nodal field variables from successive iterations, that is, , where or or . Here the superscript refers to the iteration level and refers to the nodal point index. To accelerate the convergence, we underrelax the results from successive iteration by a factor of 0.25.

#### 4. Results and Discussion

The parameters which influence the double diffusive natural convection process in the present study areβbuoyancy ratio ββ, Lewis number ββ, mass stratification parameters (), injection/suction () velocity parameter ββ, window width (), and Rayleigh number ββ. Grid selection tests have been carried out on five different grid systems consisting of , , , , and elements. On these five grid systems simulations have been carried out for various combinations of above-mentioned parameters, and obtained stream function, temperature, and concentration values are compared. Few of these comparison plots are presented in Figure 3. In Figures 3(a)β3(c) ββ, , values are plotted alongββ-axis at fixedββββforββ,ββ, () = = 0.1,ββ,ββ, and the corresponding plots along -axis at fixedββββare in Figures 3(d)β3(f), respectively. It is observed from these plots and the corresponding data (Table 1) that as one moves from to higher grid systems, there is only a marginal change in the solutions. Hence grid system has been chosen for extensive numerical simulations. Details of the study showing the influence of various parameters on heat and mass transfer process in the presence of concentration stratification are given in the following subsections.

##### 4.1. Influence of Parameter ββ

The influence of velocity parameter ββ on the heat transfer process in the presence of concentration stratification has been analyzed by comparing the cumulative heat fluxes (CHFLXs) along the left vertical wall of the porous enclosure. In Figure 4 CHFLXs along the left vertical wall are presented for , , ,ββ, ,ββ. At all values of ββ, CHFLXs are seen to increase as one moves along left vertical wall, that is, . This increase is relatively sharp along the leading half portion of the left vertical wall, that is, . In this region CHFLXs are seen to increase with increasing values of ββ. However along , CHFLXs are seen to decrease with increasing ββ. Consequently the global heat flux (GHFLX) or Nusselt number () is seen to decrease with increasing values of ββ. To explain these observations one has to investigate the temperature and flow fields. In Figures 5(a)β5(d) and Figures 6(a)β6(d) isotherms and streamlines corresponding to the above set of parameters are traced.

From the isotherm plots in Figures 5(a)β5(d) one can notice the manifestation of sharp thermal boundary layer (TBL) along the left vertical wall. While, for the caseββ, the TBLs continue to remain sharp all along the left wall, those corresponding to tend to get increasing blunt for as ββ increases. From the isotherm values one can clearly notice that the introduction of forced convection has led to greater transfer into the core of the domain. Clearly it indicates a raise in average temperature in the porous enclosure with increasing values of velocity. With increasing ββ one can also notice an increase in the transverse temperature gradients in the core of the domain.

The streamline plots in Figures 6(a)β6(d) depict the manifestation of a dominant primary circulation zone in the core of the domain and secondary circulation patterns centered on the windows. On increasing the velocity while the primary circulation zone undergoes a dominant longitudinal stretch, which nearly touches both the vertical walls, the secondary zones are seen to increasingly swell and dominantly impinge onto the primary circulation zone. Clearly the longitudinal stretch in primary circulation zone leads to better mixing of hot and cold fluids thereby leading to better heat transfer from the hot left vertical wall into the porous enclosure. Clearly the increase in the forced convection component leads to better heating in the porous enclosure and an overall raise in temperature of the porous enclosure. Hence the observed fall in the global cumulative heat flux (GCHLX) with increasing ββ turns out to be natural consequence.

In Figure 7 the cumulative mass fluxes (CMFLXs) along the left vertical wall are presented. At all values ofββββ similar to CHFLXs, the CMFLXs increase along the left vertical wall. With increasing ββ along the leading part of the left vertical wall a raise in CMFLXs are noticed. However, owing to the presence of concentration stratification, the extent of the stretch along the left vertical wall wherein such a raise is noticed decreases from to .

Again owing to the presence of concentration stratification, the variation in CMFLXs along the upper portion of the left vertical wall (i.e., ) gets complex. An initial raise in ββ, 0.0 to 0.1 leads to a fall in CMFLXs, corresponding to , and thereby a fall in global cumulative mass flux (GCMFLX) or Sherwood number (). A further increase in ββ from 0.1 to 0.25 leads to shoot-up in CMFLXs and GCMFLX and on further increasing ββ to 0.5 a significant fall in CMFLXs and GCMFLX is noticed. To understand such a trend in mass fluxes the concentration distribution in the whole domain for different values ofβ ββ is traced and presented in Figures 8(a)β8(d).

The isoconcentration contour pattern depicted by plots clearly indicates the additional dimension to the complexity of interaction between various forces responsible for concentration transport in the porous enclosure. One can notice the presence of concentration boundary layers, whose characteristics are sensitive to velocity, along the left vertical wall. The plots clearly depict an overall increase in the concentration transport with increasing values of ββ that is, velocity. With the increase in ββ circular packets of concentric concentration are seen to manifest on the windows.

##### 4.2. Influence of Window Width ()

The influence of size of window on the flow, temperature, and concentration fields is analyzed and presented in Figures 9(a) and 9(b) and Figures 10(a)β10(i). In Figures 9(a) and 9(b) CHFLX and CMFLX along the left vertical wall are presented for , , , , ,ββ. At all values of both CHFLX and CMFLX increase as one gradually moves along the vertical wall. All along the left vertical wall both local heat and mass fluxes increase with increasing values of . Both Nusselt number and Sherwood number are seen to increase with increasing values of . In Figures 10(a)β10(i) streamlines, isotherms, isoconcentration contours corresponding to = 0.1, 0.2, 0.3, , , ββand ,ββ. From the streamlines in Figures 10(a)β10(c) one can notice the manifestation of a prominent primary circulation zone in the core of the domain and two secondary circulations centered on the window slits at all values of . Increasing the window slit size leads to overall enlargement of both the primary and secondary circulation zones. From the isotherm plots in Figures 10(d)β10(f) one can notice the presence of TBLs along the left vertical walls. These TBLs get sharper with increasing values of . The change in the isotherm pattern and magnitudes of isotherms suggests an overall enhancement in heat transfer in the enclosure with increasing . From Figures 10(g)β10(i) similar variation is noticed in concentration fields, with the increase in . Concentration boundary layers (CBLs) noticed along the left vertical wall get sharper with increasing values of .

##### 4.3. Influence of Buoyancy Ratio ββ

To analyze the influence of buoyancy ratio ββ, which denotes the relative strengths of two buoyancy forces, namely, mass and thermal buoyancy forces, simulations are carried out for a wide range of ββ, covering . Here it can be noted that the algebraic sign of ββ denotes the effective direction of these buoyancy forces. Thermal buoyancy forces act vertically upward, and the species buoyancy may act in either direction depending on the molecular weight relative to the fluid. A heavier species contributes to a buoyant force that acts vertically downward, thereby opposing the vertically upward thermal buoyancy, and is negative in sign. Ifββ, the flow is driven by buoyancy alone. Simulations have been carried out with different degrees of mass stratifications. In Figure 11(a), CHFLXs along the left vertical wall are plotted for different values of ββ when , , , ,ββ. From Figure 11(a) one can see that CHFLXs increase with increasing values of ββ. All along the left vertical wall local heat fluxes (LHFLXs) are in a raise as ββ is raised. At all values of ββ a sharp raise in CHFLXs is noticed along the leading edge of the left vertical wall, that is, . To understand this behavior, flow and temperature fields are traced and presented in Figure 12.

The streamlines in Figures 12(a)β12(g) clearly depict that the flow field is highly sensitive both to the sign and magnitude of ββ. When the downward oriented species buoyancy forces are dominantly opposing the upward oriented thermal buoyancy forces (i.e., when ), the flow is primarily covered by the elongated semicircular pattern, emanating from the window located on the top horizontal wall, between two vertical flow fields which are adjacent to the vertical walls of the enclosure. As the ββ is raised to β1, where in the downward facing species buoyancy forces are opposing but are equal to the upward facing thermal buoyancy forces, the semicircular flow fields manifest on windows located on bottom horizontal wall and thereby reduce the dominance of downward protruding semicircular pattern that emanates from top window slit. A further increase in ββ gradually leads to the dominance of the thermal buoyancy forces. As thermal buoyancy forces get dominant, a primary circular zone manifests between the two semicircular patterns centered on the window slits. With a further raise in ββ to 0.5, where in both the buoyancy forces are upward oriented, the primary circulation zone prominently falls major part of the enclosure, except the region close to the window slits. As the ββ is raised to 1.5 one can notice that semicircular flow patterns again tend to get severe and the primary circulation zone begins to shrink in size. With a further increase in ββ, the primary circulation zone completely vanishes allowing the two semicircular patterns to elongate and cover whole of the enclosure.

The corresponding isotherms and isoconcentration contours are presented in Figures 13(a)β13(g) and Figures 14(a)β14(g). The majority of isotherms which are vertically oriented whenββββgradually undergo a clockwise diagonal twist as ββ is raised to β1.0. When ββ is raised to β0.5, the diagonally oriented isotherms which emanate from the lower left corner begin to turn towards the horizontal walls, especially while they are being traced in the core of the domain. On increasing ββ in steps to 1.5 while the isotherms close to the vertical wall tend to get vertical, they take a horizontal path with increasingly sharp bends near the vertical walls. The raise in the sharpness of TBLS clearly supports the observed increase in with increasing values of ββ.

Features observed in the isoconcentration contours are nearly similar to those observed in isotherm plots especially when . The deviation in the iso-concentration contours from the corresponding isotherms when can be attributed to the presence of species stratification in the porous enclosure. Like in isotherm plots, here too one can notice the manifestation of CBLs along the vertical walls. The increase in Sherwood numbers with increasing ββ (Figure 11(b)) is clearly supported by the increasingly sharpening CBLs with increasing ββ. Both the isotherm and isoconcentration patterns depict an enhanced heat and species transportation in the enclosure with increasing values of ββ.

##### 4.4. Influence of Lewis Number ββ

In Figure 15(a) influence of on CHFLX along the left vertical wall is presented. Along the leading edge of the left vertical wall (i.e., for ) CHFLXs remain insensitive to the differences in thermal and mass diffusivities. However along one can notice a fall in CHFLX all along the wall with increasing values of . Consequently GCHFLX or decreases with increasing ββ. It is also to be noticed that at all values of , there is a gradual increase in CHFLX as one traverses along the left vertical wall (i.e., ). When the mass diffusivity is larger than thermal diffusivity, the corresponding heat fluxes and ββare larger. This indicates the possibility of greater mixing of hot and cold fluids with increasing ββ and thereby leading to an overall raise in the temperature in the porous enclosure. To further investigate such a happening, flow and temperature fields are traced in the form of streamlines and isotherms in Figures 16 and 17.

From the streamline plots in Figures 16(a)β16(d) firstly one can notice the manifestation of prominent primary circulation centered in the core of the domain, also the presence of secondary semi-circular zones centered on the window slits on the horizontal walls. The increase in the diagonal stretch in the primary zone and the swelling of the secondary zones leading to the corner elongation/stretching of the primary circulation zones clearly depict the increase in the mixing of hot and cold fluids with increase in ββ. Also from the isotherm plots in Figures 17(a)β17(d) one can see the presence of TBLS along the left vertical wall. The effect of increasing thermal diffusivities is clear from the increase in the slope of isotherms, especially in the core of the domain. The semi-circular isotherm pattern seen in the leading portion of the lower horizontal wall is due to injection of fluid at ambient temperature through the window. In Figure 15(b) CMFLXs, corresponding to above set of parameters, along the left vertical wall are presented for different values of ββ. All along the wall CMFLXs increase with increasing values of ββ. Consequently GCMFLX or Sherwood number () increases with increasing ββ. It is also to be noted that while the magnitudes of Sh are lower than , their relative raise with increasing ββ is considerably larger than the relative fall in . From Figure 15(b) one can also notice that there is sharp raise in CMFLXs along the leading edge of the left vertical wall, that is, , and thereby CMFLXs are seen to increase gradually. To further investigate the concentration contours are traced in Figures 18(a)β18(d). From the iso-concentration plots in Figures 18(a)β18(d) one can notice the manifestation of sharp concentration boundary layer (CBL) along the left vertical wall. These CBLs tend to get blunt as one moves beyondββ. From the contour magnitudes one can also notice overall increase in the concentration level in the core of the enclosure with increasing ββ. Clearly increase in thermal diffusivities favors the concentration transport.

##### 4.5. Influence of Rayleigh Number ββ

In Figure 19(a) influence of ββ on CHFLX along the left vertical wall is presented for , ,ββ, , ,ββ. Clearly CHFLX increases as one traverses the left vertical wall starting from the lower left corner. Also all along the left vertical wall LHFLXs are seen to increase with increasing values of ββ. GHFLX or increases with ββ. In Figure 19(b) CMFLXs along the left vertical wall are presented for the above set of parameters. Clearly CMFLX increases as one marches along the left vertical wall, with a relatively sharp raise when . At any point along the left vertical wall LHFLX is seen to increase with increasing ββ. Owing to the presence of mass stratification in the porous enclosure, while the GMFLX or Sherwood number undergoes a considerable enhancement for , it changes marginally when . To further analyze the situation the flow, temperature, and concentration fields are analyzed as streamlines, isotherms, and iso-concentration contours in Figures 20, 21 and 22.

In Figures 20(a)β20(f) streamlines corresponding to ,ββ, , ,ββ, are provided. In the range there is a considerable variation in streamline pattern as ββ is increasing. At all values of ββ one can notice the presence of both the primary circulation zone centered in the middle of the enclosure and two secondary circulation zones centered on the window slits. The primary circulation zone which is vertically oriented and is squeezed between the vertical flows due to injection and suction effects gradually undergoes diagonal bulging as ββ is increased. With the increase in ββ, the breath of the two vertical flow patterns adjacent to the vertical wall tends to decrease. For relatively a small variation is noticed in streamline pattern. The enlargement in the primary circulation for suggests a better mixing of hot and cold fluid suggesting a possibility of sustained increase in the sharpness of TBLs in corresponding isotherms.

Now in Figures 21(a)β21(f) isotherms corresponding to the above set of parameters are provided. From the isotherm pattern and the magnitude one can notice the presence of TBLs with increasing sharpness as ββ is gradually increased from 10 to 150. On the isotherm magnitudes one can clearly notice an overall enhancement in the temperature distribution in the enclosure. In Figures 22(a)β22(f) one can notice the presence of CBLs along the left vertical wall. While the CBLs corresponding get sharpened, those corresponding to undergo only a little variation especially when , that is, towards at the end of the left vertical wall. This can be attributed to the presence of concentration stratification in the enclosure. The overall concentration distribution in the enclosure is on the raise with increasing ββ.

#### 5. Conclusions

Numerical computations have been carried out on Darcian double diffusive mixed convection process in a in a fluid saturated porous enclosure under effect with mass stratification by finite element method. Flow, temperature, and concentration fields are analyzed through streamline, isotherm, iso-concentration tracing and local and cumulative heat flux, cumulative concentration heat flux calculations. Following are the conclusions based on the results obtained from the present study.(i)While decreases with increasing velocity, is found to initially decrease but later increase. Prominent primary circulation pattern with secondary circulation zones that stretch on increasing velocity manifests in flow field. TBLs and CBLs, thermal and mass gradients highly sensitive to the velocity are noticed in the temperature and concentration fields.(ii)Both and increase with increasing size of window (). Increasing () enlarges both primary and secondary circulation zones and sharpens the TBLs and CBLs.(iii)Flow pattern undergoes significant pattern variation as the opposing thermal and mass buoyancy forces turn to favor each other. Semicircular flow pattern is transformed to a multicellular flow with a prominent primary circulation zone. Both and increase with increasing values of ββ.(iv)While decreases with increasing , it increases with increasing ββ. is found to increase both with increasing and . Flow, temperature, and concentration fields are sensitive both to ββ and ββ.(v)Darcy case of an increase in the levels of concentration stratification leads to a fall in and values.

#### Nomenclature

: | Suction/injection velocity |

: | Buoyancy ratio |

: | Dimensionless species concentration |

: | Nondimensional species concentration |

: | Mass diffusivity |

: | Width of the inlet |

: | Dimensionless width of the inlet |

: | Typical element in finite element formulation |

: | Gravitational acceleration |

: | Height of the cavity |

: | Thermal conductivity |

: | Permeability of the porous medium |

: | Outward unit normal to the surface |

: | Quadratic interpolation function |

: | Nusselt number |

: | Rayleigh number () |

: | Dimensional mass stratification parameter |

: | Nondimensional mass stratification parameter |

: | Nondimensional temperature |

: | Temperature of the flow through at the inlet |

: | Temperature of the left vertical wall |

: | Dimensional velocity components in and directions |

: | Nondimensional velocity components in and directions |

: | Convective velocity |

: | Absolute value of the velocity of the forced flow at the inlet |

: | Galerkin weight function used in the finite element formulation |

: | Dimensional cartesian coordinates |

: | Nondimensional cartesian coordinates |

: | Coordinates of the suction/injection starting point on the bottom wall |

: | Coordinates of the suction/injection starting point on the top wall. |

#### Greek Symbols

: | Thermal diffusivity |

: | Thermal expansion coefficientββ () |

: | Concentration expansion coefficientββ() |

: | Boundary of the domain |

: | Nondimensional stream function |

: | Kinematic viscosity of fluid |

: | Domain considered in the problem |

: | Fluid densityββ |

: | Dimensionless temperature. |

#### Subscripts

: | Ambient points |

: | Pressure |

: | Evaluated at wall temperature |

: | Evaluated at point . |

#### References

- D. A. Nield and A. Bejan,
*Convection in Porous Media*, Springer Verlag, New York, NY, USA, 1999. - D. Ingham and I. Pop,
*Transport Phenomena in Porous Media*, Elsevier, Oxford, UK, 2005. - K. Vafai,
*Hand Book of Porous Media*, Taylor & Francis, New York, NY, USA, 2nd edition, 2005. - I. Pop and D. Ingham,
*Convective Heat Tranfer: Mathematical and Computational Modelling of Viscous Fluid and Porous Media*, Pergamon, Oxford, UK, 2001. - A. Bejan and A. D. Krus,
*Heat Transfer Hand Book*, Wiley, New York, NY, USA, 2003. - D. Ingham, A. Bejan, and A. Mamut,
*Emerging Technologies and Techniques in Porous Media*, Kluwer, Dordrecht, The Netherlands, 2004. - D. Angirasa, βMixed convection in a vented enclosure with an isothermal vertical surface,β
*Fluid Dynamics Research*, vol. 26, no. 4, pp. 219β233, 2000. View at Publisher Β· View at Google Scholar Β· View at Scopus - V. J. Bansod and R. K. Jadhav, βEffect of double stratification on mixed convection
heat and mass transfer from a vertical surface in a fluid-saturated porous medium,β
*Heat TransferAsian Research*, vol. 39, no. 6, pp. 378β395, 2010. View at Google Scholar - A. J. Chamkha, βNon-Darcy fully developed mixed convection in a porous medium
chanel with heat generation/absorption and hydro magnetic effects,β
*Numerical Heat Transfer A*, vol. 32, no. 6, pp. 653β675, 1997. View at Google Scholar - A. J. Chamkha, βMixed convection flow along a vertical permeable plate embedded in a porous medium in the presence of a transverse magnetic field,β
*Numerical Heat Transfer A*, vol. 34, no. 1, pp. 93β103, 1998. View at Publisher Β· View at Google Scholar Β· View at Scopus - W. J. Chang and W. L. Chang, βMixed convection in a vertical tube partially filled with porous medium,β
*Numerical Heat Transfer A*, vol. 28, no. 6, pp. 739β754, 1995. View at Publisher Β· View at Google Scholar Β· View at Scopus - C. H. Chen, βTechnical note: non-Darcy mixed convection from a horizontal surface with variable surface heat flux in a porous medium,β
*Numerical Heat Transfer A*, vol. 30, no. 8, pp. 859β869, 1996. View at Publisher Β· View at Google Scholar Β· View at Scopus - F. C. Chou and P. Y. Chung, βEffect of stagnant conductivity on non-darcian mixed convection in horizontal square packed channels,β
*Numerical Heat Transfer A*, vol. 27, no. 2, pp. 195β209, 1995. View at Publisher Β· View at Google Scholar Β· View at Scopus - V. Prasad, F. C. Lai, and F. A. Kulacki, βMixed convection in horizontal porous layers heated from below,β
*Journal of Heat Transfer*, vol. 110, no. 2, pp. 395β402, 1988. View at Publisher Β· View at Google Scholar Β· View at Scopus - C.-K. Chen and C.-H. Chen, βNonuniform porosity and non-Darcian effects on conjugate mixed convection heat transfer from a plate fin in porous media,β
*International Journal of Heat and Fluid Flow*, vol. 11, no. 1, pp. 65β71, 1990. View at Publisher Β· View at Google Scholar Β· View at Scopus - M. F. El-Amin, N. A. Ebrahiem, A. Salama, and S. Sun, βRadiative mixed convection over an isothermal cone embedded in a porous medium with variable permeability,β
*Journal of Applied Mathematics*, vol. 2011, Article ID 124590, 2011. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at Scopus - E. M. A. Elbashbeshy, βThe mixed convection along a vertical plate embedded in non-Darcian porous medium with suction and injection,β
*Applied Mathematics and Computation*, vol. 136, no. 1, pp. 139β149, 2003. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at Scopus - I. A. Hassanien and T. H. Al-arabi, βNon-Darcy unsteady mixed convection flow near the stagnation point on a heated vertical surface embedded in a porous medium with thermal radiation and variable viscosity,β
*Communications in Nonlinear Science and Numerical Simulation*, vol. 14, no. 4, pp. 1366β1376, 2009. View at Publisher Β· View at Google Scholar Β· View at Scopus - S. Jayanthi and M. Kumari, βEffect of variable viscosity on non-Darcy free or mixed convection flow on a vertical surface in a non-Newtonian fluid saturated porous medium,β
*Applied Mathematics and Computation*, vol. 186, no. 2, pp. 1643β1659, 2007. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at Scopus - K. Khanafer and K. Vafai, βDouble-diffusive mixed convection in a lid-driven enclosure filled with a fluid-saturated porous medium,β
*Numerical Heat Transfer A*, vol. 42, no. 5, pp. 465β486, 2002. View at Publisher Β· View at Google Scholar Β· View at Scopus - K. C. Wong and N. H. Saeid, βNumerical study of non-Darcian effects on jet impingement cooling in a horizontal porous layer in the mixed convection regime,β
*International Communications in Heat and Mass Transfer*, vol. 36, no. 1, pp. 45β50, 2009. View at Publisher Β· View at Google Scholar Β· View at Scopus - F. C. Lai, βMixed convection in saturated porous media,β in
*Hand Book of Porous Media*, K. Vafai, Ed., pp. 605β661, Marcel Dekker, New York, NY, USA, 2000. View at Google Scholar - B. V. R. Kumar, S. V. S. S. N. V. G. K. Murthy, V. Sangwan, M. Nigam, and P. Chandra, βNon-darcy mixed convection in a fluid-saturated square porous enclosure under suction effect: part I,β
*Journal of Porous Media*, vol. 13, no. 6, pp. 537β554, 2010. View at Publisher Β· View at Google Scholar Β· View at Scopus - R. Nazar, N. Amin, D. Filip, and I. Pop, βThe Brinkman model for the mixed convection boundary layer flow past a horizontal circular cylinder in a porous medium,β
*International Journal of Heat and Mass Transfer*, vol. 46, no. 17, pp. 3167β3178, 2003. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at Scopus - N. P. Singh, A. K. Singh, and H. Sigh, βMixed convection flow past a porous vertical plate bounded by a porous medium in a rotating system in the presence of magnetic filed,β
*Journal of Porous Media*, vol. 13, no. 7, pp. 623β633, 2010. View at Google Scholar - S. Mahmud and I. Pop, βMixed convection in a square vented enclosure filled with a porous medium,β
*International Journal of Heat and Mass Transfer*, vol. 49, no. 13-14, pp. 2190β2206, 2006. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at Scopus - D. Srinivasacharya and Ch. RamReddy, βMixed convection in a doubly stratified micropolar fluid saturated non-darcy porous medium,β
*The Canadian Journal of Chemical Engineering*, vol. 9999, pp. 1β10, 2011. View at Google Scholar - V. Kurdyumov and A. Linan, βFree and forced convection around line sources of heat and heated cylinders in porous media,β
*Journal of Fluid Mechanics*, vol. 427, pp. 389β409, 2001. View at Google Scholar - K. Vafai and H. Hadim, βOverview of current computational studies of heat transfer in
porous media and their applications-natural and mixed convection,β in
*Advances in Numerical Heat Transfer*, W. J. Mikowycz and E. M. Sparrow, Eds., vol. 2, pp. 331β369, Taylor & Francis, New York, NY, USA, 2000. View at Google Scholar - M. A. Waheed, G. A. Odewole, and S. O. Algbe, βMixed convective heat transfer in rectangular enclosures filled with porous media,β
*ARPN Journal of Engineering and Applied Sciences*, vol. 6, no. 8, pp. 47β60, 2011. View at Google Scholar - N. H. Saeid and I. Pop, βNon-Darcy natural convection in a square cavity filled with a porous medium,β
*Fluid Dynamics Research*, vol. 36, no. 1, pp. 35β43, 2005. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at Scopus - D. A. S. Rees and D. S. Riley, βFree convection above a near horizontal semi-infinite heated surface embedded in a saturated porous medium,β
*International Journal of Heat and Mass Transfer*, vol. 28, no. 1, pp. 183β190, 1985. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at Scopus - D. A. S. Rees and I. Pop, βNon-Darcy natural convection from a vertical wavy surface in a porous medium,β
*Transport in Porous Media*, vol. 20, no. 3, pp. 223β234, 1995. View at Publisher Β· View at Google Scholar Β· View at Scopus - D. A. S. Rees, βThe effect of inertia on free convection from a horizontal surface embedded in a porous medium,β
*International Journal of Heat and Mass Transfer*, vol. 39, no. 16, pp. 3425β3430, 1996. View at Publisher Β· View at Google Scholar Β· View at Scopus - B. V. R. Kumar and Shalini, βDouble diffusive natural convection in a doubly stratified wavy porous enclosure,β
*Applied Mathematics and Computation*, vol. 171, no. 1, pp. 180β202, 2005. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at Scopus