表征体元尺度渗流的离散统一动理学格式

结合表征体元尺度的通用渗流模型，提出离散统一动理学格式（DUGKS）渗流方法，分别用均匀网格和非均匀网格计算二维Poiseuille、Couette、方腔流等经典渗流问题，检验DUGKS渗流方法的有效性和非均匀网格应用的优势，将DUGKS渗流方法应用到裂缝系统中.     

基于格子Boltzmann方法的多孔介质流动模拟GPU加速

利用NVIDIA CUDA平台,在GPU上结合稀疏存贮算法实现基于格子Boltzmann方法的孔隙尺度多孔介质流动模拟加速,测试该算法相对基本算法的性能.比较该算法在不同GPU上使用LBGK和MRT两种碰撞模型及单、双精度计算时的性能差异.测试结果表明在GPU环境下采用稀疏存贮算法相对基本算法能大幅提高计算速度并节省显存,相对于串行CPU程序加速比达到两个量级.使用较新构架的GPU时,MRT和LBGK碰撞模型在单、双浮点数精度下计算速度相同.而在较上一代的GPU上,计算精度对MRT碰撞模型计算速度影响较大.     

基于区域分裂的非均匀Lattice Boltzmann方法

基于区域分裂技术,设计了一种使用非均匀网格的LatticeBoltzmann（LB）模型,其基本思想是：将流场分解为若干个规划的子区域,在每个子区域使用均匀网格上的LB模型求解,子区域的边界条件利用插值方法决定,用该模型对二维圆柱绕流问题进行了模拟,并与其它方法的结果作了比较。     

热对流条件下颗粒沉降的格子Boltzmann方法模拟

采用格子Boltzmann方法模拟了在热对流条件下的颗粒沉降问题, 在研究单颗粒在等温流体、热流体和冷流体中运动的基础上, 进一步模拟了两个不同温度的颗粒在流体中的沉降.结果表明:两等温颗粒的沉降方式与雷诺数Re以及格拉晓夫数Gr密切相关, 而两不同温度的颗粒与两等温颗粒的沉降规律有显著不同.无论初始位置如何, 冷颗粒最终总位于热颗粒下方运动, Re较大时, 发生连续的拖曳、接触现象, 而Re较小时, 冷颗粒会以较大的沉降速度远离热颗粒. 关键词： 格子Boltzmann方法 颗粒沉降 热对流     

多孔介质中流体流动及扩散的耦合格子Boltzmann模型

多孔介质中高Péclet数和大黏性比下混溶流体的流动和扩散广泛存在于二氧化碳驱油、化工生产等工业过程中.用数值方法对该问题进行研究时,关键在于如何正确描述高Péclet数和大黏性比下多孔介质内流体的行为.为此,提出了一种基于多松弛模型和格子动理模型的耦合格子Boltzmann模型.通过Chapman-Enskog分析,证明该模型能有效求解不可压Navier-Stokes方程和对流扩散方程.数值结果表明,该模型不仅具有二阶精度和良好的稳健性,而且对于高Péclet数和大黏性比的问题具有良好的数值稳定性,为模拟此类问题提供了有效工具.     

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

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

多孔介质中化学反应对非等粘流体混合过程影响的格子Boltzmann研究

利用格子Boltzmann方法和GPU计算技术,在孔隙尺度上模拟多孔介质中包含界面化学反应的粘性指进现象,定量分析化学反应对流体混合的影响.采用单浓度变量的双稳态模型来描述界面反应,而各向同性的多孔介质则通过四参数法生成.研究发现化学反应能减小指进界面厚度,抑制流体的混合,甚至会出现反混合现象,并且随着反应速率的增加,影响越明显.     

Non—equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method

In this paper,we propose a new approach to implementing boundary conditions in the lattice Boltzmann method (LBM).The basic idea is to decompose the distribution function at the boundary node into its equilibrium and non-equilibrium parts,and then to approximate the non-equilibrium part with a first-order extrapolation of the nonequilibrium part of the distribution at the neighbouring fluid node.Schemes for velocity and pressure boundary conditions are constructed based on this method.The resulting schemes are of second-order accuracy.Numerical tests show that the numerical solutions of the LBM together with the present boundary schemes are in excellent agreement with the analytical solutions.Second-order convergence is also verified from the results.It is also found that the numerical stability of the present schemes is much better than that of the original extrapolation schemes proposed by Chen et al.(1996 Phys.Fluids 8 2527).     

多孔介质内混溶流体间粘性指进现象的格子Boltzmann模型

提出一个模拟多孔介质内混溶流体间粘性指进现象的格子Boltzmann模型.采用双分布函数分别求解压力场和浓度场.在浓度场平衡态分布函数中引入与浓度扩散相关的参数.通过调节参数,使粒子碰撞过程中的松弛时间保持恒定.模拟了粘度相同的流体间的混相驱替问题.不同网格下的模拟结果均与解析解吻合良好,验证了模型的可行性.进一步研究粘度比和贝克莱数(Pe)对粘性指进现象的影响.结果显示,增大粘度比会促进"手指"的增长.当粘度比不变时,存在Pe的临界值.当Pe超过临界值时,"手指"前缘会出现分裂现象.对横向平均浓度场的研究显示,混合区域的长度随时间的变化分为两个阶段,它首先随着t1/2成线性增长,然后随着t成线性增长.     

格子Boltzmann方法中三种流固耦合格式的对比

详细对比分析模拟流固两相流的格子Boltzmann方法中三类不同的流固耦合格式(半反弹方法、亚网格方法和曲边界方法)的特点,发现半反弹方法所获得的速度、角速度存在较大的扰动,而亚网格颗粒方法则能够获得平滑的发展曲线.