CFD中文网

今天收到一封邮件，邀请我审一本书。

@李东岳 哇，牛杯plus，感谢大佬。我学习一下

大佬，矩方法有什么经典教科书吗？

在思考一个审稿人的意见，太难了

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

@东岳 我这里太阳很大

今天工作室讨论到了求解间断的问题，一个马上要毕业的博士师兄竟然说中心差分格式耗散大。更加刷新我三观的是，他的导师竟然竟然竟然附和他，我表示反对后，给他们讲不用中心差分求解间断是因为数值振荡的问题，但他们坚持说是因为耗散大的原因。真是醉了！！！
真想把我看过的所有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.我也打算自己编程算，最近在动手编，可能编的过程会加深理解吧

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

@东岳 好的，谢谢大佬

@东岳 one dimensional homogeneous equilibrium model

如题，想做一个低温流体动态充填的数值研究。找了很多文献，都是列基本的控制方程，如何实现方程组封闭以及求解的顺序，一直想不通。

@东岳 感谢大佬回答，那意思是液体密度为定值，气体的密度可以根据状态方程得到。那混合物密度怎么求，在含气率未知的情况下。

最近在看两相流，想不出方程组如何封闭，以及如何求的含气率。
含相变的两相流动，用一维均相模型，除了三大控制方程外，看了好多文献，还用到了两相混合密度=气相密度乘含气率+液相密度乘（1-含气率），类似地，还有混合焓、单相焓与质量含气率的关系，含气率又和质量含气率有关系。以上感觉方程无法封闭呀，难道是假设每个单相的密度已知？