Document Type: Full Lenght Research Article
Authors
^{1} Faculty of Mechanical Engineering, Semnan University, Semnan, Iran
^{2} Faculty of Chemical Engineering, Semnan University, Semnan, Iran
Abstract
Keywords
1.
Introduction
Such condition in which the flow is driven by buoyancy force produced by combined act of thermal and concentration gradients is known as a doublediffusive convection. This phenomenon occurs in many industrial engineering and environmental applications such as drying processes, heat exchangers, chemical processes, energy storages, etc. There are noticeable studies reporting doublediffusive convection in nonporous enclosures, while only a few studies have investigated this phenomenon in a porous enclosure in the literature. Authors draw readers' attention to some more pertinent researches published in this field in the following paragraph.
Researches carried out by Hu and ELWakil [1], Ostrach [2 4], and Lee et al. [5] are primary works related to blend of heat and mass transfer in enclosures. Nithiarasu et al. [6] studied doublediffusive natural convection in a porous enclosure with a freely convecting wall. Their effort exhibited a direct bearing of the buoyancy ratio on flow pattern and thermosolutal behavior. Doublediffusive natural convection in parallelogrammic enclosures filled with fluidsaturated porous media using the Darcy model has been investigated by Costa [7]. He has found that strong changes occur in the parallelogrammic enclosure when the Rayleigh number, the inclination angle, and the aspect ratio of the enclosure change. Goyeau et al. [8] reported a work on doublediffusive natural convection in a porous cavity using the DarcyBrinkmann formulation. They indicate that the impact of the Darcy number on doublediffusion process is more complex than in thermal convection, and also the behavior has been already examined for fluids. Chamkha and AlMudhaf [9] carried out a work about the doublediffusive natural convection in a tilted porous cavity considering heat source or sink. They conclude that as the cavity angle increases the average Nusselt and Sherwood numbers decreases. They also find a critical angle for which the Sh and Nu numbers attain their maximum values. Eyden et al. [10] explored turbulent doublediffusive natural convection in trapezoidal enclosures. They had two main sections; temperature driven natural convection and combined temperature and mass driven natural convection. For the second case, they injected Argon and Nitrogen in the enclosure and studied the effects of injected rates on flow pattern and thermosolutal behaviors. Bourich et al. [11] conducted a work on doublediffusive natural convection within a porous cavity partially heated from below and differentially salted. They find that the effect of buoyancy ratio on the heat, mass, and flow pattern depends on the kind of solution, the Lewis and Rayleigh numbers, and the location of the heat source. The range of buoyancy ratio for which a certain solution is maintained, depends on both Lewis number and Rayleigh number. Doublediffusive natural convection in a square enclosure with discrete heat sources has been reported by Teamah et al [12]. The results reveal that the average Nusselt number and average Sherwood number are amplified with an increment in the Rayleigh number, the dimensionless heater length, and Prandtl number. El Ayachi et al. [13], Bourich et al. [14], Teamah [15], Belazizia et al. [16] and Nikbakhti et al. [17] also reported related studies in this manner.
To the best knowledge of the authors, there is no study to examine doublediffusive natural convection in an arced cavity. This type of geometry is widely employed in the industrial and engineering applications such as underground electric transmission cables, thermal insulation, cryogenics, thermal energy storages, porous heat exchangers, and nuclear waste repositories. This gap in the literature motivates the authors to scrutinize the doublediffusive free convection in a fluid saturated arcshaped porous enclosure.
2. Mathematical Modeling
2.1 Conceptual model
Figure 1 indicates the geometry used in this study. The enclosure has been filled with porous matrix. The fluid
flow is considered to be incompressible, steady, laminar, and twodimensional with constant properties. The density varies in accordance with the Boussinesq model. The outer semicircular wall acts as a cold source and the inner wall acts as a hot source. Two horizontal walls are insulated.
2.2 Governing equations
The two dimensional continuity, momentum, energy, and mass equations for laminar natural convection in polar coordinate system are as following [18]:
(1) 
(2) 
(3) 
(4) 
(5) 
Some equations should be defined as following to evaluate the heat transfer rate, local and average Nusselt number, and Sherwood number:
(6) 

(7) 
Esq. 16 are in dimensionless form. The dimensionless variables may be defined as following:
(8) 
The variable with superscript “” denotes dimensional variables.
2.3 Boundary conditions
According to Fig. 1 the governing equations (14) are subjected to the following conditions:

(9) 

(10) 

(11) 
3. Numerical Technique
The coupled governing continuity, momentum, energy, and concentration equations are discretized using the finite volume method (FVM). The coupling between velocity and pressure in momentum equation is accomplished by using the SIMPLE technique [19]. To discretize the equations, a combination of central difference and upwind methods is employed [20]. To make sure that the convergence occurs in the sequential iterative solution, the relative variations of definite variables between two successive iterations have been required to be less than the specified accuracy level of .
A grid independence study has been examined for
Fig.1 Schematic diagram of the arcshaped porous enclosure 
, , , , and to obtain average Sherwood number. Table 1 shows the effect of the number of grids on the average Sherwood number. The results reveal that the homogeneous grid system of 10171 is acceptable to achieve precise result in the present study.
In order to validate our numerical solution, local results obtained by our code has been compared by the published results of Iwatsu et al. [21]. The engineering parameters for such validation are mimicked with the conditions as invoked by Iwatsu et al. [21]. They investigate heat and flow pattern in liddriven square cavity with two insulated vertical walls while top wall is hot and bottom wall is cold. Comparison includes one graph and one table. In figure 2 Nusselt number at the top wall for Re=100 and Gr=10^{6} has been compared between two works and a good agreement is seen. In table 2, the average Nusselt number at the top wall for Gr=100 at different Re has been compared and the errors has been calculated. These results also show good agreement between these two works.
4. Results and Discussion
In the present section, results are illustrated in the form of streamlines, isotherms, and isoconcentrations. Besides, graphs showing thermal and solutal behaviors of flow are exhibited in this section. All results have been obtained based on pertinent dimensionless parameters of Grashof number, Lewis number, Darcy number, and Buoyancy ratio. The ranges of varying parameters are , , , whereas Darcy number and Prandtl number are fixed at 10^{4} and 1, respectively. All findings are presented through the
Table 1: Effect of grid number on at , , , and 


Number of Grids 
Error 


4121 
3.0108 



6141 
2.9506 
2% 


9161 
2.9329 
0.6% 


10171 
2.9285 
0.15% 


151101 
2.9259 
0.09% 


Fig.2 Comparison of the local Nusselt number between the present code and the results of Iwatsu et al. [21]


Table 2: Comparison of the Average Nu at the top wall between the present work and that of Iwatsu et al. [21] for Gr=100 

Parameter 
Present Work 
Iwatsu et al. [21] 
Error 

Re=100 
1.99 
1.94 
2.57% 

Re=400 
3.90 
3.84 
1.56% 

Re=1000 
6.42 
6.33 
1.42% 

following paragraphs.
Streamlines, temperature and concentration contours have been plotted in figure 3 for N=5 and N=10 at Le=10 and Ra=105. For these two buoyancy ratios, we meet aiding doublediffusive flow, because solutal buoyancy force amplifies thermal buoyancy force. Vortices are evidently strong as the maximum stream functions for N=5 and 10 are 2.57 and 1.52, respectively. For aiding flow case, left hand vortices are clockwise and right hand vortices are counterclockwise. Temperature contours in the second line of the figure indicate that for N=10 the thermal boundary layer on the two sides of hot curved wall are thinner than that for N=5. As a result, the Nusselt number is higher in those areas for N=10. On the contrary, the wall thermal boundary layer for N=10 is thicker than that for N=5 in the middle of hot curved. Concentration contours in the third line prove intense gradients of concentration near the corners of hot wall portraying narrow solutal boundary layer. In the middle of the wall, dense solutal boundary layers extending towards two sides are seen.
Figure 4 indicates pathlines, isotherms, and isoconcentrations for N=5 and N=10 at fixed Lewis number (Le=10) and fixed Rayleigh number (Ra=105). The figure shows opposing flow regime in which solutal buoyancy force opposes the thermal buoyancy force. In contrast to aiding flow, left hand vortices are counterclockwise and right hand vortices are clockwise in this case. As it is expected, the strength of the vortices is lower than that of aiding flow case, so that the maximum stream function for N=5 and 10 reaches 1.38 and 0.82, respectively. Temperature field exemplifies gradients more consistent than aiding flow. It is evident from the figure that the temperature field has been stratified in radial direction, especially for N=10. Concentration field can be seen in the tird line of the figure. Isoconcentration lines are more uniformed than those of aiding flow regime. Concentration gradients in the middle of the hot wall are much denser than the corners of it. This is because of the conflicting effect of solutal buoyancy force. It is worth to mention that N=0 ascribes to negligible species impacts expressing pure thermal convection and N=∞ refers to predominant solutal impacts.
Figure 5 represents the flow pattern, isotherms, and
isoconcentrations for various Lewis numbers at fixed buoyancy ratio, Darcy number, and Rayleigh number. Flow intensity, which is characterized by the maximum stream function, is increased by decreasing the Lewis number. As shown in the figure occurs in Le=10 and occurs in Le=0.1. The figure also displays the vortices' core approaches horizontal walls by increasing Lewis number. Isotherms show uniform temperature distribution in the bulk of the cavity and increasing the Lewis numbers results in more uniform temperature distribution. The third line represents isoconcentrations which have an opposite behavior compared to the isotherms.
Figure 6 displays the contours for temperature and concentration for different values of Da, Le, and Ra numbers and for a fixed value of N. For Da=10^{4} (low Darcy number) the microscopic drag force is high and the fluid experiences a pronounced large resistance as it flows through the porous matrix, causing the flow to cease in the most of the cavity, which suppresses the buoyancy forces. In this case, an increase in Rayleigh number has an insignificant effect on heat and mass transfer because in higher Ra with natural convection dominated regime, porous matrix has a high resistance to the natural flow. Accordingly, the changes in isotherms and isoconcentrations are not remarkable in low Darcy number (Da=10^{4}) as it can be seen in the figure. As the Darcy number increases (Da=10^{2}), the resistance of porous medium decreases and the fluid easily flows in the cavity. In this case, by increasing the Rayleigh number the natural convection regime will be dominant and due to of low resistance of porous matrix, the alterations from conduction dominated regime to convection dominated one are totally noticeable in the figure.
Figure 6 also shows that at high Lewis number (Le=10) the Isotherms are almost parallel to the curved walls of cavity for nearly all values of Rayleigh numbers, indicating that most of the heat transfer is carried out by heat conduction. This is due to an increase in the thermal boundarylayer thickness. On the contrary, the isoconcentration contours are enhanced for higher Lewis number (Le=10). This happens because the solutal boundarylayer thickness decreases and consequently enhances thermosolutal activities in the enclosure when the Lewis number increases.
Sherwood number and Nusselt number along the hot curved wall for various Lewis numbers and N=10 are depicted in figures 7 and 8. For both cases, the maximum value occurs in the middle of the cavity, whereas the minimum value occurs in angles 0^{ͦ} and 180^{ͦ}. This happens because the denser gradients and accordingly thinner boundary layers take place in the middle of the cavity. The distinct difference between the two graphs is that the maximum value for Sh number increases with increasing Le number, whereas the maximum value for
Nu decreases with increasing Le. Actually, an increase in the Lewis number decreases the depth of solutal
boundary layer leading to amplifying Sherwood number, while it results in a reduction of the flow intensity, which leads to decreasing Nusselt number. Furthermore, it can be seen thet the greater mass transfer rate occurs at the corners of the cavity at Le=0.1, because the mass boundary layer is thinner for this case at the pointed areas.
Figures 9 and 10 exhibit Sherwood number and Nusselt number along the hot curved wall for different values of Le at N=10. As we know N<0 represents the aiding flow regime. Hence, opposite behavior of the previous graph is expected. In this case, the minimum value occurs in the middle of the cavity, whereas the maximum value occurs in angles 0^{ͦ} and 180^{ͦ}. The minimum value of Sh for Le=0.1 is the greatest one, because the boundary layer is the thinnest for this Le number. For other values of Le, the minimum values are almost equal. Thermal behavior is different, so that from angle 0^{ͦ} to angle 45^{ͦ} Nusselt number augments with an increment in Lewis number, but from angle 45^{ͦ} to angle 90^{ͦ} it reduces with an increment in Lewis number. This behavior recurs for angles between 90^{ͦ} to 180^{ͦ}, because the symmetry line locates at the angleof 90^{ͦ}.
The effect of buoyancy ratio on average Sherwood number and average Nusselt number for different Lewis numbers are presented in Figures 11 and 12. For aiding flow case (N<0), Nu and Sh decrease with increasing N, but Nu and Sh increase with increasing N at opposing flow case (N>0). The values of Nu and Sh for aiding flows are greater than those of opposing flows at equal absolute value of N, because the magnitude of the velocity and the flow intensity are greater for aiding flows. For instance, at Le=10, for N=10 (aiding flow)
Fig.9 Sherwood number along the hot curved wall for N=10, Ra=1e5, Da=1e4, and for various Le numbers 
Fig.10 Nusselt number along the hot curved wall for N=10, Ra=1e5, Da=1e4, and for various Le numbers 
Fig.11 Average Sherwood number along the hot curved wall versus of buoyancy ratio for Ra=1e5, Da=1e4, and for various Le numbers 
Fig.12 Average Nusselt number along the hot curved wall versus buoyancy ratio for Ra=1e5, Da=1e4, and for various Le numbers 
average Sherwood and Nusselt numbers are 8.93 and 2.84, respectively, while for N=10 (opposing flow) average Sherwood and Nusselt numbers are 5.66 and. 2.28, respectively. Moreover, the figures show that the minimum values in both average Nusselt and Sherwood numbers come about at a buoyancy ratio of 1 (N=1). The last point deduced from the graphs is that for a certain value of N, an increment in the Lewis number, which means an increase of mass diffusivity regarding to thermal diffusivity, results in an increase/decrease of concentration/temperature gradients
5. Conclusions
In the current study, doublediffusive natural convection in a porous arcshaped enclosure has been numerically studied using DarcyBrinkman formulation. The impacts of Darcy number, Rayleigh number, Buoyancy ratio, and Lewis number on flow pattern and heat and mass transfer process have been investigated. It was shown for aiding flow (N<0) the strength of the flow is greater than that of opposing flow (N>0). It leads to greater heat and mass transfers for aiding flow case. Any increment in the Darcy number varies heat and mass profiles visibly, as isotherms and isoconcentrations are distributed in the bulk of the cavity for higher Darcy numbers. An increment in the Lewis number results in an improvement/deterioration of mass/heat transfer process.
Nomenclature 

Concentration (Kg m3) 

Concentration at cold wall (Kg m3) 

Concentration at hot wall (Kg m3) 

D 
Binary diffusion coefficient 

Da 
Darcy number, 

Le 
Lewis number, 

N 
Buoyancy ratio, 

Local Nusselt number () 

Average Nusselt number () 

Dimensionless Pressure () 

Prandtl number, 

R 
Radial direction 

Rayleigh number, 

Reynolds number, 

Sh 
Sherwood number () 

Temperature (K) 

Temperature of cold wall (K) 

Temperature of hot wall (K) 

Velocity (m/s) 

Greek 

Thermal diffusivity (m2s1), 

Coefficient of Thermal expansion () 

Coefficient of Solutal expansion () 

Porosity 

Thermal conductivity (W/mK) 

Dynamic viscosity (Nsm2) 

Kinematic viscosity (m2s1), 

Angular direction 

Density (kg/m3) 

Stream function 

Superscript 

* 
Dimensional parameters 

References
[1]. C.Y. Hu, M.M. ELWakil, Simultaneous heat and mass transfer in a rectangular cavity, Proceeding of International Heat Transfer Conference 5 (1974) 2428.
[2]. S. Ostrach, Natural convection with combined driving forces, Physico Chemical Hydrodynamics 1 (1980) 233247.
[3]. S.Ostrach, Natural convection heat transfer in cavities and cells, 7th International Heat Transfer Conference 1 (1983) 365379.
[4]. S. Ostrach, Fluid mechanics in crystal growththe 1982 freeman scholar lecture, Journal of Fluid Engineering 105 (1983) 520.
[5]. T.S. Lee, P.G. Parikh, A. Acrivos, D. Bershader, Natural convection in a vertical channel with opposing buoyancy forces, International Journal of Heat and Mass Transfer 25 (1982) 499511.
[6]. P. Nithiarasu, T. Sundararajan, K.N. Seetharamu, Doublediffusive natural convection in a fluid saturated porous cavity with a freely convecting wall, International Communications in Heat and Mass Transfer 24 (8) (1997) 1121–1130.
[7]. V.A.F. Costa,Doublediffusive natural convection in parallelogrammic enclosures filled with fluid saturated porous media, International Journal of Heat and Mass Transfer 47 (2004) 2699–2714
[8]. B.Goyeau, J.P. Songbe and D. Gobin, Numerical study of doublediffusive natural convection in a porous cavity using the Darcybrinkman formulation, International Journal of Heat and Mass Transfer, 39 (7) (1996), 13631378.
[9]. A.J. Chamkha, A. AlMudhaf, Doublediffusive natural convection in inclined porous cavities with various aspect ratios and temperaturedependent heat source or sink, Heat and Mass Transfer/Waerme und Stoffuebertragung 44 (6) (2008) 679–693.
[10]. J.T. Vander Der Eyden, TH. H. Van Der Meer, K. Hanjalic, E. Biezen, J. Bruining, Doublediffusive natural convection in trapezoidal enclosuresm International Journal of Heat and Mass transfer, 41 (13) (1998), 18851898.
[11]. M. Bourich, M. Hasnaoui, A. Amahmid, Doublediffusive natural convection in a porous enclosure partially heated from below and differentially salted, International Journal of Heat and Fluid Flow 25 (2004) 1034–1046.
[12]. M.A. Teamah, M.M.Khairat Dawood, M. ElMaghlany, Doublediffusive natural convection in a square cavity with segmental heat sources, European Journal of Scientific Research, 54 (2) (2011), 287301.
[13]. R. El Ayachi, A. Raji, M. Hasnaoui, A. Abdelbaki, M. Naimi, Resonance of doublediffusive convection in a porous medium heated with a sinusoidal exciting temperature, Journal of Applied fluid Mechanics, 3 (2) (2010), 4352.
[14]. M. Bourich, A. Amahmid, M. Hasnaoui, Doublediffusive convection in a porous enclosure submitted to cross gradients of temperature and concentration, Energy Conversion and Management, 45 (2004), 1655–1670.
[15]. M.A. Teamah, Doublediffusive laminar natural convection in a symmetrical trapezoidal enclosure, Alexandria Engineering Journal, 45 (3) (2006), 251263.
[16]. A. Belazizia, S. Benissaad, S. Abboudi, Doublediffusion natural convection of binary fluid in a square enclosure with top Active vertical wall, Advances in Theoretical and Applied Mechanics, 5 (3) (2012), 119131.
[17]. R. Nikbakhti, A. rahimi, Doublediffusive natural convection in a rectangular cavity with partially thermally active side walls, Journal of the Taiwan Institute of Chemical Engineers 43 (2012) 535–541.
[18]. M.S.Valipour, A.Z. Ghadi, Numerical investigation of fluid flow and heat transfer around a solid circular cylinder utilizing nanofluid, International communications in heat and mass transfer, 38 (2011), 12961304.
[19]. S.V. Patankar, Numerical Heat Transfer and Fluid Flow, Hemisphere Publishing Corporation, Washington, 1980.
[20]. H.K. Versteeg, W. Malalasekera, An introduction to computational fluid dynamics. The finite volume method, John Wiley & Sons Inc, New York, 1995.
[21]. R. Iwatsu, J.M. Hyun, K. Kuwahara, Mixed convection in a driven cavity with a stable vertical temperature gradient, International Journal of Heat and Mass Transfer 36 (1993) 16011608.