圆柱周围的流动 — CAE术语解说

有几个原因。①形状简单，具有丰富的解析解和实验数据——Williamson（1996年）和Papaioannou（2006年）等权威实验·数值数据已公开。②根据雷诺数Re，覆盖Stokes流动（Re<1）→定常分离（Re=10〜40）→卡门涡街（Re=40〜200）→乱流（Re>1000）等多样的流动体制。③卡门涡街的Strouhal数（St = f*D/U）可与实验值St≈0.2进行比较，用于验证求解器是否能够正确捕捉涡的发生频率。这是查看CFD动态精度的最佳问题。

圆柱后方周期性地交替出现涡的现象，1911年由Theodore von Karman进行理论分析而得名。当Re ≈ 40超过时，圆柱后方形成的对称分离泡变得不稳定，上下交替地开始释放涡。这对圆柱产生周期性的升力变动——当这种振动与结构物的固有振动频率一致时，就会发生共振「涡激振动（VIV: Vortex Induced Vibration）」。这是桥梁和烟囱在风中发出声音的原因，1940年的塔科马海峡大桥崩塌也是这一类型。

要准确捕捉Re=100的卡门涡街，需要：①对于2D分析，用结构化六面体网格使圆柱表面的y+ < 1（即使是层流分析，边界层分辨率也很重要），②保持时间步长Co < 0.5左右（CFL条件），③计算时间足够长（直到稳定的周期振动建立）——这些是基本要点。验证指标是阻力系数Cd、升力系数Cl的振幅，以及Strouhal数St。St = 0.166〜0.195（Re相关）与实验值的一致性在±5%以内即可认为合格。

当Re > 200左右时，3D效应开始出现。在2D分析中，涡结构在z方向无限延伸，但在实际3D流动中，z方向也会发展不稳定性，涡结构会扭转·分裂——这被称为A/B模式。Re = 190〜260时出现A模式（波长≈3〜4D），Re = 250〜300时出现B模式（波长≈0.8〜1D）。3D LES分析可以精确捕捉这些，但计算成本大幅增加。从工程目的考虑，如果只是「想求涡释放周频率」，2D非定常RANS通常已足够。

烟囱·桥梁·海洋导管的涡激振动设计是典型应用。用圆柱模型计算涡释放频率 f = St * U / D，与结构物的固有振动频率进行比较。如果一致就会发生VIV，因此采取以下措施——①改变形状（把截面改为椭圆或D型，削弱涡释放），②向结构物添加制振阻尼器，③安装螺旋形翅片（Strakes）来破坏涡的规律性——。石油平台的海底导管在潮流作用下的VIV是疲劳设计的决定性因素，实际工程中进行CFD与FEM疲劳的耦合分析。