paper

Lattice Boltzmann study on Kelvin-Helmholtz instability: the roles of velocity and density gradients

arXiv:1010.4124 · doi:10.1103/PhysRevE.83.056704

Abstract

A two-dimensional lattice Boltzmann model with 19 discrete velocities for compressible Euler equations is proposed (D2V19-LBM). The fifth-order Weighted Essentially Non-Oscillatory (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. 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 for the KHI decreases with increasing the width of velocity transition layer but increases with increasing the width of density transition layer . After the initial transient period and before the vortex has been well formed, the linear growth rates, and , vary with and approximately in the following way, and , where , , and are fitting parameters and is the effective interaction width of density transition layer. When the linear growth rate 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, \emph{et al.}, Phys. Plasma \textbf{17}, 042103 (2010)].

Accepted for publication in PRE

References in corpus (4)

Cited by in corpus (36)