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

Gordonaero

一维原始控制方程：
连续方程：
\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算法多的东西

Gordonaero

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

李东岳

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