“…15,16 It is a mesoscopic approach that solves the Boltzmann equation (BE) with the Bhatnagar-Gross-Krook (BGK) approximation. 17 In comparison with the conventional computational fluid dynamics methods, based on the Navier-Stokes equations, the LBM offers some interesting advantages such as easy calculation procedure, a simple implementation for parallel computation, and simulation of complex fluid problems such as nanofluids, [18][19][20][21] MHD flow, 22 porous media, [23][24][25] multiphase flow, 26 melting PCM, 27 and turbulent flow. 22,28,29 Generally, there are three categories of thermal lattice Boltzmann models: the multispeed model, 30 the passive scalar model, 31 and the double-distribution model.…”