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.     

Lattice Boltzmann Solid Particles in a Lattice Boltzmann Fluid

We define a lattice Boltzmann model of solid, deformable suspensions immersed in a fluid itself described in terms of the lattice Boltzmann method. We discuss the rules governing the internal dynamics of the solid object as well as the rules specifying the interaction between solid and fluid particle. We perform a numerical drag experiment to validate the model. Finally we consider the case of a population of flexible chains in suspension in a shear stress flow and study the influence on the velocity profile.     

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     

Lattice Boltzmann simulations of a dumbbell moving in a Poiseuille flow

In this paper, the lattice Boltzmann method is applied to simulate a dumbbell moving in a pressure-driven flow in a planar channel with the stress-integration method for the evaluation of hydrodynamic force acting on the cylinders. The simulation results show that the dumbbell also has the important feature of the Segr\'e--Silberberg effect like a particle in a Poiseuille flow. The dumbbell trajectories, orientations, the cylinders vertical velocities and angular velocities all reach their equilibrium values separately independent of their initial positions. It is also found that the dumbbell equilibrium positions depend on the flow Reynolds number, blockage ratio and elastic coefficient. This study is expected to be helpful to understand the dynamics of polymer solutions, polymer synthesis and reaction, etc.     

Extended Lattice Boltzmann Model

Conventional lattice Boltzmann models for the simulation of fluid dynamics are restricted by an error in the stress tensor that is negligible only for small flow velocity and at a singular value of the temperature. To that end, we propose a unified formulation that restores Galilean invariance and the isotropy of the stress tensor by introducing an extended equilibrium. This modification extends lattice Boltzmann models to simulations with higher values of the flow velocity and can be used at temperatures that are higher than the lattice reference temperature, which enhances computational efficiency by decreasing the number of required time steps. Furthermore, the extended model also remains valid for stretched lattices, which are useful when flow gradients are predominant in one direction. The model is validated by simulations of two- and three-dimensional benchmark problems, including the double shear layer flow, the decay of homogeneous isotropic turbulence, the laminar boundary layer over a flat plate and the turbulent channel flow.     

Lattice Boltzmann simulations of flow past a cylindrical obstacle

We present lattice Boltzmann simulations of flow past a cylindrical obstacle. Our study is based on the Lévy walk model of turbulence in a lattice Boltzmann model. We discuss pressure around the cylinder with laminar and turbulent incident flows, as well as the dependence of the von Karman street on the analog of integral scale in our model.     

Towards a Renormalized Lattice Boltzmann Equation for Fluid Turbulence

A coarse-grained Lattice Boltzmann equation is examined in which the effects of unresolved (subgrid) scales are formally incorporated within a renormalized relaxation time of the collision operator. Actual values of the renormalized relaxation time are analyzed for the practical case of high-Reynolds flows past slant bodies (airfoils).     

用晶格玻尔兹曼方法研究颗粒在涡流中的运动

用晶格玻尔兹曼方法研究小颗粒在涡流中的运动.涡流由流经空腔的流体产生.用动量交换法和压力张量积分法计算颗粒在涡流中的运动轨迹、速度和角速度.最后用张量积分法计算两个不同半径的颗粒在涡流中的运动. 关键词： 晶格玻尔兹曼方法 涡流 颗粒     

振动纤维捕集颗粒的格子Boltzmann模拟

采用格子Boltzmann方法对振动纤维捕集颗粒进行数值研究.纤维附近采用多块网格加细技术计算,颗粒采用Lagrange跟踪方法模拟.研究雷诺数为200的流向振动纤维绕流的AⅡ、AⅢ、AIV、S四种涡结构下的亚微米煤粉颗粒的捕集问题.结果表明纤维的流向强迫振动能够显著提高颗粒的捕集效率.且迎风面的捕集效率提升不大,背风面的提升则更为显著.颗粒撞击角度的统计反映了背风面捕集效率提升的细节.另外流动处于AⅢ模态时,每个周期内脱落两个正涡一个负涡,被捕集颗粒的初始位置分布不对称.而其它模态基本关于流场中心对称.     

方腔内微细颗粒物运动特性的格子Boltzmann方法模拟

对Masselot and Chopard提出的模拟气固两相流动的格子Boltzmann-格子气(LBE-LGA)方法进行推广,能够反映两相间拖曳作用.利用该方法研究封闭方腔内的气固两相流运动特性,分析斯托克斯数St和模拟颗粒数目对颗粒群运动的影响,并与文献结果进行比较,表明LBE-LGA方法模拟颗粒运动是可行的.     

Lattice Boltzmann method for three-dimensional moving particles in a Newtonian fluid

A lattice Boltzmann method is developed to simulate three-dimensional solid particle motions in fluids. In the present model, a uniform grid is used and the exact spatial location of the physical boundary of the suspended particles is determined using an interpolation scheme. The numerical accuracy and efficiency of the proposed lattice Boltzmann method is demonstrated by simulating the sedimentation of a single sphere in a square cylinder. Highly accurate simulation results can be achieved with few meshes, compared with the previous lattice Boltzmann methods. The present method is expected to find applications on the flow systems with moving boundaries, such as the blood flow in distensible vessels, the particle－flow interaction and the solidification of alloys.     

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.     

椭圆柱体在牛顿流体中运动的格子Boltzmann方法模拟

用格子Boltzmann方法建立了椭圆柱体的二维动力学模型，利用所建立的模型，数值模拟了 牛顿流体中不同形状的椭圆柱体在相同初始条件下的运动和同一椭圆柱体在不同初始条件下 的运动，并通过比较相同条件下圆柱体的运动，讨论了椭圆柱体二维运动的特征，得到了一 些有意义的结果. 关键词： 格子Boltzmann方法 椭圆柱体 牛顿流体     

辐射传输的格子Boltzmann模拟

辐射动力学理论是描述辐射传输是重要手段,基于此,本文建立了辐射能和辐射动量守恒方程,并基于Chapman-Enskog多尺度展开方法实现了从辐射传输Boltzmann方程到宏观方程的推导,进而建立了适于一维辐射传输的2分量格子Boltzmann模型。数值结果与精确解吻合较好,表明本文提出的LBM方法具有很好的准确性和稳...     

微通道内气体流动的格子Boltzmann方法研究

采用格子Boltzmann方法模拟了微通道在滑移区内不同Knudsen数下的微气体Poiseuille流,分析了微气体流动的速度分布以及流量与压降的关系,并给出了相对滑移长度和Poiseuille数随Knudsen数的变化特性。研究结果表明,微气体Poiseuille流的速度轮廓呈抛物线分布,但是边界速度大于0,出现速...     

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.     

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

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

非等温流向振荡单圆柱绕流的格子Boltzmann模拟

本文采用格子Boltzmann方法模拟对非等温横向振荡圆柱绕流(Re=100)进行模拟,详细考察振荡频率(F=f_o/f_s=0.2～8.0)、振幅(A/D=0.1～0.4,A为圆柱流向最大位移,D为圆柱直径)对绕流换热的影响。结果表明:低频率、低振幅流向振荡运动将使圆柱换热增强,高频率、高振幅振荡使圆柱换热效果减弱;圆柱表面平均N_u脉动均方根值随着振荡频率、振幅的增大而增大;在锁定区域内,平均N_u及其脉动均方根值均出现极大值。     

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.