并列双椭圆柱绕流的格子Boltzmann-虚拟区域方法的模拟研究

采用格子Boltzmann-虚拟区域方法对并列双椭圆柱绕流进行了模拟研究。首先,通过与并列双圆柱的结果进行对比,验证了数值方法的有效性。其次对雷诺数为200时两种间距(g=0.5a和2.0a,g为柱体表面间距,a为椭圆长轴)的情况进行了研究,考察了椭圆长轴与短轴之比,即α值对柱体升阻力系数及涡结构的影响。研究发现,与圆柱相比,对于g=0.5a椭圆柱的升阻力系数可能出现两种变化,一是升阻力随时间演化较规则,接近周期性;二是流场可能长时间偏向于其中一个椭圆柱,这些变化与α的值有关。对于g=2.0a,两个椭圆柱后的某一区域内会出现四列涡街,经过一段时间,四列涡街又会演化成两列向两侧扩张的涡街。     

双圆柱绕流特性的模拟研究

采用格子Boltzmann方法对低雷诺数下气体绕流圆柱的规律进行了研究. 对比计算了双圆柱在不同圆心距、不同Re数、不同来流速度与双圆柱圆心连线角度的情况下,各个圆柱的受力大小和曳力系数. 结果表明,若Re数为20, 改变圆柱间距,圆柱间距在1.2d和1.4d之间时,下游圆柱所受曳力有极小值;双圆柱间距为1.6d时,双圆柱受到总曳力最小;圆柱间距大于2d时,上游颗粒受到的曳力不再受到下游颗粒的影响. 若圆柱间距为1.2d, 改变雷诺数,Re数在30和40之间,下游圆柱所受曳力有极小值. 另外,来流速度角度对圆柱的受力影响很大. 上述规律为低Re数下圆柱绕流的深入研究与应用打下基础.      

并列圆柱绕流的格子Boltzmann数值模拟

采用非均匀不可压格子Boltzmann模型对低雷诺数下并列圆柱绕流进行了数值模拟,给出了数值计算结果,分析了间距g对圆柱尾流及升力、阻力的影响,并在此基础上得到了4种尾迹模式.此外,研究了流场的初始扰动对流动分岔现象的影响,发现在适当的扰动下可以很快得到同步同相的尾流.对Re=160和200下圆柱的升、阻力进行了对比,结果表明升力和阻力受间距g的影响大于雷诺数.     

格子Boltzmann方法模拟多孔介质惯性流的边界条件改进

格子Boltzmann方法可以有效地模拟水动力学问题,边界处理方法的选择对于可靠的模拟计算至关重要.本文基于多松弛时间格子Boltzmann模型开展了不同边界条件下,周期对称性结构和不规则结构中流体流动模拟,阐述了不同边界条件的精度和适用范围. 此外,引入一种混合式边界处理方法来模拟多孔介质惯性流, 结果表明:对于周期性对称结构流动模拟,体力格式边界条件和压力边界处理方法是等效的,两者都能精确地捕捉流体流动特点; 而对于非周期性不规则结构,两种边界处理方法并不等价,体力格式边界条件只适用于周期性结构;由于广义化周期性边界条件忽略了垂直主流方向上流体与固体格点的碰撞作用,同样不适合处理不规则模型;体力-压力混合式边界格式能够用来模拟周期性或非周期性结构流体流动,在模拟多孔介质流体惯性流时,比压力边界条件有更大的应用优势,可以获得更大的雷诺数且能保证计算的准确性.      

利用格子Boltzmann方法数值模拟Rayleigh-Benard对流

采用格子Boltzmann方法对较大Rayleigh数范围下的二维Rayleigh-Benard对流进行了模拟研究.引入能量分布函数,利用该能量分布函数与粒子速度分布函数耦合来求解一个热流场,能量分布函数与粒子速度分布函数和Boltzmann方程构成了一个新的双分布格子Boltzmann模型.在考虑密度随温度变化的情况下,进行数值模拟,得到了Rayleigh-Benard对流速度、温度随时间的变化规律、系统的流线和等温线分布及平均Nusselt数与Rayleigh数的之间的关系,与相关文献数据进行了对比,模拟结果非常吻合,证明了改进的双分布格子Boltzmann模型的有效性.     

倾斜壁面浸润性梯度驱动液滴运动过程的格子Boltzmann模型研究

基于介观模型的多组分伪势格子Boltzmann方法,模拟了倾斜壁面浸润性梯度驱动液滴的运动过程,研究了壁面浸润性梯度、壁面倾斜角度对液滴运动过程的影响．结果表明,对于一定倾斜角度的壁面,当壁面上浸润性梯度足够大时,液滴能够克服重力的作用实现“爬坡”;液滴在运动过程中,其前进及后退接触角与当地静态接触角间存在差值;增大壁面浸润性梯度时,液滴能够获得更快的加速,并且前进及后退接触角与当地静态接触角之间的差值也随之增大;增大壁面倾斜角度时,液滴的运动受到阻碍,前进及后退接触角与当地静态接触角的差值小幅减小．     

幂律流体圆柱绕流的格子波尔兹曼模拟

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明:插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.      

利用格子Boltzmann方法模拟单个气泡在复杂流道内的运动特性

利用格子Boltzmann方法模拟了单个气泡在具有三个半圆形喉部的复杂流道内的上升过程.通过分析气泡运动过程中的形态、运动轨迹及运动速度的变化,研究复杂流道对气泡运动特性的影响.研究结果表明,在上升过程中,由于壁面的影响,气泡的形状发生严重的变形,运动轨迹也发生相应的偏转.通过实验结果的对比,验证了模拟结果的正确性.结果表明格子Boltzmann方法可以用于模拟具有复杂边界的两相流问题.     

并列双方柱绕流的Lattice Boltzmann模拟分析

采用Lattice Boltzmann方法对并列双方柱绕流问题进行数值模拟.对方柱间距比s/D=0.2、0.3、0.4、0.5、1.0、1.4、1.5、1.6、1.8、2.0、2.5共11种情况下的流场分别进行计算,给出了对应的流线图、方柱的升力图及阻力图,同时计算了方柱各边中点的压力,给出了压力与间距比的变化曲线,并对各个流线图、升力图、阻力图及压力图进行了分析讨论.结果表明:当s/D≤1.5时,流动呈偏流型,形成周期漩涡;当s/D＞1.5时,流动呈对称型,形成对称漩涡;s/D=1.5为流动从偏流型向对称型转换的临界间距比.该结果对工程中建筑物位置的安排及其安全系数的提高有一定的参考价值.     

建筑物周围流场格子Boltzmann仿真

基于格子BGK方程的三维15点格子Boltz-mann模型,对建筑物周围流场进行了数值仿真,对其流场的划分进行了分析,并对不同雷诺数下建筑物周围的流场进行了数值模拟计算,得到了比较满意的结果,为环境质量的评估提供了新的途径.     

复杂微通道气体流动的格子Boltzmann模拟

构建了一个模拟复杂微通道内气体流动的多松弛格子Boltzmann模型。该模型采用动力学曲面滑移边界,考虑了微尺度效应和努森层影响。此外,为了更准确地描述微通道内气体的滑移速度,在模型中引入孔隙局部Kn数来代替平均Kn数。之后采用Poiseuille流对模型进行验证,模拟结果与用直接模拟蒙特卡洛方法和分子模拟结果吻合较好,证明了该模型模拟微通道内处于滑移区和过渡区气体流动的有效性。最后,采用该模型模拟多孔介质内气体渗流过程。结果表明,随着孔隙平均Kn数的增加,多孔介质内的高渗区域增加,且优先从小孔隙中开始增加,这是由于小孔隙中微尺度效应更加明显,相对大孔隙流动阻力更小所致。     

Lattice Boltzmann simulations of turbulent shear flow between parallel porous walls

The effects of two parallel porous walls are investigated, consisting of the Darcy number and the porosity of a porous medium, on the behavior of turbulent shear flows as well as skin-friction drag. The turbulent channel flow with a porous surface is directly simulated by the lattice Boltzmann method （LBM）. The Darcy-Brinkman- Forcheimer （DBF） acting force term is added in the lattice Boltzmann equation to simu- late the turbulent flow bounded by porous walls. It is found that there are two opposite trends （enhancement or reduction） for the porous medium to modify the intensities of the velocity fluctuations and the Reynolds stresses in the near wall region. The parametric study shows that flow modification depends on the Darcy number and the porosity of the porous medium. The results show that, with respect to the conventional impermeable wall, the degree of turbulence modification does not depend on any simple set of param- eters obviously. Moreover, the drag in porous wall-bounded turbulent flow decreases if the Darcy number is smaller than the order of O（10-4） and the porosity of porous walls is up to 0.4.     

The simulation of turbulent particle‐laden channel flow by the Lattice Boltzmann method

We perform direct numerical simulation of three‐dimensional turbulent flows in a rectangular channel, with a lattice Boltzmann method, efficiently implemented on heavily parallel general purpose graphical processor units. After validating the method for a single fluid, for standard boundary layer problems, we study changes in mean and turbulent properties of particle‐laden flows, as a function of particle size and concentration. The problem of physical interest for this application is the effect of water droplets on the turbulent properties of a high‐speed air flow, near a solid surface. To do so, we use a Lagrangian tracking approach for a large number of rigid spherical point particles, whose motion is forced by drag forces caused by the fluid flow; particle effects on the latter are in turn represented by distributed volume forces in the lattice Boltzmann method. Results suggest that, while mean flow properties are only slightly affected, unless a very large concentration of particles is used, the turbulent vortices present near the boundary are significantly damped and broken down by the turbulent motion of the heavy particles, and both turbulent Reynolds stresses and the production of turbulent kinetic energy are decreased because of the particle effects. We also find that the streamwise component of turbulent velocity fluctuations is increased, while the spanwise and wall‐normal components are decreased, as compared with the single fluid channel case. Additionally, the streamwise velocity of the carrier (air) phase is slightly reduced in the logarithmic boundary layer near the solid walls. Copyright © 2015 John Wiley & Sons, Ltd.     

基于格子Boltzmann两相模型的粘弹流体挤出胀大分析

挤出胀大的数值模拟是非牛顿流体研究中具有挑战性的问题．本文运用格子Boltzmann方法(LBM)分析Oldroyd-B和多阶松弛谱PTT粘弹流体的挤出胀大现象,采用颜色模型模拟出口处粘弹流体和空气的两相流动,通过重新标色获得两种流体的界面,并最终获得胀大的形状．Navier-Stokes方程和本构方程的求解采用双分布函数模型．将胀大的结果与解析解、实验解和单相自由面LBM结果进行了比较,发现格子Boltzmann两相模型结果与解析解和实验结果相吻合,相比于单相模型,收敛速度更快,解的稳定性更高．研究了流道尺寸对胀大率的影响,并对挤出胀大的内在机理进行了分析．     

粘弹流体轴对称流动的格子Boltzmann方法

基于BGK碰撞模型，通过在迁移方程中引入作用力项，建立了粘弹流体的轴对称格子Boltzmann模型．通过Chapman-Enskog展开，获得了准确的柱坐标下轴对称宏观流动方程．采用双分布函数对运动方程和本构方程进行迭代求解，模拟分析了粘弹流体管道流动，获得了流场中的速度和构型张量的分布，通过与解析解进行比较，验证了模型的准确性．研究了作为粘弹流体流动基准问题的收敛流动，对涡旋位置进行了定量分析，将回转长度的计算结果与有限体积法进行了比较，两种数值结果十分吻合．研究结果表明，模型能够准确表征粘弹流体的轴对称流动，具有较广阔的应用前景．     

微通道人工泡状流的Lattice Boltzmann数值模拟

采用高频电控热激发汽泡的方式构造微通道人工泡状流,可以有效抑制微通道沸腾流动的不稳定性和强化传热。本文基于Lattice Boltzmann大密度比多相流复合模型,数值研究了通道内人工泡状流的流动和传热,通过比较分析不同发泡频率的泡状流,量化分析了汽泡运动和增长对微通道流动与传热的相互影响。一方面着重分析了汽泡运动对微通道运动边界层以及汽泡相变增长对热边界层的影响,另一方面也研究了边界层对汽泡动力行为的影响,所得结论对研究抑制微通道沸腾流动不稳定性和强化传热有参考意义。     

Lattice Boltzmann simulation of two-sided lid-driven flow in a staggered cavity

In this article, the lattice Boltzmann method is employed in order to explore incompressible fluid flow inside a two-sided lid-driven staggered cavity. Results of the lattice Boltzmann simulation for antiparallel motion of lids are compared with the data from existing literature. For parallel motion of lids, the characteristics of flow pattern for a variety of Re numbers (50–3200) are presented. An asymmetric steady-state flow pattern for parallel motion of lids is obtained.     

Numerical simulation of laminar jet-forced flow using lattice Boltzmann method

In the paper, a numerical study on symmetrical and asymmetrical laminar jet-forced flows is carried out by using a lattice Boltzmann method （LBM） with a special boundary treatment. The simulation results are in very good agreement with the available numerical prediction. It is shown that the LBM is a competitive method for the laminar jet-forced flow in terms of computational efficiency and stability.     

Lattice Boltzmann simulation of flow past a circular cylinder near a moving wall

The two‐dimensional flows past a circular cylinder near a moving wall are simulated by lattice Boltzmann method. The wall moves at the inlet velocity and the Reynolds number ranges from 300 to 500. The influence of the moving wall on the flow patterns is demonstrated and the corresponding mechanism is illustrated by using instability theory. The correlations among flow features based on gap ratio are interpreted. Force coefficients, velocity profile and vortex structure are analyzed to determine the critical gap ratio. Copyright © 2011 John Wiley & Sons, Ltd.     

Lattice Boltzmann methods for shallow water flow applications

We apply the lattice Boltzmann (LB) method for solving the shallow water equations with source terms such as the bed slope and bed friction. Our aim is to use a simple and accurate representation of the source terms in order to simulate practical shallow water flows without relying on upwind discretization or Riemann problem solvers. We validate the algorithm on problems where analytical solutions are available. The numerical results are in good agreement with analytical solutions. Furthermore, we test the method on a practical problem by simulating mean flow in the Strait of Gibraltar. The main focus is to examine the performance of the LB method for complex geometries with irregular bathymetry. The results demonstrate its ability to capture the main flow features. Copyright © 2007 John Wiley & Sons, Ltd.