› 势能流模拟器 返回
流体分析

势能流模拟器

二维势能流圆柱绕流可视化,含流线与压力系数彩色图。体验Magnus效应(旋转圆柱产生的升力)及d'Alembert悖论。

流动模式
参数
来流速度 U∞1.0 m/s
环量 Γ0.0 m²/s
特征长度 a1.00 m
显示
实时数值
1.00
来流速度 U∞ [m/s]
0.0
环量 Γ [m²/s]
2
驻点数
2.00
瞬时最大流速 |V| [m/s]
0.00
升力系数 CL
0.0
升力 L=ρU∞Γ [N/m]
−3.00
Cp 上端
−3.00
Cp 下端
势流·粒子动画
高速·低压 (Cp<0)
低速·高压 (Cp→1)

示踪粒子随速度场平流。颜色=压力系数 Cp,黄色=驻点,绿色=升力矢量。用预设切换均匀流、源、圆柱、自由涡和兰金卵形。

理论与主要公式
复势的叠加(均匀流+偶极子+涡):
$w(z)=U_\infty z+\dfrac{\mu}{2\pi z}-\dfrac{i\Gamma}{2\pi}\ln z,\quad u-iv=\dfrac{dw}{dz}$

圆柱半径(偶极子强度 μ):$a=\sqrt{\dfrac{\mu}{2\pi U_\infty}}$

流函数(圆柱+环量):$\Psi = U_\infty r\!\left(1-\frac{a^2}{r^2}\right)\!\sin\theta - \frac{\Gamma}{2\pi}\ln r$

升力(Kutta–Joukowski):$L = \rho\, U_\infty\, \Gamma$

压力系数:$C_p = 1 - \left(\frac{V}{U_\infty}\right)^2$

什么是势流模拟器

🙋
在这个模拟器中增大“环量 Γ”,圆柱周围的流动就开始旋转起来。这有什么用?
🎓
大致来说,这就是产生“升力”的机理。例如旋转的棒球会拐弯的“马格努斯效应”,就是这种环量造成的压力差引起的。用上方的滑块增大 Γ,可以看到圆柱上方流速升高(蓝色部分)、压力下降,于是产生向上的力。
🙋
原来如此!不过我听说如果流动完全对称,流过圆柱的流体前后相互抵消,阻力为零。这就是“达朗贝尔佯谬”吧?
🎓
没错!在理想势流中,压力分布前后对称,所以阻力为零。真实流体有黏性,所以不会这样。在模拟器中把“环量 Γ”调回 0,看看圆柱表面的颜色(压力系数 $C_p$),前后几乎是同样的颜色(压力)吧?这就是佯谬的可视化。
🙋
明白了!那“Q”滑块是做什么用的?调成负值好像会被吸进去。
🎓
它表示“源”或“汇”。用“源”预设把 Q 设为正值,从一点向四周流出的源与均匀流合成,得到绕前端圆钝的半无限物体(兰金半体)的流动。驻点出现在源的上游 x=−Q/(2πU∞) 处,Q=8 m²/s、U∞=1 m/s 时约在前方 1.27 m。Q 为负则成为汇。实务中,先用这种简单模型来理解发动机进气口或排风扇吸入口附近的流动。

常见问题

指在理想流体(无黏、不可压)的势流中,均匀流中的物体不受阻力这一理论结果。在本模拟器的“圆柱(偶极子)”中把环量 Γ 设为 0,圆柱前后的压力系数分布对称(前后驻点 Cp=1,上下端 Cp=−3),可以确认压力阻力为零。
把环量强度 Γ(gamma)设为正值或负值。增大 Γ,圆柱周围流线的不对称性增强,产生升力。同时可以在压力系数色图中看到圆柱上下产生压力差的样子。
Q 用于“源”和“兰金卵形”预设。Q 为正时是流体从一点涌出的“源”,为负时是流体被吸入的“汇”。均匀流与源叠加得到绕兰金半体的流动,源与汇组合得到绕兰金卵形(封闭物体)的流动,可以确认物体大小和驻点位置由 Q 和 U∞ 决定。
圆柱内部按速度为零处理,流线、示踪粒子和压力色只画在圆柱外侧。复势公式本身在圆柱内侧也有值,但没有物理意义,因此不显示。

现实世界中的应用

航空工程(翼型理论):这是薄翼理论的基础思想。通过加入环量来解释作用在机翼上的升力,库塔–茹科夫斯基定理是机翼设计的根本。在模拟器中改变 Γ 观察升力的变化,就是对它的直接体验。

运动球类的流体分析:足球的香蕉球、棒球的曲线球、乒乓球的上旋,都是由“马格努斯效应”即旋转产生环量引起的。势流模型是以理想化方式理解这一效应的第一步。

风力发电·涡轮设计的初步研究:先把叶片周围复杂的流动建模为无黏势流,用来大致把握压力分布和流线。环量的概念对理解涡流动尤为重要。

建筑·环境风工程的基础研究:在大致把握建筑物和结构物周围的风流时,有时会使用势流简化模型。源·汇模型有助于想象通风口附近的流动。

常见误解与注意事项

刚开始使用本模拟器时,有几个初学者特别容易掉进去的陷阱。第一,别忘了“势流是‘完美’流体的模型”。由于忽略了黏性,它与真实流动有根本区别。例如有人看到加入环量 Γ 产生升力,就草率地认为“这样就能解释机翼的全部升力了!”,但在真实机翼上,黏性产生的“环量”与“后缘条件”错综复杂地交织在一起。请把本工具看作提取了核心物理概念的“理想化第一步”。

第二,注意环量大小与驻点位置。环量较小时,圆柱的驻点沿表面向下(或向上)移动;Γ 超过 4πU∞a 后,驻点离开表面,流场中只出现 1 个驻点,并形成绕圆柱的封闭流线。U∞=1 m/s、a=1 m 时分界为 Γ≈12.6 m²/s,Γ=15 m²/s 时驻点移到离中心约 1.8 m 处。这是正确的势流解,但真实的旋转圆柱要获得这么大的环量需要非常快的转速。先从与 U∞a 同量级的 Γ(例如 U∞=5 m/s、a=1 m 时 5 m²/s 左右)开始试,更容易看清表面驻点移动与升力的关系。

第三,要理解“点源·汇(Q)”是数学上的抽象工具。实务中说“1 m²/s 的汇”很难有直观感受吧?可以把它看作用单个“点”来代表例如从直径 10 cm 的圆形风管吸入一定量空气的效果。改变 Q 会使流线图案剧烈变化,这是“奇点”所具有的数学强度的体现,诀窍是把它解释为有实际尺寸的吸入口的“一阶近似”。

使用指南

  1. 在 0.2~5 m/s 范围内设置来流速度 U∞
  2. 在 0.4~2.0 m 调整特征长度 a(圆柱为半径,兰金卵形为源与汇间距的一半)
  3. 在“圆柱(偶极子)”中把环量 Γ 在 −18~+18 m²/s 范围内变化,观察马格努斯效应引起的升力变化
  4. 切换流动模式(均匀流/源/圆柱(偶极子)/自由涡/兰金卵形),比较流线和 Cp 色图的变化
  5. 读取上端·下端压力系数、升力系数 CL、升力 L 和驻点数(由于达朗贝尔佯谬阻力总为 0,因此不显示)

具体计算示例

半径 R=1.0 m(直径 D=2.0 m)的圆柱受流速 U=5 m/s(大气中 ρ=1.225 kg/m³)的气流作用时:雷诺数 Re=ρUD/μ≈6.8×10⁵。无环量(Γ=0)时按理想流体的达朗贝尔佯谬 CD=0,但加入 Γ=+5 m²/s 后,按库塔–茹科夫斯基定理产生 L'=ρUΓ=1.225·5·5≈30.6 N/m 的升力,CL=Γ/(U·R)=1.0。同时 Cp(上端)≈−3.65、Cp(下端)≈−2.38,可以确认流线变得不对称。

实务注意事项