3. 边界条件处理 边界条件的处理在LBM中不难但是十分重要,大多时候报错都是因为边界处理的不正确或者不合适导致的。这里程序中在左右使用的是周期边界条件,上下壁面所使用的是Zou-He[2]无滑移边界条件,Zou-He边界条件是计算宏观量前进行处理的。计算公式如下: 下壁面 f_2 = f_4 f_5=f_7+\frac{\rho \le...
由于LBM MATLAB参考较少,为了方便大家学习,少走我走过的弯路,陆陆续续将一些自己写过的简单算法分享一下。这篇知乎分享的代码略有阉割,主要是有一些数据输出的部分被删除,各位看官斟酌参考。 本程序可模拟不同亲疏水性壁面接触角,也可进行Laplace方程验证,主要参考默罕默德和KIM的两本书,由于版权问题,就不放链接了...
首先，设定程序中的基本参数。这些参数包括壁面性质（obst=1表示固体，0表示流体），以及与当前格子相邻的9个速度分量方向格子状态（obstk值描述边界类型，0为流体，1为固体边界，2表示左右周期性边界）。接着，初始化液滴状态，此步骤为程序开始前的必要准备。计算初始压力，此操作有助于模拟环境压力对...