一种滑移区气体流动的格子Boltzmann曲边界处理新格式

针对滑移区复杂气-固边界存在速度滑移现象,提出了一种基于格子Boltzmann方法的非平衡态外推与有限差分相结合的曲边界处理新格式.该格式具有可考虑实际物理边界与网格线偏移量的优势,较传统half-way DBB(diffusive bounce-back)格式更能准确反映实际边界情况,同时还可获取壁面处气体宏观量及其法向梯度等信息.采用本文所提曲边界处理格式模拟分析了滑移区气体平直/倾斜微通道Poiseuille流、微圆柱绕流和同心微圆柱面旋转Couette流问题.研究结果表明,采用曲边界处理新格式所得结果与理论值以及文献结果符合良好,适用于滑移区气体流动的复杂边界处理,且比half-way DBB格式具有更高的精度,较修正DBB格式具有更好的适应性.     

Local second-order boundary methods for lattice Boltzmann models

A new way to implement solid obstacles in lattice Boltzmann models is presented. The unknown populations at the boundary nodes are derived from the locally known populations with the help of a second-order Chapman-Enskog expansion and Dirichlet boundary conditions with a given momentum. Steady flows near a flat wall, arbitrarily inclined with respect to the lattice links, are then obtained with a third-order error. In particular, Couette and Poiseuille flows are exactly recovered without the Knudsen layers produced for inclined walls by the bounce back condition.     

疏水表面滑移流动及减阻特性的格子Boltzmann方法模拟

采用格子Boltzmann方法研究了固体壁面对流体的作用强度与其润湿性的关系,在此基础上进一步模拟了疏水表面微通道内的流体流动,获得了润湿性对疏水表面滑移流动及减阻特性的影响规律,证实了疏水表面表观滑移的存在性并揭示了其产生机理.结果表明,疏水性作用在疏水表面的近壁区诱导了一个低密度层,而表观滑移则发生在低密度层上.表观滑移是疏水表面具有减阻作用的直接原因,减阻效果随滑移长度的增大而增大.对于特定的流体系统,滑移长度是疏水表面的固有属性,仅是壁面润湿性的单一函数,而与流动本身的性质无关.     

The lattice Boltzmann phononic lattice solid

I present a Boltzmann lattice gas-like approach for modeling compressional waves in an inhomogeneous medium as a first step toward developing a method to simulate seismic waves in complex solids. The method is based on modeling particles in a discrete lattice with wavelike characteristics of partial reflection and transmission when passing between links with different properties as well as phononlike interactions (i.e., collisions), with particle speed dependent on link properties. In the macroscopic limit, this approach theoretically yields compressional waves in an inhomogeneous acoustic medium. Numerical experiments verify the method and demonstrate its convergence properties. The lattice Boltzmann phononic lattice solid could be used to study how seismic wave anisotropy and attenuation are related to microfractures, the complex geometry of rock matrices, and their couplings to pore fluids. However, additional particles related to the two transverse phonons must be incorporated to correctly simulate wave phenomena in solids.     

Gas Flow in Microchannels – A Lattice Boltzmann Method Approach

Gas flow in microchannels can often encounter tangential slip motion at the solid surface even under creeping flow conditions. To simulate low speed gas flows with Knudsen numbers extending into the transition regime, alternative methods to both the Navier–Stokes and direct simulation Monte Carlo approaches are needed that balance computational efficiency and simulation accuracy. The lattice Boltzmann method offers an approach that is particularly suitable for mesoscopic simulation where details of the molecular motion are not required. In this paper, the lattice Boltzmann method has been applied to gas flows with finite Knudsen number and the tangential momentum accommodation coefficient has been implemented to describe the gas-surface interactions. For fully-developed channel flows, the results of the present method are in excellent agreement with the analytical slip-flow solution of the Navier–Stokes equations, which are valid for Knudsen numbers less than 0.1. The present paper demonstrates that the lattice Boltzmann approach is a promising alternative simulation tool for the design of microfluidic devices.     

The effect of surface roughness on rarefied gas flows by lattice Boltzmann method

This paper studies the roughness effect combining with effects of rarefaction and compressibility by a lattice Boltzmann model for rarefied gas flows at high Knudsen numbers. By discussing the effect of the tangential momentum accommodation coefficient on the rough boundary condition, the lattice Boltzmann simulations of nitrogen and helium flows are performed in a two-dimensional microchannel with rough boundaries. The surface roughness effects in the microchannel on the velocity field, the mass flow rate and the friction coefficient are studied and analysed. Numerical results for the two gases in micro scale show different characteristics from macroscopic flows and demonstrate the feasibility of the lattice Boltzmann model in rarefied gas dynamics.     

一个新的用于不可压磁流体动力学的格子波尔兹曼模型

绝大多数现有的格子波尔兹曼磁流体动力学模型其实是用可压缩方法来模拟不可压磁流体。而这些可压缩效应在数值模拟中往往会带来意想不到的误差。在这篇文章中，我们提出了一个全新的可用于的不可压格子波尔兹曼磁流体动力学模型，并且进行了哈特曼流的数值模拟。模拟结果与哈特曼流的解析解非常吻合。这个方法需要一个假设条件来消除误差。我们做了大量的数值试验，并且与Dellar教授的模型进行了详细的分析与比较。     

Boundary conditions for lattice Boltzmann simulations

A heuristic interpretation of no-slip boundary conditions for lattice Boltzmann and lattice gas simulations is developed. An improvement is suggested which consists of including the wall nodes in the collision operation.     

Numerical analysis of fluid flow through a cylinder array using a lattice Boltzmann model

In this paper we present a detailed computational study of an incompressible Newtonian fluid flow across a periodic array of two-dimensional cylinders which is a simplest non-trivial representation of a porous media. A two-dimensional Lattice Boltzmann Method is used to solve the governing Navier－Stokes equation taking into account of viscous dissipation effects and influence of nonlinear fluid drag. Both the flow fields and the Darcy－Forchheimer drag coefficient as a function of the solid volume fraction are calculated for a wide range of flow Reynolds numbers. The predictions were compared with the results from conventional numerical and empirical models for verification. Apart from confirming that inertial effects can cause a significant deviation from Darcy's law for large velocities the results also show that the characteristics of the vorticity field vary considerably as the Reynolds number increases, which will have major implications to the transport of passive particulate substances within the pores and their removal rate.     

On the generalized Hermite-based lattice Boltzmann construction,lattice sets,weights, moments,distribution functions and high-order models

The influence of the use of the generalized Hermite polynomial on the Hermite-based lattice Bofiz- mann （LB） construction approach, lattice sets, the thermal weights, moments and the equilibrium distribution function （EDF） are addressed. A new moment system is proposed. The theoretical possibility to obtain a unique high-order Hermite-based singel relaxation time LB model capable to exactly inatch some first hydrodynamic inoments thermally i） on-Cartesian lattice, ii） with thermal weights in the EDF, iii） whilst the highest possible hydrodynamic moments that are exactly reatched are obtained with the shortest on-Cartesian lattiee sets with some fixed real-valued temperatures, is also analyzed.     

声悬浮过程的格子Boltzmann方法研究

采用轴对称多弛豫时间格子Boltzmann(LB)方法,研究了圆柱形封闭谐振腔中圆盘形样品的声悬浮过程.模拟结果表明,(001)模式下谐振腔的共振长度L=0.499λ,在谐振腔中心引入样品后共振漂移量δL≈-0.9,这与线性声学理论计算结果基本相符.声悬浮力的LB模拟过程包含了黏滞性效应和共振漂移效应,所获得的模拟结果与理论公式计算值在量值上一致,而且其在细节上更符合实验现象.此外,LB模拟还揭示出了声悬浮过程中的声压波形畸变、声流和声辐射压等非线性声学效应.     

Boltzmann equation on a lattice: Existence and uniqueness of solutions

Cercignani, Greenberg, and Zweifel proved the existence and uniqueness of solutions of the Boltzmann equation on a toroidal lattice under the assumption that the collision kernel is bounded. We give an alternative, considerably simpler, proof which is based on a fixed point argument.     

用格子Boltzmann方法数值模拟三维空化现象

空化是一种微观、瞬时、随机、多相的复杂现象,其过程中所产生的极端条件以及伴随的一系列空化效应,将对液流系统产生破坏性和建设性两方面的作用.采用基于Shan-Chen模型的单组分多相流格子Boltzmann方法对水体中的三维空化现象进行了数值模拟,研究了低压下水体中气核半径与空化现象的相互关系,成功再现了低压下水体中微小气核发展成气泡的过程,并进一步研究了水体依次流经低压区、高压区时空化产生、发展、溃灭的全过程.数值模拟结果和理论预测结果符合良好. 关键词： 单组分多相流 格子Boltzmann方法 三维空化     

A lattice Boltzmann representation of fluid closures for Landau damping

A lattice Boltzmann representation of the zeroth moment density closure of Hammett et al. is shown to reproduce Landau exponential decay of a perturbed density profile provided five streaming velocities are included in the lattice model. By contrast, the three streaming velocity model displays numerical instabilities when the dissipation coefficient is reduced to a level that Landau damping dominates the asymptotic density decay.     

微尺度振荡Couette流的格子Boltzmann模拟

采用有效多松弛时间-格子Boltzmann方法(Effective MRT-LBM)数值模拟了微尺度条件下的振荡Couette和Poiseuille流动. 在微流动LBM中引入Knudsen边界层模型,对松弛时间进行修正. 模拟时平板或外力以正弦周期振动,Couette流中考虑了单平板振动、上下板同相振动这两类情况. 研究结果表明,修正后的MRT-LBM模型能有效用于这类非平衡的微尺度流动模拟;对于Couette流,随着Kn数的增大,壁面滑移效应变得越明显. St越大,板间速度剖面的非线性特性越剧烈;两板同相振荡时,若Kn,St均较小,板间流体受到平板拖动剪切的影响很小,板间速度几乎重叠在一起;在振荡Poiseuille流动中,St数增大到一定值时,相位滞后现象减弱;相对于Kn数,St数对振荡Couette 和Poiseuille流中不同位置处速度相位差的产生有较大影响. 关键词： 格子Boltzmann方法 有效MRT模型 Knudsen层 振荡流     

Simulating high Reynolds number flow in two-dimensional lid-driven cavity by multi-relaxation-time lattice Boltzmann method

By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50,000 in the literature; however, this paper can readily extend the maximum to 1000,000 with the above combination.     

Derivation of Slip boundary conditions for the Navier-Stokes system from the Boltzmann equation

Rarefied gas flow behavior is usually described by the Boltzmann equation, the Navier-Stokes system being valid when the gas is less rarefied. Slip boundary conditions for the Navier-Stokes equations are derived in a rigorous and systematic way from the boundary condition at the kinetic level (Boltzmann equation). These slip conditions are explicitly written in terms of asymptotic behavior of some linear half-space problems. The validity of this analysis is established in the simple case of the Couette flow, for which it is proved that the right boundary conditions are obtained.     

Thermal equation of state for lattice Boltzmann gases

The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar--Gross--Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invariance, recovering the Galilean invariance for the isothermal lattice Boltzmann Bhatnagar--Gross--Krook equation (LBGKE) induces a new natural thermal-dynamical system, which is compatible with the elementary statistical thermodynamics.     

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.     

纳米结构超疏水表面冷凝现象的三维格子玻尔兹曼方法模拟

采用三维多相流格子玻尔兹曼方法 (lattice Boltzmann method, LBM),对纳米结构超疏水表面液滴的冷凝行为进行模拟研究.通过Laplace定律和光滑表面的本征接触角理论对三维LBM模型进行定量验证.模拟分析了超疏水表面纳米阵列的几何尺寸和润湿性的局部不均匀性对冷凝液滴形核位置和最终润湿状态的影响规律.结果表明,较高的纳米阵列使液滴在纳米结构间隙的上部侧面和底部优先形核长大,通过采用上下不均匀的间隙可避免液滴在底部形核长大,而在上部侧面形核的冷凝液滴在生长过程中向上运动,其润湿状态由Wenzel态转变为Cassie态;较低的纳米阵列使液滴在纳米结构底部优先形核长大,液滴的最终润湿状态为Wenzel态;润湿性不均匀的纳米结构表面使液滴在阵列顶端亲水位置处优先形核长大,成为Cassie态.冷凝液滴在不同几何尺寸的纳米结构表面上的最终润湿状态的模拟结果与文献报道的实验结果符合良好.通过模拟还发现,冷凝液滴在生长过程中的运动行为与液滴统计平均作用力的变化有关.本文的LBM模拟再现了三维空间中液滴的形核、长大和润湿状态转变等物理现象.