基于LBM-IBM-DEM的圆形颗粒在粘性流体中沉降的耦合模型及数值模拟[D]. 湘潭大学, 2013. 的解读.

感谢李浩同学的文章。在这里记录几个问题。

一、这个公式中，第二个fβ是上一步中的分布函数，而fα与fβ的区别是，以BounceBack边界为例：

%   7 3 6  ^y   
%    \|/   |   x   
%   4-1-2   --->   
%    /|\   
%   8 5 9
cx= [  0,   1,  0, -1,  0,    1,  -1,  -1,   1];     
cy= [  0,   0,  1,  0, -1,    1,   1,  -1,  -1];

opp = [ 1,   4,  5,  2,  3,    8,   9,   6,   7]; %分布函数在中心点是从1开始编号。

然后for i=1,9

fα(i)=fβ(opp(i))。

二、在做沉降的例子时候，似乎需要让小球的沉降时间等于流体时间的百分之一。