Burnett-level discrete Boltzmann modeling of compressible flows under force

In this paper, a Burnett-level discrete Boltzmann method (DBM) is proposed for the compressible flow in a force field, and a discrete velocity set with 25 velocities is constructed for the DBM, featuring good spatial symmetry. In the discrete Boltzmann equation, both the discrete equilibrium distribution function and the force term satisfy 25 independent moment relations and are computed with the matrix inversion method. This approach ensures high physical accuracy, computational efficiency, and ease of implementation. Through the Chapman–Enskog expansion analysis, it is demonstrated that the current DBM can recover the Burnett equations for the compressible system under force in the continuum limit. Moreover, the DBM has the capability of capturing essential thermodynamic nonequilibrium behaviors. Finally, the model is validated through five typical benchmarks, including the free falling, Sod shock tube, sound wave, thermal Couette flow, and Rayleigh–Taylor instability.