流体方程为何按连续性、运动、能量顺序排列

一、核心问题与学习目标

在教材、论文和仿真文档中，流体基本方程几乎总按连续性方程→运动方程→能量方程的顺序书写。这并非随意约定，而是有深层物理原因。本文面向有微积分与矢量运算基础的开发者和计算流体力学从业者，帮助读者：

• 理解三个方程对应质量、动量、能量的守恒层级；
• 从速度矩展开视角推导压强演化与绝热关系；
• 明确不同气体自由度下绝热指数的差异；
• 掌握数值仿真中变量更新的正确顺序与稳定性要点。

二、排序的物理逻辑

2.1 守恒律层级

方程	守恒量	物理含义
连续性方程	质量	密度变化由速度散度决定
运动方程	动量	速度变化由压强梯度和外力决定
能量方程	能量	温度或内能变化由做功和热传导决定

先确定”有多少物质”，再确定”物质如何运动”，最后确定”运动中能量如何变化”——这正是从低阶矩到高阶矩的自然展开。

2.2 速度矩视角

在统计力学中，设粒子完整速度为 $\mathbf{w}=\mathbf{v}+\mathbf{u}$，其中 $\mathbf{v}$ 为平均流速，$\mathbf{u}$ 为热运动速度（$\langle\mathbf{u}\rangle=0$），$\rho$ 为密度。定义 $S{ij}=\rho\langle u_i u_j\rangle$ 为二阶速度矩张量，$M{ijk}=\rho\langle u_i u_j u_k\rangle$ 为三阶矩。零阶矩对应质量守恒，一阶矩对应动量守恒，二阶矩对应压强与能量演化。方程自然按阶数从低到高排列。

三、从二阶矩到压强方程的推导

3.1 二阶矩演化方程

在 $k=2$ 且不取迹时，可推出：

$\frac{d}{dt}S<em>{ij}=-3\partial_k(v</em>{(i}S<em>{jk)})-\partial_k M</em>{ijk}+\dot{S}^c_{ij}$

第一项描述平均速度对速度弥散的输运，第二项含三阶矩体现更高阶统计影响，$\dot{S}^c_{ij}$ 为碰撞项。若不截断，方程将无限展开，形成层级耦合结构。

3.2 引入各向同性压强假设

假设流体压强各向同性：

$S<em>{ij}=P\delta</em>{ij}$

对指标取迹（$i=j$ 求和），左边得 $\frac{d}{dt}(3P)$；右边第一项缩并后为 $-5P\,\partial_i v_i$，碰撞项在各向同性假设下仅贡献保持形式的部分。于是：

$\frac{d}{dt}(3P)=-5P\,\partial_i v_i$

即：

$\frac{d}{dt}P=-\frac{5}{3}P\,\partial_i v_i$

此结果对应单原子理想气体（自由度 $f=3$，$\gamma=1+2/f=5/3$）。若为双原子气体 $f=5$，则指数为 $\gamma=7/5$；多原子气体 $f=6$ 时 $\gamma=4/3$。读者不可将 $5/3$ 视为通用结论。

3.3 结合连续性方程

连续性方程为：

$\frac{d}{dt}\rho+\rho\,\partial_i v_i=0\quad\Rightarrow\quad\partial_i v_i=-\frac{1}{\rho}\frac{d}{dt}\rho$

代入压强方程：

$\frac{d}{dt}P=-\frac{5}{3}P\left(-\frac{1}{\rho}\frac{d}{dt}\rho\right)=\frac{5}{3}\frac{P}{\rho}\frac{d}{dt}\rho$

整理得：

$\frac{1}{P}\frac{d}{dt}P=\frac{5}{3}\frac{1}{\rho}\frac{d}{dt}\rho$

两边积分：

$\ln P=\frac{5}{3}\ln\rho+C\quad\Rightarrow\quad P\propto\rho^{5/3}$

此即单原子理想气体的绝热关系。若考虑热传导或黏性耗散，指数将偏离 $5/3$；等温过程则指数为 $1$。

3.4 不同自由度气体绝热指数对比

气体类型	自由度 $f$	$\gamma=1+2/f$	绝热关系
单原子（He、Ar）	3	5/3 ≈ 1.667	$P\propto\rho^{5/3}$
双原子（N₂、O₂）	5	7/5 = 1.4	$P\propto\rho^{7/5}$
多原子（CO₂、H₂O）	6	4/3 ≈ 1.333	$P\propto\rho^{4/3}$

四、数值实现：变量更新顺序与稳定性

在编写流体仿真代码时，通常按以下顺序更新：

1. 根据当前速度场，更新密度（连续性方程）
2. 根据密度和压强，更新速度（运动方程）
3. 根据速度和密度，更新压强或能量（能量方程）

不能随意调换的原因： 运动方程依赖密度与压强，能量方程依赖速度散度。若先更新能量再更新速度，将导致压强与速度间出现时间不一致，引发数值不稳定。

稳定性注意事项

• 时间积分方法：压强项建议采用隐式或半隐式格式，避免显式格式在高马赫数下的刚性问题；
• CFL条件：时间步长需满足 $\Delta t \leq C \cdot \Delta x / (|\mathbf{v}|+c_s)$，其中 $c_s$ 为声速；
• 不可压缩流：当 $\partial_i v_i=0$ 时，压强方程退化为泊松方程 $\nabla^2 P = \rho\,\partial_i\partial_j(v_i v_j)$，需用迭代法求解，逻辑仍是从低阶约束到高阶演化。

五、结果验证方法

1. 量纲检查：确认 $P\propto\rho^{\gamma}$ 两边量纲一致；
2. 极限测试：令 $\partial_i v_i=0$，压强应保持不变；
3. 守恒检查：周期性边界下总质量应守恒；
4. 对比解析解：一维绝热膨胀满足 $P/P_0=(\rho/\rho_0)^{\gamma}$，取 $\rho_0=1$、$P_0=1$、$\gamma=5/3$，当 $\rho=0.5$ 时 $P\approx0.315$，可用于代码对比。

六、常见问题

Q：为什么指数是5/3而非3/5？
A：因 $S{ij}=P\delta{ij}$ 取迹得 $3P$，结合连续性方程消去速度散度后，最终指数为 $5/3$，而非其倒数。

Q：碰撞项为何不影响绝热指数？
A：碰撞项被假设为仅维持各向同性、不引入额外能量耗散，故只影响弛豫时间尺度，不改变 $P$ 与 $\rho$ 的幂律关系。

Q：忽略黏性会有什么误差？
A：实际流体中黏性会导致动能转化为热能，使绝热指数偏离理论值。在高雷诺数流动中影响较小，但在边界层或低速流动中不可忽略。

七、总结

流体方程按连续性→运动→能量排序，本质是质量、动量、能量守恒的层级体现，也是速度矩从低阶到高阶的展开。通过推导 $\frac{d}{dt}P=-\frac{5}{3}P\,\partial_i v_i$ 并结合连续性方程，得到单原子气体 $P\propto\rho^{5/3}$ 的绝热关系。掌握这一逻辑，既能正确书写方程，也能指导仿真中的变量更新与稳定性设计。后续可进一步研究热传导、黏性和湍流模型对压强演化的修正。