OpenFOAM动网格案例中相分数突然增大并发散

[Stan同学]

在计算多相流动网格案例时，计算一段时间后库朗数突然增大并导致计算发散。起初以为是网格变形过大导致的发散，但查看日志文件和结果发现，计算并非是结构物受到波浪冲击后达到最大变形时终止的，而是在逐渐回到原位置时终止的，此时相分数alpha.water超过了1并在急剧增大。这时结构物附近有一定的越流（如下图所示，截图截得不太整齐请大家见谅），且计算停止时（44.2s）的气液交界面显然是有很大问题的。因此我认为应该是相方程计算不收敛进而导致的整个计算终止。

日志文件中相分数在44.076s开始大于1：
最后计算崩溃：

因此，想请问大家是否在多相流、动网格的计算中是否遇到过类似的问题，有什么办法可以解决这种相分数越界的问题吗？

[lyj_]

您好，请问您的这个问题解决了吗？我目前也遇到同样的问题

[李东岳]

这种就是发散了。没有一步到位的debug方法

[Amadeus]

你的这个算例是并行计算吗? 如果是并行，那么建议把网格划分方式先改成

numberOfSubdomains 6;

method          simple;

simpleCoeffs
{
    n               (3 2 1);
    delta           0.001;
}

这种比较少的核心数跑一次.

如果串行也会出现这种情况就确实很麻烦....