写在前面的话:其实有的时候我也很混乱角点用不用特殊处理,但是我在github上下的一份PPT里说:当两条边界上的值不同时,角点处理非常重要!所以当程序运行出问题却找不到原因的时候,可以把角点单独拿出来写一下,看是否是这方面的原因。
先放一张图,这是在Timm 那本书上截下来的。上次已经说过Timm将边界条件分成两种,link-wise(左)和wet-node(右)。在角点这主要介绍2种方法,反弹以及非平衡态反弹。
1、反弹
角点的反弹同直线边界一样,未知分布函数是通过迁移过来已知分布函数获得的。
步骤:
第一步:确定待求的角点分布函数。如左图凹角处待求的分布函数是f1 f2 f5,凸角处待求的分布函数是f5。
第二步:直接使用反弹法获得未知的分布函数
代码: