The Rayleigh number of the systems simulated in this paper are within the turbulent regime (Ra > 1010). Hence, the RNG k−ε model is employed to simulate the turbulent flow. The turbulence kinetic energy (k) and dissipation rate (ε) are calculated from transport equations (5) and (6)), respectively