同位网格下的压力计算问题


  • OpenFOAM讲师

    最近做了一个一维蒸发的测试算例,发现of算出来的速度和相分数分布是对的,但是压力分布却是错的,记得之前跟老师交流过,同位网格很容易出现这个问题,但是想不太通为什么会这样,哪位可以指导一下?~



  • 什么样的物理过程?发云图看看?显示cell value



  • @yhdthu 现在不都是使用同位网格吗,用rhie-chow插值方法可以避免棋盘型网格出现。压力分布出错是不是和压力源项的添加有关呢?


  • OpenFOAM讲师

    @李东岳 前辈好~模型大概如下:
    0_1510880579824_WechatIMG16443.jpeg

    初始状态静止,vapor侧的壁面加热,liquid蒸发界面开始向上移动,界面速度和水相速度都有解析解,另外还有间断压力。

    界面速度
    0_1510880898762_WechatIMG16444.jpeg

    liquid速度
    0_1510880905880_WechatIMG16445.jpeg

    压力间断
    0_1510880911743_WechatIMG16446.jpeg

    我用的是interPhaseChangeFoam算的,自己改了改源项,发现除了压力其他能对上解析解。

    2_1510881232808_WechatIMG16449.jpeg 1_1510881232807_WechatIMG16448.jpeg 0_1510881232807_WechatIMG16447.jpeg

    我想,由于解的是one-field方程,速度和压力都是mixed的,应该是光滑的;另外没做界面重构,interface无法keep sharp,这会导致压力方程源项配置出问题,导致压力计算错误。


  • OpenFOAM讲师

    @mohui 同位网格下的cell velocity和surface velocity更新不一致,尤其在多相流计算中体现更明显,我想这是因为间断边界处理不当导致的



  • @yhdthu 简短边界是指气液界面处吗?另外看你的压力图哪里不对呢?


  • OpenFOAM讲师

    @mohui 是的,interface处是有间断的,但是one-field形式的控制方程体现不出来间断
    压力分布不对,看一下理论解就明白了~


Log in to reply