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

@yu_tian 我去翻了一下代码，CofG是指重力中心，不是很明白这个参数的意义，如果这个值不能代指旋转中心，也就是6里的origin的话，那就是旋转轴必须通过原点。需要你再确认一下。

@yu_tian CofG是指旋转中心吗？如果是的话，下边那个radialVelocity指的转速需要给定绕三个轴旋转的速度，（0 0 360）就是指绕x轴转速为0，绕y轴转速为0，绕z轴转速为360，归一化之后旋转方向就是（0 0 1），就是。如果你要改旋转轴，就要同时改旋转中心和转速，转速归一化之后就是旋转方向，加上旋转中心就是旋转轴。

@李东岳 不只是从4.1到8，5和6我也试过，感觉就是alpha的上界越界越来越严重，但也不是所有类似算例都跑不了。和波浪没关系，测试的时候是没有波浪的。现在我倒是找到办法了，就是继续降低pcorr的tolerance，从-7降到-10或者-11，然后nCorrector继续加，总是能收敛的。只要是收敛了，结果差的就很有限。只是计算资源消耗更厉害。

目前我主要是用4.1来模拟潮流能水轮机，考虑自由表面，总体上稳定性还行，不太需要专门考虑收敛性的问题。少数情况下有问题的话，调一下pcorr的收敛条件就行，还不行的话加一下nCorrectors，或者略微改一下网格，重新生成一下就行。

但是最近要考虑波浪，于是换成了8，但是在4.1上可以收敛的算例放在8上就收敛不了，而且不能通过过去在4.1上有效的手段改善。然后我又试了一下5和6，发现在一些情况下可以收敛，但是需要调整，比如使用更小的时间步。

我查了一下OpenFOAM.org上的更新说明，发现了几乎每个版本更新都有对MULES的修正或者改进。然后我就试了一下调一下MULES的设置，发现计算过程中确实有差别，但是最终还是会发散。

另外还有两个很奇怪的问题，或者也可能是一个问题的两种反映。一个是基本无论怎么调参数，发散的时间都差不多，基本上是转子转完一圈就开始出现明显不正常的交界面库朗数上升，或者alpha的下界向负数快速增大。另一个问题是，我把4.1中算了一定时间的数据用mapFields导进8的算例里作为初始值，在降低时间步长之后是可以稳定计算的，而且alpha下界越界不大，上界基本不会越界。

我很好奇的是，OpenFOAM中多相流求解器的稳定性是随版本更新逐渐降低的吗，或者至少是在存在逆压梯度和分离流的情况下是否是逐渐降低的？

@oitocfd 不是有小问题，而是我没有试，一般我要解决的问题用SST就行了，不用考虑其他的

@oitocfd 可以直接上SST算

@李东岳 应该不是，关闭所有CFD中文网网页后在重新打开，此时没有手动操作但显示在线（也就是上次勾选了保持登录状态），试图发帖，无论试几次，无论第几次，都会显示这个Forbidden。

发帖的时候右下角会出现Forbidden，然后发不了，可通过重新登陆解决，希望改进一下。

这个帖子其实包括两个问题，首先是OpenFOAM自带的的波浪边界是否支持纯波浪。目前我基本没有头绪，可以确定的是波可以造出来，但波形完全不对，Airy波，二阶和五阶Stokes都不行。

然后是波浪衰减，我已经基本搞明白了改哪些参数有用了。第一个是fvSchemes里的ddt，一般默认是Euler，波浪这里最好改成Crank-Nicolson，也就是算例waves里的设置。这样的话时间步长在1/250波浪周期就有不错的结果（看趋势或许更长也可以，但我没试），如果用Euler可能要1/1000以下才行。我上边贴了一篇论文，里边有详细的比较。Value的话0.9就行，再高提升有限而且不稳定。

第二个是我波浪参数有问题，要么是波高太大，要么是波陡太大，总之增大波长或者降低波高（交界处网格尺寸不变）都能有效改善衰减，可能是因为这个波浪参数非线性太强。同时，将Airy波（一阶波）改成Stokes波也会有明显改善，但二阶波和五阶波差别不大。

现在的话，水深和波长不变，波高降到2.4，改用Stokes二阶波和Crank-Nicolson，连波流耦合导致的波浪变形一起算上，波高衰减不超过3%，已经能用了。

而且，如果用outletPhaseMeanVelocity做出口边界，用不用消波好像没区别。

@oitocfd 非常感谢，这个问题我已经解决了波高衰减的那一部分，但是纯波浪那部分还没头绪。