OpenFOAM中多相流求解器稳定性的问题
-
目前我主要是用4.1来模拟潮流能水轮机，考虑自由表面，总体上稳定性还行，不太需要专门考虑收敛性的问题。少数情况下有问题的话，调一下pcorr的收敛条件就行，还不行的话加一下nCorrectors，或者略微改一下网格，重新生成一下就行。
但是最近要考虑波浪，于是换成了8，但是在4.1上可以收敛的算例放在8上就收敛不了，而且不能通过过去在4.1上有效的手段改善。然后我又试了一下5和6，发现在一些情况下可以收敛，但是需要调整，比如使用更小的时间步。
我查了一下OpenFOAM.org上的更新说明，发现了几乎每个版本更新都有对MULES的修正或者改进。然后我就试了一下调一下MULES的设置，发现计算过程中确实有差别，但是最终还是会发散。
另外还有两个很奇怪的问题，或者也可能是一个问题的两种反映。一个是基本无论怎么调参数，发散的时间都差不多，基本上是转子转完一圈就开始出现明显不正常的交界面库朗数上升，或者alpha的下界向负数快速增大。另一个问题是，我把4.1中算了一定时间的数据用mapFields导进8的算例里作为初始值，在降低时间步长之后是可以稳定计算的，而且alpha下界越界不大，上界基本不会越界。
我很好奇的是，OpenFOAM中多相流求解器的稳定性是随版本更新逐渐降低的吗，或者至少是在存在逆压梯度和分离流的情况下是否是逐渐降低的？
-
