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

[yhdthu]

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

[李东岳]

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

[mohui]

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

[yhdthu]

@李东岳 前辈好～模型大概如下：

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

界面速度

liquid速度

压力间断

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

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

[yhdthu]

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

[mohui]

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

[yhdthu]

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