Next Article in Journal
A Novel Low-Query-Budget Active Learner with Pseudo-Labels for Imbalanced Data
Next Article in Special Issue
MHD Natural Convection and Radiation over a Flame in a Partially Heated Semicircular Cavity Filled with a Nanofluid
Previous Article in Journal
Security Risks to Petroleum Industry: An Innovative Modeling Technique Based on Novel Concepts of Complex Bipolar Fuzzy Information
Previous Article in Special Issue
Generalized Exp-Function Method to Find Closed Form Solutions of Nonlinear Dispersive Modified Benjamin–Bona–Mahony Equation Defined by Seismic Sea Waves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Double Diffusive Natural Convection in a Square Cavity Filled with a Porous Media and a Power Law Fluid Separated by a Wavy Interface

1
Mechanical Engineering Department, College of Engineering, University of Ha’il, Ha’il 81451, Saudi Arabia
2
Department of Mathematics, Capital University of Science and Technology, Islamabad 44000, Pakistan
3
Department of Industrial Engineering and Systems, College of Engineering, Princess Nourah bint Abdulrahman University, Riyadh 11671, Saudi Arabia
4
Department of Mathematics & Statistics, PMAS-Arid Agriculture University, Rawalpindi 46300, Pakistan
5
Department of Mechanical Engineering, College of Engineering, Imam Mohammad Ibn Saud Islamic University, Riyadh 11432, Saudi Arabia
6
Laboratory of Metrology and Energy Systems, University of Monastir, Monastir 5000, Tunisia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(7), 1060; https://doi.org/10.3390/math10071060
Submission received: 11 January 2022 / Revised: 25 February 2022 / Accepted: 20 March 2022 / Published: 25 March 2022

Abstract

:
This study deals with the influence of a wavy interface separating two layers filled with power law fluid and porous media, respectively. The governing equations are solved using the Finite Element Method (FEM) and the numerical model is validated by comparing with experimental findings. The parameters governing the studied configuration are varied as: Rayleigh number (103Ra ≤ 106), power law index (0.6 ≤ n ≤ 1.4), Darcy number (10−2Da ≤ 10−6), buoyancy ratio (0.1 ≤ N ≤ 10) and Lewis number (1 ≤ Le ≤ 10). It is inferred that the temperature gradient increases by augmenting the Rayleigh number, as the flow is observed from the vertical to horizontal direction in both layers. Constant enhancement in the heat and mass transfer is also observed by enriching the buoyancy effect. Moreover, the average Nusselt and Sherwood numbers decline by increasing the width of the porous layer.

1. Introduction

The double diffusion is a phenomenon that occurs in several engineering and industrial applications such as pulp paper, oil drilling, heat removal, heat storage, and food processing [1,2,3,4,5,6]. There are relatively few works that have been conducted on non-Newtonian fluid, instead of Newtonian fluid. Al-Amir et al. [7] reported the impact of the Prandtl number on the natural convection in cavity-containing non-Newtonian and nanofluid porous mediums, which are separated by the sinusoidal interface. It has been realized that the average Nusselt number rises by enhancing the Darcy and Prandlt number, and reduces by enriching the power law index. Alsabery et al. [8] examined the trapezoidal-shaped cavity including two layers (porous and non-Newtonian). The results confirmed that the flow rises remarkably when using silver nanofluid, and the effectiveness of heat transfer is perceived by varying the angle of inclination. The impact of wavy interfaces on the natural convection in non-Darcy porous cavities is studied by Nguyen et al. [9], using the incompressible smoothed particle hydrodynamics (ISPH) method. The results indicated that the average Nusselt number decreases by increasing the amplitude and undulation number of the interface separating the layers. Alsabery et al. [10] performed a numerical study to investigate the effect of the inclination on natural convection in a cavity filled with a porous media and a non-Newtonian fluid. The findings showed that convection is more intense for lower values of power law index. Power law fluids have attracted the researchers due to their importance in several engi-neering applications [11]. More related studies can be found in References [12,13,14,15,16]. The natural convection in porous layers has been studied Al-Srayyih et al. [17] using Galerkin Finite Element Method (GFEM). It was found that the enhancement of the heat transfer occurs using nanofluids. Barnoon et al. [18] investigated the coupled radiation-convection in a cavity filled with a non-Newtonian fluid and equipped with internal obstacles. The authors concluded that the tilt angle has an important effect of the values of Nusselt number. Jabber et al. [19] investigated the 2D natural convection in an enclosure filled with a porous medium saturated with a anofluid and equipped by wavy interfaces. The authors mentioned that power law index causes a reduction of the heat transfer rate. The effect of a power law fluid on natural convection in a cavity having wavy wall has been investigated by Chen et al. [20]. It was shown that the rate of heat transfer of pseudoplastic fluid is better compared to Newtonian fluid. Kefayati et al. [21] applied LBM to realize the behavior of non-Newtonian fluid under natural convection; a uniform magnetic field is also applied in this study. The results indicated that the heat transfer increases with the power law index and decreases with the Hartman number. Saleh et al. [22] studied a differentially heated cavity equipped with rotating obstacles. The results showed that the values of Nusselt number remain constant at L/D > 0.77. Other research papers related to the effect of separating interfaces in cavities are reported [23,24,25,26,27,28]. Turan et al. [29] studied the natural convection of a non-Newtonian fluid. They investigated the influence of the Prandtl number, and proposed new correlations for the Nusselt number for Newtonian and non-Newtonian fluids.
Based on the above-described literature review, and to the best of the authors’ knowledge, there is currently no published work on wavy interfaces in square cavities filled with two layers (non-Newtonian and porous (Newtonian)). The applications of the current study are in engineering and industry. We can further study this by increasing the layers. The effect of various controlled parameters, such as the Rayleigh number, power law index, Darcy number, buoyancy ratio, and Lewis number, is examined in detail.

2. Problem Formulation

The problem configuration is schematically illustrated in (Figure 1). It consists of 2D double-diffusive convection in a cavity containing two fluid layers (non-Newtonian and porous Newtonian fluid), which are separated by a sinusoidal wall. L and H represent the length and height of the cavity, respectively. θ h and C h are the high temperature and concentration at the left wall, respectively. Similarly, θ c and C c are the low temperature and concentration at the right wall, respectively. The remaining horizontal walls are considered to be adiabatic. The flow is considered to be laminar, steady, and incompressible, with the application of the Darcy–Brinkman–Forchheimer model.
The wavy interface is derived from the following equation:
X = H p + A   S i n 2 π K L   Y + φ
where Hp and A represent the width of the porous layer and the amplitude, respectively. φ denotes the phase shift, which is taken as π 2 .

3. The Governing Equations

Under the above-considered assumptions, the dimensionless systems of the equations for the fluid region and porous layer are written as follows [30,31,32,33,34,35,36]:
For the non-Newtonian fluid layer:
U X + V Y = 0
U U X + V U Y = P X + P r R a 2 X μ a M U X + Y μ a M U Y + V X
U V X + V V Y = P Y + P r R a 2 Y μ a M V Y + X μ a M U Y + V X + P r θ + N C
U θ X + V θ Y = 1 R a 2 θ X 2 + 2 θ Y 2
U C X + V C Y = 1 L e   R a 2 C X 2 + 2 C Y 2
μ a = M 2 U X 2 + V Y 2 + V X + U Y 2 n 1 2
where N is the buoyancy ratio and n is the power law index.
The following also applies:
P r = μ a ρ α
R a = ρ β g y L 3 T H T C μ a α
L e = a e D
For the porous fluid layer:
U ¯ X + V ¯ Y = 0
U ¯ U ¯ X + V ¯ U ¯ Y = P ¯ X + P r m R a m 2 U ¯ X 2 + 2 U ¯ Y 2 P r m   U ¯ D a m   R a m U ¯   U ¯ D a m   1.75 150
U ¯ V ¯ X + V ¯ V ¯ Y = P ¯ Y + P r m R a m 2 V ¯ X 2 + 2 V ¯ Y 2 + P r m θ ¯ + N C ¯ P r m   V ¯ D a m   R a m V ¯   U ¯ D a m   1.75 150
U ¯ θ ¯ X + V ¯ θ ¯ Y = 1 R a m 2 θ ¯ X 2 + 2 θ ¯ Y 2
U ¯ C ¯ X + V ¯ C ¯ Y = 1 L e   R a m 2 C ¯ X 2 + 2 C ¯ Y 2
where the following parameters are entered in the above problem:
D a = K L 2
The relationship among the actual and modified Prandtl, Rayleigh and Darcy numbers can be written as follows: P r m = P r ,   R a m = R a and D a m = D a .
The boundary conditions of the proposed problem in the corresponding regions are as follows:
  • U ¯ = 0, θ ¯ = 1, V ¯ = 0, C ¯ = 1 (left wall)
  • U = 0, θ = 0, V = 0, C = 0 (right wall)
  • θ ¯ Y = C ¯ Y = 0 , U ¯ = V ¯ = 0 and θ Y = C Y = 0 , U = V = 0 (top/bottom walls)
The numbers on the vertical hot wall are given as follows:
N u = θ ¯ X X = 0 :   Nusselt   number   ( local )
S h = C ¯ X X = 0 :   Sherwood   number   ( local )
N u a v g = 0 1 N u   d Y :   Nusselt   number   ( average )
S h a v g = 0 1 S h   d Y   :   Sherwood   number   ( average )

4. Solution Methodology

The mathematical models presented in Equations (2)–(7) and (11)–(15) is solved using the higher order Galerkin finite element method. As a first step, a weak formulation is developed by choosing a suitable test space. Afterwards, a hybrid mesh, consisting of both triangular and quadrilateral elements, is generated to cover the computational domain. A finite element method involving the cubic polynomials (P3) is implemented to compute the velocity, temperature, and concentration fields, while the pressure is approximated by the quadratic (P2) finite element space of functions. The stability and robustness of this higher-order pair of FEMs has been tested in [37]. The system of discretized equations is simplified using the adaptive Newton’s method. For further details regarding the solver, the reader is referred to [38].

4.1. Grid Convergence

As shown in Table 1, eight number of elements are compared to check the mesh independency at Gr = 105, Da = 10−3, ϵ = 0.75 and Hp = 0.1. From this comparison it is clear that the deviations of Nuav and Shav between the grids G7 and G8 are very small. Thus, for time economy and results accuracy the grid G7 was retained for all the performed simulations. The grid for the proposed computational model is displayed in Figure 2.

4.2. Code Validation

The current results are checked by comparing them to those of Gibanov et al. [39], by evaluating the average Nusselt number under the same conditions. As is shown in Table 2, the error did not exceed 0.36%, which represents good agreement among the two results. Moreover, the numerical code is checked by comparing the isotherms and streamlines with those found by Gibanov et al. [39]. The results of the current numerical model are in concordance with the results of Gibanov et al. [39] (Figure 3).
In addition to the above validation, to strengthen the reliability of the implemented FEM, we also confirmed the validation by comparing it with the experimental study of Corvaro and Paroncini [40], as demonstrated in Figure 4.

5. Results and Discussion

The results of the different numerical simulation cases are presented by streamlines, isotherms and isoconcentration contours, as well as variations in the local Nu and Sh, and variations in the average Nu and Sh. All the simulations are performed for the fixed values of n = 0.6, Da = 10−3, ϵ = 1.0, Le = 2.5, Pr = 1, N = 0.1, M = 2, Hp = 0.5, A = 0.05, and Ra = 105.
For each case, one parameter is varied. In fact, the results will consecutively present the effect of the variation in the Rayleigh number (Ra), power law index (n), Darcy number (Da), buoyancy number (N), and Lewis number (Le). Furthermore, the effects of the width of the porous media (Hp), the undulation number of interface (K), and its amplitude (A) on the average Nu and Sh are presented for a range of the power law index (n), varied from 0.6 to 1.6.

5.1. Impact of Rayleigh Number

Figure 5 presents the effect of the Rayleigh number on the streamlines (left side), isotherms (center), and isoconcentrations (right side). Due to its great effect on the heat transfer rate, the Rayleigh number is varied from 103 to 106. For a low Rayleigh number (Ra = 103), the center of the main vortex is located on the non-Newtonian side of the cavity and rotates counterclockwise. By increasing the Ra number, the rotation of the main vortex is intensified. Remarkable penetration in the direction of the porous layer is detected, due to the elongation of the main vortex.
As can be observed for Ra = 106, the intensity of the streamlines occurs at the core of the cavity, indicating an increase in the heat transport, due to the natural convection. Here, it is noticed that the center of the vortex is still on the non-Newtonian side. By increasing the Ra number, the behavior of the isotherms is modified, due to the enhancement of the temperature gradient. When the Ra number is equal to 103, the isotherms are parallel to each other, indicating the domination of the conductive heat transfer mode, but as its value rises, an increase in the temperature gradient is observed. The thermal flow passes from a conductive regime to a convective regime by increasing the Ra number. Indeed, the stratification of the isotherms is pronounced when the Ra number passes from 103 to 106. The thermal gradient is more densely packed close to the downside of the hot wall, indicating an important thermal boundary layer. This fact is confirmed by the coordinates of the height of the local Nusselt numbers in Figure 4. The isoconcentration behavior is greatly affected by the increase in the Ra numbers, which indicates a progressive move from a diffusive to convective mass transfer. When Ra = 106, the solutal gradient near the hot wall is important, which indicates the importance of the convective transfer from the non-Newtonian fluid layer to the porous layer.
Figure 6 presents the variation in the local Nu and Sh for different Ra numbers. The curves prove the high increase in the Nu and Sh values when the Rayleigh number passes from 105 to 106. Table 3 presents a comparison between the average Nu and average Sh for various Rayleigh numbers, under the same conditions. For Ra = 103 and Ra = 104, the Nu values are low compared to the other Rayleigh values. This is due to the lower heat exchange between the porous medium layer and the non-Newtonian layer. The maximum values of the average Nu and Sh numbers confirm the importance of the mass transfer effect on behalf of the heat transfer, especially when the Ra numbers increase. This is due to the current case condition, in which the buoyancy number is equal to 0.1.

5.2. Impact of Power Law Index

Figure 7 presents the streamlines (left side), the isotherms (center), and the isoconcentration (right side). By definition, the power law index affects the heat generation and, consequently, the viscosity. The expected results need to demonstrate the regression of the heat transfer and show a great impact on the mass transfer between the two layers. As the power law index rises, the intensity and movement of the streamlines are slightly affected in the cavity, which means that the energy required to rotate the vortex decreases. This fact is explained further by plotting the variation in the local Nu and Sh, and calculating their mean values at the hot wall (Table 4). Regardless of the fact that they exhibit the same behavior, it is clear in Figure 8 that the variation in the local Nu and Sh drops by increasing the n values.

5.3. Impact of Darcy Number

Figure 9 presents the effect of the Darcy number (Da) variation on the streamlines, isotherms, and isoconcentrations, while the other parameters are constant. When the Da is equal to 10−6 and 10−5, the streamlines are localized in the non-Newtonian zone, which indicates that no exchange is performed between the porous and non-Newtonian layers. The isothermal and isoconcentration lines are parallel in the porous zone, compared to the non-Newtonian zone. By increasing the Darcy numbers, the spread of the streamline through the wavy interface is increasingly developed. The thermal and solutal gradient near the downside of the hot wall is noticed, due to the rise in the Darcy numbers. The results show that for Da = 10−2, thermal stratification takes place and solutal distortion appears in the core cavity, since the value of the Darcy number increases to 10−3.
The variation in the local Nu and Sh for different Darcy numbers (Da) is shown in Figure 10. For low values of Darcy numbers (10−6 and 10−5), the heat and mass transfer are almost constant and at low levels. A slight increase is detected when the Darcy number is equal to 10−4. A noticeable increase in the heat and mass transfer effect takes place for Da = 10−3 and 10−2, which is described by the rise in the vortex rotation intensity through the two sides of the wavy interface. Table 5 presents the evolution of the average Nu and Sh for the various values of Da, when Ra = 105, Da = 10−3, ϵ = 0.75, and Hp = 0.1. The great increase in the heat and mass transfer values can be deduced by the increase in the diffusivity of the fluid flow through the interface to the porous medium.

5.4. Impact of Buoyancy Ratio

The effect of the variation in the buoyancy number N from 0.1 to 10, when the other values remain constant, is presented in Figure 11. The behavior of the streamlines is similar to the Ra variation cases. The main cell is lengthened to the whole cavity space by increasing the N values, although the center of the vortex is still localized in the non-Newtonian zone. The thermal gradient is densely packed near the downside of the hot wall in the porous media zone by the increase in the N values, which results in intensification of the heat transfer through the interface. Moreover, the solutal gradient is tightened near the downside of the hot wall and the topside of the cold wall, and stratification of the isoconcentrations takes place in the core region of the cavity.
These behaviors could be explained by the profiles of the local Nu and Sh in Figure 12, and their average values in Table 6. The enhancement of the heat and mass transfer through the interface is measured by increasing the buoyancy effect.

5.5. Impact of Lewis Number

Figure 13 indicates the streamlines at the left side, the isotherms at the center, and the isoconcentration at the right side. It can be observed that there is not a significant change in the streamline (left side), isotherm (center), and isoconcentration (right side) contours by enhancing the Lewis number for both layers. This is confirmed by Figure 14, where the behavior of the local Nu is almost the same for different values of Lewis number. On the other hand, the local Sh is enriched by increasing the Lewis number. The particular reason for this minor change in these results is the small value of the buoyancy ratio (N = 0.1). Furthermore, it is observed that there is a minor effect of concentration in the momentum equation. Table 7 depicts the trend of the average Nu and Sh for different Lewis numbers (Le), which reflects the minor decline in the local Nusselt number as the Lewis number rises, but the average Sherwood numbers increase considerably.
In the next portion of this discussion, some plots of the average Nu and Sh are drawn on the power law index. Figure 15 indicates that the average Nu and Sh decline with an increase in the width of the porous layer with power law indices. The behavior of the average Nu and Sh with the undulation parameter and amplitude of the interface is shown in Figure 16 and Figure 17, respectively. It is visualized, in both figures, that the results switch towards the decline direction, by increasing the values of K and A, correspondingly.

6. Conclusions

Double-diffusive natural convection has been investigated in a 2D cavity, composed of two layers (non-Newtonian and porous), separated by a wavy interface, using FEM. The effects of the Rayleigh number, power law index, Darcy number, buoyancy ratio, and Lewis number are studied. The main conclusions can be summarized as follows:
  • The temperature gradient rises by augmenting the Rayleigh number, as the flow is observed from the vertical to horizontal direction in both layers of the cavity, which seems to be a remarkable change in heat transfer. Similar behavior is observed in the concentration case.
  • By increasing the values of the power law indices, a decline in the average Nu and Sh is observed.
  • When Da  10 6 to 10 5 , the streamlines are localized in the non-Newtonian zone, which may indicate that no exchange is observed between the porous and non-Newtonian layers. By increasing the value of Da, the wavy interface is changed.
  • Constant enhancement in heat and mass transfer is noticed in a wavy interface by enriching the buoyancy effect.
  • The local Nu is almost the same for different values of Lewis number. On the other hand, the local Sh is enriched by increasing the Lewis number.
  • The average Nusselt and Sherwood numbers decline by increasing the width of the porous layer with power law indices.
  • It is visualized, in both figures, that the average Nu and Sh move towards the decline direction, by increasing the values of K and A, correspondingly.

Author Contributions

Conceptualization, L.K., S.H. and C.M.; methodology, K.G., C.M. and M.J.; software, L.K. and S.H.; validation, K.G. and M.J.; formal analysis, L.K. and S.H.; investigation, K.G. and M.J.; resources, K.G.; data curation, K.G., C.M. and M.J.; writing—original draft preparation, L.K., S.H., C.M., K.G. and M.J.; writing—review and editing, K.G. and M.J.; project administration, K.G.; funding acquisition, K.G. All authors have read and agreed to the published version of the manuscript.

Funding

Princess Nourah bint Abdulrahman University Researchers Supporting Project number PNURSP2022R41, Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kefayati, G.; Tang, H. Three-dimensional Lattice Boltzmann simulation on thermosolutal convection and entropy generation of Carreau-Yasuda fluids. Int. J. Heat Mass Transf. 2019, 131, 346–364. [Google Scholar] [CrossRef]
  2. Xin, S.; Le Quéré, P.; Tuckerman, L.S. Bifurcation analysis of double-diffusive convection with opposing horizontal thermal and solutal gradients. Phys. Fluids 1998, 10, 850–858. [Google Scholar] [CrossRef] [Green Version]
  3. Sezai, I.; Mohamad, A.A. Double diffusive convection in a cubic enclosure with opposing temperature and concentration gradients. Phys. Fluids 2000, 12, 2210–2223. [Google Scholar] [CrossRef]
  4. Kuznetsov, G.V.; Sheremet, M.A. A numerical simulation of double diffusive conjugate natural convection in an enclosure. Int. J. Therm. Sci. 2011, 50, 1878–1886. [Google Scholar] [CrossRef]
  5. Khezzar, L.; Siginer, D.; Vinogradov, I. Natural convection of power law fluids in inclined cavities. Int. J. Therm. Sci. 2012, 53, 8–17. [Google Scholar] [CrossRef]
  6. Ali, L.; Ali, B.; Liu, X.; Ahmed, S.; Shah, M.A. Analysis of bio-convective MHD Blasius and Sakiadis flow with Cattaneo-Christov heat flux model and chemical reaction. Chin. J. Phys. 2021; in press. [Google Scholar] [CrossRef]
  7. Al-Amir, Q.; Ahmed, S.Y.; Hamzah, H.K.; Ali, F.H. Effects of Prandtl Number on Natural Convection in a Cavity Filled with Silver/Water Nanofluid-Saturated Porous Medium and Non-Newtonian Fluid Layers Separated by Sinusoidal Vertical Interface. Arab. J. Sci. Eng. 2019, 44, 10339–10354. [Google Scholar] [CrossRef]
  8. Alsabery, A.; Chamkha, A.; Hussain, S.; Saleh, H.; Hashim, I. Heatline visualization of natural convection in a trapezoidal cavity partly filled with nanofluid porous layer and partly with non-Newtonian fluid layer. Adv. Powder Technol. 2015, 26, 1230–1244. [Google Scholar] [CrossRef]
  9. Nguyen, M.T.; Aly, A.; Lee, S.-W. Effect of a wavy interface on the natural convection of a nanofluid in a cavity with a partially layered porous medium using the ISPH method. Numer. Heat Transf. Part A Appl. 2017, 72, 68–88. [Google Scholar] [CrossRef]
  10. Alsabery, A.; Hussain, S.; Saleh, H.; Hashim, I. Inclination angle effect on natural convection in a square cavity partially filled with non-Newtonian fluids layer. AIP Conf. Proc. 2015, 1678, 060007. [Google Scholar]
  11. Bin Kim, G.; Hyun, J.M.; Kwak, H.S. Transient buoyant convection of a power-law non-Newtonian fluid in an enclosure. Int. J. Heat Mass Transf. 2003, 46, 3605–3617. [Google Scholar] [CrossRef]
  12. Wu, W.-T.; Massoudi, M. Recent Advances in Mechanics of Non-Newtonian Fluids. Fluids 2020, 5, 10. [Google Scholar] [CrossRef] [Green Version]
  13. Lamsaadi, M.; Naimi, M.; Hasnaoui, M.; Mamou, M. Natural Convection in a Vertical Rectangular Cavity Filled with a Non-Newtonian Power Law Fluid and Subjected to a Horizontal Temperature Gradient. Numer. Heat Transf. Part A Appl. 2006, 49, 969–990. [Google Scholar] [CrossRef]
  14. Lamsaadi, M.; Hasnaoui, M. Natural convection of non-Newtonian power law fluids in a shallow horizontal rectangular cavity uniformly heated from below. Heat Mass Transf. 2005, 41, 239–249. [Google Scholar] [CrossRef]
  15. Alsabery, A.I.; Chamkha, A.J.; Saleh, H.; Hashim, I. Natural Convection Flow of a Nanofluid in an Inclined Square Enclosure Partially Filled with a Porous Medium. Sci. Rep. 2017, 7, 2357. [Google Scholar] [CrossRef]
  16. Chamkha, A.J.; Ismael, M. Natural Convection in Differentially Heated Partially Porous Layered Cavities Filled with a Nanofluid. Numer. Heat Transf. Part A Appl. 2014, 65, 1089–1113. [Google Scholar] [CrossRef]
  17. Al-Srayyih, B.M.; Gao, S.; Hussain, S.H. Effects of linearly heated left wall on natural convection within a superposed cavity filled with composite nanofluid-porous layers. Adv. Powder Technol. 2019, 30, 55–72. [Google Scholar] [CrossRef]
  18. Barnoon, P.; Toghraie, D.; Dehkordi, R.B.; Afrand, M. Two phase natural convection and thermal radiation of Non-Newtonian nanofluid in a porous cavity considering inclined cavity and size of inside cylinders. Int. Commun. Heat Mass Transf. 2019, 108, 104285. [Google Scholar] [CrossRef]
  19. Jabbar, M.Y.; Ahmed, S.Y.; Hamazh, H.; Ali, F.H.; Khafaji, S.O.W. Effect of Layer Thickness on Natural Convection in a Square Enclosure Superposed by Nano-Porous and Non-Newtonian Fluid Layers Divided by a Wavy Permeable Wall. Int. J. Mech. Mechatron. Eng. 2019, 19, 3. [Google Scholar]
  20. Chen, C.K.; Cho, C.C. Natural Convection in Wavy-Wall Cavities Filled with Power-Law Fluid. Int. J. Mech. Mechatron. Eng. 2014, 8, 1387–1391. [Google Scholar]
  21. Kefayati, G.; Tang, H. MHD thermosolutal natural convection and entropy generation of Carreau fluid in a heated enclosure with two inner circular cold cylinders, using LBM. Int. J. Heat Mass Transf. 2018, 126, 508–530. [Google Scholar] [CrossRef]
  22. Saleh, H.; Alsabery, A.I.; Hashim, I. Natural convection in polygonal enclosures with inner circular cylinder. Adv. Mech. Eng. 2015, 7, 1687814015622899. [Google Scholar] [CrossRef] [Green Version]
  23. Massarotti, N.; Nithiarasu, P.; Zienkiewicz, O. Natural Convection in Porous Medium-Fluid Interface Problems: A Finite Element Analysis by Using the CBS Procedure. Int. J. Numer. Methods Heat Fluid Flow 2001, 11, 473–490. [Google Scholar] [CrossRef]
  24. Gobin, D.; Goyeau, B.; Neculae, A.P. Convective Heat and Solute Transfer in Partially Porous Cavities. Int. J. Heat Mass Transf. 2005, 48, 1898–1908. [Google Scholar] [CrossRef]
  25. Hirata, S.C.; Goyeau, B.; Gobin, D.; Carr, M.; Cotta, R. Linear Stability of Natural Convection in Superposed Fluid and Porous Layers: Influence of the Interfacial Modelling. Int. J. Heat Mass Transf. 2007, 50, 1356–1367. [Google Scholar] [CrossRef]
  26. Aguilar-Madera, C.G.; Valdés-Parada, F.J.; Goyeau, B.; Ochoa-Tapia, J.A. One-Domain Approach for Heat Transfer Between a Porous Medium and a Fluid. Int. J. Heat Mass Transf. 2011, 54, 2089–2099. [Google Scholar] [CrossRef]
  27. Aguilar-Madera, C.G.; Valdés-Parada, F.J.; Goyeau, B.; Ochoa-Tapia, J.A. Ochoa-Tapia, Convective Heat Transfer in a Channel Partially Filled with a Porous Medium. Int. J. Therm. Sci. 2011, 50, 1355–1368. [Google Scholar] [CrossRef]
  28. Abu-Hijleh, B.; Al-Nimr, M. The Effect of the Local Inertial Term on the Fluid Flow in Channels Partially Filled with Porous Material. Int. J. Heat Mass Transf. 2001, 44, 1565–1572. [Google Scholar] [CrossRef]
  29. Turan, O.; Sachdeva, A.; Chakraborty, N.; Poole, R.J. Laminar natural convection of power-law fluids in a square enclosure with differentially heated side walls subjected to constant temperatures. J. Non-Newton. Fluid Mech. 2011, 166, 1049–1063. [Google Scholar] [CrossRef]
  30. Kefayati, G.; Tang, H.; Chan, A.; Wang, X. A lattice Boltzmann model for thermal non-Newtonian fluid flows through porous media. Comput. Fluids 2018, 176, 226–244. [Google Scholar] [CrossRef]
  31. Akbal, S.; Baytaş, F. Effects of non-uniform porosity on double diffusive natural convection in a porous cavity with partially permeable wall. Int. J. Therm. Sci. 2008, 47, 875. [Google Scholar] [CrossRef]
  32. Mchirgui, A.; Hidouri, N.; Magherbi, M.; Ben Brahim, A. Second law analysis in double diffusive convection through an inclined porous cavity. Comput. Fluids 2014, 96, 105–115. [Google Scholar] [CrossRef]
  33. Bera, P.; Pippal, S.; Sharma, A.K. A thermal non-equilibrium approach on double-diffusive natural convection in a square porous-medium cavity. Int. J. Heat Mass Transf. 2014, 78, 1080–1094. [Google Scholar] [CrossRef]
  34. Zhu, Q.; Zhuang, Y.; Yu, H. Three-dimensional numerical investigation on thermosolutal convection of power-law fluids in anisotropic porous media. Int. J. Heat Mass Transf. 2017, 104, 897–917. [Google Scholar] [CrossRef]
  35. Kefayati, G. Simulation of double diffusive natural convection and entropy generation of power-law fluids in an inclined porous cavity with Soret and Dufour effects (Part I: Study of fluid flow, heat and mass transfer). Int. J. Heat Mass Transf. 2016, 94, 539–581. [Google Scholar] [CrossRef]
  36. Kefayati, G. Simulation of double diffusive natural convection and entropy generation of power-law fluids in an inclined porous cavity with Soret and Dufour effects (Part II: Entropy generation). Int. J. Heat Mass Transf. 2016, 94, 582–624. [Google Scholar] [CrossRef]
  37. Hussain, S.; Schieweck, F.; Turek, S. Effcient Newton multigrid solution techniques for higher order space time Galerkin discretizations of incompressible flow. Appl. Numer. Math. 2014, 83, 51–71. [Google Scholar] [CrossRef]
  38. SHussain, S.; Mehmood, K.; Sagheer, M.; Farooq, A. Entropy generation analysis of mixed convective flow in an inclined channel with cavity with Al2O3-water nanofluid in porous medium. Int. Commun. Heat Mass Transf. 2017, 89, 198–210. [Google Scholar] [CrossRef]
  39. Gibanov, N.S.; Sheremet, M.A.; Ismael, M.A.; Chamkha, A.J. Mixed convection in a ventilated cavity filled with a triangular porous layer. Transp. Porous Media 2017, 89, 198–210. [Google Scholar] [CrossRef]
  40. Corvaro, F.; Paroncini, M. Experimental analysis of natural convection in square cavities heated from below with 2D-PIV and holographic interferometry techniques. Exp. Therm. Fluid Sci. 2007, 31, 721–739. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram.
Figure 1. Schematic diagram.
Mathematics 10 01060 g001
Figure 2. The grid for the proposed model.
Figure 2. The grid for the proposed model.
Mathematics 10 01060 g002
Figure 3. Comparison of the flow structure and temperature field with those of Gibanov et al. [39]. Reprinted with permission from ref. [39]. Copyright 2017 Springer Nature.
Figure 3. Comparison of the flow structure and temperature field with those of Gibanov et al. [39]. Reprinted with permission from ref. [39]. Copyright 2017 Springer Nature.
Mathematics 10 01060 g003
Figure 4. Comparison with the results of Corvaro and Paroncini [40] (Reprinted with permission from ref. [40]. Copyright 2007 Elsevier): flow field (experimental in (a)), isotherms (numerical in (b)), and isotherms (experimental in (c)). The corresponding present numerical study is given in (df).
Figure 4. Comparison with the results of Corvaro and Paroncini [40] (Reprinted with permission from ref. [40]. Copyright 2007 Elsevier): flow field (experimental in (a)), isotherms (numerical in (b)), and isotherms (experimental in (c)). The corresponding present numerical study is given in (df).
Mathematics 10 01060 g004
Figure 5. Streamlines (column 1), isotherms (column 2), and isoconcentration (column 3).
Figure 5. Streamlines (column 1), isotherms (column 2), and isoconcentration (column 3).
Mathematics 10 01060 g005
Figure 6. Variation in local Nu and Sh.
Figure 6. Variation in local Nu and Sh.
Mathematics 10 01060 g006
Figure 7. Streamlines (column 1), isotherms (column 2), and isoconcentration (column 3).
Figure 7. Streamlines (column 1), isotherms (column 2), and isoconcentration (column 3).
Mathematics 10 01060 g007
Figure 8. Variation in local Nu and Sh.
Figure 8. Variation in local Nu and Sh.
Mathematics 10 01060 g008
Figure 9. Streamlines (column 1), isotherms (column 2), and isoconcentration (column 3).
Figure 9. Streamlines (column 1), isotherms (column 2), and isoconcentration (column 3).
Mathematics 10 01060 g009
Figure 10. Variation in local Nu and Sh.
Figure 10. Variation in local Nu and Sh.
Mathematics 10 01060 g010
Figure 11. Streamlines (column 1), isotherms (column 2), and isoconcentration (column 3).
Figure 11. Streamlines (column 1), isotherms (column 2), and isoconcentration (column 3).
Mathematics 10 01060 g011
Figure 12. Variation in local Nu and Sh.
Figure 12. Variation in local Nu and Sh.
Mathematics 10 01060 g012
Figure 13. Streamlines (column 1), isotherms (column 2), and isoconcentration (column 3).
Figure 13. Streamlines (column 1), isotherms (column 2), and isoconcentration (column 3).
Mathematics 10 01060 g013
Figure 14. Variation in local Nu and Sh.
Figure 14. Variation in local Nu and Sh.
Mathematics 10 01060 g014
Figure 15. Impact of average Nu and Sh.
Figure 15. Impact of average Nu and Sh.
Mathematics 10 01060 g015
Figure 16. Impact of average Nu and Sh.
Figure 16. Impact of average Nu and Sh.
Mathematics 10 01060 g016
Figure 17. Impact of average Nu and Sh.
Figure 17. Impact of average Nu and Sh.
Mathematics 10 01060 g017
Table 1. Convergence of grid analysis.
Table 1. Convergence of grid analysis.
GridNELDOFs N u a v g S h a v g
G123033223.4283555.767679
G235249743.4180985.727896
G354073193.4201595.715686
G41006130073.4222585.703224
G51520190403.4240145.701316
G62470297573.4246015.700223
G76544769333.4246135.698674
G8167901921483.4249155.698709
Table 2. Validation of the code of average Nu with [39]. Adapted with permission from ref. [39]. Copyright 2017 Springer Nature.
Table 2. Validation of the code of average Nu with [39]. Adapted with permission from ref. [39]. Copyright 2017 Springer Nature.
Present Study[39]Error
12.5787912.53350.36%
Present ModelGibanov et al. [39].
Table 3. Average Nu and Sh comparison.
Table 3. Average Nu and Sh comparison.
R a   =   10 3 R a   =   10 4 R a   =   10 5 R a   =   10 6
N u a v g 1.0035441.2941513.4250278.530151
S h a v g 1.0218812.0175485.69900213.07502
Table 4. Average Nu and Sh comparison.
Table 4. Average Nu and Sh comparison.
n   =   0.6 n   =   0.8 n   =   1 n   =   1.2 n   =   1.4
N u a v g 3.4249153.3385053.2728343.2218583.181187
S h a v g 5.6987095.5193095.3834655.2782735.194299
Table 5. Average Nu and Sh comparison.
Table 5. Average Nu and Sh comparison.
D a   =   10 6 D a   =   10 5 D a   =   10 4 D a   =   10 3 D a   =   10 2
N u a v g 1.35561.3737421.6528943.4251394.703467
S h a v g 1.5378051.5828512.509855.6993067.181993
Table 6. Average Nu and Sh comparison.
Table 6. Average Nu and Sh comparison.
N = 0.1N = 1N = 5N = 10
N u a v g 3.4249154.1029395.8312497.130527
S h a v g 5.6987097.08079110.4767512.8777
Table 7. Average Nu and Sh comparison.
Table 7. Average Nu and Sh comparison.
Le = 1Le = 2.5 Le = 5Le = 10
N u a v g 3.4837563.4251363.3899143.366655
S h a v g 3.483955.6993837.70360910.18238
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kolsi, L.; Hussain, S.; Ghachem, K.; Jamal, M.; Maatki, C. Double Diffusive Natural Convection in a Square Cavity Filled with a Porous Media and a Power Law Fluid Separated by a Wavy Interface. Mathematics 2022, 10, 1060. https://doi.org/10.3390/math10071060

AMA Style

Kolsi L, Hussain S, Ghachem K, Jamal M, Maatki C. Double Diffusive Natural Convection in a Square Cavity Filled with a Porous Media and a Power Law Fluid Separated by a Wavy Interface. Mathematics. 2022; 10(7):1060. https://doi.org/10.3390/math10071060

Chicago/Turabian Style

Kolsi, Lioua, Shafqat Hussain, Kaouther Ghachem, Muhammad Jamal, and Chemseddine Maatki. 2022. "Double Diffusive Natural Convection in a Square Cavity Filled with a Porous Media and a Power Law Fluid Separated by a Wavy Interface" Mathematics 10, no. 7: 1060. https://doi.org/10.3390/math10071060

APA Style

Kolsi, L., Hussain, S., Ghachem, K., Jamal, M., & Maatki, C. (2022). Double Diffusive Natural Convection in a Square Cavity Filled with a Porous Media and a Power Law Fluid Separated by a Wavy Interface. Mathematics, 10(7), 1060. https://doi.org/10.3390/math10071060

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop