涡旋形成模拟器 返回
流体力学 / CFD

涡旋形成模拟器

点击/拖动"搅拌"流体,生成涡旋结构。基于朗肯涡模型实时计算速度场,粒子追踪可视化,内置卡门涡街与涡对预设。

控制面板

🖱 拖动:正转涡 | 右键+拖动:反转涡
0.008
200
30
0.60

显示选项

预设场景

操作

统计信息

活动涡数
0
总环量 Γ
0.0
粒子数
0
FPS
0

理论说明

朗肯涡模型
涡元 i 位于 (xi,yi),强度 Γi,在点 (x,y) 处诱导速度:
u = −Γ(y−yi) / 2π(r²+ε²)
v = +Γ(x−xi) / 2π(r²+ε²)

开尔文环量定理
无黏、不可压缩 → 有效环量 Γ = ∮v·ds 守恒

CAE应用
涡方法(Vortex Method)、涡激振动(VIV)、湍流建模
拖动以生成涡旋

什么是涡旋形成模拟器

🧑‍🎓
这个模拟器里,我点一下屏幕,怎么就出现一个转圈圈的东西了?那是什么呀?
🎓
简单来说,你点一下,就是在流体里“戳”出了一个涡旋。它就像你搅动咖啡时出现的漩涡。在实际工程中,比如飞机机翼尾部脱落的涡旋,或者风吹过电线时产生的涡旋,原理都类似。你试着在模拟器里多点几个地方,或者拖动鼠标画一条线,就能看到多个涡旋是怎么生成和相互影响的。
🧑‍🎓
诶,真的吗?那旁边那个“黏性”滑块是干嘛的?我调大了它,涡旋好像很快就散开了。
🎓
没错!你观察得很准。“黏性”在这里代表流体的“粘稠度”或者说内摩擦。简单来说,黏性越大,流体越“粘”,涡旋受到的阻力就越大,所以衰减得越快。比如,蜂蜜里的涡旋就比水里的消失得快得多。你试着把黏性滑块拉到最大,再生成一个涡旋,看看它是不是几乎瞬间就消失了?这就是模拟流体阻尼的效果。
🧑‍🎓
哦!那“涡核半径”这个参数又是什么意思?我调小它,感觉涡旋中心那个转得特别快的区域变小了。
🎓
完全正确!涡核半径 $\epsilon$ 定义了涡旋“眼睛”的大小。在涡核内部,流体像刚体一样旋转;外部则像势流。如果没有这个核心半径,理论上涡心速度会变成无穷大,这显然不现实。工程现场常见的是,龙卷风的涡核半径很小,破坏力就集中在中心;而大气中的气旋涡核很大。你改变这个参数后,会看到速度分布和粒子旋转的紧密程度都跟着变化,非常直观。

物理模型与关键公式

本模拟器基于朗肯涡模型来计算速度场。该模型巧妙地避免了涡心速度奇点问题,将涡旋分为内部刚体旋转区和外部势流区。

$$v_{\theta}(r) = \begin{cases}\frac{\Gamma r}{2\pi \epsilon^2}& \text{for }r \le \epsilon \quad \text{(刚体旋转区)}\\ \frac{\Gamma}{2\pi r}& \text{for }r > \epsilon \quad \text{(势流区)}\end{cases}$$

其中,$v_{\theta}$ 是切向速度,$r$ 是到涡心的距离,$\Gamma$ 是涡旋的环量(强度),$\epsilon$ 就是你刚才调整的涡核半径。

为了模拟涡旋的耗散(即你看到的“散开”),我们引入了黏性衰减项。这通过一个随时间衰减的涡量强度来实现。

$$\Gamma(t) = \Gamma_0 \cdot e^{-\nu t}$$

这里,$\Gamma_0$ 是初始环量,$\nu$ 是与你操作的“黏性(衰减率)”参数相关的衰减系数,$t$ 是时间。这个公式描述了涡旋强度如何像钟摆一样,因摩擦而逐渐减弱。

现实世界中的应用

航空航天工程:飞机机翼在产生升力的同时,翼尖会拖出强烈的翼尖涡旋。这些涡旋对后续飞机的飞行安全构成威胁(尾流湍流),因此在机场起降时需要严格间隔。通过模拟,可以优化机翼形状以减弱涡旋强度。

土木与桥梁工程:风吹过桥梁的桥墩或缆索时,会形成交替脱落的卡门涡街(模拟器中有预设)。这会产生周期性的作用力,可能引发桥梁的涡激振动,历史上曾导致美国塔科马海峡大桥的垮塌。模拟有助于评估和设计抗风结构。

能源与动力工程:在风力发电机的叶片设计、水轮机转轮设计中,控制涡旋的生成和脱落至关重要。不合理的涡旋会导致效率下降、噪音增大和结构疲劳。通过CAE模拟可以优化叶片型线,提升能量捕获效率。

气象与环境科学:龙卷风、热带气旋(台风/飓风)本质上都是大规模的 atmospheric vortex。虽然尺度天差地别,但其核心的涡旋动力学原理是相通的。研究其形成和运动规律,对于灾害预报和防灾减灾有重要意义。

常见误解与注意事项

开始使用这个工具时,有几个需要注意的地方。首先,人们常常认为“鼠标移动得越快就能产生越强的涡旋”,但实际上并非如此简单。这个模拟器是沿着鼠标轨迹赋予一种称为“环量Γ”的涡旋强度。也就是说,相比移动的“速度”,移动“轨迹所包围的面积”影响更大。缓慢但画一个大圆,有时比快速小幅度移动能产生更强的涡旋。这在实际工作中也提醒我们,在设置边界条件时要注意不要混淆“速度”和“涡度”。

其次是参数“粘性”与“涡核半径ε”的关系。两者都涉及涡旋的扩散方式,看起来效果相似,但物理意义完全不同。“粘性”是流体本身的性质,表示涡旋能量耗散为热(衰减)的过程。而“涡核半径ε”是计算模型上的参数,决定以多大尺度平滑地表现涡核中心区域的结构。例如,将ε极端调小(如0.01)并将粘性设为零,涡旋会变得非常尖锐且几乎不衰减,容易导致数值不稳定(涡旋强度爆发性增长)。在实际的CFD工作中,区分理解用于建模的参数与物性值至关重要。

最后,别忘了这是一个“二维”模拟的根本限制。屏幕是平的,对吧?现实中的流动是三维的,因此这里看起来漂亮的涡丝,实际上会作为涡管复杂地伸展或缠绕。即使在这个工具中再现了卡门涡街,实际圆柱绕流中也会迅速产生三维的紊乱(展向波动)。二维计算对于抓住现象本质非常有效,但在寻求定量值的实际工程分析中,必须始终考虑三维效应,这是一个主要的陷阱。

相关的工程领域

这种基于“涡方法”的思想,被应用于CAE流体分析的各个领域。首先要提到的是航空航天工程。从飞机机翼分离的涡旋(翼尖涡)是影响后续飞机的危险乱流的原因。为了高效追踪这类大型涡结构,历史上一直在研究作为拉格朗日粒子法一种的“涡方法”。此外,直升机旋翼周围的流动也是复杂的涡旋纠缠。

另一个是汽车工程,特别是空气动力学设计。车体后方产生的湍流(尾流)是空气阻力的主要原因。为了把握这里产生的大型涡结构,会使用DES(分离涡模拟)之类的高级湍流模型。DES是一种混合方法,在物体附近使用RANS(雷诺平均),在分离涡发展的区域使用LES(大涡模拟)求解,其思想接近于直接计算这个工具中所见的“清晰涡旋”行为。

更令人意外的是,它还与工厂管道设计相关。从泵排出的流体通过管道弯头时,会产生称为二次流的涡旋。这会扰乱下游流量计的测量精度,或促进特定部位的腐蚀。预测这类内部流动的涡旋也是CFD的重要工作,其基础的涡度概念是相通的。

为了深入学习

如果你熟悉了这个工具并“想了解更多”,强烈建议下一步学习“涡量方程”。这个从纳维-斯托克斯方程推导出的方程,直接描述了涡旋“生成·输运·扩散·消亡”的机制。在二维情况下,对于有粘性的流体,形式如下:

$$ \frac{\partial \omega}{\partial t} + (\vec{v} \cdot \nabla) \omega = \nu \nabla^2 \omega $$

左边第二项对应涡的输运(本模拟器中用粒子可视化的部分),右边对应粘性引起的涡扩散(通过参数“粘性”调整的部分)。从这个方程可以看出,三维情况下会加入更复杂的“涡拉伸”项。教材方面,《流体力学(上)》(今井功)中关于涡运动的章节非常清晰明了。

实践学习方面,可以一边使用这个模拟器,一边以假设→验证的循环记录改变参数时的行为。例如,“增大涡核半径ε,两个接近的涡旋合并会变快还是变慢?”等等。培养这种感觉后,尝试用工程级的CFD软件(也有如OpenFOAM等免费软件)运行“二维空腔流动”等基础教程,虽然会出现网格和求解器等新要素,但你应该能体会到背后涡的物理本质并未改变。接下来的主题,可以进阶学习湍流的能量级串,或者流固耦合(FSI)基础,这样就能理解在这个工具中观察到的现象,如何发展成实际的“结构物振动”问题。