恒定流与非恒定流

在流体力学中，流动状态是描述流体随时间和空间变化的关键属性。依据流体的宏观属性随时间是否变化，可将流动划分为**恒定流（Steady Flow）和非恒定流（Unsteady Flow）**两大类。

2. 控制方程的时间项区别

2.1 连续方程

• 恒定流
  [
  \nabla\cdot(\rho \mathbf{V}) = 0
  ]
  只保留空间散度项，适用于密度常数或可压缩流的稳态情形。

• 非恒定流
  [
  \frac{\partial \rho}{\partial t} + \nabla\cdot(\rho \mathbf{V}) = 0
  ]
  时间导数项描述了局部密度随时间的积累或稀释。

2.2 动量方程（Navier–Stokes）

• 恒定流
  [
  \rho (\mathbf{V}\cdot\nabla)\mathbf{V}= -\nabla p + \mu \nabla^{2}\mathbf{V}+ \rho \mathbf{g}
  ]
  左侧仅为对流加速度。

• 非恒定流
  [
  \rho \left(\frac{\partial \mathbf{V}}{\partial t}+(\mathbf{V}\cdot\nabla)\mathbf{V}\right)= -\nabla p + \mu \nabla^{2}\mathbf{V}+ \rho \mathbf{g}
  ]
  额外的局部加速度 \(\partial \mathbf{V}/\partial t\) 捕捉了启动、停机等瞬态效应。

2.3 能量方程

• 恒定流
  [
  \rho \mathbf{V}\cdot\nabla h = \frac{Dp}{Dt} + \Phi + \dot{q}
  ]
  只含对流输运。

• 非恒定流
  [
  \rho \frac{\partial h}{\partial t}+ \rho \mathbf{V}\cdot\nabla h = \frac{Dp}{Dt} + \Phi + \dot{q}
  ]
  时间项反映了内部能量的瞬时变化。

3. 判别与分析方法

1. 实验观测：使用压力传感器、流速计等在固定点记录随时间的信号；若信号波动幅度在测量误差范围内，可视为恒定。
2. 数值模拟：在 CFD 中设置 steady‑state（稳态）或 transient（瞬态）求解器；若残差在迭代过程中趋于零且不随时间步长变化，则为恒定流。
3. 特征时间尺度：比较流动特征时间 \(t_c\) 与外部激励时间 \(t_e\)。若 \(t_c \ll t_e\)，流动可近似为恒定；反之则必须采用非恒定模型。

4. 典型案例

4.1 恒定管道流

• 问题描述：水在水平圆管中以恒定入口流速 \(V_{in}\) 进入，管长远大于入口扰动的衰减长度。
• 分析步骤
  1. 采用 泊肃叶（Poiseuille） 公式：
     [
     \Delta p = \frac{32 \mu L V_{avg}}{D^{2}}
     ]
  2. 计算雷诺数 \(Re = \frac{\rho V_{avg} D}{\mu}\)，判断是否为层流。
  3. 若 \(Re<2300\)，使用层流解析解；若 \(Re>4000\)，采用经验摩阻系数 \(f\)。
• 结果：压力降随管长线性增长，流速分布在横截面上呈抛物线形，整个过程满足 \(\partial/\partial t = 0\)。

4.2 非恒定储罐排放

• 问题描述：一圆形储罐底部开孔，液体在重力作用下排出，排放过程从满罐到空罐。
• 控制方程（质量守恒）
  [
  \frac{dV(t)}{dt}= -A_o \sqrt{2g,h(t)}
  ]
  其中 \(V(t)\) 为罐内液体体积，\(A_o\) 为孔口面积，\(h(t)\) 为液面高度。
• 求解步骤
  1. 将体积与液面高度关联： \(V(t)=\frac{\pi D^{2}}{4}h(t)\)。
  2. 代入得到一阶非线性常微分方程：
     [
     \frac{dh}{dt}= -\frac{2A_o}{\pi D^{2}}\sqrt{2g,h}
     ]
  3. 分离变量并积分，得到排放时间公式：
     [
     t = \frac{\pi D^{2}}{4A_o\sqrt{2g}},\bigl(\sqrt{h_0}-\sqrt{h}\bigr)
     ]
• 特征：液面高度随时间呈平方根衰减，明显的时间导数项表明该过程是非恒定流。

4.3 脉冲喷嘴

• 情形：在航空发动机喷嘴中，燃料以脉冲方式喷出，频率可达数千赫兹。
• 分析要点
  • 必须保留 \(\partial \mathbf{V}/\partial t\) 项，使用 瞬态 Navier–Stokes。
  • 采用 LES（大涡模拟） 或 DNS（直接数值模拟） 捕捉高频涡旋的产生与衰减。
  • 通过 Fast Fourier Transform (FFT) 分析压力信号的频谱，验证数值模型的时域精度。

5. 实际工程中的选择原则

1. 时间尺度远大于流动特征时间 → 可采用恒定流假设，简化计算。
2. 启动/停机、调节阀门、周期性激励 → 必须使用非恒定流模型。
3. 数值资源受限 → 先做稳态求解获取基准解，再在关键部位进行局部瞬态细化。
4. 安全与可靠性要求高（如核电冷却、航空发动机） → 全程采用非恒定流仿真，以捕捉潜在的瞬态失稳。

6. 小结

• 恒定流与非恒定流的根本区别在于 是否保留时间导数项。
• 在控制方程层面，这直接导致连续方程、动量方程和能量方程的形式差异。
• 判别流动是否为恒定，需结合实验观测、特征时间尺度以及数值模拟的残差行为。
• 典型工程案例（管道稳态流、储罐排放、脉冲喷嘴）展示了两类流动在理论推导和实际计算中的不同处理方式。
• 合理选择模型能够在保证精度的前提下显著降低计算成本，是流体力学工程实践的关键技巧。