3D DNS of laminar Rayleigh-Bénard convection in a cylinder for incompressible fluid flow
Pan-Xin Li,Ben-Wen Li,Lu Chen,Xue-Yuan Leng,Xiao-Hong Luo,Yang Yu,Jing-Kui Zhang
DOI: https://doi.org/10.1016/j.cjph.2022.07.011
IF: 3.957
2022-07-18
Chinese Journal of Physics
Abstract:A numerical investigation for three-dimensional laminar Rayleigh-Bénard (RB) convection in a cylinder with Γ=1 and 2, Pr=1 is conducted for incompressible fluid flow. For the original governing equations, two sets of dimensionless forms due to different definition of dimensionless variables are presented. The relation between these two sets is deduced theoretically and confirmed numerically. To solve these two sets of dimensionless governing equations (Set I & Set II), the Chebyshev-Fourier collocation spectral method (CSM) is developed for space discretization. For decoupling velocity and pressure, two methods, i.e., the improved projection scheme (IPS) and the artificial compressibility method (ACM), are adopted and compared for 2×103≤Ra≤2×104 at Γ=1 . Then, four group codes are designed corresponding to the combinations of equations (Set I & Set II) and methods (IPS & ACM), and compared. Results show that: all the four codes can capture the first critical Rayleigh number with 2% uncertainty; for the flow pattern, after the first critical Rayleigh number, two opposite flow patterns appear randomly for the same Set (Set I or Set II) but different methods; four heat transfer, correlations between the Nusselt number and the Rayleigh number are constructed, and the relative errors between them are within 5%; dimensionless governing equations Set II can provide higher accuracy than Set I dose, and the IPS can give higher efficiency than the ACM dose. Thereafter, the synthesis of Set II and the IPS is used to further investigate the RB convection for 2×104 at Γ=1 and 3×103≤Ra≤3×104 at Γ=2 . The results show that, the flow patterns are completely symmetric for both Γ=1 and 2. The increase of Γ can enhance heat transfer at a fixed Ra number. Besides, the correlations between the Nusselt number and Ra number for the specified range of Ra number at Γ=1 and 2 are given. Graphical abstract The laminar Rayleigh-Bénard (RB) convection in a cylinder is numerically investigated by DNS using collocation spectral method. We summarized two sets of dimensionless governing equations for the same physical problem and deduced their relations. Two combined numerical methods are developed and compared for the laminar RB convection, and concluded that, the improved projection scheme (IPS) can achieve better performance than the artificial compressibility method (ACM). The evolution of RB convection fluid flow can be observed in detail and the Bénard cells may appear conversely for the same conditions, as reported in references. Finally, the scaling law, i.e., the correlation of Nu with respect to Ra is revealed. Contours of Θ in vertical plane for different Ra with Γ=1 . Download : Download high-res image (130KB) Download : Download full-size image
physics, multidisciplinary