针对岩石物理试验中出现的孔隙流体(油水)两相分离现象，应用格子Boltzmann(LB)方法中的两相不相溶流体的伪势模型，对油水界面动力学行为进行微观数值模拟，分析多孔介质中两相流动的微观特征，并从理论上给出两相不相溶流体界面张力因子Gf值的确定方法。模拟由于表面张力造成的油水两相分离现象，在此基础上研究润湿性对真实储层岩心孔隙流体两相分离的影响，并实现全程动态可视化。研究表明，用LB方法进行储层岩石油水两相分离简便易行、形象直观，是研究流体分离规律和特点的重要评价方法。

首次用格子Boltzmann方法中的伪势模型对液滴撞击固壁的动力学行为进行了数值模拟。详细研究了液滴在壁面上的流动状态以及各种因素对撞击过程的影响。通过数值模拟得到:壁面的可润湿性越小,液滴越容易发生反弹,液滴的回缩速度越快;液滴的撞击速度越大,所得到的相对直径越大,回缩速度越快;液滴的粘性越小,所得到的相对直径越大;液滴的表面张力越大,液滴越容易发生反弹现象。另外,液滴的最大相对直径与We数满足一定的线性关系,这些结果与前人的理论预测和实验结果完全吻合。

格子Boltzmann方法的理论及应用

LBM方法，格子Boltzmann方法，计算流体力学，何雅玲编著

对格子Boltzmann方法求解含第三类边界条件的扩散方程进行了理论和数值研究,构造了一种新的基于bounce-back的边界处理数值格式,用来处理复杂边界问题。借助渐近分析,证明了新方法的数值相容性。用数值算例从不同角度分析了算法的精度和稳定性等,与已有算法相比,新方法在精度、稳定性和效率方面均有较大提高。最后通过一个复杂边界反应扩散的示例演示了新方法应用于复杂多孔介质内多物理化学输运模拟的可行性和有效性。

针对二维对流扩散方程，基于D2Q4格子速度，用Chapman-Enskog多尺度分析技术，将时间尺度取为二阶，空间尺度取为一阶，推导了各个速度方向上的平衡态分布函数所满足的条件，给出了简单且对称的平衡态分布函数表达式，所得到的平衡态分布函数能正确地恢复出二维对流扩散方程，从而构建了一种新的求解二维对流扩散方程的D2Q4格子Boltzmann（LB）模型。用所给LB模型对扩散方程和两个不同初边界条件的对流扩散方程进行了数值求解，数值实验结果表明数值解与精确解吻合较好，与相关文献结果比较边界误差要小得多，验证了模型的有效性。

基于格子Boltzmann 方法的页岩气微观流动模拟

Boltzmann方法的理论及应用-何雅玲&格子Boltzmann方法的原理及应用-郭照立 二合一版本

matlab仿真曲线生成代码格子-玻尔兹曼方法-GPU 该程序是GPU CUDA版本D3Q19 BGK格子Boltzmann方法计算流体动力学求解器，用于模拟稳态/非稳态3D单相牛顿流，其中不考虑运动边界和体力（即重力）。 与CPU串行代码相比，此GPU代码在NVIDIA Geforce 2080ti中的速度快250倍以上，在NVIDIA Geforce 1050ti中的速度快140倍以上，并且具有相同的精度。 要运行它，必须具有CUDA Toolkit的NVIDIA GPU。 它包含：A. LBM讲义； B。 作者的博士学位论文（第4章详细介绍了LBM的实现）； C． 关于统一笛卡尔网格生成的论文（CartGen：鲁棒，高效且易于实现的统一/八叉树/嵌入式边界笛卡尔网格生成器）； D． 三个用于表面重建/平滑的Matlab工具（MyCrustOpen，fitNormal和smoothpatch）； E. Matlab函数geo_preprocess生成统一的笛卡尔非人体拟合网格； F.四个模拟案例： Lid_driven_cavity：稳定的层流 Poiseulle_flow：稳定的

为研究平板间方柱绕流上下平板对置于其中的方柱绕流所产生的影响，采用格子Boltzmann方法对二维平板间低雷诺数(Re=100)方柱绕流问题进行了数值模拟研究．分析了3种不同阻塞比下，平板边壁对方柱的升、阻力系数、Strouhal数和尾涡流场的影响．结果表明：平板对方柱绕流特性有明显的影响，随着阻塞比的增加，阻力系数和Strouhal数均增大，与无边壁相比阻力系数可增加达30％，而升力系数却随之减小．计算结果与相关实验数据相吻合，验证了格子Bohzmann方法对钝体绕流非定常问题模拟的有效性．