“…In order to alleviate the heavy computational burden of directly solving Boltzmann equation, a variety of Boltzmann equationbased methods, such as the unified gas-kinetic scheme (UGKS) [10,11], the discrete velocity method (DVM) [12], the discrete unified gas-kinetic scheme (DUGKS) [13,14], the lattice Boltzmann method (LBM) [2,[15][16][17][18][19][20][21][22][23][24][25][26], have been presented and well developed. Recently, Discrete Boltzmann Method (DBM) has also been developed and widely used in various complex flow simulations [27][28][29], such as multiphase flows [30], flow instability [31,32], combustion and detonation [33,34], etc. From the viewpoint of numerical algorithm, similar to finite-different LBM, the velocity space is substituted by a limited number of particle velocities in DBM.…”