“…From a computational resource perspective, the remarkable merits are brevity of programming, numerical potency, inherent parallelism, and ease treatment of intricate boundary conditions. This kind of method has comprehensive capacities in quite several fields, from phonon transport [13] to approximate incompressible flows [14][15][16][17][18][19][20][21][22][23][24][25], full compressible flows [26][27][28][29][30][31][32][33][34][35][36][37], dendrite growth [38,39] and thermal multiphase flows [40]. Recently, the mesoscopic kinetics method is also becoming increasingly popular in computational mathematics and engineering science for solving certain NPDEs, including Burgers' equations [41,42], Korteweg-de Vries equation [43], Gross-Pitaevskii equation [44], convection-diffusion equation [45][46][47][48][49][50][51], Kuramoto-Sivashinsky equation [52], wave equation [53,54], Dirac equation [55], Poisson equation…”