该博客介绍了格子玻尔兹曼方法（LBM）的基本算法，用于计算流体运动，避免了组装线性方程组的复杂性。博主详细阐述了LBM的时间步迭代过程、边界条件处理，包括固壁边界和速度边界，并给出了顶盖驱动流的计算结果。文章最后提供了相关代码供读者参考和运行。
上一篇文章计算了低雷诺数时的方舱蠕动流，但简化过度，只计算了稳定解，且没有考虑粘性——可以看到速度、压强分布完全对称。
吴桐与：有限差分法——蠕动流的数值模拟zhuanlan.zhihu.com偶然得知格子玻尔兹曼方法（Lattice Boltzmann Method, LBM）无需求解线性方程组即可计算流体运动，让我顿时兴奋了起来。写过计算力学程序的伙伴们都知道，组装系数阵是一件多么令人头痛的事情——
LBM的简单算法
LBM将离散了粒子的速度空间和真实空间，并同时离散分布函数，用一个碰撞算子模拟一个时间步内的分布函数的演变，继而计算包括密度、压强、内能等诸多流场信息。我也是这两天拜读的LBM经典著作，才略入此门懂点皮毛。详细的理论推导请参考下面的教材，本文只介绍我所参考的算法。
The Lattice Boltzmann Method Principles and Practicedoi.orgLBM每前进一个时间步完成一次循环，因此是一种迭代算法，一个循环内的主要步骤为：
- 由分布函数
计算密度
和速度
,
- 计算平衡态分布函数
- 由碰撞算子更新
最低0.47元/天 解锁文章
扫一扫
4万+
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
