这种情况大家怎么改进的？不会需要找人润色吧，好像得一大笔钱

@东岳 我这里太阳很大

今天工作室讨论到了求解间断的问题，一个马上要毕业的博士师兄竟然说中心差分格式耗散大。更加刷新我三观的是，他的导师竟然竟然竟然附和他，我表示反对后，给他们讲不用中心差分求解间断是因为数值振荡的问题，但他们坚持说是因为耗散大的原因。真是醉了！！！
真想把我看过的所有CFD书拍他们脸上！
整天说谁谁谁听不进别人的意见，这个人师才真的是。
以前以为这位导师水平很高，现在感觉不过尔尔。
吐槽一下，不吐不快！

@lllwonderliquid 非常感谢

问题好，回答也好，受教了

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

@lllwonderliquid
你好，
“做的模拟本质上确实为不可压缩，流体的压缩性是通过气相与液相的相互传质体现的，对于独立的两相而言，不考虑压缩性。”
请问你这个方法流体力学上叫什么？有没有相关的介绍资料？感觉这个方法也可以用在沸腾传热上。想请教一下你。

一维原始控制方程：
连续方程：
\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算法求压力了，那该怎么求压力呢？也没有多余的压力方程来求啊。
求大神解答，感激不尽。

@cheng_xu 文献挺多的，最早的有Numerical simulation of transients in feed systems for cryogenic rocket engines，近期的有Homogeneous Equilibrium Mixture Model for Simulation of Multiphase/Multicomponent Flows.我也打算自己编程算，最近在动手编，可能编的过程会加深理解吧

超临界情况下，教材上讲流体是处于气液不分状态的，例如甲烷，那这个时候还能算两相问题吗？还是直接按单相处理就可以，特别是物性参数这一块。