【LBM】开源软件(7)dat文件导入

本次参考:Yulan Fang。这纯属我个人的自学笔记,因为palabos的参考资料数时有限。Yulan Fang大佬真的很强,参考他(她)的资料,一步一步来。没有任何商业用途,如果介意,联系我,立删!!

这是一个非常有意思的操作,外部流场障碍物流动模型。之前我在手动编程时偶尔也会做一些类似的模型,这次用palabos试试。

思路和手动编程的思路是一样的,就是二值化,并将1赋值给壁面。


png图像

利用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);//设置为反弹边界。

同时还得依照画布的尺寸将长度设置为等价长度。

结果:


初始


1


2


3


4

同时,参考【LBM】开源软件(6),也可以转化成VTK格式,用paraview进行后处理。


paraview(1)


paraview(2)
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。