让我们打开tutorial_1文件夹(测试环境:Ubuntu 18.04)
打开Makefile文本文件
入门Palabos的话,Makefile里面有两个主要的注意点
- palabos的路径,此处…/…/…表示返回上级三次,如果在其他文件夹运行程序,可以改成palabosRoot = /home/yulan/Palabos/palabos-v2.0r0。
- projectfiles = ****.cpp则是需要运行的程序。
调出终端
此处tutorial_1文件夹具备程序运行的所有条件,在终端先输入make编译一下,结束后可以输入./tutorial_1_1可以单核运行。如果你已经安装了并行计算软件,比如openmpi,则可以输入mpirun -np 4 ./tutorial_1_1.
结果
运行结果如图,讲述的是某个大区域,设置一个方形小区域的密度稍微大于其他处,流体的流动情况。
在保存图片的文件夹中打开终端,输入convert -delay 5 -loop 0 u*.gif output.gif可以获得动态图片。
代码解读
第37行,格子表示器,表明了使用D2Q9模型。
上方框定义了的正方形中心点的坐标,即图上所说的半径(nx/6),实际上表示的是正方形中心到各个边的距离。
下方框Box2D表示形成区域,四个数是x轴坐标起始点和y轴坐标起始点。
此处也很明显地表示了方形区域的密度是稍微大一点的,方框内密度是 rho0+deltaRho,其他区域密度是 rho0。
从上往下,依次为:
输出图片保存位置,此处为tmp文件夹;
程序迭代的次数,此处表示1000次就结束;
nx和ny,整个计算流域的尺寸,此处为600×600;
边界循环条件,true则表示流体触碰到一个边界则会从另一侧边界回来,改成false,这样就没有边界了。如下图,同样的流体模拟,下图是不是和最上面的结果很不同呢?
Line 92:iT%5==0 这里的5表示每迭代5次保存一次图片,数字越小生成的图像越精密。(而总迭代次数在上文提到过,可以自行设置)
这个算法是不是很巧妙? _iT_的数值每次涨1,以取步数余数为0时写出图片。