涡流模拟器 返回
流体力学

涡流 模拟器

探索自由涡、强迫涡与兰金涡的流场特性。调节环量 Γ、涡核半径和角速度,实时查看速度分布、压力场及空化风险预测。

涡模型与参数
涡流类型
环量 Γ (m²/s)
m²/s
涡核半径 rᶜ (m)
m
角速度 ω (rad/s)
rad/s
参考压力 p∞ (Pa)
Pa
流体密度 ρ (kg/m³)
kg/m³
Warning:存在空化风险
最低压力低于水的蒸气压(20°C时约2338 Pa),可能发生空化(汽蚀)现象。
计算结果
vθ_max (m/s)
p_min (Pa)
p at r=rᶜ
空化判断
流线与速度矢量图
切向速度 vθ(r) 分布
压力分布 p(r)
Pres
理论与主要公式

自由涡

$$v_\theta = \frac{\Gamma}{2\pi r},\quad p = p_\infty - \frac{\rho\Gamma^2}{8\pi^2 r^2}$$

强迫涡

$$v_\theta = \omega r,\quad p = p_\infty - \tfrac{1}{2}\rho\omega^2 r^2$$

兰金涡(组合)

$$v_\theta = \begin{cases}\omega r & r \lt r_c \\ \dfrac{\omega r_c^2}{r}& r \ge r_c\end{cases}$$

什么是涡流

🙋
“自由涡”和“强迫涡”是什么?听起来好复杂。
🎓
简单来说,你可以把它们想象成两种不同的“旋转”方式。自由涡就像浴缸排水时形成的漩涡,中心转得快,外面转得慢。强迫涡则像你搅拌一杯咖啡,整个液体像一块固体一样旋转,离中心越远转得越快。在实际工程中,比如水泵的入口或龙卷风,就是这两种模式的混合体。你可以在模拟器里选择“自由涡”或“强迫涡”类型,马上就能看到它们的速度分布完全不同。
🙋
诶,真的吗?那“环量Γ”这个参数是干嘛的?单位是m²/s,感觉好奇怪。
🎓
你可以把环量Γ理解成这个漩涡的“旋转强度”。它越大,意味着漩涡转得越“猛”。在自由涡里,它直接决定了速度大小:$v_\theta = \frac{\Gamma}{2\pi r}$。工程现场常见的是,飞机机翼产生的升力,就和它周围空气的环量大小直接相关。试着在模拟器里拖动“环量Γ”的滑块,你会看到整个速度曲线和压力曲线都跟着升高或降低,中心压力变化特别明显。
🙋
压力变化?对了,下面那个“汽蚀风险”的Warning是怎么来的?压力低了就会汽蚀吗?
🎓
没错!汽蚀(也叫空化)就是当液体局部压力太低,低到它的饱和蒸气压时,液体就会“沸腾”产生小气泡,这些气泡破裂时会损坏设备。比如在船舶螺旋桨的叶尖,高速旋转容易形成强涡,中心压力极低,就会发生汽蚀。模拟器会根据你设置的流体密度和参考压力,实时计算压力分布,并判断是否有区域的压力低于水的蒸气压(约2338 Pa)。你试着把环量Γ调得非常大,或者把参考压力p∞调小,很可能就会触发汽蚀Warning了!

物理模型与关键公式

自由涡(势涡)模型:这是一种无旋流动,其切向速度与半径成反比,中心速度理论上为无穷大(实际中受粘性等因素限制)。压力随着向中心靠近而急剧下降。

$$v_\theta = \frac{\Gamma}{2\pi r},\quad p = p_\infty - \frac{\rho\Gamma^2}{8\pi^2 r^2}$$

变量含义:$v_\theta$是切向速度(m/s),$r$是到涡心的半径(m),$\Gamma$是环量(m²/s),$p$是当地压力(Pa),$p_\infty$是无穷远处的参考压力(Pa),$\rho$是流体密度(kg/m³)。

强迫涡(刚体旋转)模型:流体像刚体一样旋转,角速度恒定。切向速度与半径成正比,压力分布为抛物线型。

$$v_\theta = \omega r,\quad p = p_\infty - \frac{1}{2}\rho\omega^2 r^2$$

变量含义:$\omega$是旋转角速度(rad/s)。其他变量同上。兰金涡是这两种模型的结合:在涡核半径$r_c$内为强迫涡,在$r_c$外为自由涡。

现实世界中的应用

航空航天工程:飞机机翼尾缘脱落的翼尖涡,其外围部分可近似为自由涡。强大的翼尖涡会对后续飞机造成影响,这也是机场起降需要间隔时间的原因之一。

流体机械(泵与涡轮):离心泵或水轮机的叶轮内部流动包含强烈的强迫涡特征。设计不良会导致叶轮中心或叶片背面压力过低,引发汽蚀,严重损坏设备。

气象学与自然灾害:龙卷风是典型的兰金涡实例。其破坏性的风场结构——内部核心高速旋转(近似强迫涡),外部风速衰减(近似自由涡)——可以用本模拟的模型进行理想化分析。

船舶工程:船舶螺旋桨在运行时,桨叶尖端会泄出螺旋形的涡旋(梢涡),这些涡旋中心压力极低,是螺旋桨发生空化腐蚀和产生噪音的主要源头。

常见误解与注意事项

开始使用本模拟器时,有几个容易误解的地方需要注意。首先,“自由涡是无粘性的理想模型”。虽然常以水槽漩涡为例,但实际中壁面摩擦会导致明显不同的流动行为。请理解模拟器中优美的曲线是“用于理解本质的理想形态”。其次是参数设置。“环量Γ”和“角速度ω”看似无关,但在兰金涡模型中通过核心半径rc相互关联。例如,若设置ω=5 rad/s、rc=0.1m,则核心边界处的速度为v_θ=ω*rc=0.5 m/s,该值将自动决定外部自由涡的Γ值(Γ=2π*rc*v_θ)。若忽略这种连续性,可能产生非物理的速度突变,需特别注意。最后关于压力。切勿直接采信计算结果显示中心压力趋于负无穷大的情况。现实中压力降至蒸汽压时会发生空化从而改变流动状态,且中心点的数学奇异性本身就是模型的局限性。在实际工程中,定量分析“低压区域的范围与程度”更为重要。例如在泵设计中,“基于所需扬程预留多少空化余量(NPSH)”正是这类分析的延伸。

使用指南

  1. 在环量输入框(valGammaNum)设定自由涡强度,单位m²/s,典型范围0.5~5.0;对应物理模型中Γ=2πrvθ
  2. 调节涡核半径参数(valRcNum),单位mm,控制兰金涡流的势流与刚体旋转过渡区域,常用值5~50mm
  3. 设定角速度(valOmegaNum),单位rad/s,驱动强迫涡中心的旋转速度,范围10~500 rad/s
  4. 点击模拟按钮实时生成速度矢量场、径向压力分布曲线及空化数Cv=2(p-pv)/(ρv²)等参数输出

具体计算示例

以水泵叶轮出口涡流为例:设Γ=2.5m²/s,Rc=15mm,ω=300rad/s。在距涡核r=50mm处,自由涡速度vθ=Γ/(2πr)=2.5/(2π×0.05)≈7.96m/s;压力降Δp≈0.5ρvθ²=0.5×1000×7.96²≈31.7kPa;若环境压力为101kPa、汽化压力为2.3kPa(20°C水),空化数Cv≈2×(101-2.3)/(1000×7.96²)≈3.08,>1.0时无空化风险

实务注意事项

  1. 涡核半径过小(<5mm)易导致数值计算奇点;过大(>80mm)则强迫涡区域占比过大,失去兰金涡特征
  2. 高速离心泵(ω>400rad/s)应重点监测压力谷值,当Cv<0.5时需提高进口压力或降低转速
  3. 在CFD前处理中,用本模拟器计算的初始涡场作为边界条件输入ANSYS-Fluent,可加快收敛速度15~20%
  4. 强迫涡中角速度ω与环量Γ非独立:ω=Γ/(2πRc²)的约束关系决定了涡核内刚体旋转的自洽性