In this paper, a double-population thermal lattice Boltzmann method (DDF) is used due to its characteristics such as higher stability and ability to simulate natural convection problems. In this method density distribution function is used to solve momentum equations and energy distribution function is used to solve energy equation. In this study, a multi relaxation time model (MRT) and a single one (SRT) are proposed to estimate the density and energy distributions of collision operators, respectively. Then, using SRT-MRT model for collision operator, the pure heat conductivity on the square plate, natural convection in the cavity and force convection in the channel are simulatedwith four boundary conditions. Then, the results of the present method are compare with the result of analytical solutions and previous studies to show the ability of the double populations lattice Boltzmann method with the SRT-MRT collision operator are validated.