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