“…The lattice Boltzmann equations (LBE) for the flow and thermal energy transports are written as follow: (23) where, the equilibrium distribution functions (f k eq and g k eq for the flow and thermal fields, respectively) are given as (24) (25) Here, Δt, Δx (=e k Δt), c (=Δx/Δt=1.0), u, T, e k and w k are lattice time step, lattice space step, lattice streaming speed, velocity vector, temperature, discrete velocity vector and equilibrium distribution of weight for direction k, respectively. The relaxation parameters (τ ν and τ g for the momentum and thermal energy transport fields, respectively) are defined as (26) (27) Under the assumptions that the buoyancy force is directly proportional to the temperature gradient and negligible radiative heat transfer, the Boussinesq approximation is accounted [5,34,47,48] by calculating the buoyancy force in the vertical direction as follow.…”