关于p_rgh,p,rho*gh的讨论

[yhdthu]

众所周知，用p_rgh计算会有较好的稳定性。如果计算结束，我们想要关注压力云图，必定是看p。
因为p_rgh = p - rho*gh
那么问题来了，rho*gh会影响p_rgh的值，所以它也叫虚假压力。
由于N-S方程求解的是静压，p_rgh和rho*gh将p分成了两个部分，两者在不同参考系下一般是不同的，它们是怎么相互作用组合成了我们要的压力p呢？

[李东岳]

rhogh会影响p_rgh的值

为什么？

[yhdthu]

东岳前辈好，我很久之前在cfd-online上看到一位外国朋友分享的p_rgh和rho*gh的理解和云图计算结果（不好意思找了很久找不见了），显示这两者可以是不连续的，但是二者加和后得的p却是连续的，这东西跟选择的参考系有关系，比如我把网格放在“水平面以下”，这时候rho*gh就有物理意义。

另外我在看代码的时候，也没搞清楚怎么解着解着就得出了连续的p？还请前辈指教

PS：我做了个二维水翼的空化模拟，但是我认为压力云图显示的有问题（是个动画，这里没发放，我就大概描述一下）。这个空化现象是有周期性的（即空化在前缘产生，到生长，最后脱落），在第一个周期压力显示正常，即低压区都是空化区，第二个及后几个周期的压力开始乱跳，有些远离水翼壁面的非空化区压力也非常低，但没有空化产生，或者局部压力突然升高？（不知道是不是破碎产生的高压，我认为这个暂时不用考虑）

还请大家能给出意见，谢谢

[yhdthu]

另外，我在查看interPhaseChangeFoam求解器代码时，一直搞不太清楚在求解压力时引入的源项的演绎过程，这个隐式求解是怎么搞的呢？为什么表面张力也不见了呢？

        fvScalarMatrix p_rghEqn
        (
            fvc::div(phiHbyA) - fvm::laplacian(rAUf, p_rgh)
          - (vDotvP - vDotcP)*(pSat - rho*gh) + fvm::Sp(vDotvP - vDotcP, p_rgh)
        );

[队长别开枪]

@yhdthu
连续方程

动量插值

[yhdthu]

非常感谢，我明白了，因为这个模型中的连续方程不是divergence(U)不为0～

[yhdthu]

@队长别开枪 你好，我还有个疑问，关于相方称中的人工压缩项，我发现在fvsolutio中一般是关闭的（cAlpha＝0），我试着打开了这项算了一下（令其值为1），结果发现，残差变得稳定很多，但是计算结果（云图）效果非常的不好，在interface处呈现锯齿状的结构，请问这是怎么回事呢？另外其他收敛曲线（p_rgh，k，epsilon）都一直在震荡，请问如何解决震荡的问题呢？

PS：我用的是rke湍流模型，sauer空化模型。

[CFD中文网]

表面张力在速度方程中

我很久之前在cfd-online上看到一位外国朋友分享的p_rgh和rhogh的理解和云图计算结果（不好意思找了很久找不见了）。 显示这两者可以是不连续的，但是二者加和后得的p却是连续的，这东西跟选择的参考系有关系，比如我把网格放在“水平面以下”，这时候rho*gh就有物理意义。另外我在看代码的时候，也没搞清楚怎么解着解着就得出了连续的p？还请前辈指教

举例：鼓泡床越往下压力越大，压力初始这么给定比较好：
p : 2 3 4 5 6 7 8
有：
rhogh: 1 2 3 4 5 6 7
我们给定p-rhogh边界条件这样就可以了
p-rhogh:1 1 1 1 1 1 1
是不是简单很多？

[yhdthu]

@cfd-china 嗯那 我看到了，但这项不是属于在interface附近的表面力项，类似于体积力g也加入到压力更新么？

另外这个参考高度z，取得是xoy面么？如果是自动取模型最低端肯定是对的，但是如果模型随意画的那岂不是有问题了？？

[队长别开枪]

@yhdthu
不可压流动中驱动流体流动的是压力梯度，一般对流场中的一个特定位置设定参考位置和参考压力值。如果包含第一类压力边界条件，使用p_rgh进行设置更为方便。而且多相流中，使用p_rgh能使得计算更加稳定和准确，非得使用p的话，就得像fluent使用诸如body-force-weighted等特殊的压力插值格式计算压力梯度了。在不同的流体交界面附近，压力p是连续的，但是压力梯度是非连续的，因为两边的流体密度不一样。
参考高度的不同只会影响p的值，不会影响其梯度的值，因为它们之间相差一个梯度为零的等值压力场。
最后，体积力项包含在p和p_rgh的转换关系里，详细参考东岳流体文档 http://dyfluid.com/interFoam.html ，公式(22)-(25)。

[队长别开枪]

@yhdthu
代数方法求解体积分数方程我不熟悉啊，我用的是自己的几何重构VOF模块+openfoam流动求解器。

[yhdthu]

@队长别开枪  好的，谢谢你的回复，我再看看