Three-dimensional finite-difference lattice Boltzmann model and its application to inviscid compressible flows with shock waves

In this paper, a three-dimensional (3D) finite-difference lattice Boltzmann model for simulating compressible flows with shock waves is developed in the framework of the double-distribution-function approach. In the model, a density distribution function is adopted to model the flow field, while a total energy distribution function is adopted to model the temperature field. The discrete equilibrium density and total energy distribution functions are derived from the Hermite expansions of the continuous equilibrium distribution functions. The discrete velocity set is obtained by choosing the abscissae of a suitable Gauss–Hermite quadrature with sufficient accuracy. In order to capture the shock waves in compressible flows and improve the numerical accuracy and stability, an implicit–explicit finite-difference numerical technique based on the total variation diminishing flux limitation is introduced to solve the discrete kinetic equations. The model is tested by numerical simulations of some typical compressible flows with shock waves ranging from 1D to 3D. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.     

Three-dimensional coupled double-distribution-function lattice Boltzmann models for compressible Navier–Stokes equations

Two three-dimensional (3D) lattice Boltzmann models in the framework of coupled double-distribution-function approach for compressible flows, in which specific-heat ratio and Prandtl number can be adjustable, are developed in this paper. The main differences between the two models are discrete equilibrium density and total energy distribution function. One is the D3Q25 model obtained from spherical function, and the other is the D3Q27 standard lattice model obtained from Hermite expansions of the corresponding continuous equilibrium distribution functions. The two models are tested by numerical simulations of some typical compressible flows, and their numerical stability and precision are also analysed. The results indicate that the two models are capable for supersonic flows, while the one from Hermite expansions is not suitable for compressible flows with shock waves.     

Development and Comparative Studies of Three Non-Free Parameter Lattice Boltzmann Models for Simulation of Compressible Flows

This paper at first shows the details of finite volume-based lattice Boltzmann method (FV-LBM) for simulation of compressible flows with shock waves. In the FV-LBM, the normal convective flux at the interface of a cell is evaluated by using one-dimensional compressible lattice Boltzmann model, while the tangential flux is calculated using the same way as used in the conventional Euler solvers. The paper then presents a platform to construct one-dimensional compressible lattice Boltzmann model for its use in FV-LBM. The platform is formed from the conservation forms of moments. Under the platform, both the equilibrium distribution functions and lattice velocities can be determined, and therefore, non-free parameter model can be developed. The paper particularly presents three typical non-free parameter models, D1Q3, D1Q4 and D1Q5. The performances of these three models for simulation of compressible flows are investigated by a brief analysis and their application to solve some one-dimensional and two-dimensional test problems. Numerical results showed that D1Q3 model costs the least computation time and D1Q4 and D1Q5 models have the wider application range of Mach number. From the results, it seems that D1Q4 model could be the best choice for the FV-LBM simulation of hypersonic flows.     

A Multi-Speed Compressible Lattice-Boltzmann Model

An extension of the lattice Boltzmann BGK method to compressible flows is presented that combines three novel additions: (1) particles move density and energy weights in multiple velocity bins (11 for 1-D flow) to nearby cell centers. (2) the equilibrium distribution remains an unexpanded Maxwellian; and (3) transport and relaxation to equilibrium are performed implicitly at each node. These advances allow for the parallel modeling of high Mach number shocks and high Reynolds number flows, while avoiding advective numerical diffusion, the need for Riemann solvers, and non-linear limiters. A 1D shock tube application is shown. Generalization to higher dimensions and multi-materials are discussed.     

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 isvalidated 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.     

A Switch Function-Based Gas-Kinetic Scheme for Simulation of Inviscid and Viscous Compressible Flows

In this paper, a switch function-based gas-kinetic scheme (SF-GKS) is presented for the simulation of inviscid and viscous compressible flows. With the finite volume discretization, Euler and Navier-Stokes equations are solved and the SF-GKS is applied to evaluate the inviscid flux at cell interface. The viscous flux is obtained by the conventional smooth function approximation. Unlike the traditional gas-kinetic scheme in the calculation of inviscid flux such as Kinetic Flux Vector Splitting (KFVS), the numerical dissipation is controlled with a switch function in the present scheme. That is, the numerical dissipation is only introduced in the region around strong shock waves. As a consequence, the present SF-GKS can well capture strong shock waves and thin boundary layers simultaneously. The present SF-GKS is firstly validated by its application to the inviscid flow problems, including 1-D Euler shock tube, regular shock reflection and double Mach reflection. Then, SF-GKS is extended to solve viscous transonic and hypersonic flow problems. Good agreement between the present results and those in the literature verifies the accuracy and robustness of SF-GKS.     

激波的格子玻尔兹曼方法模拟

格子玻尔兹曼方法已被广泛应用于模拟流体,但主要应用于模拟不可压流体。应用二维正方格阵上的3速9点模型模拟激波的传播。考察了激波关系式、激波的宽度和结构。模拟结果与由Chapman-Enskog展开导出的流体方程的解相符。这说明格子玻尔兹曼方法能够模拟可压缩流体。     

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.     

Simple and efficient relaxation methods for interfaces separating compressible fluids,cavitating flows and shocks in multiphase mixtures

Numerical approximation of the five-equation two-phase flow of Kapila et al. [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10) (2001) 3002–3024] is examined. This model has shown excellent capabilities for the numerical resolution of interfaces separating compressible fluids as well as wave propagation in compressible mixtures [A. Murrone, H. Guillard, A five equation reduced model for compressible two phase flow problems, Journal of Computational Physics 202(2) (2005) 664–698; R. Abgrall, V. Perrier, Asymptotic expansion of a multiscale numerical scheme for compressible multiphase flows, SIAM Journal of Multiscale and Modeling and Simulation (5) (2006) 84–115; F. Petitpas, E. Franquet, R. Saurel, O. Le Metayer, A relaxation-projection method for compressible flows. Part II. The artificial heat exchange for multiphase shocks, Journal of Computational Physics 225(2) (2007) 2214–2248]. However, its numerical approximation poses some serious difficulties. Among them, the non-monotonic behavior of the sound speed causes inaccuracies in wave’s transmission across interfaces. Moreover, volume fraction variation across acoustic waves results in difficulties for the Riemann problem resolution, and in particular for the derivation of approximate solvers. Volume fraction positivity in the presence of shocks or strong expansion waves is another issue resulting in lack of robustness. To circumvent these difficulties, the pressure equilibrium assumption is relaxed and a pressure non-equilibrium model is developed. It results in a single velocity, non-conservative hyperbolic model with two energy equations involving relaxation terms. It fulfills the equation of state and energy conservation on both sides of interfaces and guarantees correct transmission of shocks across them. This formulation considerably simplifies numerical resolution. Following a strategy developed previously for another flow model [R. Saurel, R. Abgrall, A multiphase Godunov method for multifluid and multiphase flows, Journal of Computational Physics 150 (1999) 425–467], the hyperbolic part is first solved without relaxation terms with a simple, fast and robust algorithm, valid for unstructured meshes. Second, stiff relaxation terms are solved with a Newton method that also guarantees positivity and robustness. The algorithm and model are compared to exact solutions of the Euler equations as well as solutions of the five-equation model under extreme flow conditions, for interface computation and cavitating flows involving dynamics appearance of interfaces. In order to deal with correct dynamic of shock waves propagating through multiphase mixtures, the artificial heat exchange method of Petitpas et al. [F. Petitpas, E. Franquet, R. Saurel, O. Le Metayer, A relaxation-projection method for compressible flows. Part II. The artificial heat exchange for multiphase shocks, Journal of Computational Physics 225(2) (2007) 2214–2248] is adapted to the present formulation.     

一个三维多松弛时间全速域格子Boltzmann模型

构建一个既适用于低速不可压流体又适用于高速可压缩流体的三维自由参数多松弛时间格子Boltzmann模型.模型中,根据SO（3）群的不可约表述基函数构造转化矩阵,根据恢复可压Navier-Stokes方程的需要选取非守恒矩平衡值.通过von Neumann稳定性分析模型参数对数值稳定性的影响,并给出建议选择范围.模型经过基准问题的验证,模拟结果与解析解及其它数值结果符合较好.     

A Hybrid Lattice Boltzmann Flux Solver for Simulation of Viscous Compressible Flows

In this paper, a hybrid lattice Boltzmann flux solver (LBFS) is proposed for simulation of viscous compressible flows. In the solver, the finite volume method is applied to solve the Navier-Stokes equations. Different from conventional Navier-Stokes solvers, in this work, the inviscid flux across the cell interface is evaluated by local reconstruction of solution using one-dimensional lattice Boltzmann model, while the viscous flux is still approximated by conventional smooth function approximation. The present work overcomes the two major drawbacks of existing LBFS [28–31], which is used for simulation of inviscid flows. The first one is its ability to simulate viscous flows by including evaluation of viscous flux. The second one is its ability to effectively capture both strong shock waves and thin boundary layers through introduction of a switch function for evaluation of inviscid flux, which takes a value close to zero in the boundary layer and one around the strong shock wave. Numerical experiments demonstrate that the present solver can accurately and effectively simulate hypersonic viscous flows.     

Direct Numerical Simulation of Three-Dimensional Richtmyer--Meshkov Instability

Direct numerical simulation （DNS） is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interracial instability is known as the Richtmyer-Meshkov （R-M） instability. The compressible Navier-Stoke equations are discretized with group velocity control （GVC） modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interracial instability on produced vorticity and mixing is discussed.     

Limited projections laser speckle tomography of complex flows

Different approaches are discussed applied to reconstruction of local flow parameters using line-of-sight laser probing. The first one is based on single projection measurement and detailed numerical simulation of flow pattern. The second one is based on multi-projectional line-of-sight measurement and reconstructing local flow parameters by solving inverse integral transforms. Holographic interferometry (HI) and speckle photography (SP) are used for laser probing of complex 3D media under study. Computer-assisted tomography (CAT) approach based on Radon integral transform is described and examples of reconstructed flowfields are given including compressible flows with shock waves. Statistical information about flow studied is extracted as well, and maps of local statistical turbulence parameters are quantitatively determined by using obtained SP-data and a new Erbeck–Merzkirch integral transform.     

基于格子Boltzmann方法的超声速非平衡流模拟

基于D1Q4可压缩格子Boltzmann模型,按照流通矢量分裂方法的思路,采用坐标旋转技术构造求解三维带化学反应Navier-Stokes方程对流通量求解器.结合有限体积法求解三维化学非平衡流Navier-Stokes方程,采用时间算子分裂算法解决化学反应刚性问题,数值模拟超声速化学非平衡流的三个经典算例.数值结果表明:在高马赫数下,采用D1Q4可压缩格子Boltzmann模型构造的三维对流通量求解器数值模拟中没有出现非物理解,同时在超声速化学非平衡流场中正确分辨激波、燃烧波等物理现象,精度和分辨率均较高,验证了本文构造的三维对流通量求解器的可靠性,拓宽了D1Q4可压缩格子Boltzmann模型的应用范围,为计算超声速化学非平衡流提供一种新方法.     

Assessment of localized artificial diffusivity scheme for large-eddy simulation of compressible turbulent flows

The localized artificial diffusivity method is investigated in the context of large-eddy simulation of compressible turbulent flows. The performance of different artificial bulk viscosity models are evaluated through detailed results from the evolution of decaying compressible isotropic turbulence with eddy shocklets and supersonic turbulent boundary layer. Effects of subgrid-scale (SGS) models and implicit time-integration scheme/time-step size are also investigated within the framework of the numerical scheme used. The use of a shock sensor along with artificial bulk viscosity significantly improves the scheme for simulating turbulent flows involving shocks while retaining the shock-capturing capability. The proposed combination of Ducros-type sensor with a negative dilatation sensor removes unnecessary bulk viscosity within expansion and weakly compressible turbulence regions without shocks and allows it to localize near the shocks. It also eliminates the need for a wall-damping function for the bulk viscosity while simulating wall-bounded turbulent flows. For the numerical schemes used, better results are obtained without adding an explicit SGS model than with SGS model at moderate Reynolds number. Inclusion of a SGS model in addition to the low-pass filtering and artificial bulk viscosity results in additional damping of the resolved turbulence. However, investigations at higher Reynolds numbers suggest the need for an explicit SGS model. The flow statistics obtained using the second-order implicit time-integration scheme with three sub-iterations closely agrees with the explicit scheme if the maximum Courant–Friedrichs–Lewy is kept near unity.     

Physical investigation of acoustic waves induced by the oscillation and collapse of the single bubble

The objective of this paper is to apply both experimental and numerical methods to investigate acoustic waves induced by the oscillation and collapse of a single bubble. In the experiments, the schlieren technique is used to capture the temporal evolution of the bubble shapes, and the corresponding acoustic waves. The results are presented for the single bubble generated by a low-voltage bubble generator in the free field of water. During the numerical simulations, a three-dimensional (3D) weakly compressible model is introduced to investigate the single bubble dynamics, including the generation and propagation of acoustic waves. The results show that (1) Compression wave, rarefaction wave and shock wave are generated during expansion stage, collapse stage and rebound stage of the bubble respectively. (2) Compression waves are induced by the rapid expansion of the bubble and eventually steepen into one shock wave propagating outward in the liquid, then another strong shock wave is emitted at the final collapse stage. The velocity and pressure of the liquid field increases after the shock wave. (3) Rarefaction waves are generated during the collapse stage due to the contraction of the bubble. The rarefaction wave reduces the liquid pressure and its spatial distribution is dispersive. The pressure of these acoustic waves and their effect on the liquid velocity attenuate with the increase of propagation distance.     

A multiple-relaxation-time lattice Boltzmann method for high-speed compressible flows

This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time(SRT) lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations,it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version.     

Lower order three-dimensional Burgers equation having non-Maxwellian ions in dusty plasmas

The dust acoustic(DA) shock wave with dust charge fluctuations, non-Maxwellian ions, and non-isothermal electrons is studied theoretically. The perturbation technique is employed to derive the lower order three-dimensional(3D) Burgers equation, and shock wave solution is explored by the tan-hyperbolic method. The effects of flat trapped and trapped electron distributions in the presence of Maxwellian and non-Maxwellian ions on characteristics shock waves are observed. The temperature ratio of non-Maxwellian ion temperature and non-isothermal electron temperature is found to play an important role in forming the shock-like structure.     

Density-PDFs and Lagrangian statistics of highly compressible turbulence

In isothermal, highly compressible turbulent flows, density fluctuations follow a log-normal distribution. We establish a connection between these density fluctuations and the probability-density-functions (PDF) of Lagrangian tracer particles advected with the flow. Our predicted particle statistics is tested against large scale numerical simulations, which were performed with 512³ collocation points and 2 million tracer particles integrated over several dynamical times.     

高密度比多介质可压缩流动的PPM方法

介绍一种基于流体体积分数模型的界面捕捉方法,数值模拟高密度比及含强激波的多介质可压缩流动.将PPM方法应用到多介质体积分数形式的Euler方程组,用双波近似求解一般刚性气体状态方程Riemann问题,可方便处理界面强剪切的滑移线问题,并给出了水下激波和气泡相互作用以及多相Richtmyer Meshkov不稳定性的计算结果.