Palabos Tutorial_1_1解读

让我们打开tutorial_1文件夹（测试环境：Ubuntu 18.04）

打开Makefile文本文件

入门Palabos的话，Makefile里面有两个主要的注意点

1. palabos的路径，此处…/…/…表示返回上级三次，如果在其他文件夹运行程序，可以改成palabosRoot = /home/yulan/Palabos/palabos-v2.0r0。
2. 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时写出图片。