质量守恒定律与连续性方程
质量守恒定律是流体力学的基石之一，它直接导出描述流体运动的连续性方程。下面从理论出发，系统阐述质量守恒定律的含义、连续性方程的推导过程以及在实际工程中的常见应用。
- 全局形式：对任意固定控制体 (V)（边界为 (\partial V)），流体质量的总量随时间的变化等于进入控制体的质量通量减去离开控制体的质量通量。用数学语言可写成
[
\frac{d}{dt}\int_{V}\rho , dV + \oint_{\partial V}\rho \mathbf{u}\cdot\mathbf{n}, dS = 0,
]
其中 (\rho) 为密度，(\mathbf{u}) 为速度向量，(\mathbf{n}) 为外向法向。
- 局部形式：对任意微小体积应用散度定理并利用任意性，可得到
[
\frac{\partial \rho}{\partial t} + \nabla!\cdot(\rho\mathbf{u}) = 0.
]
该式即为 连续性方程 的最一般形式，适用于任意可压缩流体。
连续性方程的推导
- 选取微元体：取坐标系中的一个无穷小体积 (dV = dx,dy,dz)。
- 质量变化：该体积在时间 (\Delta t) 内的质量增量为
[
\Delta m = \rho(x,y,z,t+\Delta t)dV - \rho(x,y,z,t)dV.
]
- 质量通量：沿 (x) 方向左、右两面的质量流出分别为
[
\rho u_x|{x},dy,dz,\Delta t,\qquad
\rho u_x|{x+dx},dy,dz,\Delta t,
]
同理可写出 (y) 与 (z) 方向的通量。
4. 守恒平衡：质量增量等于进入流体的通量减去离开流体的通量，整理后并除以 (\Delta t,dV)，取极限 (\Delta t\to0) 得到局部守恒式。
最终得到的局部连续性方程与上文全局形式等价。
连续性方程的常见形式
- 不可压缩流体（(\rho=\text{常数})）
[
\nabla!\cdot\mathbf{u}=0.
]
- 一维稳态可压缩流（仅随 (x) 变化，(\partial/\partial t=0)）
[
\frac{d}{dx}(\rho u)=0 ;\Longrightarrow; \rho u = \text{常数}.
]
- 轴对称流（使用圆柱坐标 ((r,\theta,z))）
[
\frac{1}{r}\frac{\partial}{\partial r}(r\rho u_r)+\frac{\partial}{\partial z}(\rho u_z)=0.
]
示例一：一维稳态管道流
设管道截面积 (A(x)) 随 (x) 变化，流体密度 (\rho) 为常数，速度为 (u(x))。连续性方程化为
[
\frac{d}{dx}\bigl(\rho u A\bigr)=0 ;\Longrightarrow; \rho u A = \text{常数}= \dot{m},
]
其中 (\dot{m}) 为质量流率。若已知入口条件 ((u_1, A_1))，则任意截面上的速度
[
u(x)=\frac{u_1 A_1}{A(x)}.
]
此结果常用于喷嘴、收缩段以及膨胀段的设计。
示例二：可压缩气体的喷管流动
考虑等熵、定常、轴对称喷管，使用一维连续性方程
[
\rho u A = \dot{m}= \text{常数}.
]
结合能量方程与状态方程 (p=\rho R T)，可得到马赫数随截面面积的关系（诺特-诺特方程）
[
\frac{A}{A^*}= \frac{1}{M}\Bigl[\frac{2}{\gamma+1}\bigl(1+\frac{\gamma-1}{2}M^2\bigr)\Bigr]^{\frac{\gamma+1}{2(\gamma-1)}},
]
其中 (\gamma) 为比热比，(M) 为马赫数，(A^*) 为临界截面积。该公式的推导本质上仍然基于质量守恒。
数值实现要点
- 离散化：在有限体积法中，连续性方程的积分形式直接用于构造质量守恒的数值通量。
- 通量计算：常用的通量求解器包括 Rusanov、HLLC 等，它们保证在数值上满足质量守恒。
- 稳态收敛：对不可压缩流体，常采用压力校正（Projection）或 SIMPLE 系列算法，以消除数值散度误差。
- 网格质量：网格的非结构化或高度扭曲会导致离散散度不为零，需通过网格光滑或高阶重构降低误差。
小结
质量守恒定律是流体力学唯一不容违背的基本原理，它在任意控制体上都必须成立。通过对该定律的局部化，可得到连续性方程——描述流体密度与速度场相互关系的核心方程。根据流体是否可压缩、是否稳态以及几何对称性，连续性方程可以简化为不同的形式，进而在管道设计、喷管分析以及数值模拟中发挥关键作用。掌握其推导过程、典型应用以及数值实现细节，是进行任何流体力学研究和工程设计的前提。