多孔介质中含溶解反应的互溶驱替过程格子Boltzmann研究

采用格子Boltzmann方法研究了孔隙尺度下多孔介质内含流固溶解反应的互溶驱替过程,重点研究了被驱替流体与驱替流体黏性差异较大的情况下,溶解反应引起的多孔介质内部结构变化对驱替过程的影响;定量分析了不同达姆科勒数及佩克莱数下多孔介质孔隙率和驱替过程驱替效率随时间的演变.研究结果表明:达姆科勒数较大时,溶解反应的发生会在多孔介质内部生成虫洞,导致一部分被驱替流体不能被波及,驱替流体沿虫洞离开多孔介质,造成驱替效率的减少.在此基础上,随着达姆科勒数的增大,孔隙率变化越大,生成的虫洞越宽,最终驱替效率变大,但仍小于无溶解反应时的驱替效率;随着佩克莱数的增大,指进增长速度越快,孔隙率变化越小,驱替效率越小.     

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

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

多相非牛顿流体驱替过程的格子Boltzmann模拟

将非牛顿流体模型引入多相格子Boltzmann方法（LBM）,研究多相非牛顿流体的驱替过程.首先以牛顿流体驱替过程作对照,比较非牛顿流体的驱替特性.然后基于简单结构模型,分析不同界面张力下,非牛顿流体的驱替能力和驱替过程中被驱替相的形态变化规律,揭示界面张力的影响机理.结果表明：界面张力在驱替过程中起阻碍作用,其它参数相同情况下,界面张力越低,流体的驱替能力越强.最后,结合四参数随机生长（QSGS）结构生成方法,模拟不同压差、不同界面张力以及不同润湿性条件下多孔复杂结构内的非牛顿多相驱替过程,分析这些参数对过程的影响.     

复杂微通道内非混相驱替过程的格子Boltzmann方法

本文采用改进的基于伪势模型的格子Boltzmann方法研究复杂微通道内的非混相驱替问题.这种方法克服了原始伪势模型中计算结果对网格步长的依赖.首先用Laplace定律验证模型的正确性,然后用该方法研究壁面润湿性、粗糙结构、黏性比以及距离对非混相驱替过程的影响.模拟结果表明:与壁面粗糙结构和黏性比相比,壁面润湿性的影响是决定性的因素.随着接触角的增加,驱替效率增加,当接触角大于某一值后,驱替效率不再变化;随着黏性比的增加,驱替效率增加;而壁面粗糙性对驱替过程的影响较复杂,只有凸起半圆的半径在一定范围内增加时,驱替效率增加;距离较小时将促进驱替过程.     

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

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

微多孔介质非牛顿流体格子Boltzmann模拟

根据修正Blake-Kozeny方程对Herschel-Bulkley模型的非牛顿流体在微多孔介质中流动的定量描述,通过增加外力项建立了非牛顿流体在宏观多孔介质中的格子Boltzmann模型,并将该模型应用到平板多孔介质中非牛顿流体的电渗驱动流模拟。对多孔介质渗透性参数和非牛顿流体流变性参数的影响进行模拟和分析,获得了各种参数对多孔介质中非牛顿流体流动速度的影响规律。     

致密多孔介质中气体渗流的格子Boltzmann模拟

为了能用格子Boltzmann方法来正确地刻画致密多孔介质中微尺度流动问题,对单通道模型进行推广,并将其应用于孔隙群里气体渗流的数值模拟.通过对部分具有代表性的多孔介质中的实际流动问题进行模拟,研究渗透率与平均压力和Knudsen数之间的相互关系.基于理论分析及相关文献中的试验结果,验证模拟结果,为用格子Boltzmann方法深入研究气体渗流问题奠定基础.     

多孔介质内溶解与沉淀过程的格子Boltzmann方法模拟

本文采用格子Boltzmann方法模拟了多孔介质内的溶解和沉淀现象, 并分析了雷诺数、施密特数、达姆科勒数对多孔介质孔隙结构及浓度分布的影响. 结果表明: 对于多孔介质内的溶解(沉淀)过程, 当雷诺数越大时, 孔隙率越大(小), 平均浓度值越小(大); 当达姆科勒数或施密特数较小时, 溶解和沉淀过程均受反应控制, 此时反应在多孔介质的固体表面较为均匀的发生; 当达姆科勒数或施密特数较大时, 溶解和沉淀过程均受扩散控制, 此时反应主要发生在上游及大孔隙区域.     

多孔介质内交变流动与换热的格子Boltzmann研究

针对制约脉冲管制冷机效率提高的蓄冷器内交变流动与换热问题,本文使用格子Boltzmann方法计算并分析了多孔介质内交变流动与换热特性。结果表明,交变流动压力波振幅对速度与压力和温度间的相位差的影响很小;通道界面处速度与压力的相位差在某一孔隙率下有一极小值;多孔介质内最大阻力系数与雷诺数的关系可以为设计蓄冷器时填料结构的...     

倾斜多孔介质方腔内纳米流体自然对流的格子Boltzmann方法模拟

利用格子玻尔兹曼方法 (lattice Boltzmann method, LBM) 对倾斜多孔介质方腔内 Al₂O₃-H₂O 纳米流体的自然对流进行数值模拟, 考虑了孔隙率 (0.3 ≤?≤ 0.9)、瑞利数 (10³ ≤ Ra ≤ 10⁶)、纳米颗粒体积分数(0 ≤ ? ≤ 0.04) 和倾斜角 (0°≤γ≤ 120°) 等因素的影响, 研究了正弦温度分布边界条件下倾斜多孔介质方腔内纳米流体的自然对流传热机理. 结果表明: 若?和γ保持不变时, 随着 Ra 数的增大, 热壁面处的平均努塞尔数 (Nu_ave 数) 呈现出先减小后增大的趋势; 对于给定的 Ra 数, 当γ = 0°时, 随着孔隙率的增大, 热壁面处Nu_ave 数逐渐增大, 当γ = 40°, 80°和 120°时, Nuave 数在?= 0.7 左右时达到最大值; 若?和 Ra 数保持不变, 当γ = 40°时, 方腔内的自然对流换热效率最强, 当γ = 80°时热壁...     

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

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

Lattice Boltzmann方法模拟多孔介质Klinkenberg效应

格子Boltzmann数值模拟方法是研究复杂的多孔介质结构特别是Klinkenberg效应的有效方法之一，对处理复杂边值问题尤其有效，用格子Boltzmann方法研究了气流穿越多孔介质问题，并将数值计算结果与实验结果进行了比较，结果表明格子Boltzmann方法是数值模拟气流穿越多孔介质问题的有效方法之一。     

Natural Convection of Temperature-Sensitive Magnetic Fluids in Porous Media

In this article, natural convection of a temperature-sensitive magnetic fluid in a porous media is studied numerically by using lattice Boltzmann method. Results show that the heat transfer decreases when the ball numbers increase. When the magnetic field is increased, the heat transfer is enhanced; however, the average wall Nusselt number increases at small ball numbers but decreases at large ball numbers due to the induced flow being more likely confined near the bottom walls with a high number of obstacles.     

Lattice Boltzmann simulation of transverse wave travelling in Maxwell viscoelastic fluid

A nine-velocity lattice Boltzmann method for Maxwell viscoelastic fluid is proposed. Travelling of transverse wavein Maxwell viscoelastic fluid is simulated. The instantaneous oscillating velocity, transverse shear speed and decay rateagree with theoretical results very well.     

Lattice Boltzmann method for three-dimensional moving particles in a Newtonian fluid

A lattice Boltzmann method is developed to simulate three-dimensional solid particle motions in fluids. In the present model, a uniform grid is used and the exact spatial location of the physical boundary of the suspended particles is determined using an interpolation scheme. The numerical accuracy and efficiency of the proposed lattice Boltzmann method is demonstrated by simulating the sedimentation of a single sphere in a square cylinder. Highly accurate simulation results can be achieved with few meshes, compared with the previous lattice Boltzmann methods. The present method is expected to find applications on the flow systems with moving boundaries, such as the blood flow in distensible vessels, the particle－flow interaction and the solidification of alloys.     

Lattice Boltzmann Solid Particles in a Lattice Boltzmann Fluid

We define a lattice Boltzmann model of solid, deformable suspensions immersed in a fluid itself described in terms of the lattice Boltzmann method. We discuss the rules governing the internal dynamics of the solid object as well as the rules specifying the interaction between solid and fluid particle. We perform a numerical drag experiment to validate the model. Finally we consider the case of a population of flexible chains in suspension in a shear stress flow and study the influence on the velocity profile.