第二个方案的源是在CFD中文网上看到的一篇帖子，其具体链接忘了。这个帖子给了一个github的链接，就是将OpenFOAM输出的颗粒位置信息转变为真实的位置信息的脚本。其链接在此。

1. 背景

我们知道，paraview之所以打不开OF输出的颗粒文件，是因为OF输出的颗粒信息的位置文件中的内容并不是真实的颗粒坐标。

而上面那个脚本就是将这个坐标转换为真实的位置坐标。

2. 具体步骤

GitHub链接里其实给的比较清楚了，说明文件写的十分良心。
主要有两种方法，一是边算边转换，而是算完之后统一转换。
说明文件中以MPPICFoam下的一个案例文件作为说明

cp -r $FOAM_TUTORIALS/lagrangian/MPPICFoam/column MPPICFoam_column
cd MPPICFoam_column

第一种方式

1、在system/controlDic文件的最后添加

functions
{writeCloudOldStyle1{type        writeCloudOldStyle;libs        ("liblagrangianExtraFunctionObjects.so");writeControl writeTime;clouds(kinematicCloud);}
}

2、运行算例
3、执行命令

find -name positions | while read line; do mv $line $line.coord; mv $line.orig $line; done
find -name positions | while read line; do sed -i -e 's=^\(.*object.*\)positions.orig;=\1positions;=' $line; done

之后就能用paraview成功打开文件进行查看了。

第二种方式

第二种方式是再计算完成之后进行。
1、需要在`system文件夹下创建一个positionConvertDict文件。

 /*--------------------------------*- C++ -*----------------------------------*\| =========                 |                                                 || \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           ||  \\    /   O peration     | Version:  5.x                                   ||   \\  /    A nd           | Web:      www.OpenFOAM.org                      ||    \\/     M anipulation  |                                                 |\*---------------------------------------------------------------------------*/FoamFile{version     2.0;format      ascii;class       dictionary;location    "system";object      positionConvertDict;}// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //functions{writeCloudOldStyle1{type        writeCloudOldStyle;libs        ("liblagrangianExtraFunctionObjects.so");writeControl writeTime;clouds(kinematicCloud);}}

2、执行命令

MPPICFoam -postProcess -dict system/positionConvertDict

这里的命令也可以换成DPMFoam开头，对于自己整的求解器可能会显示没有该命令的情况，但是可以把自己求解器的文件格式改成DPMFoam一样，在用DPMFoam命令转换，亲测有效。
3、执行命令

find -name positions | while read line; do mv $line $line.coord; mv $line.orig $line; done
find -name positions | while read line; do sed -i -e 's=^\(.*object.*\)positions.orig;=\1positions;=' $line; done

bug

这两个方法在并行计算时会存在一点小bug，因为他这个流程是我先生成一个真实颗粒位置的文件positions.org（再合并后的文件夹里），然后再把原始的颗粒位置文件（paraview无法识别的）的内容替换成positions.org。第三步的命令就是这个操作会检索所有可能的position名称的文件，但是在proc文件中并不会出现positions.org，因此会报错。
解决办法是把proc文件给删掉，然后再执行第三部命令