参数设置
流动控制
拖拽提示：在画布上拖拽圆柱改变位置。播放中涡动画运行，实时时间显示。可保存多个流场快照，最大5个，可在后续对比中查看流场的时间演变。
流场基本方程
有旋转的势流（Kutta-Joukowski）：
$$\psi = U r\!\left(1-\frac{R^2}{r^2}\right)\!\sin\theta - \frac{\Gamma}{2\pi}\ln r$$
Kutta-Joukowski升力定理（单位长度升力）：
$$L = \rho\, U\, \Gamma$$
卡尔曼涡频率公式：
$$f_{vs}= St \cdot \frac{U}{D}, \quad St \approx 0.2 \;(Re = 10^3 \text{～}10^5)$$
Newton-Weise阻力方程：
$$F_D = \frac{1}{2}\,C_D\,\rho\,U^2\,D$$
圆柱周围流动、升力、阻力计算是什么？
👤
在这个模拟器中，当我改变旋转参数时，流场会怎样变化呢？
🔬
大概来说，流体在圆柱周围有两种力在作用：一是"回旋流"，二是"向下拉力"。旋转参数，也就是圆柱的"旋转速度"来表示有多强，改变它的滑块从0开始逐渐变大，你会看到流体会在圆柱的上下不对称流动，中间有力流过去。这就是马格努斯效应，也就是曲球效应。这个模拟器会自动计算出相关的升力、阻力和卡尔曼涡的特征周期。
👤
嗯，这样的话，那「流场」是怎么进行计算的呢？仅仅看到流动的样子就能明白原理吗？
🔬
好问题！这里有几个基本原理。首先，流体流动用"势流+涡"来表现。这个模拟器是，目前设定的圆柱和旋转参数根据势流方程组求解计算，利用"Kutta-Joukowski升力定理"来估计升力。旋转的强度越大，上方的流速越快、下方越慢，这样就产生压力差，向上的升力就越大。随着旋转增强，流体被圆柱推向某一侧，流动模式就反向改变。这就叫"马格努斯力"。同时，该模拟器还根据 Reynolds 数计算阻力系数，用 $F_D = \frac{1}{2}\rho U^2 D C_D$ 式计算阻力。参数改变后，这些数值实时改变，你能看出流动如何对应这些力学参数。
常见问题
问：改变圆柱直径和流速时，什么会发生变化？
答：改变直径会改变升力、阻力和流场绕柱体的周期。改变流速会改变马格努斯效应强弱。流体旋转参数越大，绕柱周围流场绕过圆柱的频率越高，升力和阻力也会增加。你可以实时看到变化。
问：卡尔曼涡频率是如何计算的？
答：斯特劳哈尔数恒定（约 St≈0.2），基于流速和圆柱直径计算频率：f = St·U/d。在本模拟器中使用势流解析的基础上，涡频率保持恒定，这是一个经验值，不会像完整 CFD 那样动态变化。
问：这个模拟器能再现真实的流动现象吗？
答：虽然本工具基于经典势流模型和经验公式，但由于势流不包括黏性和边界层效应，实际流动中的细节差异会有。该工具用于理解势流基本原理和快速参数估算，实际应用需要 CFD 仿真验证。
问：如果我想用马格努斯效应来说明具体的工程案例，你能给个建议吗？
答：经典的棒球曲球，还有从事自旋体旋转（如网球发球旋转）都是例子。还有旋转气缸用于防风（如 Flettner 转子）或降低拖力的技术，本模拟器中通过调节旋转参数，可以观察升力与旋转速度的关系，为实际工程设计提供定性认识。
实务应用中的推广
通过本模拟器，用户能快速了解到流体力学的核心概念。在掌握本模拟器中看到的流场形态、升力与阻力的时间变动后，建议用户接触更高级的 CFD 工具，比如 ANSYS、OpenFOAM 等，去预测真实流动过程，比对网格依赖性，进行更精密的设计优化。同时，势流理论远非全部真实，应当指出。本工具的流体流场，无黏粘性，这是当然的，仅当流体流动的粘性以及流体与物体表面相互作用无视不计时，本方法才给出接近真实的答案。对于圆柱直径较小、流速较低或极高的雷诺兹数区间中，黏性效应或涡分离会改变流场样貌。这个模拟显示的是，某一参数域的简化表示，仅供初学者学习或定性估算。
其次，参数设置的现实性。假设，直径d=0.05m、流速u=2.0m/s、流体密度ρ=1.225kg/m³的时候：雷诺兹数 Re≈6667、斯特劳哈尔数 St≈0.2，可得涡频率 f=8Hz，阻力系数 CD≈1.2，单位长度阻力 FD=0.37N/m。实际配管束中，由于相邻管之间有干涉，CD 值会增加到 1.5～2.0，本模拟实现的单管孤立配置不能直接用于成束管解析。
最后，阻力计算的注意事项。这里使用的阻力系数 CD 依赖雷诺兹数。即使在同一流体、同一流速条件下，CD 仍随 Re 而变。当 Re>1000 时涡放出现象显著，此时 St≈0.2 是圆柱标准值，但实装置中因表面粗糙度、端部效应可能偏离 ±15%。FD 是单位长（[N/m]）阻力，计算实管总阻力时需乘以实际管长，并结合支持节点间的共振频率来判断是否需要加强或减振。当流动激励的振幅增长到直径的 0.5 倍以上时，会进入自激振动（流动诱发振动，FIV）甚至 galloping 不稳定，此时升力 L 时间变动剧烈，需要考虑流速限制或减衰措施。本工具是理想化、定性的参考，工程应用需谨慎。
具体计算例子
圆柱直径 d=0.05 m，空气物性取本工具数值 ρ=1.204 kg/m³、ν=1.511×10⁻⁵ m²/s，来流速度 U=5 m/s，循环量 Γ=0。雷诺数 Re = U·D/ν = (5×0.05)/(1.511×10⁻⁵) ≈ 16,500。St≈0.20 时，涡脱落频率 f=(0.20×5)/0.05≈20 Hz。该 Re 范围内 C_D≈1.2，因此 F_D=0.5×C_D×ρ×U²×D≈0.90 N/m。若加入 Γ=2.0 m²/s，则升力 L=ρUΓ=1.204×5×2.0≈12.0 N/m。