interPhaseChangeFoam空化求解，负绝对压力问题

[yhdthu]

@李东岳 东岳前辈，目前有个问题还是无法解决，当发生cavitation collapse 的时候，即相变由低密度到高密度时，压力会发生一个很大的波动，这会影响绝对压力的计算，我现在还是没想通数值上到底是怎么变化的过程，会导致这种情况的发生，请前辈指教～

[队长别开枪]

@yhdthu 物理上空穴（气相）破裂的地方就会出现压力波，这也是气蚀现象的主要来源。

[yhdthu]

@队长别开枪 你好，这个物理现象我知道，但我不明白数值是怎么变化的出现了这种现象

[yhdthu]

附上目前问题压力云图：

t1时刻：

t2时刻：

t3时刻：

我又从新撸了一遍理论，目前的理解是这样：
压力泊松方程求解的是压力值，而不是求解的压力修正方程，但之后更新通量phi和体心值U时，用的是压力梯度，这也就是为什么压力虽然算不对，但并不影响U和alpha的计算结果。
由于压力泊松方程是椭圆形方程，跟给定的边界条件有很大影响，所以目前把出口边界改成了advective（nuTlida，U，p_rgh，alpha.water，除了nut）试试看算的如何。
除此之外，请问大家还有其他建议么？

[赵一铭]

这个问题我在CFD界发过一个类似的文章，也是一个类似背压升高的情况，不过不知道你这个是不是这个原因，结果怎么样？

[dyhe]

你好，你是在算naca0015吗？

[yhdthu]

@赵一铭 前辈好，那篇文章我拜读了，所以才想通的这些问题，目前正在计算，我不确定的是advective边界的使用方法是否正确，我是将除了nut外的所有边界均设为advective，请问这样对不对呢？

OUTLET
{
    type            advective;
    phi             phi;
}

[yhdthu]

@dyhe 我算的是naca66

[dyhe]

@yhdthu
你好，你那边空泡的生长，发展，与溃灭过程模拟的咋样，脱落频率和试验比的上吗?Schner模型会产生压力的强烈振荡吗？求指教:cheeky:

[yhdthu]

@dyhe 还没算完，周期表较长了

[yhdthu]

还有个问题，为啥这个湍流粘性最小值是负的呢？请问这个残差曲线该怎么看？

[dyhe]

@yhdthu
最小值为负是数值误差吧，我看后来收敛到0了
这个图不是残差图

[yhdthu]

@dyhe 我不太清楚他这个怎么定义的？

[dyhe]

@dyhe 
当空泡溃灭时，你的升力会出现这样的突变吗

[yhdthu]

@dyhe 可以把你的压力变化云图贴上来，我目前正在解决这个问题，没有看升阻力系数

[dyhe]

@yhdthu



针对此问题，你有没有什么想法，关于空泡流模拟，你有没有成功的算例呢？

[yuan_neu]

可以通过调整空化模型的参数来避免负压问题
具体可以参考引用的文献

[dyhe]

这个压力突变是不是空泡破裂产生的压力波造成的呢，你为什么一直在改边界条件，是不是openfoam自带的算例的边界条件有问题？

[yhdthu]

@dyhe 空泡破碎是物理解释，我要寻找数值原因

[zwl]

请问为什么会有两个质量源项和两个体积源项？