LBM中的角点处理及部分代码（以D2Q9为例）

写在前面的话：其实有的时候我也很混乱角点用不用特殊处理，但是我在github上下的一份PPT里说：当两条边界上的值不同时，角点处理非常重要！所以当程序运行出问题却找不到原因的时候，可以把角点单独拿出来写一下，看是否是这方面的原因。

先放一张图，这是在Timm 那本书上截下来的。上次已经说过Timm将边界条件分成两种，link-wise（左）和wet-node（右）。在角点这主要介绍2种方法，反弹以及非平衡态反弹。

 

1、反弹

角点的反弹同直线边界一样，未知分布函数是通过迁移过来已知分布函数获得的。

步骤：

第一步：确定待求的角点分布函数。如左图凹角处待求的分布函数是f1 f2 f5，凸角处待求的分布函数是f5。