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



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



  • rhogh会影响p_rgh的值

    为什么?



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

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

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

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



  • 另外,我在查看interPhaseChangeFoam求解器代码时,一直搞不太清楚在求解压力时引入的源项的演绎过程,这个隐式求解是怎么搞的呢?为什么表面张力也不见了呢?
    fvScalarMatrix p_rghEqn
    (
    fvc::div(phiHbyA) - fvm::laplacian(rAUf, p_rgh)
    - (vDotvP - vDotcP)(pSat - rhogh) + fvm::Sp(vDotvP - vDotcP, p_rgh)
    );



  • @yhdthu
    连续方程
    0_1482041760699_1.JPG
    动量插值
    0_1482041777348_2.JPG
    0_1482041785920_3.JPG



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



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

    PS:我用的是rke湍流模型,sauer空化模型。



  • 表面张力在速度方程中哦。

    我很久之前在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
    是不是简单很多?



  • @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流动求解器。



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


登录后回复
 

与 CFD 中国 的连接断开,我们正在尝试重连,请耐心等待