渦度方程式
理论与物理
什么是涡量
老师，涡量方程这名字听起来就很难，首先“涡量”到底是什么？
涡量 $\boldsymbol{\omega}$ 是表示流体局部旋转的矢量量，定义为速度场的旋度（curl）。
在二维情况下，它变成一个标量 $\omega = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}$。这是量化流体微元旋转程度的量。
原来如此，并不需要流体整体旋转，只要局部有旋转就有涡量，是这个意思吗？
没错。例如，即使是直线剪切流 $u = ky$、$v = 0$，也有 $\omega = -k \neq 0$，所以存在涡量。也就是说，涡量是与“涡”的存在不同的概念。
涡量方程的推导
那么，描述涡量随时间变化的方程是怎么推导出来的呢？
对Navier-Stokes方程两边作用 $\nabla \times$ 即可得到涡量方程。对于不可压缩流体，结果如下。
左边是涡量的物质导数（拉格朗日导数），表示跟随流体粒子观察到的涡量变化率。
右边的两项分别有什么含义呢？
第一项 $(\boldsymbol{\omega} \cdot \nabla)\mathbf{u}$ 是涡管拉伸/倾斜项（vortex stretching）。在三维流动中，当涡管被拉伸时其横截面变细，根据角动量守恒，涡度会增大。龙卷风或浴缸排水涡变细时旋转加快，就是这个效应。
第二项 $\nu \nabla^2 \boldsymbol{\omega}$ 是粘性扩散项，表示涡度通过分子粘性扩散的效果。运动粘性系数 $\nu$ 越大，涡度扩散消散得越快。
二维情况下涡管拉伸项会消失，对吧？
观察得很仔细。在二维情况下，$\boldsymbol{\omega}$ 只有 $z$ 方向分量，所以 $(\boldsymbol{\omega} \cdot \nabla)\mathbf{u} = 0$，涡量方程就变成了简单的对流扩散方程。
这意味着涡量只能在边界面上产生。
涡量的生成机制
涡量不是凭空产生的吗？
对于不可压缩、正压流体，涡量不会在流体内部产生。涡量的主要来源如下。
- 固体壁面: 无滑移条件导致壁面产生速度梯度，从而产生涡量。Lighthill（1963）提出的壁面涡量通量表示为 $\nu \frac{\partial \omega}{\partial n}\big|_{wall}$
- 斜压效应: 当密度梯度与压力梯度不平行时，由项 $\frac{1}{\rho^2}(\nabla \rho \times \nabla p)$ 产生涡量。海洋的热盐环流就是一个很好的例子
- 体积力的非均匀分布: 例如旋转坐标系中的科里奥利力等
在CFD中，要准确捕捉壁面的涡量该怎么做？
壁面附近的网格分辨率至关重要。需要选择将壁面第一层网格的 $y^+$ 控制在1以下，或者使用壁面函数从壁面剪切应力近似涡量通量。
与Kelvin环量定理的关系
涡量方程和Kelvin环量定理是怎么联系起来的？
对于无粘性、正压流体且仅受有势力作用的情况，沿闭合曲线的环量 $\Gamma = \oint \mathbf{u} \cdot d\mathbf{l}$ 是守恒的。根据Stokes定理，$\Gamma = \int_S \boldsymbol{\omega} \cdot d\mathbf{S}$，这对应于涡量方程中粘性项和斜压项为零的情况。
实际的CFD中有粘性，所以环量不守恒，对吧。涡量方程就像是环量定理的推广吗？
这个理解是正确的。涡量方程是包含了粘性和非正压效应的、更一般的描述。
涡量方程的历史——亥姆霍兹涡定理（1858年）与涡管守恒
涡量（Vorticity）的概念和支配涡管（Vortex Tube）行为的“亥姆霍兹涡定理”由Hermann von Helmholtz于1858年发表。该定理的本质：①涡管作为流体粒子线移动（涡强度随对流而守恒），②涡管不能终止于端点——要么形成闭合环，要么终止于流体边界。这也解释了为什么龙卷风无法脱离地面，以及海面可以“吸起”涡旋。开尔文（Kelvin勋爵）将其发展为“开尔文环量定理（在正压完全流体中，物质曲线的环量Γ=const.）”。该定理是螺旋桨和机翼升力生成机制（附着涡）的根本原理，也是现代涡格法（VLM）和涡面元法等CFD代码的数学基础。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。一开始水流会不稳定地喷溅，过一会儿才变成稳定水流，对吧？描述这个“变化过程中”的就是时间项。心脏搏动导致血流脉动，发动机阀门开闭引起流动变化，这些都是非定常现象。那么定常分析是什么？就是只看“经过足够时间、流动稳定之后”——也就是令此项为零。计算成本会大幅降低，因此先用定常求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶扔进河里会怎样？会被水流带着往下游漂，对吧？这就是“对流”——流体的运动搬运物体的效果。暖风的暖气能送到房间另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——这项包含“速度×速度”，因此是非线性的。也就是说，流速越快，这项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置不管的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水，对吧？因为蜂蜜的粘性（$\mu$）高，所以不易流动。粘性越大，扩散项越强，流体的运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，在Re数大的流动中，对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：推注射器的活塞，液体就会从针头有力地射出，对吧？为什么？因为活塞侧压力高，针头侧压力低——这个压力差就是推动流体的力。大坝泄洪也是同样原理。天气图上等压线密集的地方会怎样？没错，会刮强风。“有压力差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD中的“压力”多为表压而非绝对压力。切换到可压缩分析时结果突然出错，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：被加热的空气会上升——为什么？因为比周围空气轻（密度低），被浮力推上去了。这个浮力作为源项添加到方程中。此外，燃气灶火焰产生化学反应热、工厂电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，都用源项表示。忘记源项会怎样？自然对流分析中忘记加入浮力，流体就完全不动——就像冬天房间里开了暖气，暖空气却不上浮，得到这种物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3 时）：将密度视为常数。马赫数0.3以上需考虑可压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要捕捉激波）、自由表面流动（需要VOF/Level Set等方法）
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 入口条件中从体积流量换算时，注意截面积单位 |
| 压力 $p$ | Pa | 区分表压与绝对压力。可压缩分析使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气: 约1.225 kg/m³@20°C、水: 约998 kg/m³@20°C |
| 粘性系数 $\mu$ | Pa·s | 注意与运动粘性系数 $\nu = \mu/\rho$ [m²/s] 混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转换的判断指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
涡量-流函数法
用CFD求解涡量方程时，最基本的方法是什么？
对于二维不可压缩流动，涡量-流函数法（vorticity-stream function method）是经典方法。联立求解涡量方程和流函数的泊松方程。
速度通过 $u = \frac{\partial \psi}{\partial y}$、$v = -\frac{\partial \psi}{\partial x}$ 恢复。这种方法的优点在于压力从未知数中消失了。
和直接求解N-S方程相比有什么优势？
在二维情况下，未知数只有 $\omega$ 和 $\psi$ 两个，避免了压力-速度耦合问题（无需SIMPLE法等迭代）。但扩展到三维时涡量变成3分量矢量，变得复杂，所以实际中不太使用。
空间离散化
空间方向的离散化怎么做？
来看看有限差分法（FDM）中典型的离散化。使用均匀网格间距 $h$ 的中心差分，泊松方程为
对流项若用中心差分，当网格雷诺数 $Re_h = |u|h/\nu > 2$ 时会出现数值振荡。为避免此问题，常使用迎风格式或QUICK法。
有限体积法（FVM）的话会有什么不同？
FVM中需要评估单元界面上的涡量通量。对流通量使用二阶精度迎风格式（Second Order Upwind）或TVD格式，扩散通量使用中心差分。OpenFOAM的scalarTransportFoam求解器是用FVM求解涡量对流扩散的参考实现。
时间积分
时间方向怎么推进？
总结一下代表性的时间积分格式。
| 格式 | 精度 | 稳定性 | 特点 |
|---|---|---|---|
| 前向欧拉 | 1阶 | 条件稳定 | 受限于 $\Delta t < h^2/(4\nu)$ |
| Crank-Nicolson | 2阶 | 无条件稳定 |