3. 二维圆柱绕流的涡街模拟与斯特劳哈尔数提取
问题定义与几何建模
二维圆柱绕流是流体力学中的经典问题。当雷诺数处于约 47 至 200 之间时，圆柱尾部会出现交替脱落的涡旋，形成卡门涡街。本例取圆柱直径 D = 0.02 m，来流速度 U∞ = 0.05 m/s，空气运动粘度 ν = 1.5×10⁻⁵ m²/s，则雷诺数 Re = U∞·D/ν ≈ 66.7，处于稳定的层流涡街区间。
计算域采用矩形区域：上游距圆心 10D，下游距圆心 25D，上下各距圆心 10D。圆柱表面设置无滑移壁面，上下边界设为对称（自由滑移），出口为压力出口，入口为速度入口。网格采用结构化四边形网格，圆柱附近加密，第一层网格高度约 0.0002 m，周向布置 120 个节点。
求解器设置与关键参数
在 Fluent 中启动 2D 双精度求解器，选择基于压力的瞬态求解器。粘性模型选择 Laminar（层流），因为 Re ≈ 67 远低于湍流临界值。压力-速度耦合采用 SIMPLE 算法，压力离散用 Second Order，动量离散用 Second Order Upwind，瞬态格式用 Bounded Second Order Implicit。
时间步长取 Δt = 0.01 s，对应每个涡脱落周期约 20 个时间步。总模拟时间设为 20 s，确保涡街充分发展并达到统计稳态。库朗数控制在 1 以下，保证时间推进的稳定性。
升力系数监测与频率提取
在圆柱表面创建升力系数监测（Report → Forces → Lift Coefficient），方向设为 Y 轴。同时监测阻力系数用于辅助判断。升力系数 Cl 随时间呈正弦振荡，其振荡频率即为涡脱落频率 f。
对升力系数时程曲线进行快速傅里叶变换（FFT）或通过峰值间隔法提取频率。峰值间隔法：记录相邻两个升力系数峰值对应的时间 t₁ 和 t₂，则 f = 1/(t₂ - t₁)。为提高精度，取多个周期的平均值。
斯特劳哈尔数计算与对比
斯特劳哈尔数定义为 St = f·D/U∞。对于圆柱绕流，在 Re = 50~200 范围内，实验关联式给出 St ≈ 0.198（Roshko 1954）或 St = 0.212 - 4.5/Re（Fey 等）。本例 Re = 66.7，代入 Fey 关联式得 St ≈ 0.212 - 4.5/66.7 ≈ 0.145。
假设模拟提取的涡脱落频率 f ≈ 0.36 Hz，则 St = 0.36 × 0.02 / 0.05 = 0.144，与关联式预测值 0.145 吻合良好，相对误差约 0.7%。
关键验证指标：
- 升力系数振荡幅值 Cl_amp ≈ ±0.3~0.5（Re=67 时）
- 阻力系数平均值 Cd_mean ≈ 1.4~1.5
- 涡脱落频率 f 与 St 关联式偏差 < 5%
- 升力系数与阻力系数相位差约 90°
常见问题与排障
| 现象 | 可能原因 | 处理方式 |
|---|---|---|
| 升力系数不振荡 | 时间步长过大或网格过粗 | 减小 Δt 至 0.005 s，加密尾流区网格 |
| 涡街不对称 | 计算域上下不对称或收敛不足 | 检查边界条件对称性，延长模拟时间 |
| 频率提取偏差大 | 采样周期不足或FFT分辨率低 | 增加模拟总时长至 30 s 以上 |
| 残差震荡不收敛 | 库朗数过大 | 减小时间步长或改用 Coupled 求解器 |
建议在圆柱后方 1.5D 处设置一个监测点，记录该点速度的时程曲线，其振荡频率应与升力系数频率一致，可交叉验证涡脱落频率的提取结果。