Li Ying-Jun

Three-Dimensional Lattice Boltzmann Model for High-Speed Compressible Flows

A highly efficient three-dimensional (3D) Lattice Boltzmann (LB) model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 (2004) 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative (NND) scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, (i) Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; (ii) reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.

Highly Efficient Lattice Boltzmann Model for Compressible Fluids: Two-Dimensional Case

We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsutahara, Phys. Rev. E 69 (2004) 035701(R)] and reasonable dispersion and dissipation terms. The dispersion term effectively reduces the oscillation at the discontinuity and enhances numerical precision. The dissipation term makes the new model more easily meet with the von Neumann stability condition. This model works for both high-speed and low-speed flows with arbitrary specific-heat-ratio. With the new model simulation results for the well-known benchmark problems get a high accuracy compared with the analytic or experimental ones. The used benchmark tests include (i) Shock tubes such as the Sod, Lax, Sjogreen, Colella explosion wave, and collision of two strong shocks, (ii) Regular and Mach shock reflections, and (iii) Shock wave reaction on cylindrical bubble problems. With a more realistic equation of state or free-energy functional, the new model has the potential tostudy the complex procedure of shock wave reaction on porous materials.

Flux Limiter Lattice Boltzmann Scheme Approach to Compressible Flows with Flexible Specific-Heat Ratio and Prandtl Number

We further develop the lattice Boltzmann (LB) model [Physica A 382 (2007) 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar—Gross—Krook (BGK) collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference (FLFD) scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including (i) The thermal Couette flow; (ii) One- and two-dimensional Riemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.

Two-Dimensional Rayleigh–Taylor Instability in Incompressible Fluids at Arbitrary Atwood Numbers

The Rayleigh–Taylor instability in two-dimensional incompressible fluids at arbitrary Atwood numbers is studied by expanding the perturbation velocity potential to third order. The second and third harmonic generation effects of single-mode perturbation are analyzed, as well as the nonlinear correction to the exponential growth of the fundamental modulation. The mode coupling coefficients are dependent on the Atwood numbers. Our simulations support the weakly nonlinear results. We find that the ratio of the nonlinear saturation amplitude ηs and the perturbation wavelength λ is dependent on the Atwood number AT and the relation is ηs/λ = (1/π) .

Multiple-Relaxation-Time Lattice Boltzmann Approach to Richtmyer—Meshkov Instability

The aims of the present paper are twofold. At first, we further study the Multiple-Relaxation-Time (MRT) Lattice Boltzmann (LB) model proposed in [Europhys. Lett. 90 (2010) 54003]. We discuss the reason why the Gram—Schmidt orthogonalization procedure is not needed in the construction of transformation matrix M; point out a reason why the Kataoka—Tsutahara model [Phys. Rev. E 69 (2004) 035701 (R)] is only valid in subsonic flows. The von Neumann stability analysis is performed. Secondly, we carry out a preliminary quantitative study on the Richtmyer—Meshkov instability using the proposed MRT LB model. When a shock wave travels from a light medium to a heavy one, the simulated growth rate is in qualitative agreement with the perturbation model by Zhang—Sohn. It is about half of the predicted value by the impulsive model and is closer to the experimental result. When the shock wave travels from a heavy medium to a light one, our simulation results are also consistent with physical analysis.

Flux Limiter Lattice Boltzmann for Compressible Flows

In this paper, a new flux limiter scheme with the splitting technique is successfully incorporated into a multiple-relaxation-time lattice Boltzmann (LB) model for shacked compressible flows. The proposed flux limiter scheme is efficient in decreasing the artificial oscillations and numerical diffusion around the interface. Due to the kinetic nature, some interface problems being difficult to handle at the macroscopic level can be modeled more naturally through the LB method. Numerical simulations for the Richtmyer—Meshkov instability show that with the new model the computed interfaces are smoother and more consistent with physical analysis. The growth rates of bubble and spike present a satisfying agreement with the theoretical predictions and other numerical simulations.

Polar Coordinate Lattice Boltzmann Kinetic Modeling of Detonation Phenomena

A novel polar coordinate lattice Boltzmann kinetic model for detonation phenomena is presented and applied to investigate typical implosion and explosion processes. In this model, the change of discrete distribution function due to local chemical reaction is dynamically coupled into the modified lattice Boltzmann equation which could recover the Navier Stokes equations, including contribution of chemical reaction, via the Chapman-Enskog expansion. For the numerical investigations, the main focuses are the none quilibrium behaviors in these processes. The system at the disc center is always in its thermodynamic equilibrium in the highly symmetric case. The internal kinetic energies in different degrees of freedom around the detonation front do not coincide. The dependence of the reaction rate on the pressure, influences of the shock strength and reaction rate on the departure amplitude of the system from its local thermodynamic equilibrium are probed.

Finite-Difference Lattice Boltzmann Scheme for High-Speed Compressible Flow: Two-Dimensional Case

Lattice Boltzmann (LB) modeling of high-speed compressible flows has long been attempted by various authors. One common weakness of most of previous models is the instability problem when the Mach number of the flow is large. In this paper we present a finite-difference LB model, which works for flows with flexible ratios of specific heats and a wide range of Mach number, from 0 to 30 or higher. Besides the discrete-velocity-model by Watari [Physica A 382 (2007) 502], a modified Lax–Wendroff finite difference scheme and an artificial viscosity are introduced. The combination of the finite-difference scheme and the adding of artificial viscosity must find a balance of numerical stability versus accuracy. The proposed model is validated by recovering results of some well-known benchmark tests: shock tubes and shock reflections. The new model may be used to track shock waves and/or to study the non-equilibrium procedure in the transition between the regular and Mach reflections of shock waves, etc.

On the Second Harmonic Generation through Bell—Plesset Effects in Cylindrical Geometry

Generation of the second harmonic initiated by Bell—Plesset effects in a cylindrical geometry is studied analytically. For an initial single-mode velocity perturbation, the second-order mode-coupling formula is obtained by expanding the perturbation displacement and velocity potential up to the second-order accuracy. It is found that the initially symmetric interface evolves into a significant bubble-spike asymmetric pattern. The second-order solutions clearly show that the amplitude of the spike grows faster than that of the bubble. The temporal evolutions of the amplitudes of the bubble and spike are dependent on the interface velocity V0. The larger interface velocity leads to the smaller amplitude of the perturbation at an arbitrary interface position in a cylindrically convergent geometry.

FFT-LB Modeling of Thermal Liquid-Vapor System

We present an improved lattice Boltzmann (LB) model for thermal liquid-vapor system. In the new model, the Windowed Fast Fourier Transform (WFFT) and its inverse are used to calculate both the convection term and the external force term of the LB equation. By adopting the WFFT scheme, Gibbs oscillations can be damped effectively in unsmooth regions while high resolution feature of the spectral method can be retained in smooth regions. As a result, spatial discretization errors are dramatically decreased, conservation of the total energy is much better preserved, and the spurious velocities near the liquid-vapor interface are significantly reduced. The high resolution, together with the low complexity of the WFFT approach, endows the proposed method with considerable potential for studying a wide class of problems in the field of multiphase flows.