Heat Flow Visualization for Natural Convection in Rhombic Enclosures: Bejan’s Heatline Approach

The phenomena of natural convection within a rhombic enclosure filled with air (Pr = 0.71) for (a) isothermally (case 1) and (b) non-isothermally (case 2) heated bottom walls with various aspect-ratios has been studied numerically. In all the cases, top horizontal wall is maintained adiabatic and side walls are maintained cold. Galerkin finite element method with penalty parameter is used to solve non-linear coupled partial differential equations for flow and temperature fields. Poisson equation of streamfunction and heatfunction is also solved using Galerkin method. Simulations are carried out over a range of Rayleigh numbers and numerical results are presented in terms of streamfunction, heatfunction and temperature contours. Streamlines are useful to visualize the fluid flow whereas heatlines are used to study the heat energy distribution within the rhombic cavity. Heatlines are further used to visualize the trajectories of heat flow and zones of high thermal mixing. At lower Ra, heatlines are smooth circular arcs with low magnitude streamfunctions and heatfunctions and thus the heat transfer is conduction dominant. Asymmetric flow is observed for all the cases due to geometrical asymmetry. As Ra increases, buoyant force starts dominating and the magnitudes of streamfunctions and heatfunctions are found to be greater due to enhanced convection effect. Heatlines are distorted greatly showing complex heat distribution inside the cavity. It is observed that primary heat circulation cell is larger for greater tilt angles and thus thermal mixing is high. Heat transfer rates are also studied via local and average Nusselt numbers as functions of Ra and Pr on bottom, left and right walls. Various quantitative and qualitative features of Nusselt numbers have also been explained based on heatlines.