LBM中基于半程反弹的统一边界条件研究

格子Boltzmann方法能够有效地模拟复杂流场的流体流动,边界处理方法的选择对于模拟结果的可靠与否起着至关重要的作用.本文基于半程反弹原理对插值格式的曲线边界条件进行改进,提出一种统一的曲线边界处理新方法,并与格子Boltzmann模拟中常用的几种曲线边界条件做了对比研究.模拟结果显示该方法的计算结果与理论解精确吻合...     

大血管中血液流动的LBM模拟

真实的血液系统一般都具有复杂的几何边界.利用曲线边界条件构建二维血管中血液流动的格子Boltzmann模型,并对Poiseuille流和血液在肺部大血管中的流动进行了数值模拟,数值结果和解析解符合得很好,充分说明了该模型的有效性.所构建的模型可望进一步推广应用于具有复杂几何边界的一般流体系统运动的数值模拟.     

有限体积法与LBM分区耦合模拟方腔自然对流

在已有的密度分布函数重构算子的基础上,推导出了温度分布函数的重构算子,解决了格子Boltzmann方法(LBM)与有限体积法耦合计算传热问题的关键难题.选二维方腔自然对流对耦合方法进行了考核.在瑞利数Ra=103~106范围内,耦合结果同商业软件FLUENT结果符合得很好,并且各物理量在耦合界面处连续且光滑过渡.通过残差曲线可以看出,耦合模型在密网格以及大瑞利数情况下,数值稳定性要好于单一LBM.     

基于GPGPU的LBM障碍绕流的实时模拟

基于显卡的通用计算(GPGPU)是近年来并行计算和快速绘制的热点.格子Boltzmann方法(LBM)作为流体动力学的新方法,其并行性好,常常用于基于物理的流体模拟,且具有适用于复杂边界障碍的特性,但计算较为复杂.利用GPGPU技术来加速LBM的流体计算模型,构建了基于图形处理器(GPU)的流体计算框架,实现了格子Boltzmann计算的D2Q9和D3Q15模型,并用于实时的障碍绕流模拟.     

致密多孔介质中气体非达西现象的LBM模拟

针对气体在致密多孔介质中低速渗流时,其渗流规律在渗流曲线的低压段表现出对达西定律线性关系的偏离,存在着非达西现象。采用格子Boltzmann方法,研究气体和多孔介质的特性对气体渗流Klinkenberg效应的影响因素。结果表明:在气体渗流曲线的低压力梯度段,随着气体黏度系数、净围压、渗透率和孔隙率的变小,渗流曲线的非线性临界点向压力梯度增大的方向移动,对达西定律线性关系的偏离更明显。说明在低渗和低压情况下Klinkenberg效应不能被忽略,气体黏度系数和孔隙率对Klinkenberg效应作用有影响;当净围压或渗透率很大时,气体渗流流量和压力梯度符合达西定律线性关系。     

液滴撞击液膜喷溅过程的LBM模拟

气液两相流动现象广泛存在于自然界和工程应用中,而液滴撞击液膜后产生喷溅的过程则是这一类问题的典型代表.采用格子Boltzmann方法（LBM）对上述过程进行了数值模拟.采用基于“单相”（single-phase）模型的LB两相流方法计算得到了3组Reynolds数（Re）和Weber数（We）组合下液滴撞击液膜后所产生的3种效应.结果表明,随着We数和Re数的不断增大,液滴撞击液膜后将产生铺展、喷溅以及喷溅并伴有小液滴脱落等不同现象,计算所得的结果与实验及理论分析的结果相吻合,表明了LBM研究气液两相流动问题的可行性.     

格子Boltzmann方法模拟磁场作用下的融化传热过程

采用格子Boltzmann方法对磁场作用下方腔内相变材料的融化传热及液相流体流动过程进行数值模拟,重点分析磁场倾角对融化相变传热过程的影响,为实现相变传热的主动控制技术提供参考.研究表明磁场可以控制融化相变传热过程,当磁场倾角大于180°时,磁场加速度方向指向热壁面方向,磁场强化换热,热壁面平均Nu_h数增加,融化速率加快;当磁场倾角小于90°、磁场加速度方向指向冷壁面侧时,磁场弱化换热,热壁面平均Nu_h数下降,融化速率降低;磁场倾角为231°时,磁场强化传热效果最佳.     

一种微型泵的设计

设计了一种可灵活应用于微量分析系统的微型泵,详细分析了其结构和工作原理,说明了主要的制造工艺,并给出了它在精密化学分析中的应用实例,最后总结了它的优缺点。该微型泵采用硅微细加工技术制作,使用了电磁制动无阀式结构的设计方案,并使用双向驱动方式提高效率。该微型泵具有可靠性高、寿命长、体积小、成本低等优点,充分满足了分析系统小型化、利于批量生产和降低成本的要求。     

低Reynolds数下单排方柱绕流尾迹涡型的LBM模拟

针对工程中遇到的周期柱体绕流问题,采用格子Boltzmann方法(LBM)对低Reynolds数下单排周期排列方柱的定常绕流进行数值模拟,详细分析了柱体后面的尾迹涡型随着Reynolds数的变化情况,并通过对不同Reynolds数下尾迹涡型结构的分析,得到了二射流、三射流和四射流发生合并现象的临界Reynolds数,其中二射流出现合并的临界Reynolds数与文献的数值模拟结果一致,三射流和四射流出现合并的临界Reynolds数则偏大,但更接近实验结果。     

外环流旋转活塞泵转子加工形线的确定

针对在外环流旋转活塞泵转子的加工中传统经验方法存在的不足之处 ,提出了对转子加工形线的一种全新的确定方法顶点曲率圆法 .通过对转子理论形线的分析 ,给出了确定转子理论参数的方法 .同时对顶点曲率圆法作了详细的推导 ,得到了由转子理论参数确定转子加工参数的方法     

New boundary treatment methods for lattice Boltzmann method

0　ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｓｉｍｐｌｅｓｔａｎｄｃｕｒｒｅｎｔｌｙｗｉｄｅｌｙｕｓｅｄｌａｔｔｉｃｅＢｏｌｔｚｍａｎｎｅｑｕａｔｉｏｎｉｓｔｈｅｏｎｅｗｉｔｈｔｈｅＢＧＫ (Ｂｈａｔｎａｇａｒ Ｇｒｏｓｓ Ｋｒｏｏｋ)ｃｏｌｌｉｓｉｏｎｏｐｅｒａｔｏｒ[1]ｆα(ｒ+ｅαδｔ,ｔ+δ1) -ｆα(ｒ ,ｔ)=1τ[ｆ( 0 )α (ｒ ,ｔ) -ｆα(ｒ ,ｔ) ],α =0 ,1,2 ,… ,ｂ(1)ｗｈｅｒｅｆα(ｒ ,ｔ)ｉｓｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｆｕｎｃｔｉｏｎａｌｏｎｇｔｈｅα ｔｈｄｉｒｅｃｔｉｏｎ;ｆ( 0 )α :ｔｈｅｃｏｒｒｅｓ…     

Simulation of mesoscale interfacial properties using the lattice Boltzmann method

We derive the mesoscopic interparticle potentials from macroscopic thermodynamics for van der Waals, Redlich-Kwong, and Redlich-Kwong-Soave equations of state and find that all these potentials are very similar to the Lennard-Jones potential. To in-vestigate the interfacial property at the mesoscale level, we incorporate free energy functions into the single-component multiphase lattice Boltzmann model and obtain the saturated density coexistence curves and interface mass density profiles across the interface using this method with different equations of state. The simulation results accurately reproduce the properties of equilibrium thermodynamics. Numerical results for single-component phase transitions indicate that a bubble-growth process is obtained and the equilibrium phase diagram is achieved at a given temperature. Bulk free energy, the interfacial energy coefficient, and other properties of nonequilibrium thermodynamic parameters, which are used to examine interfacial properties, are obtained in these simulations, and all these parameters are found to obey irreversible thermodynamics.     

用lattice Boltzmann方法模拟矩形竖腔内混合对流

用lattice Boltzmann方法(LBM)模拟了纵横比(高度与宽度之比)大于1的矩形竖腔内同时存在自然对流和强制对流的混合流动时流场与温度场分布.比较了不同雷诺数(Re)、瑞利数(Ra)下矩形腔的流场分布在不同纵横比时的特点,研究了这些参数值对混合对流换热的影响,提出混合对流存在时Re和Ra对流动和换热的影响可由Richardson数(Ri)综合反映,并得到在不同纵横比时平均Nusselt数(Nu)随Ra和Re的变化关系式.与常用的LBGK模型不同,用了LB力矩模型.     

Simulation of phase transition process using lattice Boltzmann method

A new lattice Boltzmann model based on SC model, is proposed to describe the liquid-vapor phase transitions. The new model is validated through the simulation of the one-component phase transition process. Compared with the simulation results of van der Waals gas and the Maxwell equal-area construction, the results of the new model have a better agreement with the analytical solutions than those of SC and Zhang models. Since the obtained temperature range and the maximum density ratio in this model are expanded, and the magnitude of maximum spurious current is only between those of SC and Zhang models, it is believed that this new model has better stability than SC and Zhang models. Consequently, the application scope of this new model is expanded compared with the existing phase transition models. According to the principle of the corresponding states in engineering thermodynamics, the simulations of ammonia and water phase transition process were implemented using this new model with different equations of state. Compared with the experimental data of ammonia and water, the results show that the Peng-Robinson is the best equation of state to describe the phase transition process of ammonia and water. Especially, the simulation results of ammonia with Peng-Robinson equation of state have an excellent agreement with its experimental data. Therefore these simulation results have a significant influence on the real engineering applications.     

利用格子Boltzmann方法计算页岩渗透率

针对页岩中孔隙主要分布在纳米尺度的特点,建立Knudsen数修正固体边界并考虑镜面反弹的格子Boltzmann模型。利用二维平板模型,考察Knudsen数对渗透率的影响,得出在一定误差条件下,须考虑Knudsen数修正的最大孔隙宽度。由页岩的SEM扫描图像重构得到三维数字岩心,利用新模型进行模拟流动,并计算得到绝对渗透率。结果表明:利用Knudsen数修正后,通道中间部分流体速度增大,固体边界处速度变小;压力不变,随着孔隙直径变小,Knudsen数增大,渗透率变小。     

基于多松弛格子Boltzmann方法插值反弹运动边界研究

采用多松弛格子Boltzmann方法,结合曲边插值边界处理方法,在固定笛卡尔坐标系均匀网格下,模拟了运动圆柱绕流问题.结合非平衡外推思想,提出了一种新的预估格点分布函数的方法:通过"圆柱表面垂直方向附近的流场内格点的密度和圆柱的速度"构造平衡分布函数,并将其作为新生成格点的分布函数.利用对绝对误差的频谱分析技术,与现有两种方法(基于平均和预估思想)对比结果显示:新提出的方法在具有更高的精度的同时,可以保证仅需要流体域内较少的流动信息,从而为流体内的颗粒碰撞等问题提供了有利条件.通过误差频谱结果,进一步对相应方法中的误差根源进行了讨论,为进一步改进方法提供了依据.对相同雷诺数下加大一倍流场区域的计算结果表明,加大网格可以减小升阻力误差在低频时频谱幅值.     

用LBM方法模拟壁面驱动粘性不可压半圆形空腔流

基于格子Boltzmann方法(LBM)数值模拟壁面驱动的粘性不可压半圆形空腔流. 采用具有二阶精度的曲线边界处理方法, 得到了不同雷诺数下的流线图、 涡线图及速度分量沿半圆形中心线的分布. 在小雷诺数的条件下, 流动状态仅由一个涡组成; 随着雷诺数的增加, 出现一个二级涡, 涡的大小与雷诺数有关. 数值结果表明, 格子Boltzmann方法简单有效, 适合处理该问题.     

Investigation on the mechanism of convective heat and mass transfer with double diffusive effect inside a complex porous medium using lattice Boltzmann method

A new lattice Boltzmann model is proposed for the coupling of multi-physics in natural convection by introducing another distri-bution function to represent the scalar transport of mass and an additional source term on the right hand side of the lattice Boltzmann equation (LBE) based on the TD2G9 model. According to the Boussinesq assumption and considering the coupled diffusive effect, the governing equations for the coupling of multi-physics in natural convection is proposed based on the non-equilibrium thermodynamic theory. Combined with the algorithm for the reconstruction of a porous medium, this model is used to investigate the coupled heat and mass transfer under multi-physics of natural convection inside a porous medium at the pore scale. The characteristic of profiles of dimensionless velocity, temperature and concentration are shown graphically for different values of Rayleigh number (Ra) and porosity. Furthermore, the influence of the temperature gradient on the mass transfer inside a porous medium has been examined numerically. Thus, the mechanism of the coupled heat and mass transfer inside a porous medium is investigated numerically at the pore scale innovatively.     

基于格子玻尔兹曼方法的致密砂岩驱替模拟

基于数字岩心建模方法和格子玻尔兹曼方法,开展致密砂岩水驱模拟研究,分析致密砂岩水驱微观过程和机制。首先利用大庆油田龙虎泡地区高台子组致密砂岩CT扫描数据,建立三维数字岩心模型;然后利用Shan-Chen格子玻尔兹曼方法,建立油水两相驱替数值模型;之后对比分析不同注入速度和注入压力条件下驱替效率、换油率与驱替参数的关系,并给出致密砂岩三维微观孔隙尺度的驱替过程。结果表明:相同累积注入量条件下,提高水驱注入速度和驱替压差可以提高最终驱替效率;水驱注入速度和驱替压差存在驱替效率拐点;恒速水驱的平均最终驱替效率为49.6%,恒压水驱的平均最终驱替效率为47.3%。     

Application of the lattice Boltzmann method to electrohydrodynamics： Deformation and instability of liquid drops in electrostatic fields

It is attractive and meaningful to apply the lattice Boltzmann method (LBM) to electrohydrodynamic (EHD) problems. As an important and fundamental problem in EHD, the behavior of freely-suspended liquid drops subjected to electrostatic fields was simulated with a 2-D LBM by means of the permittiv- ity-density relation proposed by the authors and Peng’s multiphase LBM which can obtain high liquid-to-vapor density ratio and low liquid-phase compressibility. The phenomena and laws of deformation, especially the instable behavior of the liquid drops were well simulated. And two important critical parameters, i.e. the critical permittivity ratio and the critical electric-field strength were obtained. These results agree with the corresponding theories and experiments. The simulations show the particular advantages and potential of the LBM in studying multiphase electrohydrodynamic problems.