圆柱绕流问题在早期的研究过程中局限性大，研究成果也都比较简单。1985年，伊杰里奇克[1]通过实验得出在均匀来流条件下，圆柱体的阻力系数Cd大小与来流雷诺数(Re) 有很大关系，并给出了相应的关系曲线。中国科学院力学研究所凌国灿[2]在1985年对圆柱绕流旋涡运动及离散涡方法进行了研究。1992年，清华大学苑明顺[3]用大涡模拟方法对亚临界条件下二维圆柱绕流问题进行了数值模拟，在二维模型中引入了涡量耗散项来模拟涡街的三维涡量耗散效应，对比分析了有关流动特性，结果表明，采用大涡模拟能够得到比较合理的主要流动参数。王亚玲等[4]利用计算流体力学软件CFX-4，对黏性不可压缩流体的圆柱绕流在定常情况下进行了三维数值模拟。国内外****通过大量的实验和数值方法研究了圆柱绕流的气动力特性及变化规律[5-7]，但是关于高雷诺数下圆柱绕流近壁区的流动特征及流动机理等方面的分析还不够深入，对圆柱绕流过程中边界层转捩的影响更是很少关注。
在圆柱绕流过程中，转捩是其中一个重要的影响因素。转捩现象是流体力学中的经典难题，是流体从层流状态到湍流状态的过渡过程。而且在实际的工程应用中，转捩现象是很多飞行器外绕流的一个重要过程，如果能够准确地预测转捩点的位置并且控制转捩的发生，就能够大幅度提高飞行器的空气动力性能。同样对于圆柱绕流问题，如果能够很好地预测并控制转捩过程，模拟出圆柱尾迹的复杂流动状况，对圆柱绕流问题的研究具有非常重要的意义。
基于工程需要，国内外研究人员根据实验结果及经验提出了一些模拟转捩过程的模型，随时间的发展主要是早期提出的低雷诺数湍流模型及其修正形式，其中以Schmidt和Pantankar[8]的工作较有代表性；之后提出的考虑间歇性的转捩模型，其中以Vancoillie和Dick[9]以及Steelant和Dick[10]的工作较有代表性；以及目前提出的基于局部变量的转捩模型，其中有代表性的是Walters和Leylek[11]的工作，以及Langtry和Menter等[12-14]提出的完全基于当地变量的γ-
Transition SST模型的优点为：①避免了一般情况下对平均场进行积分的过程，计算周期短，计算要求较低；②对自由来流湍流度、分离和压力梯度等影响转捩的因素敏感；③与SST k-ω湍流模型相容，便于与现有的CFD程序结合。
本文采用Transition SST模型对亚临界区(Re=1.4×105)、临界区(Re=5.0×105和Re=8.5×105)、超临界区(Re=1×106) 和过临界区(Re=3×106) 的圆柱绕流问题进行了数值模拟，实现了对圆柱绕流的转捩预测。在亚临界区(Re=1.4×105) 和临界区(Re=8.5×105) 的条件下，将其结果与实验结果及采用SST k-ω两方程湍流模型得到的结果进行了对比分析。结果表明，通过Transition SST模型较为准确地预测出了圆柱绕流边界层的转捩过程和圆柱绕流近壁区的流动特征，验证了Transition SST模型在模拟圆柱绕流方面的优越性及该模型预测转捩位置的能力。在此基础上，应用Transition SST模型通过流场结构及摩擦力系数等曲线分析了不同雷诺数下圆柱绕流边界层转捩位置和圆柱绕流近壁区流动特征的变化规律。
1 数值计算方法 1.1 基本控制方程 对于不可压缩黏性流体，其流动控制方程可表示为
 | (1) |
 | (2) |
式中：ui为流体沿xi方向的速度分量，在二维情况下i=1，2，u1=u和u2=v分别为水平和垂直方向的速度分量；t为时间；p为压力；Re=ρ∞U∞D/μ为雷诺数，ρ∞为来流密度，U∞为来流速度，D为圆柱直径，μ为黏性系数。
1.2 湍流模型 本文采用的SST k-ω湍流模型中关于湍动能k以及比耗散率ω的输运方程[15]为
 | (3) |
和
 | (4) |
式中：Gk为湍动能的速度梯度；Gω为比耗散率的速度梯度；Γk和Γω分别为k和ω的有效扩散系数；Yk和Yω分别为关于k和ω的湍流耗散项；Dω为交叉扩散项；ρ为来流密度。
SST k-ω湍流模型定义的有效扩散系数[15]Γk和Γω分别为
 | (5) |
 | (6) |
式中：σk和σω分别为关于k和ω的湍流普朗特数；μt为湍流黏性系数[15]，可以通过下式得到：
 | (7) |
其中：S为应变率；a1为模型参数，a1=0.31。
σk和σω[15]分别被定义为
 | (8) |
 | (9) |
式中：F1和F2为混合函数。
湍流黏性系数μt中a*的定义[15]如下：
 | (10) |
式中：Ret=ρk/μω；Rk=6；a0*=βi/3，βi=0.072；一般对于高雷诺数的SST k-ω两方程湍流模型，a*=a∞*=1。
1.3 转捩模型 将关于间歇因子γ和当地边界层动量厚度雷诺数
 | (11) |
式中：转捩源项[12]定义如下：
 | (12) |
 | (13) |
 | (14) |
 | (15) |
其中：Flength为控制转捩区长度的参数；Ω为旋涡强度；Fonset为涡量雷诺数ReV的函数[12]：
 | (16) |
 | (17) |
 | (18) |
 | (19) |
 | (20) |
 | (21) |
 | (22) |
式中：y为沿外法线方向与壁面之间的距离；Fonset用来启动Pγ，ReV/(2.193Reθc)为当地转捩的判据，当该比值超过1时，Pγ启动，即γ开始增长，Reθc为临界雷诺数，转捩雷诺数
基于动量厚度雷诺数
 | (23) |
式中：
 | (24) |
 | (25) |
 | (26) |
 | (27) |
 | (28) |
 | (29) |
 | (30) |
 | (31) |
式中：θBL为边界层动量厚度；Reθt为当地转捩雷诺数；Fθt为开关函数，从边界层内到边界层外由1逐渐变为0；Fwake的作用是确保Fθt在圆柱下游的尾流区域中为0；cθt和σθt分别为常数[12]：cθt=0.03，σθt=2.0。
2 数值模拟条件及网格 2.1 数值计算条件及物理模型 本文研究的是来流速度U∞沿x轴正方向流向直径为D的二维圆柱绕流问题，如图 1所示。图中：θ为圆柱表面的方位角。分别取来流雷诺数Re=1.4×105(亚临界区)、Re=5.0×105、8.5×105(临界区)、Re=1×106(超临界区) 和Re=3×106(过临界区)，根据实验条件选择来流的湍流强度为0.8%。
 |
| 图 1 流动示意图 Fig. 1 Schematic of flow |
| 图选项 |
本文时均阻力系数Cd、斯特劳哈尔数St、时均摩擦力系数Cf分别由下列各式确定：
 | (32) |
 | (33) |
 | (34) |
式中：Fd为圆柱所受到的时均阻力；f为涡脱落的频率；τ0为壁面剪切应力。
2.2 计算网格及边界条件 计算域取40D×55D的矩形域，其中计算域左边界和计算域上下边界与圆柱中心的距离均为20D，计算域右边界与圆柱中心的距离为35D。物面法向第1层网格高度由y+~1确定。图 2和图 3分别为整个计算域网格和圆柱表面附近网格的示意图。
 |
| 图 2 整个计算域网格示意图 Fig. 2 Schematic of grids of entire computational domain |
| 图选项 |
 |
| 图 3 圆柱表面附近网格示意图 Fig. 3 Schematic of grids of near cylindrical surface |
| 图选项 |
计算域的左边界采用速度入口边界条件，在此边界上沿x、y方向速度分量分别为u=U∞，v=0；圆柱表面采用无滑移壁面边界条件，即u=v=0；计算域的上下边界采用对称边界条件；计算域的右边界采用流动出口边界条件。
3 数值模拟结果 3.1 亚临界区数值模拟结果分析 表 1给出了亚临界区(Re=1.4×105) 时数值模拟计算得到的二维圆柱时均阻力系数Cd和斯特劳哈尔数St的值，其中也给出了Cantwell[16]和Schewe[17]通过实验测得的结果及苑明顺[3]通过数值模拟计算得到的结果。通过表 1中Cd和St的对比可以看出，在亚临界区采用Transition SST模型和SST k-ω模型都能够很好地模拟圆柱绕流问题，但与采用SST k-ω模型模拟所得到的结果相比，采用Transition SST模型模拟得到的结果与实验结果符合得更好。
表 1 在亚临界区时二维圆柱时均阻力系数Cd和St数值模拟计算结果 Table 1 Numerical predicted results of time average drag coefficient Cd and St of 2D cylinder in subcritical region
| 模型 | Cd | St | 备注 |
| SST k-ω模型 | 1.05 | 0.24 | 数值模拟结果 |
| Transition SST模型 | 1.20 | 0.21 | 数值模拟结果 |
| Cantwell[16] | 1.23 | 0.18 | 实验结果 |
| Schewe[17] | 1.18 | 0.20 | 实验结果 |
| 苑明顺[3] | 1.19 | 0.21 | 数值模拟结果 |
表选项
图 4给出了亚临界区(Re=1.4×105) 时通过数值模拟得到的圆柱表面时均压力系数Cp分布曲线，同时也给出了Cantwell[16]实验测得的结果。
 |
| 图 4 在亚临界区圆柱表面时均压力系数Cp分布曲线 Fig. 4 Distribution curves of time-averaged pressure coeffient Cp on cylindrical surface in subcritical region |
| 图选项 |
图 5给出了亚临界区(Re=1.0×105、1.4×105和2.6×105) 时通过数值模拟得到的圆柱表面时均摩擦力系数Cf分布曲线，还给出了Krish-nan等[18]通过DDES方法得到的结果以及Achenbach[19]通过实验测得的结果Exp.1(Re=1.0×105) 和Exp.2(Re=2.6×105)。
 |
| 图 5 在亚临界区圆柱表面时均摩擦力系数Cf分布曲线 Fig. 5 Distribution curves of time-averaged friction coefficent Cf on cylindrical surface in subcritical region |
| 图选项 |
由图 4和图 5可知，由于Transition SST模型对分离和压力梯度等因素比SST k-ω模型更为敏感，所以采用Transition SST模型更为准确地模拟出了圆柱绕流背风面压力和摩擦力的变化情况，使得采用Transition SST模型模拟得到的压力和摩擦力系数与实验结果更为相符。
从图 5可以进一步看出，采用Transition SST模型模拟得到的结果(Re=1.4×105) 在圆柱上表面θ≈0°~50°范围内Cf随θ增大逐渐增大，并在θ≈50°处达到最大值；在θ≈50°~88°范围内Cf随θ增大而逐渐减小，并在圆柱上表面θ≈88°处Cf的正负符号发生变化，这表明流动在这个位置发生了分离，即分离点的位置为θ≈88°。同样从图 5可以看出，采用SST k-ω模型和Krishnan等[18]通过DDES方法得到分离点的位置分别为θ≈85°和θ≈87°，采用Transition SST模型模拟得到的结果与Krishnan等[18]通过DDES方法得到的结果更加相符。通过3种模型模拟得到的分离点位置介于Re=1.0×105和Re=2.6×105时通过实验方法得到的分离点位置θ≈78°和θ≈94°之间，符合在亚临界区分离点位置随雷诺数增大逐渐后移的变化规律。而在圆柱上表面θ≈88°~120°范围内采用Transition SST模型模拟得到的结果与SST k-ω模型及DDES方法得到不同变化的原因，可能是因为Re=1.4×105时处在亚临界区向临界区的过渡阶段，流动分离后有一个再附的趋势，即有一个形成分离泡的趋势。
由于Transition SST模型对分离和压力梯度等因素的敏感性，采用Transition SST模型可以较为准确地模拟出亚临界区(Re=1.4×105) 时圆柱绕流尾迹区在一个周期内的流动特征，如图 6所示。对于圆柱绕流问题来说，其尾迹区域的近壁区流动特征主要是分离与局部二次分离产生的不同尺度旋涡的合并、成对、分叉等强烈的相互作用，旋涡脱落的过程实际上是多个涡间以及与剪切层相互作用的结果[17]。在大雷诺数下，只要旋涡强度足够大或者距离壁面足够近，二维或三维旋涡均可在近壁区诱导产生局部分离，形成新的旋涡[20-21]。
 |
| 图 6 在亚临界区时采用Transition SST模型模拟一个周期内的流线 Fig. 6 Streamline in a period based on transition SST model in subcritical region |
| 图选项 |
从图 6可以看出圆柱背风面尾迹区流动在一个周期内的变化过程。由于t1时刻尾迹区圆柱下表面附近的旋涡距离壁面足够近，在其诱导作用下t2时刻在尾迹区圆柱下表面附近发生了局部流动分离，产生了一个新的旋涡；尾迹区圆柱下表面附近的旋涡随时间沿流向向后发展逐渐变大，与尾迹区圆柱上表面附近的剪切层相互作用，在t3时刻尾迹区圆柱上表面附近产生了2个正向旋转的小涡，且此2个小涡与之前尾迹区圆柱上表面附近产生的反向旋转的小涡由于旋向以及来流的影响逐渐开始合并，在t6时刻完成合并，成为一个正向旋转的旋涡。同样由于涡间以及来流的影响，在t6时刻发展的正向旋涡和在t2时刻由诱导作用发生局部分离产生的旋涡开始逐渐合并，在t8时刻完成合并发展为一个大的分离涡，之后发生旋涡脱落。这就是尾迹区圆柱上表面附近旋涡产生及发展脱落的过程，同样可以看到尾迹区圆柱下表面附近旋涡产生及发展脱落的过程。
综上所述，由于Transition SST模型对流动分离和压力梯度等因素较为敏感，能够很好地模拟圆柱表面的压力、摩擦力变化及圆柱绕流尾迹近壁区的分离等流动特征，从而使得Transition SST模型更加适用于模拟亚临界区的圆柱绕流问题。
3.2 临界区数值模拟结果分析 表 2给出了临界区(Re=8.5×105) 时对二维圆柱绕流数值模拟得到的时均阻力系数Cd和斯特劳哈数St的值，表中还给出了Achenbach[19]和Jones等[22]的实验结果。可以看出，与采用SST k-ω模型模拟得到的结果相比，采用Transition SST模型模拟得到的Cd和St与实验结果符合得更好。
表 2 在临界区时二维圆柱时均阻力系数Cd和St数值模拟计算结果 Table 2 Numerical predicted results of time average drag coefficient Cd and St of 2D cylinder in critical region
| 模型 | Cd | St | 备注 |
| SST k-ω模型 | 0.59 | 0.41 | 数值模拟结果 |
| Transition SST模型 | 0.71 | 0.48 | 数值模拟结果 |
| Achenbach[19] | 0.67 | Exp.1 | |
| Jones等[22] | 0.46 | Exp.2 |
表选项
图 7为在临界区(Re=8.5×105) 时数值模拟计算所得圆柱表面时均压力系数Cp分布曲线，还分别给出了Achenbach[19] (Exp.1) 和Jones等[22] (Exp.2) 实验测得的结果。
 |
| 图 7 在临界区圆柱表面时均压力系数Cp分布曲线 Fig. 7 Distribution curves of time-averaged pressure coefficient Cp on cylindrical surface in critical region |
| 图选项 |
从图 7可以看出，相比于采用SST k-ω模型模拟得到的圆柱表面压力系数Cp，采用Transition SST模型模拟计算得到的圆柱表面压力系数Cp与Jones等[22]得到的实验结果(Exp.2) 更加相符。而采用SST k-ω模型和Transition SST模型模拟得到的圆柱表面压力系数Cp与Achenbach[19]得到的实验结果(Exp.1) 在圆柱迎风面都符合得很好，但是在圆柱背风面与实验结果差异都比较大。Achenbach[19]提到由于受到实验条件和方法的局限性，在实验过程中没能够很好地模拟圆柱背风面的压力分布，这是在背风面与2种模型模拟的结果和Jones等[22]实验得到的结果产生较大差异的原因。
图 8为临界区(Re=8.5×105) 时数值模拟计算所得圆柱表面时均摩擦力系数Cf分布曲线，还
 |
| 图 8 在临界区圆柱表面时均摩擦力系数Cf分布曲线 Fig. 8 Distribution curves of time-averaged friction coefficient Cf on cylindrical surface in critical region |
| 图选项 |
给出了Achenbach[19]实验测得的结果(Exp.1)。可以看出，相比于采用SST k-ω模型，采用Transition SST模型模拟得到的圆柱表面时均摩擦力系数Cf曲线与实验结果更加相符，而且模拟出了与实验结果一致的摩擦力系数Cf突然增大又回落减小的过程，也就是圆柱边界层的转捩过程。
同样，由于Transition SST模型对流动分离和压力梯度等因素的敏感性，使得图 7和图 8中采用Transition SST模型模拟得到的结果与实验结果更为相符。
而且从图 8可以看出，采用Transition SST模型模拟得到的圆柱表面时均摩擦力系数Cf在圆柱上表面θ≈99°处时符号由正变负，之后在θ≈103°处Cf的符号又由负变正，表明在θ≈99°~103°这段位置上边界层流动发生了分离和再附的过程，即产生了分离泡；并且时均摩擦力系数Cf在θ≈103°处变为正值后突然开始变大，之后又回落减小，这说明边界层内流动在这里发生了从层流到湍流的转捩过程；在θ≈120°处Cf的符号再次由正变负，说明流动在这里发生了湍流分离。而实验得到的Cf值在θ≈105°处突然变大，之后又回落减小，在θ≈147°处Cf的符号由正变负发生了流动分离。实验结果中并没有捕捉到边界层分离再附产生的分离泡现象，Achenbach[19]提出理论上应该会产生分离泡，没有产生的原因可能是因为实验条件等问题的影响。
图 8中采用SST k-ω模型与实验结果符合程度最差，且采用SST k-ω模型模拟得到的Cf在θ≈113°处由正变负，圆柱表面在此处流动发生了分离，但并没有模拟出边界层分离再附产生的分离泡现象和边界层转捩现象，不能很好地模拟二维圆柱绕流问题。
从图 8还可以看出，采用Transition SST模型模拟得到的圆柱边界层转捩区是在圆柱上表面θ≈103°~120°的位置，而实验得到的转捩区则是在圆柱上表面θ≈105°~147°的位置。采用Transition SST模型模拟得到的圆柱边界层转捩区长度明显要比实验得到的转捩区长度短。而且可以看出，采用Transition SST模型和SST k-ω模型模拟得到的圆柱边界层分离点位置分别在θ≈120°和θ≈113°处，但是实验得到的分离点位置是在θ≈147°处，采用2种模型模拟得到的分离点位置与实验结果都有一定差异，但采用Transition SST模型得到的模拟结果与实验结果更加接近。
图 9为临界区(Re=8.5×105) 条件下采用Transition SST模型模拟得到的二维圆柱表面附近流场区域的速度矢量图。图 10为图 9中Ⅰ~Ⅴ处的局部放大矢量图。
 |
| 图 9 在临界区二维圆柱速度矢量图 Fig. 9 Velocity vector diagram of 2D cylinder in critical region |
| 图选项 |
 |
| 图 10 圆柱上表面在图 9的局部放大速度矢量图 Fig. 10 Local amplification velocity vector diagram in Fig.9 on cylindrical upper surface |
| 图选项 |
从图 10(a)可以看到，在此区域内沿流向圆柱上表面外法线方向的速度梯度在逐渐增大，这使得沿流向圆柱上表面剪切应力逐渐变大，从而导致沿流向圆柱上表面时均摩擦力系数Cf逐渐变大，这也是图 8中θ从0°增大到60°时沿流向圆柱上表面时均摩擦力系数Cf逐渐增大的原因。
从图 10(b)可以看到，在此区域内沿流向圆柱上表面外法线方向的速度梯度在逐渐变小，这使得沿流向圆柱上表面剪切应力逐渐减小，从而导致沿流向圆柱上表面时均摩擦力系数Cf是逐渐减小的，这也是图 8中θ从60°增大到99°时沿流向圆柱上表面时均摩擦力系数Cf逐渐减小的原因。
从图 10(c)和图 10(d)可以看到，此区域圆柱上表面附近气流的速度矢量方向先是与主流方向相同，在θ≈99°时圆柱表面附近气流的速度矢量方向变为零，之后速度矢量与主流方向相反，在θ≈103°之后又从与主流方向相反变为与主流方向相同。这说明在θ≈99°~103°范围内流动发生了分离和再附的过程，即形成了一个分离泡。这也是图 8中在圆柱上表面θ≈99°~103°范围内时均摩擦力系数Cf为负值的原因，即Cf的正负符号2次发生变化的原因。
从图 10(e)可以看到，在θ≈120°之后圆柱表面附近流动方向与主流方向相反，说明在此处之后流动发生了分离。同样这也是图 8中在圆柱上表面θ≈120°~180°范围内Cf变为负值的原因。
图 11为圆柱上表面附近流场间歇因子γ分布云图。图中：γ=0代表层流流动，γ=1代表完全湍流，γ从0到1的变化过程即发生了转捩。通过图 11可以看出，在圆柱上表面θ≈103°处边界层发生了转捩，边界层的流动从层流流动变成了湍流流动，导致圆柱表面摩擦力突然增大。这即是图 8中在θ≈103°处圆柱上表面时均摩擦力系数Cf突然急剧增大，并在增大到最大值后回落的原因。
 |
| 图 11 圆柱上表面间歇因子云图 Fig. 11 Contour of intermittency factor on cylindrical upper surface |
| 图选项 |
图 12为圆柱上表面过顶点后局部速度云图和流线图，为了展示圆柱壁面附近的流动规律，对壁面附近的流线进行了加密处理。从图 12也可以看出，过圆柱上表面顶点之后在θ≈99°~103°的区域内边界层出现了流动分离并再附的过程，即产生了一个分离泡，与图 10(c)和图 10(d)的结果对应一致。分离再附之后，沿圆柱表面向后流动，在θ≈120°处流动最终发生湍流分离，形成尾迹流。
 |
| 图 12 圆柱上表面速度云图和流线图 Fig. 12 Contours of velocity and streamlines on cylindrical upper surface |
| 图选项 |
与亚临界区圆柱绕流问题的情况相同，由于Transition SST模型对分离和压力梯度等影响转捩的因素较为敏感，使得采用Transition SST模型能够很好地捕捉到圆柱绕流边界层内的转捩过程以及圆柱绕流尾迹近壁区的流动特征，从而可以更好地模拟临界区的圆柱绕流问题。
3.3 不同雷诺数下的模拟结果对比分析 图 13为采用Transition SST模型在不同雷诺数下对二维圆柱绕流模拟计算得到的圆柱表面时均摩擦力系数Cf和圆柱边界层的间歇因子γ分布曲线。图 14为模拟所得到的不同雷诺数条件下圆柱上表面附近局部流场的速度分布云图和流线图，为了展示圆柱壁面附近的流动变化，对圆柱壁面附近的流线进行了加密处理。
 |
| 图 13 不同雷诺数下圆柱表面时均摩擦力系数Cf和间歇因子γ分布曲线 Fig. 13 Distribution curves of time-averaged friction coefficient Cf and intermittency factor γ on cylindrical surface under different Reynolds number |
| 图选项 |
 |
| 图 14 不同雷诺数下圆柱上表面附近局部流场速度分布云图和流线图 Fig. 14 Contours of local flow velocity and streamlines near cylinder under different Reynolds number on cylindrical upper surface |
| 图选项 |
由图 13(a)可以看出，当雷诺数Re=1.4×105时(亚临界区)，在θ≈85°处圆柱表面时均摩擦力系数Cf直接从正值变为负值。从图 14(a)的流场结构可以看到，这是因为流体在绕圆柱流动过程中，在θ≈85°处边界层分离后直接形成了背风面底部分离涡流动区。
从图 13(a)可以看到，对应Re=5.0×105、8.5×105(临界区) 和Re=1×106(超临界区)，圆柱表面时均摩擦力系数Cf的符号分别在θ≈94°、99°和90°处从正值变成负值，之后分别在θ≈107°、103°和97°处Cf的符号又由负值变成正值，从图 14(b)、图 14(c)和图 12的流场结构中可以看出这是因为流体在绕圆柱流动过程中发生了一个分离到再附的过程，即形成了一个分离泡。
从图 13(a)可以看到，当Re=5.0×105时，在θ≈107°处圆柱上表面的时均摩擦力系数Cf发生了突然增大之后回落减小的过程，而且由图 13(b)可见，在此处边界层间歇因子γ发生了从0到1的变化过程，说明流动在θ≈107°处发生转捩；同样可以看出，当Re=8.5×105和1×106时，圆柱表面边界层内的流动发生了从层流到湍流的转捩过程，且转捩点分别在θ≈103°和θ≈97°附近。从以上结果可知，对于二维圆柱绕流问题，当雷诺数在临界区和超临界区时，随着雷诺数的增大，圆柱边界层转捩点的位置逐渐向来流方向移动。
从图 13(a)还可以看到，当雷诺数Re=3×106(过临界区) 时，在圆柱上表面(或下表面) 时均摩擦力系数Cf仅在θ≈120°时发生了一次从正变为负过程，结合图 14(d)的流场结构可知，在此雷诺数条件下，圆柱绕流直接发生了底部流动分离，在底部流动分离之前没有出现分离泡。从图 13(b)可以看出，当Re=3×106时，在圆柱上表面θ≈63°处边界层间歇因子γ发生了从0到1的变化过程，说明在θ≈63°时边界层流动发生了从层流到湍流的转捩，但图 13(a)中在Re=3×106时圆柱表面时均摩擦力系数Cf的值并未出现突增的现象，这是因为转捩前Cf的值就很大，而转捩引起的Cf的变化较小，导致转捩时Cf的变化不明显。
4 结论 本文采用SST k-ω两方程湍流模型与γ-
1) 采用Transition SST模型能够很好地模拟圆柱绕流问题。这是由于Transition SST模型对流动分离和压力梯度等因素较为敏感，所以可以很好地捕捉到不同雷诺数下圆柱绕流尾迹近壁区的流动特征，并且能够较为准确地预测边界层转捩位置及转捩过程。
2) 当自由来流雷诺数处于临界区、超临界区和过临界区时，圆柱绕流边界层会出现从层流到湍流的转捩过程，且圆柱绕流边界层转捩位置随雷诺数增大逐渐靠近圆柱前驻点(即θ=0°)。
3) 当自由来流雷诺数处于亚临界区时，圆柱绕流边界层发生层流分离，无分离泡及转捩现象发生；当自由来流雷诺数处于临界区和超临界区时，圆柱绕流边界层产生分离泡后发生转捩现象，并在转捩后圆柱绕流边界层发生湍流分离；当自由来流雷诺数处于过临界区时，圆柱绕流边界层无分离泡现象发生，且在转捩之后发生湍流分离。
参考文献
| [1] | 伊杰里奇克N E. 水力摩阻手册[M].北京: 航空发动机编辑部, 1985: 256-300.YI JIE RITCHIE N E. Hydraulic friction handbook[M].Beijing: Aviation Engine Newsroom, 1985: 256-300.(in Chinese) |
| [2] | 凌国灿. 圆柱绕流旋涡运动及离散涡方法[J].力学进展, 1985, 15(4): 458–470.LING G C. Vortex flow around a circular cylinder motion and the discrete vortex method[J].Journal of Mechanics, 1985, 15(4): 458–470.(in Chinese) |
| [3] | 苑明顺. 高雷诺数圆柱绕流的二维大涡模拟[J].水动力学研究与进展:A辑, 1992, 7(Suppl.): 614–622.YUAN M S. High Reynolds number flow around a cylinder of 2-D large eddy simulation[J].Water Dynamics Research and Development:Series A, 1992, 7(Suppl.): 614–622.(in Chinese) |
| [4] | 王亚玲, 刘应中, 缪国平. 圆柱绕流的三维数值模拟[J].上海交通大学学报, 2001, 35(10): 1464–1469.WANG Y L, LIU Y Z, MIU G P. Three dimensional numerical simulation of viscous flow a round circular cylinder[J].Journal of Shanghai Jiaotong University, 2001, 35(10): 1464–1469.(in Chinese) |
| [5] | 刘应中, 缪国平. 高等流体力学[M].上海: 上海交通大学出版社, 2000.LIU Y Z, MIU G P. Higher fluid mechanics[M].Shanghai: Shanghai Jiao Tong University Press, 2000.(in Chinese) |
| [6] | 易家训. 流体力学[M].北京: 高等教育学出版社, 1982.YI J X. Fluid mechanics[M].Beijing: Higher Education Press, 1982.(in Chinese) |
| [7] | ESKINAZI S. Principles of fluid mechanics[M].Boston: Allyn and Bacon Inc., 1968. |
| [8] | SCHMIDT R, PANTANKAR S. Simulation boundary layer transition with low-reynolds number k-ε turbulence models:Part 2-An approach to improving the predictions[J].Journal of Tuebomachinery, 1991, 113(1): 18–26.DOI:10.1115/1.2927732 |
| [9] | VANCOILLIE G, DICK E. A turbulence model for the numerical simulation of the transition zone in a boundary layer[J].International Journal of Engineering Fluid Mechanics, 1988, 1(1): 28–49. |
| [10] | STEELANT J, DICK E. Modeling of laminar-turbulent transition for high freestream turbulence[J].Journal of Fluids Engineering, 2001, 123(1): 22–30.DOI:10.1115/1.1340623 |
| [11] | WALTERS K D, LEYLEK J H. Impact of film-cooling jets on turbine areogynamic losses[J].Journal of Turbomachinery, 2000, 122(3): 537–545.DOI:10.1115/1.1303818 |
| [12] | LANGTRY R B, MENTER F R.Transition modeling for general CFD applications in aeronautics[C]//43rd AIAA Aerospace Sciences Meeting and Exhibit.Reston:AIAA, 2005. |
| [13] | MENTER F R, LANGTRY R B, LIKKI S R, et al.A correlation-based transition model using local variables-Part I:Model formulation[C]//49th International Gas Turbine and Aeroengine Congress and Exhibition.New York:ASME, 2006:413-422. |
| [14] | MENTER F R, LANGTRY R B, LIKKI S R, et al.A correlation-based transition model using local variables-Part II:Test cases and industrial applications[C]//49th International Gas Turbine and Aeroengine Congress and Exhibition.New York:ASME, 2006:423-434. |
| [15] | MENTER F R. Two-equation eddy-viscosity turbulence models for engineering application[J].AIAA Journal, 1994, 32(8): 1598–1605.DOI:10.2514/3.12149 |
| [16] | CANTWELL B J.A flying hot wire study of the turbulent near wake of a circular cylinder at a Reynolds number of 14, 000[D].Pasadena:California Institute of Technology, 1976. |
| [17] | SCHEWE G. On the force fluctuations acting on a circular cylinder in cross flow from subcritical up to transcritical Reynolds numbers[J].Journal of Fluid Mechanics, 1983, 133(1): 265–285. |
| [18] | SQUIRES K D, KRISHNAN V, FORSYTHE J R. Prediction of the flow over a circular cylinder at high Reynolds number using detached-eddy simulation[J].Journal of Wind Engineering and Industrial Aerodynamics, 2008, 96(10-11): 1528–1536.DOI:10.1016/j.jweia.2008.02.053 |
| [19] | ACHENBACH E. Distribution of local pressure and skin friction around a circular cylinder in cross-flow up to Re=5×106[J].Journal of Fluid Mechanics, 1968, 34(4): 625–639.DOI:10.1017/S0022112068002120 |
| [20] | 凌国灿, 牛家玉, CHWANGA T. 圆柱近尾迹流动转捩特征的数值研究[J].空气动力学报, 1997, 15(1): 130–138.LING G C, ZHU J Y, CHWANG A T. Numerical studies of cylindrical near wake flow transition characteristics[J].Acta Aerodynamica Sinica, 1997, 15(1): 130–138.(in Chinese) |
| [21] | WALKER J D A, SMITH C R, CERRA A W, et al. The impact of a vortex ring on a wall[J].Journal of Fluid Mechanics, 1987, 181: 99–140.DOI:10.1017/S0022112087002027 |
| [22] | JONES G W, CINCOTTA J J, WALKER R W.Aerodynamic forces on a stationary and oscillating circular cylinder at high Reynolds numbers:NASA TR R-300[R].Washington, D.C.:NASA, 1969. |
