A two-dimensional lattice Boltzmann model with 19 discrete velocities for compressible fluids is proposed. The fifth-order weighted essentially nonoscillatory (5th-WENO) finite difference scheme is employed to calculate the convection term of the lattice Boltzmann equation. The validity of the model is verified by comparing simulation results of the Sod shock tube with its corresponding analytical solutions [G. A. Sod, J. Comput. Phys. 27, 1 (1978).]. The velocity and density gradient effects on the Kelvin-Helmholtz instability (KHI) are investigated using the proposed model. Sharp density contours are obtained in our simulations. It is found that the linear growth rate $\gamma${} for the KHI decreases by increasing the width of velocity transition layer Dv but increases by increasing the width of density transition layer D$\rho${}. After the initial transient period and before the vortex has been well formed, the linear growth rates $\gamma${}v and $\gamma${}$\rho${}, vary with Dv and D$\rho${} approximately in the following way, ln$\gamma${}v=a-bDv and $\gamma${}$\rho${}=c+elnD$\rho${}(D$\rho${}D$\rho${}E the linear growth rate $\gamma${}$\rho${} does not vary significantly any more. One can use the hybrid effects of velocity and density transition layers to stabilize the KHI. Our numerical simulation results are in general agreement with the analytical results [L. F. Wang et al., Phys. Plasma 17, 042103 (2010)].

To probe both the hydrodynamic nonequilibrium (HNE) and thermodynamic nonequilibrium (TNE) in the combustion process, a two-dimensional multiple-relaxation-time (MRT) version of lattice Boltzmann kinetic model (LBKM) for combustion phenomena is presented. The chemical energy released in the progress of combustion is dynamically coupled into the system by adding a chemical term to the LB kinetic equation. Aside from describing the evolutions of the conserved quantities, the density, momentum, and energy, which are what the Navier-Stokes model describes, the MRT-LBKM presents also a coarse-grained description on the evolutions of some nonconserved quantities. The current model works for both subsonic and supersonic flows with or without chemical reaction. In this model, both the specific-heat ratio and the Prandtl number are flexible, the TNE effects are naturally presented in each simulation step. The model is verified and validated via well-known benchmark tests. As an initial application, various nonequilibrium behaviors, including the complex interplays between various HNEs, between various TNEs, and between the HNE and TNE, around the detonation wave in the unsteady and steady one-dimensional detonation processes are preliminarily probed. It is found that the system viscosity (or heat conductivity) decreases the local TNE, but increases the global TNE around the detonation wave, that even locally, the system viscosity (or heat conductivity) results in two kinds of competing trends, to increase and to decrease the TNE effects. The physical reason is that the viscosity (or heat conductivity) takes part in both the thermodynamic and hydrodynamic responses.

We present a polar coordinate lattice Boltzmann kinetic model for compressible flows. A method to recover the continuum distribution function from the discrete distribution function is indicated. Within the model, a hybrid scheme being similar to, but different from, the operator splitting is proposed. The temporal evolution is calculated analytically, and the convection term is solved via a modified Warming-Beam (MWB) scheme. Within the MWB scheme a suitable switch function is introduced. The current model works not only for subsonic flows but also for supersonic flows. It is validated and verified via the following well-known benchmark tests: (i) the rotational flow, (ii) the stable shock tube problem, (iii) the Richtmyer-Meshkov (RM) instability, and (iv) the Kelvin-Helmholtz instability. As an original application, we studied the nonequilibrium characteristics of the system around three kinds of interfaces, the shock wave, the rarefaction wave, and the material interface, for two specific cases. In one of the two cases, the material interface is initially perturbed, and consequently the RM instability occurs. It is found that the macroscopic effects due to deviating from thermodynamic equilibrium around the material interface differ significantly from those around the mechanical interfaces. The initial perturbation at the material interface enhances the coupling of molecular motions in different degrees of freedom. The amplitude of deviation from thermodynamic equilibrium around the shock wave is much higher than those around the rarefaction wave and material interface. By comparing each component of the high-order moments and its value in equilibrium, we can draw qualitatively the main behavior of the actual distribution function. These results deepen our understanding of the mechanical and material interfaces from a more fundamental level, which is indicative for constructing macroscopic models and other kinds of kinetic models.

A discrete Boltzmann model (DBM) is proposed to probe the Rayleigh-Taylor instability (RTI) in two-component compressible flows. Each species has a flexible specific-heat ratio and is described by one discrete Boltzmann equation (DBE). Independent discrete velocities are adopted for the two DBEs. The collision and force terms in the DBE account for the molecular collision and external force, respectively. Two types of force terms are exploited. In addition to recovering the modified Navier-Stokes equations in the hydrodynamic limit, the DBM has the capability of capturing detailed nonequilibrium effects. Furthermore, we use the DBM to investigate the dynamic process of the RTI. The invariants of tensors for nonequilibrium effects are presented and studied. For low Reynolds numbers, both global nonequilibrium manifestations and the growth rate of the entropy of mixing show three stages (i.e., the reducing, increasing, and then decreasing trends) in the evolution of the RTI. On the other hand, the early reducing tendency is suppressed and even eliminated for high Reynolds numbers. Relevant physical mechanisms are analyzed and discussed.

We investigate thermal and isothermal symmetric liquid-vapor separations via a fast Fourier transform thermal lattice Boltzmann (FFT-TLB) model. Structure factor, domain size, and Minkowski functionals are employed to characterize the density and velocity fields, as well as to understand the configurations and the kinetic processes. Compared with the isothermal phase separation, the freedom in temperature prolongs the spinodal decomposition (SD) stage and induces different rheological and morphological behaviors in the thermal system. After the transient procedure, both the thermal and isothermal separations show power-law scalings in domain growth, while the exponent for thermal system is lower than that for isothermal system. With respect to the density field, the isothermal system presents more likely bicontinuous configurations with narrower interfaces, while the thermal system presents more likely configurations with scattered bubbles. Heat creation, conduction, and lower interfacial stresses are the main reasons for the differences in thermal system. Different from the isothermal case, the release of latent heat causes the changing of local temperature, which results in new local mechanical balance. When the Prandtl number becomes smaller, the system approaches thermodynamical equilibrium much more quickly. The increasing of mean temperature makes the interfacial stress lower in the following way: $\sigma${}=$\sigma${}0[(Tc$-${}T)/(Tc$-${}T0)]3/2, where Tc is the critical temperature and $\sigma${}0 is the interfacial stress at a reference temperature T0, which is the main reason for the prolonged SD stage and the lower growth exponent in the thermal case. Besides thermodynamics, we probe how the local viscosities influence the morphology of the phase separating system. We find that, for both the isothermal and thermal cases, the growth exponents and local flow velocities are inversely proportional to the corresponding viscosities. Compared with the isothermal case, the local flow velocity depends not only on viscosity but also on temperature.

The aims of the present paper are threefold. First, we further study the fast Fourier transform thermal lattice Boltzmann (FFT-TLB) model for van der Waals (VDW) fluids proposed in Phys. Rev. E, 2011, 84(4): 046715. We analyze the merits of the FFT approach over the traditional finite difference scheme and investigate the effects of smoothing factors on accuracy and stability in detail. Second, we incorporate the VDW equation of state with flexible parameters into the FFT-TLB model. As a result, the revised model may be used to handle multiphase flows with various critical densities and temperatures. Third, we design appropriate boundary conditions for systems with solid walls. These improvements, from the views of numerics and physics, significantly extend the application scope of the model in science and engineering.