Rayleigh-Taylor (RT) instability widely exists in nature and engineering fields. How to better understand the physical mechanism of RT instability is of great theoretical significance and practical value. At present, abundant results of RT instability have been obtained by traditional macroscopic methods. However, research on the thermodynamic non-equilibrium (TNE) effects in the process of system evolution is relatively scarce. In this paper, the discrete Boltzmann method based on non-equilibrium statistical physics is utilized to study the effects of the specific heat ratio on compressible RT instability. The evolution process of the compressible RT system with different specific heat ratios can be analyzed by the temperature gradient and the proportion of the non-equilibrium region. Firstly, as a result of the competition between the macroscopic magnitude gradient and the non-equilibrium region, the average TNE intensity first increases and then reduces, and it increases with the specific heat ratio decreasing; the specific heat ratio has the same effect on the global strength of the viscous stress tensor. Secondly, the moment when the total temperature gradient in y direction deviates from the fixed value can be regarded as a physical criterion for judging the formation of the vortex structure. Thirdly, under the competition between the temperature gradients and the contact area of the two fluids, the average intensity of the non-equilibrium quantity related to the heat flux shows diversity, and the influence of the specific heat ratio is also quite remarkable.

Based on the framework of our previous work [H.L. Lai et al., Phys. Rev. E, 94, 023106 (2016)], we continue to study the effects of Knudsen number on two-dimensional Rayleigh–Taylor (RT) instability in compressible fluid via the discrete Boltzmann method. It is found that the Knudsen number effects strongly inhibit the RT instability but always enormously strengthen both the global hydrodynamic non-equilibrium (HNE) and thermodynamic non-equilibrium (TNE) effects. Moreover, when Knudsen number increases, the Kelvin–Helmholtz instability induced by the development of the RT instability is difficult to sufficiently develop in the later stage. Different from the traditional computational fluid dynamics, the discrete Boltzmann method further presents a wealth of non-equilibrium information. Specifically, the two-dimensional TNE quantities demonstrate that, far from the disturbance interface, the value of TNE strength is basically zero; the TNE effects are mainly concentrated on both sides of the interface, which is closely related to the gradient of macroscopic quantities. The global TNE first decreases then increases with evolution. The relevant physical mechanisms are analyzed and discussed.

The Richtmyer–Meshkov instability in a two-component system during the reshock process for various density ratios is studied through the discrete Boltzmann method. Detailed investigations are conducted on both hydrodynamic and thermodynamic non-equilibrium behaviors. Specifically, the analysis focuses on the density gradient, viscous stress tensor, heat flux strength, thermodynamic non-equilibrium intensity, and thermodynamic non-equilibrium area. It is interesting to observe the complex variations to non-equilibrium quantities with the changing shock front, rarefaction wave, transverse wave, and material interface. Physically, the non-equilibrium area is extended as the perturbed material interface grows after the passing of the shock wave or secondary impact. Moreover, the global non-equilibrium manifestation decreases when the transmitted shock front and transverse waves leave or when the reflected rarefaction wave weakens. Additionally, the global thermodynamic non-equilibrium effect is enhanced as the physical gradients or non-equilibrium area increase. Finally, the local non-equilibrium effect decreases when the fluid structure gradually disappears under the action of dissipation/diffusion.

In this paper, a new lattice Boltzmann model for the two-component system of coupled sine-Gordon equations is presented by using the coupled mesoscopic Boltzmann equations. Via the Chapman-Enskog multiscale expansion, the macroscopical governing evolution system can be recovered correctly by selecting suitable discrete equilibrium distribution functions and the amending functions. The mesoscopic model has been validated by several related issues where analytic solutions are available. The experimental results show that the numerical results are consistent with the analytic solutions. From the mesoscopic point of view, the present approach provides a new way for studying the complex nonlinear partial differential equations arising in natural nonlinear phenomena of engineering and science.

This study investigates the complex kinetics of thermodynamic nonequilibrium effects (TNEs) and their relative importance during the development of the Kelvin–Helmholtz instability (KHI) using high-order discrete Boltzmann models (DBMs). First, the capabilities and differences among various discrete velocity sets in capturing TNEs and distribution functions are assessed. This analysis proposes practical guidelines for constructing discrete velocity stencils to enhance phase-space discretization and improve the robustness of high-order DBM simulation. At different stages of KHI and under varying initial conditions, multiscale TNEs, such as viscous stresses of different orders, emerge with distinct dominant roles. Specifically, three scenarios are identified: (i) regimes dominated by first-order TNEs, (ii) alternation between first- and second-order TNEs, and (iii) states where second-order TNEs govern the system's behavior. To quantitatively capture these transitions, criteria for TNE dominance at different orders in KHI evolution are established based on the relative thermodynamic nonequilibrium intensity (RTNE). In scenarios dominated by second-order TNEs, differences between first-order and second-order models are compared in terms of macroscopic quantities, nonequilibrium effects, and kinetic moments, revealing the physical limitations of low-order models in capturing TNEs. Furthermore, the effectiveness, extensibility, and limitations of a representative high-order model are examined under second-order TNE-dominated conditions. To encapsulate these findings, a nonequilibrium phase diagram that visually maps the multiscale characteristics of KHI is constructed. This diagram not only provides intuitive insights into the dynamic interplay of different nonequilibrium effects but also serves as a kinetic roadmap for selecting suitable models under diverse nonequilibrium conditions.