用格子Boltzmann模型模拟非等温流场

根据微观和宏观之间的质量、动量、能量守恒准则和在原格子Boltzmann模型基础上,建立了几个新的格子Boltzmann模型,使得在外力场中的格子Boltzmann模型得到进一步完善.通过还原宏观流体力学方程,捕捉到了浮力强迫系数与Grashof数之间的关系.所得动量方程和Navier Stokes方程相比,在黏性输运项上有明显的改进,说明黏性应力不但与流体的速度梯度和流体的压缩性有关,而且还与非定常的内能梯度和动量通量有关.该模型对非等温流场的数值结果证明了其具有很好的数值稳定性和适用性. 关键词： Boltzmann模型 平衡分布函数 流体力学方程     

A Conceptual Model of Somali Jet Based on the Biot--Savart Law

We study the relationship between thermal conditions on the earth bottom boundary and the formation of Somali jet based on the Biot-Savart law and the data from National Centres for Environmental Prediction （NCEP）. As the radiation from the Sun gradually moves from the southern meridian, the temperature on the surface of Somali Peninsular and Arabic Peninsular gradually increases. During the same period the surface temperature of the Northern Indian Ocean increases much slower. It is shown that this increase of the temperature difference between the land and sea is inductive to the formation and development of Rayleigh-Benard convection and leads to the increasing relative vorticity strength between positive and negative vertical vortices over the land and sea. According to the Biot-Savart law, increase of vorticity strength will correspondingly induce the horizontal velocity. A pair of positive and negative vorticity fields over the two Peninsulars and the sea surface is effective in forming and maintaining this current. This mechanism is referred to as the ＇Somali suction pump＇. It draws air continually from the Southern hemisphere and releases it at the coastal area of Somali.     

格子Boltzmann模型平衡分布边界条件和多棱柱排列渗流流场的数值模拟

介绍了适用于多种流场数值模拟的无滑动格子Boltzmann平衡分布边界条件，这一边界条件是以Bounce-Back方法为基础并满足质量、动量守恒的准则.数值计算结果表明平衡分布边界条件克服了Bounce-Back方法在边界上所产生的滑动速度误差效应.利用平衡分布边界条件数值模拟了由棱柱形充填粒子构成的微尺度渗流流场中的Darcy-Forchheimer方程，通过与Lee 和Yang的数值结果比较，该预测结果是足够可靠的. 关键词： 平衡分布边界条件 渗流介质 Darcy-Forchheimer阻力     

SOME PROGRESS IN THE LATTICE BOLTZMANN MODEL

A lattice Boltzmann equation model has been developed by using the equilibrium distribution function of the Maxwell－Boltzmann-like form, which is third order in fluid velocity u_α. The criteria of energy conservation between the macroscopic physical quantities and the microscopic particles are introduced into the model, thus the thermal hydrodynamic equations containing the effect of buoyancy force can be recovered in terms of the Taylor and Chapman－Enskog asymptotic expansion methods. The two-dimensional thermal convection phenomena in a square cavity and between two concentric cylinders have been calculated by implementing a heat flux boundary condition. Both numerical results are in good agreement with the conventional numerical results.     

LATTICE BOLTZMANN EQUATION MODEL IN THE CORIOLIS FIELD

In a large-scale field of rotational fluid, various unintelligible and surprising dynamic phenomena are produced due to the effect of the Coriolis force. The lattice Boltzmann equation (LBE) model in the Coriolis field is developed based on previous works.^[1-4] Geophysical fluid dynamics equations are derived from the model. Numerical simulations have been made on an ideal atmospheric circulation of the Northern Hemisphere by using the model and they reproduce the Rossby wave motion well. Hence the applicability of the model is verified in both theory and experiment.     

一个完全守恒的格子Boltzmann模型

利用六角形格子离散的方法,使得每个六角形Cell里含有 3种不同运动速度的粒子微团,宏观物理量用这些粒子微团的矩来定义.根据微观和宏观之间的质量、动量、能量守恒准则,建立了一个二维的D2Q19格子Boltzmann模型,从该D2Q19模型出发可推导出宏观的流体力学方程组.用该模型对冲击波在障碍物组表面上的折射现象进行数值模拟,并将计算结果与实验进行了比较.     

CORRECTIONS TO THE COLLISION TERM IN THE BGK BOLTZMANN EQUATION

With the discrete method of the hexagonal cell and three different velocities of particle population in each cell, a two-dimensional lattice Boltzmann model is developed in this paper.^[1,2] The collision operator in the Boltzmann equation is expanded to fourth order using the Taylor expansion.^[3,4] With this model, good results have been obtained from the numerical simulation of the reflection phenomenon of the shock wave on the surface of an obstacle, and the numerical stability is also good. Thus the applicability of the D2Q 19 model is verified.     

Lattice Boltzmann Model and Geophysical Hydrodynamic Equation

A lattice Boltzmann equation model in a rotating system is developed by introducing the Coriolis force effect.The geophysical hydrodynamic equation can be derived from this model.Numerical computations are performed to simulate the cylindrical annulus experiment and Benard convection.The numerical results have shown the flow behaviour of large-scale geostrophic current and Benard convection cells,which verifies the applicability of this model to both theory and experiment.     

Lattice Boltzmann Numerical Simulation of a Circular Cylinder

The lattice Boltzmann equation (LBE) model based on the Boltzmann equation is suitable for the numerical simulation of various flow fields.The fluid dynamics equation can be recovered from the LEB model.Howeverl,compared to the Navier-Stokes transport equation,the fluid dynamics equation derived from the LBE model is somewhat different in the viscosity transport term,which contains not only the Navier-Stokes transport equation but also nonsteady pressure and momentum flux terms.The two nonsteady terms can produce the same function as the random stirring force term introduced in the direct numerical or large-eddy vortex simulation of turbulence.Through computation of a circular cylinder,it is verified that the influence of the two nonsteady terms on flow field stability cannot be ignored,which is helpful for the study of turbulence.     

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.