Rayleigh–Bénard convection of a viscoplastic liquid in a trapezoidal enclosure

• Published in: International journal of mechanical sciences 2020-08, Vol.180, p.105630, Article 105630
• Main Authors: Aghighi, M.S., Ammar, A., Masoumi, H., Lanjabi, A.
• Format: Article
• Language: English
• Subjects: Bingham model; Engineering Sciences; Finite element method; Rayleigh–Bénard; Trapezoidal enclosure; Yield stress
• ISSN: 0020-7403; 1879-2162
• DOI: 10.1016/j.ijmecsci.2020.105630

•Natural convection of viscoplastic fluids in a trapezoidal cavity is studied.•The effects of the sloping walls on the flow and temperature fields and morphology of yielded/unyielded regions are presented.•Both the heat transfer rate and plug regions increase with increasing the inclination angle of the side wall.•The critical yield number Yc increases significantly with the inclination angle.•The inclination angle reduces the effect of Rayleigh number on critical yield stress. The objective of this paper is to clarify the role of sloping walls on convective heat transport in Rayleigh–Bénard convection within a trapezoidal enclosure filled with viscoplastic fluid. The rheology of the viscoplastic fluid has been modeled with Bingham fluid model. The system of coupled nonlinear differential equations was solved numerically by Galerkin's weighted residuals scheme of finite element method. The numerical experiments are carried out for a range of parameter values, namely, Rayleigh number (5.103 ≤ Ra ≤ 105), yield number (0 ≤ Y ≤ Yc), and sidewall inclination angle (ϕ = 0,  π/6,  π/4,  π/3) at a fixed Prandtl number (Pr = 500). Effects of the inclination angle on the flow and temperature fields are presented. The results reveal that inclination angle causes a multicellular flow and appears as the main parameter to govern heat transfer in the cavity. The heat transfer rate is found to increase with the increasing angle of the sloping wall for both Newtonian and yield stress fluids. On the other hand, the plug regions also found to increase with increasing φ, which is unusual but perhaps not unexpected behavior. In the yield stress fluids, the flow becomes motionless above a critical yield number Yc because the plug regions invade the whole cavity. The critical yield number Yc is also affected by the change of inclination angle and increases significantly with the increase of φ. [Display omitted]