动量守恒定律在流体中的应用
在流体力学中，动量守恒定律是描述流体运动的基本原理之一。它揭示了在没有外力作用或外力已知的情况下，流体单元的动量变化与作用在其上的力之间的关系。本文将系统阐述动量守恒定律的理论基础、数学表达式以及在实际流体问题中的典型应用，并通过具体示例帮助读者深入理解和灵活运用该定律。
动量守恒定律在连续介质力学中的表述为 控制体动量方程（也称为 欧拉方程 的动量形式）：
[
\frac{d}{dt}\int_{V(t)} \rho \mathbf{v}, dV = \int_{V(t)} \rho \mathbf{g}, dV + \int_{S(t)} \mathbf{T}, dS
]
其中：
- ( \rho ) 为流体密度（kg·m⁻³），在可压缩流体中随位置和时间变化。
- ( \mathbf{v} ) 为流体速度矢量（m·s⁻¹）。
- ( \mathbf{g} ) 为体积力（如重力）加速度矢量。
- ( \mathbf{T} ) 为表面力张量，通常由 应力张量 ( \boldsymbol{\sigma} ) 与单位法向量 ( \mathbf{n} ) 的乘积给出：( \mathbf{T}= \boldsymbol{\sigma}\cdot\mathbf{n} )。
- ( V(t) ) 为随时间移动的控制体，( S(t) ) 为其边界面。
若取固定控制体（欧拉控制体），则可将上式写成 控制体动量守恒方程：
[
\frac{\partial}{\partial t}\int_{V}\rho \mathbf{v}, dV + \int_{S}\rho \mathbf{v}(\mathbf{v}\cdot\mathbf{n}), dS = \int_{V}\rho \mathbf{g}, dV + \int_{S}\boldsymbol{\sigma}\cdot\mathbf{n}, dS
]
该方程表明：流体动量的瞬时变化率 + 动量通量的净输出 = 体力 + 表面力。
2. 常用简化形式
在实际工程计算中，常根据流动特征对上述方程进行简化：
| 条件 | 结果 | 说明 |
|---|---|---|
| 不可压缩、稳态、无体力 | (\displaystyle \int_{S}\rho \mathbf{v}(\mathbf{v}\cdot\mathbf{n}), dS = \int_{S}\boldsymbol{\sigma}\cdot\mathbf{n}, dS) | 动量通量平衡等于表面应力。 |
| 粘性流体、忽略粘性 | (\displaystyle \int_{S}\rho \mathbf{v}(\mathbf{v}\cdot\mathbf{n}), dS = -\int_{S}p\mathbf{n}, dS) | 只保留压力作用，得到 伯努利方程 的动量解释。 |
| 一维管道流 | (\displaystyle \dot{m}(V_2-V_1) = p_1A_1-p_2A_2 + \sum F_{壁}) | (\dot{m}=\rho AV) 为质量流率，(F_{壁}) 为壁面摩擦阻力。 |
3. 动量守恒在典型流体问题中的应用
3.1 喉部流速的计算（喷嘴/收缩管）
问题描述：在理想不可压缩流体通过收缩段（喉部）时，需要求得喉部的流速 (V_2)。
解题步骤：
选取 控制体 包含入口截面 1、喉部截面 2 以及两端的假想闭合面。
对于稳态、不可压缩、无粘性、无外力的情况，动量守恒化为：
[
\rho A_1 V_1^2 - \rho A_2 V_2^2 = p_1 A_1 - p_2 A_2
]结合 连续方程 (\rho A_1 V_1 = \rho A_2 V_2) 消去密度，得到：
[
V_2 = \sqrt{\frac{2(p_1-p_2)}{\rho\left(1-\frac{A_2^2}{A_1^2}\right)}}
]
此公式直接体现了动量守恒对压力差与截面变化的关联。
3.2 涡轮或泵的轴向推力
问题描述：求转子在给定流量、入口/出口压力下的轴向推力 (F_x)。
解题要点：
采用 控制体 包含转子前后两截面。
动量守恒在轴向方向上写为：
[
F_x = \dot{m}(V_{x,,\text{out}}-V_{x,,\text{in}}) + (p_{\text{out}}A_{\text{out}}-p_{\text{in}}A_{\text{in}})
]若已知转子转速、叶片几何，可进一步利用 叶片动能理论（Euler 方程）校核结果。
3.3 船舶推进与阻力
问题描述：在船体后方的螺旋桨产生的推力如何通过动量守恒计算？
步骤概述：
选取 螺旋桨控制体，包括进、出流截面。
设进流速度为 (V_0)，出流速度为 (V_0 + \Delta V)。
质量流率 (\dot{m}= \rho A V_0)。
动量守恒给出推力：
[
T = \dot{m},\Delta V = \rho A V_0 \Delta V
]若考虑压差，则加入压强项 ( (p_{\text{out}}-p_{\text{in}})A)。
此模型是 动量理论（Momentum Theory）在推进器设计中的经典应用。
4. 示例：二维平板冲击流的动量分析
假设有一块宽度为 (b) 的平板垂直放置在均匀流场中，流体密度 (\rho)，入口速度 (V)。平板后方形成分离区，出口速度近似为零。求平板受到的冲击力。
解答：
取控制体为平板前后两截面及包围分离区的闭合面。
由于出口速度为零，动量通量仅来自入口：
[
\int_{S_{\text{in}}}\rho V^2, dA = \rho V^2 b, h
]其中 (h) 为截面高度（若为二维单元可取单位长度）。
表面力仅为压强差，若远场压强相同，则压强贡献为零。
因此冲击力 (F) 为：
[
F = \rho V^2 b h
]
该结果表明冲击力与流体动能密度 (\frac{1}{2}\rho V^2) 成正比，是动量守恒在冲击流中的直接体现。
5. 实际工程计算的注意事项
- 控制体的选取：应尽量使控制面与已知或易于测量的几何、流场相吻合，避免出现未知的剪切应力或压力分布。
- 粘性效应：在高雷诺数流动中可忽略，但在低雷诺数或壁面层分析时必须加入粘性剪切力 (\tau_w)。
- 可压缩性：当马赫数 (Ma>0.3) 时，密度随压强变化显著，需要使用 可压缩动量方程 并配合能量方程求解。
- 数值实现：在 CFD（计算流体力学）中，动量守恒方程是求解 Navier‑Stokes 方程组的核心，离散化时常采用 有限体积法，保证离散后仍满足全局动量守恒。
6. 小结
动量守恒定律是流体力学的基石之一，它通过质量流率、速度、压力以及表面应力的相互作用，提供了分析和设计流体系统的强大工具。本文从理论推导、简化形式到典型工程案例，系统展示了动量守恒在喷嘴、转子、推进器以及冲击流等多类问题中的具体应用。掌握这些基本思路后，读者能够在实际工程中快速构建控制体、列写动量平衡方程，并结合连续方程或能量方程得到所需的关键参数，为后续的详细数值模拟或实验验证奠定坚实基础。