A spectral difference lattice Boltzmann method for solution of inviscid compressible flows on structured grids

摘要

In this work, a spectral difference lattice Boltzmann method (SDLBM) is developed and applied for an accurate simulation of two,dimensional inviscid compressible flows on structured grids. The compressible form of the discrete Boltzmann-BGK equation is used in which multiple particle speeds have to be employed to correctly model the compressibility in a thermal fluid. Here, the 2D compressible Lattice Boltzmann (LB) model proposed by Watari is applied. The spectral difference (SD) method is implemented for the solution of the LB equation in which the particle distribution function is stored at the solution points while the fluxes are computed at the flux points for calculating the particle distribution function. For time accurate solutions, the fourth-order Runge-Kutta scheme is used to discretize the temporal term in the LB equation. Note that the procedure of implementing the SD method to solve the LB equation is nearly the same as that procedure developed in the literature for the solution of the Euler equations. The accuracy and robustness of the present solution methodology are demonstrated by simulating different benchmark compressible flow problems. A sensitivity study is also conducted to evaluate the effects of the numerical parameters and the grid size/distribution on the accuracy and performance of the solution. At first, three inviscid compressible flow problems, namely, the stationary isentropic vortex, the shock tube and the shock vortex interaction are solved by using the SDLBM to demonstrate the accuracy and robustness of the present solution methodology. Results computed for these test cases are in good agreement with the analytical and the available numerical solutions. To more assess the accuracy and robustness of the SDLBM, a third-order finite volume LBM (FVLBM) is also developed and the solutions obtained by these two methods for the test cases simulated are thoroughly compared with each other. It is demonstrated that the SDLBM is more accurate and robust compressible inviscid flow solver that can be used for evaluating the other compressible LB-based flow solvers. As the final problem, the supersonic flow past a bump is simulated to demonstrate the accuracy and robustness of the SDLBM in reaching to a steady-state solution. (C) 2016 Elsevier Ltd. All rights reserved.