请教关于基于压力的求解方法用于可压缩问题时,公式如何推导



  • 像PISO算法,一般是用于求解不可压问题的。但经过一定的改造,是可以求解可压问题的,比如openfoam中的rhoPimpleFoam。
    但我现在不是很明确速度预测方程,以及两次修正步的压力方程应该怎么写。大家能否为我解惑,或者提供一下这方面的权威文献?


  • 网格教授 OpenFOAM教授 管理员

    这个不是在icoFoam里面已经有了标准代码了么?
    要看PISO,找找关键词R I Issa 1987 PISO



  • @李东岳 感谢东岳老师的回复
    但我还是不太明白。我再把疑问写得具体一点,请您提出批评指正

    按照Issa所述,可压时的速度预测方程为:
    0_1504229607432_图片1.jpg
    可以整理为:
    0_1504229710957_图片3.jpg
    openfoam中通常令:
    0_1504229930771_图片5.jpg

    而经过第一次修正后的速度,密度应当满足:
    0_1504230048386_图片6.jpg
    整理为:
    0_1504230330752_图片10.jpg
    右侧的相邻单元修正未知,u_n两星 用预测步得到的 u_n星 代替:
    0_1504230710175_图片11.jpg

    现在眼看着就能将它代入连续性方程:
    0_1504230803970_傲游截图20170901095312.png
    但还有一个问题:
    0_1504230935192_图片12.jpg
    该项中还存在着一个未知的密度rho星,以至于和预测步公式的形式上产生了不统一,不能凑出HbyA了,

    Issa 在 Solution of the Implicitly Discretised Fluid Flow Equations by Operator-Splitting中的公式(46)似乎直接用n时间步的rho代替rho星了?还是我上面的推导思路不对?


  • 网格教授 OpenFOAM教授 管理员

    目前我在出差,涉及到公式的以及翻文献的还讨论不了,很抱歉!等我出差回北京我看看。
    希望论坛其他人有思路。



  • @李东岳 好的,谢谢!