Three-dimensional numerical simulations of free convection in a layered porous enclosure

Three-dimensional numerical simulations are carried out for the study of free convection in a layered porous enclosure heated from below and cooled from the top. The system is defined as a cubic porous enclosure comprising three layers, of which the external ones share constant physical properties and the internal layer is allowed to vary in both permeability and thermal conductivity. The model is based on Darcy's law and the Boussinesq approximation. A parametric study to evaluate the sensitivity of the Nusselt number to a decrease in the permeability of the internal layer shows that strong permeability contrasts are required to observe an appreciable drop in the Nusselt number. If additionally the thickness of the internal layer is increased, a further decrease in the Nusselt number is observed as long as the convective modes remain the same, if the convective modes change the Nusselt number may increase. Decreasing the thermal conductivity of the middle layer causes first an increment in the Nusselt number and then a drop. On the other hand, the Nusselt number decreases in an approximately linear trend when the thermal conductivity of the layer is increased.