一维不可压缩流体如何计算压力变化?



  • 一维原始控制方程:
    连续方程:
    \begin{equation}
    \frac{\partial{\rho}}{\partial{t}} + \frac{\partial{\rho u}}{\partial{x}} = 0
    \end{equation}
    动量方程:
    \begin{equation}
    \frac{\partial{\rho u}}{\partial{t}} + \frac{\partial{\rho u u}}{\partial{x}}=-\frac{\partial{p}}{\partial{x}}+f
    \end{equation}
    由于流体不可压,则上述方程可简化为:
    连续方程:
    \begin{equation}
    \frac{\partial{u}}{\partial{x}} = 0
    \end{equation}
    动量方程:
    \begin{equation}
    \rho \frac{\partial{u}}{\partial{t}} + \rho u \frac{\partial{u}}{\partial{x}} = - \frac{\partial{p}}{\partial{x}} + f
    \end{equation}

    方程(3)(4)合并,则就剩一个方程:
    \begin{equation}
    \rho \frac{d u}{d t} = - \frac{\partial{p}}{\partial{x}} + f
    \end{equation}
    那么这样一来,好像就不能用SIMPLE算法求压力了,那该怎么求压力呢?也没有多余的压力方程来求啊。
    求大神解答,感激不尽。



  • 方程3 4的问题就是不存在压力方程。所以才会有SIMPLE迭代算法,需要组建压力方程。



  • @东岳 东岳老师,openFoam最新的版本里面simpleFoam代码是不是就是采用的simplec算法,多出的rauUt似乎就是相比较simple算法多的东西



  • @东岳
    那意思是先用假设的压力在每个网格上解方程(5),然后再看是否满足方程(3),怎么感觉方程(5)天然满足方程(3)啊



  • SIMPLEC在这里:http://dyfluid.com/simplefoam.html
    另外,@Gordonaero 可以参考这个 http://www.cfd-china.com/topic/3395