Finite element simulations on heatline trajectories for mixed convection in porous square enclosures: Effects of various moving walls

•Mixed convection within porous square cavity involving various moving walls is relevant in many energy related systems.•Galerkin finite element method is used to solve the non-linear partial differential equations.•Energy flow via thermal management is analyzed based on heatline concept.•Effect of Darcy number, Reynolds number, Grashof number and Prandtl number on the trends of flow fields is examined.•Heat flow rates are found to be highly influenced by various motion of walls.

Finite element simulation of the mixed convection within porous square cavities for Darcy–Brinkman–Forchheimer model has been carried out in the present work. The penalty optimization based Galerkin finite element method is used to solve the partial differential equations of heat and fluid flow. Bejan’s heatline concept has been employed to visualize the heat flow within the closed cavities based on the motion of the horizontal wall(s) (cases 1a–1d) or vertical wall(s) (cases 2a–2c) involving isothermally hot bottom wall, cold side walls and insulated top wall for various fluids with Prandtl number, Prm=0.026Prm=0.026, 0.70.7 and 7.27.2, Reynolds number, Re=10–100Re=10–100 and Grashof number, Gr=103–105Gr=103–105. The higher permeability at Dam≥10−3Dam≥10−3 leads to the enhanced buoyancy convection for all the cases. Although the direction of the motion of wall(s) significantly influences the fluid flow field within the enclosure, due to the decoupling between the fluid and thermal fields at the low PemPem (Pem=0.26Pem=0.26 and 2.62.6), conductive heat transfer occurs as seen from the end-to-end heatlines. It is also found that the overall heat transfer rates at the bottom wall (Nub¯) are identical for the cases 1a–1d and cases 2a–2c at Prm=0.026Prm=0.026, irrespective of GrGr and ReRe at Dam=10−2Dam=10−2. At Prm=0.7Prm=0.7 and 7.27.2, the convection dominant heat transfer occurs for all the cases for Gr=105Gr=105, Re=10Re=10 and 100100 and Dam=10−2Dam=10−2. The strong convective circulation cells are observed at Prm=0.7Prm=0.7 and 7.27.2 for all the cases. The plume shaped isotherms are also observed along the centerline at Prm=7.2Prm=7.2, Re=10Re=10, Gr=105Gr=105 and Dam=10−2Dam=10−2 for all the cases. At Prm=7.2Prm=7.2, Gr=105Gr=105, Dam=10−2Dam=10−2 and Re=100Re=100, the multiple convective heatline cells are observed for the cases 1a–1d. It is observed that, the strengths of fluid and heat circulation cells are less at Re=100Re=100 compared to Re=10Re=10 for all the cases due to weak buoyancy force at the high ReRe. In order to achieve the high heat transfer rate at the bottom wall (Nub¯) for the mixed convection involving various moving walls, case 2b (a case of the vertically moving wall) is preferred at Prm=0.7Prm=0.7, Re=100Re=100, Gr=105Gr=105 and Dam=10−2Dam=10−2. At the high PrmPrm (Prm=7.2Prm=7.2, Re=100Re=100, Gr=105Gr=105 and Dam=10−2Dam=10−2), case 2a (a case of the vertically moving wall) is preferred based on the maximum heat transfer rate at the bottom wall (Nub¯).