Numerical investigation of natural convection in a square enclosure partially filled with horizontal layers of a porous medium

Two-dimensional buoyancy-induced flow and heat transfer inside a square enclosure partially occupied by copper metallic foam subjected to a symmetric side cooling and constant heat flux bottom heating was tested numerically. Finite Element Method was employed to solve the governing partial differential equations of the flow field and the Local Thermal Equilibrium model was used for the energy equation. The system boundaries were defined as lower heated wall by constant heat flux, cooled lateral walls, and insulated top wall. The three parameters elected to conduct the study are heater length (7 ≤ ζ ≤ 20 cm), constant heat flux (150 ≤ q″ ≤ 600 W m), and porous material thickness (5 ≤ H ≤ 20 cm). The porous material used was the copper metal foam of porosity = 0.9 and pore density PPI = 10, and saturated with a fluid of Prandtl number = 0.7. On the basis of the results obtained, it was concluded that at the porous layer thickness = 5 cm, the rate of heat transferred was (74.6%) higher than when the layer height was 20 cm (the cavity is fully filled) and at the same thickness it was found that the heat rate is (51.4%) higher than when using the half filling (H = 10 cm). Further, the local and mean Nusselt number is maximum when using the largest heater size and smallest porous layer thickness. Finally, better circulation and convective modes were observed at high values of heat flux.