The aim of this paper is twofold: the first aim is to formulate and validate a multi-scale discrete Boltzmann method (DBM) based on density functional kinetic theory for thermal multiphase flow systems, ranging from continuum to transition flow regime; the second aim is to present some new insights into the thermo-hydrodynamic non-equilibrium (THNE) effects in the phase separation process. Methodologically, for bulk flow, DBM includes three main pillars: (i) the determination of the fewest kinetic moment relations, which are required by the description of significant THNE effects beyond the realm of continuum fluid mechanics; (ii) the construction of an appropriate discrete equilibrium distribution function recovering all the desired kinetic moments; (iii) the detection, description, presentation and analysis of THNE based on the moments of the non-equilibrium distribution ( $f-f^{(eq)}$ ). The incorporation of appropriate additional higher-order thermodynamic kinetic moments considerably extends the DBM's capability of handling larger values of the liquid–vapour density ratio, curbing spurious currents, and ensuring mass/momentum/energy conservation. Compared with the DBM with only first-order THNE (Ganet al.,Soft Matt., vol. 11 (26), 2015, pp. 5336–5345), the model retrieves kinetic moments beyond the third-order super-Burnett level, and is accurate for weak, moderate and strong THNE cases even when the local Knudsen number exceeds $1/3$ . Physically, the ending point of the linear relation between THNE and the concerned physical parameter provides a distinct criterion to identify whether the system is near or far from equilibrium. Besides, the surface tension suppresses the local THNE around the interface, but expands the THNE range and strengthens the THNE intensity away from the interface through interface smoothing and widening.

A multiple-relaxation-time discrete Boltzmann model (DBM) is proposed for multicomponent mixtures, where compressible, hydrodynamic, and thermodynamic nonequilibrium effects are taken into account. It allows the specific heat ratio and the Prandtl number to be adjustable, and is suitable for both low and high speed fluid flows. From the physical side, besides being consistent with the multicomponent Navier-Stokes equations, Fick's law, and Stefan-Maxwell diffusion equation in the hydrodynamic limit, the DBM provides more kinetic information about the nonequilibrium effects. The physical capability of DBM to describe the nonequilibrium flows, beyond the Navier-Stokes representation, enables the study of the entropy production mechanism in complex flows, especially in multicomponent mixtures. Moreover, the current kinetic model is employed to investigate nonequilibrium behaviors of the compressible Kelvin-Helmholtz instability (KHI). The entropy of mixing, the mixing area, the mixing width, the kinetic and internal energies, and the maximum and minimum temperatures are investigated during the dynamic KHI process. It is found that the mixing degree and fluid flow are similar in the KHI process for cases with various thermal conductivity and initial temperature configurations, while the maximum and minimum temperatures show different trends in cases with or without initial temperature gradients. Physically, both heat conduction and temperature exert slight influences on the formation and evolution of the KHI morphological structure.

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.

Kelvin-Helmholtz (KH) instability is a fundamental fluid instability that widely exists in nature and engineering. To better understand the dynamic process of the KH instability, the influence of the tangential velocity on the compressible KH instability is investigated by using the discrete Boltzmann method based on the nonequilibrium statistical physics. Both hydrodynamic and thermodynamic nonequilibrium (TNE) effects are probed and analyzed. It is found that, on the whole, the global density gradients, the TNE strength and area firstly increase and decrease afterwards. Both the global density gradient and heat flux intensity in the vertical direction are almost constant in the initial stage before a vortex forms. Moreover, with the increase of the tangential velocity, the KH instability evolves faster, hence the global density gradients, the TNE strength and area increase in the initial stage and achieve their peak earlier, and their maxima are higher for a larger tangential velocity. Physically, there are several competitive mechanisms in the evolution of the KH instability. (i) The physical gradients increase and the TNE effects are strengthened as the interface is elongated. The local physical gradients decrease and the local TNE intensity is weakened on account of the dissipation and/or diffusion. (ii) The global heat flux intensity is promoted when the physical gradients increase. As the contact area expands, the heat exchange is enhanced and the global heat flux intensity increases. (iii) The global TNE intensity reduces with the decreasing of physical gradients and increase with the increasing of TNE area. (iv) The nonequilibrium area increases as the fluid interface is elongated and is widened because of the dissipation and/or diffusion.