OpenFOAM中多相流求解器稳定性的问题

[tidedrinker]

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

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

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

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

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

[李东岳]

OpenFOAM一般发布出来的版本都是算法上有增强的。但是你这种4.1到8，这种版本胯的太大，说不准会有什么问题。跨度太大了。类似你这种

但是在4.1上可以收敛的算例放在8上就收敛不了

如果算例小的话，如果你确定二者的设置是完全相同的话。可以去org上面去提bug。否则的话，提bug也会被驳回来。我看你好像4.1的算例没有波，8有波浪，没准是这方面原因

[tidedrinker]

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

[火山口玩泥巴]

收藏一下，感觉有用