“…The LBE method has become an alternative numerical method in fluid flow simulations including multiphase and multicomponent flows [8,9,[13][14][15][16][17][18][19][20], turbulent flows [13][14][15][20][21][22], microflows [4,5,15,20,23], and flows through porous media [2,8,9,11,12,14,24,25], Various TLBE models have also been developed [26][27][28][29][30][31][32][33][34][35] and some of them have been successfully applied in complex thermal simulations such as turbulent thermal flows [36,37] and thermal flows in porous media [38,39], The LBE method has its inherent advantages such as explicitness in formulation, convenience in implementation and parallelization. It is also well known for its capability to treat boundary conditions involving interfacial dynamics and complex geometries as the LBE method deals with the distribution functions, Eor fluid flow simulations, the most commonly used boundary treatment is the "bounce-back" method [40-^2]…”