Rayleigh-Bénard simulation using the gas-kinetic Bhatnagar-Gross-Krook scheme in the incompressible limit

In this paper, a gas-kinetic Bhatnagar-Gross-Krook (BGK) model is constructed for the Rayleigh-Bénard thermal convection in the incompressible flow limit, where the flow field and temperature field are described by two coupled BGK models. Since the collision times in the corresponding BGK models can be different, the Prandtl number can be changed to any value instead of a fixed Pr=1 in the original BGK model [P. L. Bhatnagar, E. P. Gross, and M. Krook, Phys. Rev. 94, 511 (1954)]. The two-dimensional Rayleigh-Bénard thermal convection is studied and numerical results are compared with theoretical ones as well as other simulation results.