流函数与速度势函数基础

在流体力学的研究中，为了简化复杂的矢量速度场分析，通常会引入标量函数来描述流体运动。对于二维、不可压缩流体，最常用的两种数学工具便是流函数 (Stream Function) 和 速度势函数 (Velocity Potential)。这两者不仅简化了连续性方程和欧拉方程的求解，还为分析势流（Potential Flow）奠定了理论基础。
流函数 $\psi$ 是专门为满足不可压缩流体的连续性方程而定义的标量函数。

1. 数学定义

对于一个二维不可压缩流动，其连续性方程为：
$$\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$$
为了满足上述方程，我们可以定义一个函数 $\psi(x, y)$，使得速度分量 $u$ 和 $v$ 分别表示为：
$$u = \frac{\partial \psi}{\partial y}, \quad v = -\frac{\partial \psi}{\partial x}$$
通过这种定义，无论 $\psi$ 如何选取，上述连续性方程恒成立。

2. 物理意义

• 流线 (Streamlines)： 流函数的等值线 $\psi = \text{constant}$ 即为流线。这意味着流体质点在运动过程中始终沿着 $\psi$ 的等值线移动，且流线与速度矢量 $\mathbf{u}$ 始终切向平行。
• 流量计算： 两条流线 $\psi_1$ 和 $\psi_2$ 之间通过的单位宽度体积流量 $Q$ 等于两线之间流函数值的差：
  $$Q = \int \mathbf{u} \cdot d\mathbf{n} = \psi_2 - \psi_1$$
  这使得流函数成为计算流量的极其便捷的工具。

速度势函数 $\phi$ 的定义与适用条件

速度势函数 $\phi$ 的引入则基于流动的另一种特性：无旋性 (Irrotationality)。

1. 数学定义

如果流场是无旋的，即速度场的旋度为零 ($\nabla \times \mathbf{u} = 0$)，那么根据矢量分析理论，该速度场可以表示为一个标量函数的梯度：
$$\mathbf{u} = \nabla \phi \implies u = \frac{\partial \phi}{\partial x}, \quad v = \frac{\partial \phi}{\partial y}$$
这里的 $\phi$ 即为速度势函数。

2. 适用条件

速度势函数并非在所有流动中都存在。它存在的必要条件是流场必须是无旋的。在实际物理场景中，远离边界层的高雷诺数流动通常可以近似为无旋流动。

3. 物理意义

速度势函数的等值线 $\phi = \text{constant}$ 被称为等势线。由于速度矢量 $\mathbf{u}$ 是 $\phi$ 的梯度，因此速度矢量始终垂直于等势线。

流函数与速度势函数的内在联系

当一个流动既是不可压缩的又是无旋的，这种流动被称为势流 (Potential Flow)。在这种情况下，$\psi$ 和 $\phi$ 同时存在，且满足特定的数学关系。

1. 柯西-黎曼方程 (Cauchy-Riemann Equations)

对比 $\psi$ 和 $\phi$ 对速度分量的定义：
$$u = \frac{\partial \phi}{\partial x} = \frac{\partial \psi}{\partial y}$$
$$v = \frac{\partial \phi}{\partial y} = -\frac{\partial \psi}{\partial x}$$
这两组等式正是复分析中的柯西-黎曼方程。这意味着 $\phi$ 和 $\psi$ 互为共轭调和函数。

2. 拉普拉斯方程

将上述关系代入连续性方程或无旋条件，可以导出 $\phi$ 和 $\psi$ 均满足拉普拉斯方程：
$$\nabla^2 \phi = 0, \quad \nabla^2 \psi = 0$$
这意味着势流的求解问题转化为了求解拉普拉斯方程的边界值问题，可以使用叠加原理（如将源、汇、涡叠加）来构建复杂的流场。

3. 几何关系

在势流中，流线 ($\psi = \text{const}$) 与等势线 ($\phi = \text{const}$) 在每一点都正交 (Orthogonal)。这种正交网格为流场的可视化提供了极大的便利。

综合示例：均匀流分析

考虑一个沿 $x$ 轴方向的二维均匀流，其速度场为 $u = U$（常数），$v = 0$。

1. 求解速度势 $\phi$：
   $$\frac{\partial \phi}{\partial x} = U \implies \phi = Ux + f(y)$$
   $$\frac{\partial \phi}{\partial y} = 0 \implies f'(y) = 0 \implies \phi = Ux + C$$
2. 求解流函数 $\psi$：
   $$\frac{\partial \psi}{\partial y} = U \implies \psi = Uy + g(x)$$
   $$-\frac{\partial \psi}{\partial x} = 0 \implies g'(x) = 0 \implies \psi = Uy + C$$
3. 结论：
   • 等势线 $\phi = Ux$ 是垂直于 $x$ 轴的直线。
   • 流线 $\psi = Uy$ 是平行于 $x$ 轴的直线。
   • 两者显然正交，且均满足 $\nabla^2 \phi = 0$ 和 $\nabla^2 \psi = 0$。

总结与对比

为了方便记忆，下表总结了两种函数的关键区别：