涡度方程的先端研究
涡法（Vortex Method）
🧑🎓涡度方程用不着网格也能解吧？
🎓涡法（Vortex Method）是拉格朗日方式追踪离散涡要素的方法。涡度场用离散点涡或涡团的集合表达：
$$ \boldsymbol{\omega}(\mathbf{x}, t) \approx \sum_{p=1}^{N} \Gamma_p \, \zeta_\sigma(\mathbf{x} - \mathbf{x}_p(t)) $$
这里 $\Gamma_p$ 是涡要素的循环，$\zeta_\sigma$ 是核直径 $\sigma$ 的正则化核（如Gauss核 $\zeta_\sigma(r) = \frac{1}{2\pi\sigma^2} e^{-r^2/(2\sigma^2)}$）。
🧑🎓没网格怎样算速度场？
🎓用Biot-Savart律从涡度恢复速度。二维里：
$$ \mathbf{u}(\mathbf{x}) = -\frac{1}{2\pi} \sum_{p=1}^{N} \Gamma_p \frac{(\mathbf{x} - \mathbf{x}_p)^\perp}{|\mathbf{x} - \mathbf{x}_p|^2 + \sigma^2} $$
朴素算法$O(N^2)$复杂度，但用FMM（Fast Multipole Method）能降到$O(N)$。Greengard & Rokhlin（1987）的成果。
🧑🎓涡法适用什么问题？
🎓涡局限的外部流动问题。翼后流、喷流混合、涡环碰撞等。反过来，有壁面的内部流或乱流壁层解析弱，需网格法-涡法混合研究中。
涡度和乱流的级联
🧑🎓涡度方程和乱流理论怎么关联？
🎓乱流能级串联是从涡角度理解的：大尺度涡通过vortex stretching生成小尺度涡，最终粘性散逸。
$$ \mathcal{E} = \frac{1}{2} \int |\boldsymbol{\omega}|^2 \, dV $$
enstrophy（涡度的平方积分）时间演化为：
$$ \frac{d\mathcal{E}}{dt} = \int \omega_i \omega_j S_{ij} \, dV - \nu \int |\nabla \boldsymbol{\omega}|^2 \, dV $$
第一项是stretching生成enstrophy，第二项是粘性散逸。3D乱流中生成可能超散逸，导致enstrophy有限时爆炸？（有限时特异问题）
🧑🎓有限时特异解决了吗？
🎓这是千年大奖问题之一。3D N-S方程光滑解的有限时奇点性未解。超高分辨率DNS（$4096^3$格子）数值调查中，但确定答案未得。
DNS、LES和涡度
🧑🎓DNS是直接求解涡度方程吗？
🎓通常DNS求原始变量（速度-压力），然后后处理算涡度。不过涡度基础定式DNS也有研究，Kim & Moin（1985）开创。
LES中网格以下涡用子网格模型表达。WALE（Wall-Adapting Local Eddy-viscosity）模型用速度梯度张量的对称无迹部分，壁近自然趋零，优点突出：
$$ \nu_t = (C_w \Delta)^2 \frac{(S_{ij}^d S_{ij}^d)^{3/2}}{(S_{ij} S_{ij})^{5/2} + (S_{ij}^d S_{ij}^d)^{5/4}} $$
🧑🎓最近的研究动向？
🎓关注的研究热点：
PINN（Physics-Informed Neural Network）：涡度输运方程纳入损失函数的深学习流场预测
涡度约束法（Vorticity Confinement）：数值扩散衰减的涡人为补强。航机后流解析已实用
格子Boltzmann法融合：LB法算速度，高效计算涡度研究
GPU加速涡法：FMM的GPU实现，$10^7$量涡要素实时计算可能
🧑🎓PINN用涡度方程的好处？
🎓实验数据只有部分速度时，涡度方程这个物理约束能填空白。PIV数据补完和超分辨已有应用。
Coffee Break 闲话
涡度方程与天气预报——100年的挑战
涡度方程从1900年代初就用于气象学。挪威气象学家Vilhelm Bjerknes在1904年提议「计算大气涡度变化就能预报天气」。但那时无计算机，手算一个「6小时先预」要6个月！现代超算每秒数十亿次算涡度方程，2周先预报能在72小时内出。涡度先端研究至今与天气预报精度进步直结。