应该就是这个问题，虽然我还不确定问题是怎么发生的，但我已经解决了。
Linux Mint 18 软件源自带的OpenMPI是1.10.2，好像整个OpenMPI 1.10版本包括1.10.2和1.10.7都有这个问题，然后换成2.0.4就好了。

用paraview也可以做：
首先是读取数据，只读那个面上的数据，或者说在mesh parts里只选那个面。可能用提取面也能做，我没试过。
然后用slice，把你要的那一圈切出来。
接着用plotonsortedline在paraview里生成曲线图。
最后用file里的save data把数据保存，一般来说包括你要的压力或者其他数据，和对应点的xyz坐标。
我用这个方法做水轮机叶片的压力分布是没问题。
手头算例都太大，放一个简单的导流罩内壁的压力分布，道理是一样的。

@李东岳 这个forces是写在controlDict/functions里的，在求解过程中每个（或每若干个）时间步输出一次。这个压力和转矩是随时间变化的，不是根据结果文件用postProcessor -func那一套程序得到的。所以修改结果文件的文件名应该没什么用处？

一般情况下forces是输出基于p的压力，也就是说在二相流里考虑了浮力或者说水的重力的影响，现在我想得到不包含浮力的结果，又不想使用其他间接的方法，请问forces能否输出基于p_rgh的压力和转矩？

凹凸不平的表面要做出流线可以参考这里，是关于机翼表面的流线的。如果想要流线更漂亮一点，可能需要结合surfaceLIC使用。
但这是对于平动的表面，对于旋转的表面可能更复杂一些。不过如果只要方向的话，可以省略掉calculator那一步，这样的话就没有太大区别了。
而且实际上这里的流线反映的是剪切应力的方向和大小，虽然对于刚体的表面，剪切应力的方向和相对流速的方向应该是一致的。但问题在于，我不确定在表面边界条件一致的情况下，剪切应力的大小是否和相对速度成正比，特别是在曲面上。

@Joann 我先抛个砖，你的第一个问题应该是因为在fvSchemes里定义ddt的时候使用了localEular，而不是平常用的Eular。我理解的localEular就是在不同的区域根据网格尺寸使用不同的时间步长，所以不会有全局的库朗数输出，但会有最大最小时间步长输出。

@bestucan 我决定直接放图上去了，不用GCI了，主要是因为GCI到底应该限制在多少以下没有定论，同时GCI也需要求解结果满足单调性，偏偏我这个满足不了……

关于第二张图的画法，关键是使用calculator把叶片对应截面的速度加进去。加进去的速度与叶片对应截面速度大小相同，方向相反。然后就是SurfaceVector+StreamTracer。
比较麻烦的是，因为叶片吸力面后方出现了分离（如果没有分离也没必要画这个），所以需要多个Line Source才能把吸力面后缘附近的流线画出来。因此需要多个StreamTracer，就会造成叶片附近一部分流线过密，调起来比较烦，不过也能看。
下图是画图时候的pipeline browser：

StreamTracer1+StreamTracer2-5是一种画法，Line Source是相互平行的，如下图，下图中白色实线是StreamTracer1的Line Source：

StreamTracer1+StreamTracer6是另一种画法，Line Source是接近垂直的，如下图，白色实线是StreamTracer6的Line Source：

严格意义上Slice截面应该与叶片轴线垂直，但这里差的不多，也不准备拿来发论文，所以直接按垂直z轴算的，会有一些误差。

@cresendo 能详细说一下吗？第二张图那种流线我已经搞明白怎么画了，但第一张还是毫无头绪。

对于旋转的叶片，如下图所示的叶片表面流线是否可以使用paraFoam生成？

对于旋转的叶片，下图这种考虑旋转速度的2D流线是否可以使用paraFoam生成？

我可以做到对2D面生成流线，但如果不考虑旋转，生成出来的流线如下图，基本看不出多少有意义的信息。

除最后一张，以上图片均取自岳一轩. 叶片翼刀对水平轴风力机气动性能的影响[D]. 兰州: 兰州理工大学, 2020.

时间步长会显著影响计算的精度，但是不知道是我看论文不仔细，我发现很多做潮流能的数值模拟的文章，包括期刊和学位论文，都没有对时间步长对模拟的影响进行说明。
对于网格的独立性，有GCI这个参数可以进行描述，我想知道的是对于时间步长是否有类似的参数可以描述其独立性或者收敛性。
同时，GCI在哪个范围内可以认为达成网格无关？我看的论文中基本都只给了计算出的GCI，然后说满足独立性、无关性或者收敛性，但没说大致在那个范围才算满足要求？