Analysis of Energy Flux Vector on Natural Convection Heat Transfer in Porous Wavy-Wall Square Cavity with Partially-Heated Surface

Abstract

The study utilizes the energy-flux-vector method to analyze the heat transfer characteristics of natural convection in a wavy-wall porous square cavity with a partially-heated bottom surface. The effects of the modified Darcy number, modified Rayleigh number, modified Prandtl number, and length of the partially-heated bottom surface on the energy-flux-vector distribution and mean Nusselt number are examined. The results show that when a low modified Darcy number with any value of modified Rayleigh number is given, the recirculation regions are not formed in the energy-flux-vector distribution within the porous cavity. Therefore, a low mean Nusselt number is presented. The recirculation regions do still not form, and thus the mean Nusselt number has a low value when a low modified Darcy number with a high modified Rayleigh number is given. However, when the values of the modified Darcy number and modified Rayleigh number are high, the energy flux vectors generate recirculation regions, and thus a high mean Nusselt number is obtained. In addition, in a convection-dominated region, the mean Nusselt number increases with an increasing modified Prandtl number. Furthermore, as the length of the partially-heated bottom surface lengthens, a higher mean Nusselt number is presented.

1. Introduction

The plot of the heat flow paths is important since it can provide physical insights into the energy transport process in detail. To achieve this purpose, Kimura and Bejan [1] have suggested a heatline visualization technique. Following the study of Kimura and Bejan [1], numerous researchers have explored the process of heat transport within thermal-fluid systems by utilizing the technique [2,3,4].

Hooman [5,6] has presented an energy-flux-vector method, which is basically similar to the heatline technique, for visualizing the heat flow paths. Comparing the two visualization methods, Hooman [5,6] has pointed out that the energy-flux-vector method is simpler than the heatline visualization technique since the algebraic equations do not require solving. In the literature, numerous researchers [7,8,9] have demonstrated that the process of heat transport can be completely explained by using the energy-flux-vector method.

Natural convection in porous cavities has numerous practical applications in engineering fields, including biomedical engineering, chemical and material processing, fluidized beds, geothermal engineering, thermal insulation, solar collection, and so on [10,11,12,13,14]. To achieve the purpose of heat transfer enhancement, wavy-surface geometries are often imposed [7,8]. Note that the practical applications for the natural convection in a wavy-wall cavity include heat exchangers, solar collectors, condensers in refrigerators, geothermal engineering, and so on [15,16,17]. In the literature, the investigation into the natural convection heat transfer behavior in a porous cavity with wavy surfaces has been widely discussed [18,19,20]. The results have demonstrated that the use of the wavy-surface geometries can improve the heat transfer effect. Recently, Biswal et al. [21] have utilized the energy-flux-vector method to explain the process of heat transport of natural convection within a porous cavity with curved side walls. Their results have shown that given suitable curved-side wall forms with appropriate flow conditions, the heat transfer performance can be enhanced.

As discussed above, the analysis of the energy-flux-vector method on natural convection heat transfer in a porous cavity with wavy surfaces has attracted little attention. Accordingly, in the current study, the energy-flux-vector method is utilized to analyze the heat transfer characteristics of natural convection within a porous square cavity bounded by a partially-heated flat bottom wall, low temperature left and right wavy-walls, and an insulated flat top wall. The simulations focus particularly on the effects of the modified Darcy number, modified Rayleigh number, modified Prandtl number, and length of the partially-heated bottom surface on the energy-flux-vector distribution and mean Nusselt number, respectively. 2. Mathematical Formulation 2.1. Governing Equations and Boundary Conditions

Figure 1 illustrates the studied porous square wavy-wall cavity with a partially-heated bottom surface and a characteristic length of_Lc. As shown, the partially-heated wall surface has a length of_LH, and it is placed on the center of the bottom wall. Meanwhile, the left and right walls are assumed to have a constant low temperature and a complex-wavy surface.

It is assumed that the working fluid is Newtonian and incompressible, and the flow and temperature fields are two-dimensional and steady state. In addition, it is also assumed that the Rayleigh number is less than¹⁰⁹ and thus, the assumption of laminar flow is valid. Furthermore, the Boussinesq approximation [7,8,22] is imposed, and the local thermal equilibrium condition is assumed to be achieved.

Defining the non-dimensional quantities are as follows [14]:

^x*=x_Lc, ^y*=y_Lc, _αeff=_keffερ_Cp, ^u*=u_Lcαeff, ^v*=v_Lcαeff, ^p*=p_Lc2ρ_αeff2, θ=T−_TLTH−_TL, _Κm=Kε, _Prm=μρ_αeff, D_am=_KmLc2, R_am=ρgβ_Lc3(_TH−_TL)μ_αeff

whereρ,_Cp, andprepresent the density, specific heat, and pressure, respectively;uandvare the velocity components along x- and y-axes, respectively;_THand_TLsignify the high temperature and low temperature, respectively; subscriptmindicates the modified value; and superscript∗ denotes the non-dimensional quantity. After the viscous dissipation and thermal radiation effects are ignored, the continuity, momentum and energy conservation equations described the flow behavior and heat transfer characteristics of natural convection within the porous cavity is written in the following non-dimensionalized forms [14]:

Continuity equation:

∂^u*∂^x*+∂^v*∂^y*=0

x-direction momentum equation:

^u*∂^u*∂^x*+^v*∂^u*∂^y*=−∂^p*∂^x*−_PrmD_am^u*+_Prm(^{∂2 u*}∂^x*2+^{∂2 u*}∂^y*2)−1.75150^u*2+^v*2D_am^u*

y-direction momentum equation:

^u*∂^v*∂^x*+^v*∂^v*∂^y*=−∂^p*∂^y*−_PrmD_am^v*+_Prm(^{∂2 v*}∂^x*2+^{∂2 v*}∂^y*2)−1.75150^u*2+^v*2D_am^v*+R_am⋅_Prm⋅θ

Energy conservation equation:

^u*∂θ∂^x*+^v*∂θ∂^y*=^∂2θ∂^x*2+^∂2θ∂^y*2

The dimensionless boundary conditions are expressed as follows:

Bottom partially-heat wall:		u*=v*=0,	θ=1
Bottom other walls:		u*=v*=0,	∂θ/∂n⇀*=0
Left and right wavy walls:		u*=v*=0,	θ=0
Top wall:		u*=v*=0,	∂θ/∂n⇀*=0

where^n⇀*is the normal vector.

2.2. Energy Flux Vectors and Nusselt Number

The energy flux vectors (E⇀ ) are defined as follows [5,6,7,8,9]:

E⇀=∂^H*∂^y*i⇀−∂^H*∂^x*j⇀

∂^H*∂^y*=^u*θ−∂θ∂^x*,

−∂^H*∂^x*=^v*θ−∂θ∂^y*.

Note that^H*is the non-dimensional heat function, andi⇀andj⇀are the unit components in x- and y-directions, respectively.

The mean Nusselt number (N_um) along the partially-heated bottom surface is defined as follows:

N_um=∫Nudξ

whereNu(=h_Lckeff=−(∂θ∂^n⇀*))is the Nusselt number, andhis the convection heat transfer coefficient.

2.3. Geometric Description and Numerical Method

In this study, the left and right walls of the porous cavity are assumed to have a complex-wavy surface, which is formulated as follows [7,8]:

^y*=_αwSin(2π^x*)+_αw2.5Sin(4π^x*)

where_αwis the amplitude of the wavy surface.

In the current study, the generalized coordinate transform technique, finite volume method, and SIMPLE algorithm are used to solve the governing equations presented in Equations (2)–(5). The numerical methods are identical to those used in our previous studies [7,8]. A detailed description of the numerical methods applied in this work can be found in [7,8].

2.4. Numerical Validation and Grid Independence Evaluation

The currently used numerical models and numerical methods were valid by comparing the current results with those presented by Singh et al. [14]. Table 1 shows these results. Note that the Rayleigh number ofRa=¹⁰⁵, porosity ofε=0.6, and Prandtl number ofPr=1 are given in the case. It is shown that the current results are identical to those presented in Singh et al. [14].

The mesh sizes of101×201,101×501,101×1001,201×1001have been examined for the variation of mean Nusselt number, respectively, and the results showed that the mesh size of101×1001has a grid-independent solution.

3. Results and Discussion

In the study, the ranges of the non-dimensional parameters were set as follows: modified Darcy number:D_am=¹⁰⁻⁵toD_am=¹⁰⁻²; modified Rayleigh number:R_am=¹⁰²toR_am=¹⁰⁵; modified Prandtl number:_Prm=0.1to_Prm=10; amplitude of wavy surface:_αw=0.25; and length of the partially-heated bottom surface:_LH*=0.3to_LH*=0.9.

Figure 2 and Figure 3 illustrate the distributions of the energy flux vectors and isotherms within the porous cavity for various modified Darcy numbers and modified Rayleigh numbers. Figure 4 shows the effects of the modified Darcy number and modified Rayleigh number on the mean Nusselt number. According to the definition of the energy flux vectors given in Equation (6), the conduction mechanism dominates if the flow of energy flux vectors is directly from the bottom partially-heated high-temperature wall to the left and right low-temperature wavy-walls, while the convection mechanism dominates if the energy flux vectors generate closed recirculation regions [5,6,7,8,9]. Therefore, as shown in Figure 2a, when the values of the modified Darcy number and modified Rayleigh number are both low, the heat transfer effect within the porous partially-heated cavity is dominated by a pure conduction mechanism since the closed recirculation regions in the energy-flux-vector distribution are not created. In addition, a low modified Darcy number represents a low permeability due to the presence of a high flow resistance [14]. Therefore, although a larger buoyancy effect induced by giving a high modified Rayleigh number is presented, a high flow resistance is also generated within the porous partially-heated cavity. As a result, the energy flux vectors do still not form closed recirculation regions, and thus the conduction mechanism continually dominates heat transfer behavior (see Figure 2b). Under conduction-domination, the high-temperature fluid heated by the partially-heated bottom surface is slowly diffused to left and right low-temperature wavy-wall surfaces to be dissipated (see the isotherms in Figure 2). Consequently, under a low modified Darcy number, a low mean Nusselt number is presented, irrespective of the value assigned to the modified Rayleigh number (see Figure 4).

When a high modified Darcy number and a low modified Rayleigh number is given, the closed recirculation regions of the energy flux vectors are still not formed within the porous partially-heated cavity (see Figure 3a). Although a high modified Darcy number has a high permeability and a low flow resistance, a low modified Rayleigh number induces a small buoyancy effect. Therefore, the convection effect is still weak, and the conduction mechanism continues to dominate the heat transfer effect. As a result, a low mean Nusselt number is presented (see Figure 4).

When the modified Darcy number and the modified Rayleigh number are both high, a low flow resistance and a high buoyancy effect occur within the porous partially-heated cavity. Therefore, the closed energy-flux-vector recirculations are formed, resulting in enhancement of the convection effect (see Figure 3b). The high-temperature fluid on the partially-heated bottom wall is promptly driven to the left and right low-temperature walls to be dissipated (see the isotherms in Figure 3b). As a result, a high mean Nusselt number is obtained (see Figure 4).

Figure 5 and Figure 6 illustrate the distributions of the energy flux vectors and isotherms within the porous partially-heated cavity for various modified Prandtl numbers and modified Rayleigh numbers. Figure 7 shows the effects of the modified Prandtl number and modified Rayleigh number on the mean Nusselt number. Note that the modified Darcy number is set asD_am=¹⁰⁻² in the cases. It shows that given a low modified Rayleigh number, the distributions of the energy flux vectors and isotherms are similar for various modified Prandtl numbers (see Figure 5a,b). The recirculation regions of energy flux vectors are not created, and thus the conduction heat transfer dominates. In other words, the effect of the modified Prandtl number on the heat transfer is insignificant. Therefore, a low and approximately constant mean Nusselt number is presented (see Figure 7). When a high modified Rayleigh number is given, it shows that the size of closed energy-flux-vector recirculation regions enlarges as the modified Prandtl number is increased (see Figure 6a,b). In other words, a higher modified Prandtl number has a stronger convection effect under high modified Rayleigh numbers. Consequently, a higher mean Nusselt number is presented (see Figure 7).

Figure 8 illustrates the distributions of the energy flux vectors within the porous cavity for various partially-heated lengths and modified Rayleigh numbers. Figure 9 shows the effects of the partially-heated length and modified Rayleigh number on the mean Nusselt number. Note that the modified Darcy number is set asD_am=¹⁰⁻² in these cases. For a low modified Rayleigh number, the conduction heat transfer dominates since the recirculation regions of the energy flux vectors are not formed within the porous cavity for various partially-heated lengths (see Figure 8a). Therefore, a low mean Nusselt number is presented. As the length of the partially-heated bottom surface is lengthened, a larger bottom heated area is presented. Consequently, a higher mean Nusselt number is obtained (see Figure 9). Given a high modified Rayleigh number, the energy flux vectors form closed recirculation regions. The size of the recirculation regions enlarges as the partially-heated length is extended (see Figure 8b). As a result, the strength of the convection effect enhances. Consequently, the mean Nusselt number rises with the increasing length of the partially-heated bottom surface (see Figure 9).

4. Conclusions This paper has analyzed the heat transfer behavior of natural convection in a porous square wavy-wall cavity with a partially-heated bottom surface using the energy-flux-vector method. The effects of the flow parameters and length of a partially-heated bottom surface on the energy-flux-vector distribution and mean Nusselt number have been discussed. The studied results are summarized as follows:

• Given a low modified Darcy number, the energy flux vectors did not generate recirculation regions within the porous cavity, irrespective of the value assigned to the modified Rayleigh number. The conduction heat transfer dominated, and thus the mean Nusselt number was low.
• Given a high modified Darcy number and a low modified Rayleigh number, the heat transfer effect was dominated by the conduction mechanism since no recirculation region was formed in the energy flux vectors. Therefore, a low mean Nusselt number was obtained.
• Given a high modified Darcy number with a high modified Rayleigh number, recirculation regions in the energy-flux-vector distribution were produced, resulting in a convection-domination. Consequently, a high mean Nusselt number was presented.
• In conduction-dominated region, the effect of the modified Prandtl number on the energy-flux-vector distribution and mean Nusselt number was insignificant. However, in convection-dominated region, the size of the closed energy-flux-vector recirculation region enlarged, and the value of the mean Nusselt number raised as the modified Prandtl number was increased.
• In conduction-dominated or convection-dominated regions, the mean Nusselt number was raised as the length of the partially-heated bottom surface lengthened.

Figure 1. Geometry of partially-heated wavy-wall porous cavity.

Figure 2. Distributions of energy flux vectors and isotherms within a partially-heated porous cavity given a modified Darcy number ofDam=10-5and modified Rayleigh numbers of: (a)Ram=102and (b)Ram=105. Note thatPrm=1andLH*=0.3. Note also that the black solid lines indicate the flow direction of the energy-flux-vector, while the red dashed lines indicate the isothermal contours.

Figure 3. Distributions of energy flux vectors and isotherms within a partially-heated porous cavity given a modified Darcy number ofDam=10-2and modified Rayleigh numbers of: (a)Ram=102and (b)Ram=105. Note thatPrm=1andLH*=0.3. Note also that the black solid lines indicate the flow direction of the energy-flux-vector, while the red dashed lines indicate the isothermal contours.

Figure 4. Variation of mean Nusselt number with modified Darcy number as a function of modified Rayleigh number. Note thatPrm=1andLH*=0.3.

Figure 5. Distributions of energy flux vectors and isotherms within a partially-heated porous cavity given a modified Rayleigh number ofRam=102and modified Prandtl numbers of: (a)Prm=0.1and (b)Prm=10. Note thatDam=10-2andLH*=0.3. Note also that the black solid lines indicate the flow direction of energy flux vectors, while the red dashed lines indicate the isothermal contours.

Figure 6. Distributions of energy flux vectors and isotherms within a partially-heated porous cavity given a modified Rayleigh number ofRam=105and modified Prandtl numbers of: (a)Prm=0.1and (b)Prm=10. Note thatDam=10-2andLH*=0.3. Note also that the black solid lines indicate the flow direction of energy flux vectors, while the red dashed lines indicate the isothermal contours.

Figure 7. Variation of mean Nusselt number with modified Prandtl number as a function of a modified Rayleigh number. Note thatDam=10-2andLH*=0.3.

Figure 8. Distributions of energy flux vectors within a porous cavity given partially-heated lengths ofLH*=0.3andLH*=0.9and modified Rayleigh numbers of: (a)Ram=102and (b)Ram=105. Note thatDam=10-2andPrm=1. Note also that the red solid lines indicate the partially-heated length ofLH*=0.9, and the green dashed lines indicate the partially-heated length ofLH*=0.3.

Figure 9. Variation of mean Nusselt number with partially-heated length as a function of a modified Rayleigh number. Note thatDam=10-2andPrm=1.

	Da=10−2	Da=10−4
Current results	3.441	1.067
Singh et al. [14]	3.461	1.067
Error (%)	0.6	0.0

Author Contributions

Conceptualization, C.-C.C.; Methodology, C.-C.C.; Software, C.-C.C. and Y.-T.L.; Validation, C.-C.C. and Y.-T.L.; Formal Analysis, C.-C.C.; Investigation, C.-C.C.; Resources, C.-C.C. and Y.-T.L.; Data Curation, C.-C.C. and Y.-T.L.; Writing-Original Draft Preparation, C.-C.C.; Writing-Review & Editing, C.-C.C.; Visualization, C.-C.C.; Supervision, C.-C.C.; Project Administration, C.-C.C.; Funding Acquisition, C.-C.C. and Y.-T.L.

Funding

The authors would like to thank the Ministry of Science and Technology, Taiwan, for the financial support of this study under Contract Nos. MOST 108-2221-E-150-011, MOST 108-3116-F-042A-006 and MOST 107-2221-E-150-035.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Dam	modified Darcy number,Dam=KmLc2
E⇀	energy flux vector
keff	effective thermal conductivity,Wm-1 K-1
Lc	characteristic length of square cavity,m
LH	length of partially-heated bottom surface,m
Nu	Nusselt number,Nu=hLckeff=-(∂θ∂n⇀*)
Num	mean Nusselt number,Num=∫Nudξ
Prm	modified Prandtl number,Prm=μραeff
Ram	modified Rayleigh number,Ram=ρgβLc3(TH-TL)μαeff
T	temperature,K
Greek symbols
αeff	effective thermal diffusivity,m2 s-1
αw	amplitude of wavy surface
β	thermal expansion coefficient,K-1
ε	porosity
Κ	permeability,m2
θ	dimensionless temperature,θ=T-TLTH-TL

References

1. Kimura, S.; Bejan, A. The "Heatline" visualization of convective heat transfer. J. Heat Transf.-Trans. ASME 1983, 105, 916-919.

2. Trevisan, O.V.; Bejan, A. Combined heat and mass transfer by natural convection in a vertical enclosure. J. Heat Transf.-Trans. ASME 1984, 109, 104-112.

3. Ramakrishna, D.; Basak, T.; Roy, S.; Pop, I. A complete heatline analysis on mixed convection within a square cavity: Effects of thermal boundary conditions via thermal aspect ratio. Int. J. Therm. Sci. 2012, 57, 98-111.

4. Biswal, P.; Nag, A.; Basak, T. Analysis of thermal management during natural convection within porous tilted square cavities via heatline and entropy generation. Int. J. Mech. Sci. 2016, 115-116, 596-615.

5. Hooman, K.; Gurgenci, H.; Dincer, I. Heatline and energy-flux-vector visualization of natural convection in a porous cavity occupied by a fluid with temperature-dependent viscosity. J. Porous Media 2009, 12, 265-275.

6. Hooman, H. Energy flux vectors as a new tool for convection visualization. Int. J. Numer. Methods Heat Fluid Flow 2010, 20, 240-249.

7. Cho, C.C. Heat transfer and entropy generation of mixed convection flow in Cu-water nanofluid-filled lid-driven cavity with wavy surface. Int. J. Heat Mass Transf. 2018, 119, 163-174.

8. Cho, C.C. Mixed convection heat transfer and entropy generation of Cu-water nanofluid in wavy-wall lid-driven cavity in presence of inclined magnetic field. Int. J. Mech. Sci. 2019, 151, 703-714.

9. Nayak, R.K.; Bhattacharyya, S.; Pop, I. Numerical study on mixed convection and entropy generation of a nanofluid in a lid-driven square enclosure. J. Heat Transf.-Trans. ASME 2016, 138, 012503.

10. Lu, D.A.; Flamant, G.; Snabre, P. Towards a generalized model for vertical walls to gas-solid fluidized beds heat transfer-I. Particle convection and gas convection. Chem. Eng. Sci. 1993, 48, 2479-2492.

11. Nield, D.A.; Bejan, A. Convection in Porous Media; Springer: New York, NY, USA, 2006.

12. Prud'homme, M.; Jasmin, S. Inverse solution for a biochemical heat source in a porous medium in the presence of natural convection. Chem. Eng. Sci. 2006, 61, 1667-1675.

13. Al-Amiri, A.; Khanafer, K.; Pop, I. Steady-state conjugate natural convection in a fluid-saturated porous cavity. Int. J. Heat Mass Transf. 2008, 51, 4260-4275.

14. Singh, A.K.; Basak, T.; Nag, A.; Roy, S. Heatlines and thermal management analysis for natural convection within inclined porous square cavities. Int. J. Heat Mass Transf. 2015, 87, 583-597.

15. Misirlioglu, A.; Baytas, A.C.; Pop, I. Free convection in a wavy cavity filled with a porous medium. Int. J. Heat Mass Transf. 2005, 48, 1840-1850.

16. Sultana, Z.; Hyder, M.N. Non-darcy free convection inside a wavy enclosure. Int. Commun. Heat Mass Transf. 2007, 34, 136-146.

17. Cho, C.C.; Chiu, C.H.; Lai, C.Y. Natural convection and entropy generation of Al2O3-water nanofluid in an inclined wavy-wall cavity. Int. J. Heat Mass Transf. 2016, 97, 511-520.

18. Kumar, B.V.R.; Singh, P.; Murthy, P.V.S.N. Effect of surface undulations on natural convection in a porous square cavity. J. Heat Transf.-Trans. ASME 1997, 119, 848-851.

19. Chen, X.B.; Yu, P.; Winoto, S.H.; Low, H.T. Free convection in a porous wavy cavity based on the Darcy-Brinkman-Forchheimer extended model. Numer. Heat Tranf. A-Appl. 2007, 52, 377-397.

20. Khanafer, K.; Al-Azmi, B.; Marafie, A.; Pop, I. Non-Darcian effects on natural convection heat transfer in a wavy porous enclosure. Int. J. Heat Mass Transf. 2009, 52, 1887-1896.

21. Biswal, P.; Basak, T. Heatlines visualization of convective heat flow during differential heating of porous enclosures with concave/convex side walls. Int. J. Numer. Methods Heat Fluid Flow 2018, 28, 1506-1538.

22. Singh, A.K.; Basak, T.; Nag, A.; Roy, S. Role of entropy generation on thermal management during natural convection in tilted porous square cavities. J. Taiwan Inst. Chem. Eng. 2015, 50, 153-172.