三相空化计算明渠流动失败

[火山口玩泥巴]

在interPhaseChangeFoam的基础上改了改弄了个三相空化求解器来算明渠流动，第一步就会出现如下所示错误
Debug发现错误出现在湍流模型参数的更新这里，用的RANS模型：

具体报错是FieldFunctions.C里的sumMag函数这里：

请问各位大佬，是啥原因造成的呢？

[火山口玩泥巴]

顶一下

[火山口玩泥巴]

大佬们给点建议

[火山口玩泥巴]

我现在尝试换了一套几何模型，重新画了更细的网格，网格质量除了几个面会有一点非正交性较大外其他没有问题，仍然会出现一样的问题，即在进行湍流模型参数更新的过程中浮点溢出。
大佬们有啥思路解决嘛？是边界条件的问题吗？

[李东岳]

kOmegaSST么，kOmega发散么

[火山口玩泥巴]

@李东岳 都试过了老师，都是会差不多的时间发散，即出现bounding k和bounding omega；
刚刚我用粗网格的时候又出了一些别的问题：当空化模型里面的n0设的较大就会出现解压力方程的时候浮点溢出，具体是出现在GAMGSolver的Scale函数里面的缩放因子的分母的计算那里，我折腾半天没弄明白那里的溢出原因（我看了溢出单元的压力方程源项以及空化模型提供的主对角线项，都不怎么会导致溢出），因为我暂时还没弄懂那个函数的作用，就先给压力方程换成PCG求解了，这个问题暂时就没出现了；

[火山口玩泥巴]

同时我发现相分数存在越界情况，尤其是在细的那一套网格中，越界情况尤其明显，蒸气相甚至最小值达到了-1.4 ，相方程的求解使用的是显式MULES（图方便），这是因为显式MULES效果不佳嘛？？我想在每次求解过后强行限定α的值：

alpha1 = min(max(alpha1, scalar(0)), scalar(1));

这样是否可行？

[火山口玩泥巴]

直接先用upwind处理相方程了。然后尝试使用LES的话可以正常计算粗网格，前提是得把空化核密度降低两个数量级，否则出口处会出现奇怪的低压区，这里不懂为什么会这样。

[李东岳]

你自己写的求解器么。标准求解器发散么

[火山口玩泥巴]

@李东岳 不是完全自己写的老师，是在interPhaseChangeFoam基础上改的，加了一个空气相进去。
在粗网格下面，用interFoam来算，使用komegaSST，在一定的k和omega的边界条件下不会出现发散的情况；但用自己的求解器去算即使把相变系数调到0，k和omega的边条初条怎么调都还是会发散。

[李东岳]

自己写的求解器不好debug，标准算例+标准求解器可以debug

[火山口玩泥巴]

@李东岳 是的老师，OF没提供三相空化的标准求解器，只能自己改着用，但又没啥经验，总怕自己改的这玩意稳定性太差。关键我对湍流模型了解的太浅了，k和omega越界这个东西我想改都无从下手。

[李东岳]

你用层流发散么

[火山口玩泥巴]

@李东岳 LES和laminar都试过，只要把空化核数降一些就不会发散。