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

We present a highly efficient lattice Boltzmann model for simulatingcompressible flows. This model is based on the combination of an appropriatefinite difference scheme, a 16-discrete-velocity model [Kataoka andTsutahara, Phys. Rev. E 69 (2004) 035701(R)] and reasonable dispersion anddissipation terms. The dispersion term effectively reduces the oscillationat the discontinuity and enhances numerical precision. The dissipation termmakes the new model more easily meet with the von Neumann stabilitycondition. This model works for both high-speed and low-speed flows witharbitrary specific-heat-ratio. With the new model simulation results for thewell-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 ofstate or free-energy functional, the new model has the potential tostudythe complex procedure of shock wave reaction on porous materials.     

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

A highly efficient three-dimensional （31）） 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.     

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.     

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.     

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.     

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.     

Multiple-relaxation-time lattice Boltzmann model for compressible fluids

We present an energy-conserving multiple-relaxation-time finite difference lattice Boltzmann model for compressible flows. The collision step is first calculated in the moment space and then mapped back to the velocity space. The moment space and corresponding transformation matrix are constructed according to the group representation theory. Equilibria of the nonconserved moments are chosen according to the need of recovering compressible Navier-Stokes equations through the Chapman-Enskog expansion. Numerical experiments showed that compressible flows with strong shocks can be well simulated by the present model. The new model works for both low and high speeds compressible flows. It contains more physical information and has better numerical stability and accuracy than its single-relaxation-time version.     

Numerical Simulation of Two-Dimensional Flow over Three Cylinders by LatticeBoltzmann Method

The numerical simulation using the multiple relaxation time lattice Boltzmann method (MRT-LBM) is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these three circular cylinders, one of the three cylinders on which a forced in-line vibrating is used to do this research and attempt to find out the effects of the movingcylinder and the other two rigid cylinders on the wake characteristics and vortex formation. As a benchmark problem to discuss the problem of lift coefficient r.m.s for these cylinders with spacing ratios T/D between other rigid side-by-side cylinders, and the calculation is carried out with two compared cases at Reynolds number of 100, two of the cylinders are rigid and the other one is an in-line vibrated cylinder lying downstream, in addition, forced vibrating amplitude and frequency are A/D= 0.5 and f_ν=0.4 (where A is the forced amplitude, D is the cylinder diameter, and f_ν stands for the vibrating frequency, respectively). The calculated results not onlyindicate that the spacing ratios T/D (T is the center-to-center spacing between the two upstream cylinders) have influence on the wake patterns and the formation of vortex shedding, but also analyze the lift coefficient r.m.s for the three cylinders with the spacing ratios S/D (where S is the center-to-center spacing between the center of upstream two side-by-side cylinders and downstream cylinder).     

Lattice Boltzmann simulation of liquid–vapor system by incorporating a surface tension term

In this study,we investigate the pseudopotential multiphase model of lattice Boltzmann method(LBM)and incorporate a surface tension term to implement the particle interaction force.By using the Carnahan–Starling(CS)equation of state(EOS)with a proper critical pressure–density ratio,a density ratio over 160000 is obtained with satisfactory numerical stability.The added surface tension term offers a flexible choice to adjust the surface tension strength.Numerical tests of the Laplace rule are conducted,proving that smaller spurious velocity and better numerical stability can be acquired as the surface tension becomes stronger.Moreover,by wall adhesion and heterogeneous cavitation tests,the surface tension term shows its practical application in dealing with problems in which the surface tension plays an important role.     

激光与靶相互作用体汽化模型的数值计算方法

利用克劳修斯-克拉珀龙(Clausius Clapeyron)方程描述相变过程,利用von Neumann差分方法构造辐射流体动力学方程组的计算方法,数值求解了激光辐照金属材料的热效应、汽化过程和激光维持爆震波,并同作者前期的工作作了比较.结果说明,该计算方法对激光与靶材料相互作用体汽化模型的数值计算是行之有效的,且该汽化模型为相关问题的二维、三维问题的数值计算提供了简便有效的途径.     

A lattice Boltzmann equation for diffusion

The formulation of lattice gas automata (LGA) for given partial differential equations is not straightforward and still requires some sort of magic. Lattice Boltzmann equation (LBE) models are much more flexible than LGA because of the freedom in choosing equilibrium distributions with free parameters which can be set after a multiscale expansion according to certain requirements. Here a LBE is presented for diffusion in an arbitrary number of dimensions. The model is probably the simplest LBE which can be formulated. It is shown that the resulting algorithm with relaxation parameter =1 is identical to an explicit finite-difference (EFD) formulation at its stability limit. Underrelaxation (0<<1) allows stable integration beyond the stability limit of EFD. The time step of the explicit LBE integration is limited by accuracy and not by stability requirements.     

微通道内反应流体粘性指进的数值研究

采用格子Boltzmann方法数值模拟化学反应中混溶流体在微通道中的粘性指进现象.模拟采用单浓度变量的双稳态化学反应模型,重点研究指进的形态位置随化学反应速率和稳态浓度参数（即无化学反应发生的界面浓度）的变化.结果表明：随着反应速率的增加,指进界面变薄;而稳态浓度参数的变化则影响反应区的分布以及反应速率,导致指进形态以及位置的改变,甚至出现指尖液滴分离.     

静态电磁场边值问题计算方法

介绍常用的电磁场分析方法——解析法、有限差分法、有限元法.以静态电磁场边值问题为例,分别用自己编写的MATLAB有限差分程序及大型有限元分析软件ANSYS求出了数值解,并与解析法得到的精确解进行比较,得出了用数值法求解电磁场问题基本满足工程需要的结论.采用数值计算时,对如何减小计算机内存,取消冗余数据,也做了进一步的讨论.     

生物过滤器中非均匀性流动的数值研究

生物过滤技术因其具有有效性、低成本和环境友好等优点引起了人们的广泛关注.该技术主要通过生物过滤器去除含有H2S等废气的有毒有害气体.运用格子Boltzmann方法对三种生物过滤器模型中多孔介质的非均匀性流动进行了数值模拟.数值模拟结果表明,多孔介质的性质和进口流动条件对临界Rayleigh数有显著影响,临界Rayleigh数随着多孔介质的孔隙度和Darcy数的增大而逐渐变小,并随着进口Reynolds数的增大而逐渐变大.所得结果可望为生物过滤器的优化设计提供一个合理的理论依据.     

Stability analysis of a predictor/multi-corrector method for staggered-grid Lagrangian shock hydrodynamics

This article presents the complete von Neumann stability analysis of a predictor/multi-corrector scheme derived from an implicit mid-point time integrator often used in shock hydrodynamics computations in combination with staggered spatial discretizations. It is shown that only even iterates of the method yield stable computations, while the odd iterates are, in the most general case, unconditionally unstable. These findings are confirmed by, and illustrated with, a number of numerical computations. Dispersion error analysis is also presented.     

Test of the Possible Application of the Half—Way Bounce—Back Boundary Condition for Lattice Boltzmann Methods in Complex Geometry

The way of handling boundary conditions with simple bounce-back rule in the lattice gas and lattice Boltzmann method had been considered as one of the advantage compared with other numerical schemes,The half-way bounce-back rule inherits the advantage of the bounce-back rule and improves the accuracy to the second-order on flat boundaries,In this paper,we test the possible application of the half-way bounce-back rule to the system with complex geometry,Our simulation results show that the half-way bounce-back rule is a good boundary condition in the problems without emphasis on accuracy.     

液滴撞击圆柱内表面的数值研究

针对液滴撞击圆柱内表面的过程,利用基于相场的格子Boltzmann方法模拟液滴以不同初速度、从不同初始高度、撞击不同大小的圆柱内表面时液滴的形态变化,分析了液滴自身物性(如密度和黏性等)和圆柱内表面润湿性等因素对撞击现象的具体影响.研究发现:撞击韦伯数、密度比及动力黏性比、圆柱半径等对液滴撞击后沿圆柱内表面的铺展均有一定影响,较高的韦伯数下液滴可能会发生分裂;液滴初始高度对大密度比和动力黏性比的撞击影响较小;液滴反弹现象可能出现在接触角较大时;重力作用会抑制撞击后液滴的振荡.     

基于格子Boltzmann方法的自驱动Janus颗粒扩散泳力

Janus颗粒的自驱动力研究对于纳微米尺度驱动力课题具有重要意义,本文针对Pt-SiO_2型Janus颗粒,基于格子Boltzmann模型及动量交换法提出了计算其扩散泳力的方法,通过与实验数据对比修正验证了模型准确性,并通过分析证明了此类Janus颗粒的扩散泳力与速度场无关,进一步模拟比较了不同形状颗粒的自驱运动.分析发现,对于体积相等形状不同的Janus颗粒,扩散泳力主要由轴线投影面积决定,此外反应面积也会对扩散泳力产生影响.     

基于MRT-LB伪势模型的液滴撞击液膜数值研究

运用改进的格子玻尔兹曼（LB）伪势多松弛多相模型,研究单/双液滴撞击液膜的流动特性.考察单液滴在不同气液相密度比时撞击液膜的发展.随着密度比的减小,冠状水花顶端开始向内弯曲,且底部半径显著减小.在大密度比情况下研究双液滴撞击液膜.结果表明：双液滴撞击液膜有中心射流的形成;液滴水平间距的增大,延缓中心射流的出现,并降低初始中心射流的高度;随Re数的增加,中心射流的高度明显增大.     

Analysis of Lattice Boltzmann Initialization Routines

LB simulations can be affected by the arising of initial layers due to an inconsistent initialization of the discrete LB populations. We present some previously proposed initialization routines built to overcome that problem; using the asymptotic expansion technique, we show how their features can be analyzed and, in some cases, how accuracy and computational efficiency can be improved