本次参考:Yulan Fang。这纯属我个人的自学笔记,因为palabos的参考资料数时有限。Yulan Fang大佬真的很强,参考他(她)的资料,一步一步来。没有任何商业用途,如果介意,联系我,立删!!
这是一个非常有意思的操作,外部流场障碍物流动模型。之前我在手动编程时偶尔也会做一些类似的模型,这次用palabos试试。
思路和手动编程的思路是一样的,就是二值化,并将1赋值给壁面。
利用MATLAB将图片转化成二进制格式的dat。
之后再将dat文件放入到palabos模型文件中。
在palabos的C程序中加上:
MultiScalarField2D<bool> boolMask(parameters.getNx(), parameters.getNy());
plb_ifstream ifile("1.dat");
ifile >> boolMask;
defineDynamics(lattice, boolMask, new BounceBack<T,DESCRIPTOR>, true);//设置为反弹边界。
同时还得依照画布的尺寸将长度设置为等价长度。
结果:
同时,参考【LBM】开源软件(6),也可以转化成VTK格式,用paraview进行后处理。