Palabos一般输出为vti格式文件,我这段时间为了处理数据的格式转化,研究了一下这些数据的构成。
基本上通过paraview可以很方便的保存整个时间段的模拟输出,格式为csv,然后通过matlab遍历读取csv,存入矩阵,再reshape即可得到流场的矩形数据。
如果是二维的,直接contourf则可以查看。如果是3D的数据,需要通过velocity(:,:,25)这样的手段来访问z轴第25行的截面。
值得注意的是,代码是按照文件命名读取csv,所以csv文件有几百个的时候,确保前100个的命名末尾是001,002,003。还有,点paraview工具栏的slice,保存数据时会有xyz坐标点。
代码已上传至我的github,找到matlab_csv_reader即可。
https://github.com/Yulan-Fang/SampleCodesofPalabos