流体力学Fluent报告——圆柱绕流

流体力学Fluent报告——圆柱绕流亚临界雷诺数下串列双圆柱与方柱绕流的数值模拟摘要：本文运用Fluent软件中的RNG k-ε模型对亚临界雷诺数下二维串列圆柱和方柱绕流问题进行了数值研究，通过结果对比，分析了雷诺数、柱体形状对柱体绕流阻力、升力以及涡脱频率的影响。

一般而言，Re数越大，方柱的阻力越大，圆柱体则不然；而Re越大，两种柱体的升力均越大。

相对于圆柱，同种条件下，方柱受到的阻力要大；相反地，方柱涡脱落频率要小。

Re越大，串列柱体的Sr数越接近于单圆柱体的Sr数。

关键字：圆柱绕流、升力系数、阻力系数、斯特劳哈尔数在工程实践中，如航空、航天、航海、体育运动、风工程及地面交通等广泛的实际领域中，绕流研究在工程实际中具有重大的意义。

当流体流过圆柱时 , 由于漩涡脱落，在圆柱体上产生交变作用力。

这种作用力引起柱体的振动及材料的疲劳，损坏结构，后果严重。

因此，近些年来，众多专家和学者对于圆柱绕流问题进行过细致的研究，特别是圆柱所受阻力、升力和涡脱落以及涡致振动问题。

沈立龙等[1]基于RNG k⁃ε模型，采用有限体积法研究了亚临界雷诺数下二维圆柱和方柱绕流数值模拟，得到了圆柱和方柱绕流阻力系数C与 Strouhal 数d随雷诺数的变化规律。

姚熊亮等[2]采用计算流体软件CFX中LES模型计算了二维不可压缩均匀流中孤立圆柱及串列双圆柱的水动力特性。

使用非结构化网格六面体单元和有限体积法对二维N- S方程进行求解。

他们着重研究了高雷诺数时串列双圆柱在不同间距比时的压力分布、阻力、升力及Sr数随Re数的变化趋势。

费宝玲等[3]用FLUENT软件对串列圆柱绕流进行了二维模拟，他们选取间距比L/D(L为两圆柱中心间的距离，D为圆柱直径)2、3、4共3个间距进行了数值分析。

计算均在 Re = 200 的非定常条件下进行。

计算了圆柱的升阻力系数、尾涡脱落频率等描述绕流问题的主要参量，分析了不同间距对圆柱间相互作用和尾流特征的影响。

圆柱绕流的一个重要特征是流动形态取决于雷诺数。

Lienhard[4]总结了大量的实验研究结果并给出了圆柱体尾流形态随雷诺数变化的规律。

当Re<5时，圆柱上下游的流线呈对称分布,流体并不脱离圆柱体，没有旋涡产生。

此时与理想流体相似，若改变流向，上下游流形仍相同。

当5<Re<40时,边界层发生分离,分离剪切层在圆柱体背后形成一对稳定的“附着涡”。

当40<Re<150时，流动保持层流状态并且流体旋涡交替地从圆柱后部作周期性的脱落并在尾流中形成两列交叉排列的涡，即卡门涡街。

从150<Re<300开始，旋涡内部开始由层流向湍流转捩，直至增加至3x105左右，此时圆柱体表面附近的边界层仍为层流，整个涡街逐渐转变为湍流，及e<3xl05称为亚临界区域。

当3xl05<Re<3.5x106时,边界层的流动也逐渐趋于湍流状态，尾流中没有明显的涡街结构，称为临界状态。

[5]圆柱绕流的另一个显著特征是斯特劳哈尔数是雷诺数的函数。

早在1878年，捷克科学家Strouhal [6]就对风吹过金属丝时发出鸣叫声作过研究，发现金属丝的风鸣音调与风速成正比,同时与弦线之粗细成反比，并提出计算涡脱落频率f 的经验公式:19.70.198(1)RefD Sr v ==- 式中即斯特劳哈尔数Sr 由Re 所唯一确定。

本文运用Fluent 软件中的RNG k-ε模型对亚临界雷诺数下二维串列圆柱和方柱绕流问题进行了数值研究，通过结果对比，分析了雷诺数、柱体形状对柱体绕流阻力、升力以及涡脱频率的影响。

1. 数学模型1.1 控制方程对于静止圆柱绕流，本文研究对象为二维不可压缩流动。

在直角坐标系下，其运动规律可用 N-S 方程来描述，连续性方程和动量方程分别为: 01()()i i i i i j ji j j u x u u P u u t x x x x νρ∂⎧=⎪∂⎪⎨∂∂∂∂∂⎪+=-+⎪∂∂∂∂∂⎩其中u i 为速度分量；p 为压力；ρ为流体的密度；ν为流体的动力黏性系数。

对于湍流情况，本文采用RNG k ⁃ε模型，RNG k ⁃ε模型是k ⁃ε模型的改进方案。

通过在大尺度运动和修正后的粘度项体现小尺度的影响， 而使这些小尺度运动有系统地从控制方程中去除。

所得到的 k 方程和ε方程，与标准k ⁃ε模型非常相似，其表达式如下：212()()()G ()()()G i k eff k i j j i eff k i j j ku k k t x x x u C C t x x x k k εεερραμρερεερεεεαμρ∂∂∂∂⎧+=++⎪∂∂∂∂⎪⎨∂∂∂∂⎪+=+-⎪∂∂∂∂⎩其中G k 为由于平均速度梯度引起的湍动能的产生项，eff t μμμ=+，2t C k ρμμε=，经验常数 1C ε=0.084 5，k α=εα=1.39，2C ε=1.68。

相对于标准 k ⁃ε模型，RNG k ⁃ε模型通过修正湍动粘度，考虑了平均流动中的旋转及旋转流动情况，RNG k ⁃ε模型可以更好的处理高应变率及流线弯曲程度较大的流动。

1.2 相关参数圆柱绕流的相关参数主要有雷诺数 Re 、斯特劳哈尔数Sr 、升力系数C l 和阻力系数 C d ，下面给出各个参数的计算公式和物理意义。

雷诺数Re 与圆柱绕流的状态和雷诺数有很大关系，雷诺数代表惯性力和粘性力之比：Re UL UL ρμν== 其中U 为来流速度；L 为特征长度，本文取圆柱直径或方柱边长；ρ为流体密度；μ、ν分别为流体介质动力粘度和运动粘度。

斯特劳哈尔数 Sr 是Strouhal 指出圆柱绕流后在圆柱后面可以出现交替脱落的旋涡，旋涡脱落频率、风速、圆柱直径之间存在一个关系：fL Sr U= 式中： Sr 为斯托罗哈数，取决于结构的形状断面； f 为旋涡脱落频率； L 为结构的特征尺寸； U 为来流速度。

阻力系数和升力系数是表征柱体阻力、升力的无量纲参数。

定义为： 212DD F C V A ρ= ，212L L F C V A ρ=式中ρ为流体密度； V 为来流速度；A 为迎流截面面积；D F 和L F 为柱体所受阻力和升力。

由于涡脱落的关系，阻力系数将产生振荡，本文选取平均脉动升力来研究，即取方均根值来研究。

2. 数值计算2.1 物理模型二维数值模拟双圆柱流场计算区域的选取如图 1所示，圆柱绕流以圆柱体直径为特征尺度 D ，选取圆柱半径为1.5 mm ，计算区域为9D × 32D 的矩形区域。

柱1距上游长度 5D ，下游长度27D ，保持两柱间距 L/D= 2. 5D 不变 (L 是两圆柱中心连线长度 )，两柱到上下边界距离相等。

对于方柱绕流，选择方柱边长为特征长度，D=30mm 。

图 1 串列圆柱和方柱的计算区域2.2网格划分计算区域采用分块结构化网格，柱体表面网格做加密处理，边界区网格相对稀疏。

具体网格划分情况见图 2。

其中串列圆柱网格31116个节点，30615个四边形面单元；串列方柱46446个节点，46550个四边形面单元。

2.3边界条件管道壁面和柱体表面均采用无滑移的静止壁面条件。

而入口选择速度入口，出口选择自由出流。

来溜速度大小根据Re来设置，雷诺数分300、3000、12000、30000四个等级，速度大小依次为0.1m/s、1m/s、4m/s、10m/s。

2.4计算模型本文湍流模型采用标准壁面函数的RNG k-ε模型。

采用有限容积法求解二维不可压缩粘性流体非定常流动控制方程，即把计算区域分成很圆柱近壁面网格多小的控制体，对每个控制体的各个变量进行积分。

控制方程的对流项采用二阶迎风格式离散，速度和压力采用SIMPLE算法耦合求解，将所有区域看成一个整体进行耦合计算。

动量、湍动能和湍动耗散率均采用二阶迎风格式。

先定常计算流场，再用定常计算的结果作为非定常迭代的初始值进行计算。

根据初略计算的涡脱频率，固定设置时间步长为0. 002s, 在每个时间步内设置迭代次数为20。

流体介质为液态水。

3.计算结果3.1网格模型验证为验证网格独立性，本文计算了网格节点数为8346，面单元为8932的粗网格、节点数为31116，面单元为30615的密网格、节点数为63432，面单元为67434的精密网格下Re=200、L/D=2的串列网格的Sr数，结果显示三套网格的计算结果分别为0.143、0.133、0.133。

故密网格可用。

而方柱绕流则采用同级别网格。

为验证本文计算思路与模型的正确性，本文计算了Re=200串列圆柱不同间距上下游圆柱的斯特劳哈尔数，将其结果与G.X.Wu[7]的计算数据相比较，比较图像如图3所示，最大误差为2.2%。

图 3 串列圆柱不同间距的Sr数计算对比3.2流线与涡量图6 Re=3000圆柱绕流流线图图图4 Re=3000方柱绕流流线图本文给出了计算过程中雷诺数Re=3000，t=1s时的流线图和涡量图。

3.3阻力系数本文给出了Re=3000时，圆柱绕流和方柱绕流的脉动阻力系数图如下。

图9 Re=3000圆柱绕流脉动阻力系数图8 Re=3000方柱绕流脉动阻力系数由图 9和错误!未找到引用源。

可以看到，经过一段时间后圆柱和方柱的阻力系数是振荡变化的，这是由于涡脱落流场压力发生剧烈变化，从而导致柱体表面受力不断振荡。

而相比较而言，圆柱绕流到达充分发展湍流要比方柱绕流需要的时间短。

从图 10可以看到，对于串列的柱体，上游受到阻力大于下游收到的阻力，甚至于在部分雷诺数下，某些时刻下游柱体收到的阻力会出现负值，这是由于在上游柱体尾部边界层分离，形成低压区，而漩涡的形成更是会导致柱体后面的压力动态地减小，从而使得下游柱体所受前后压差可能为负，导致下游柱体受到向前的推力。

随着雷诺数地增大，柱体收到的阻力变化也具有一定的规律。

对于圆形柱体，上游圆柱体所受阻力随着雷诺数的增大而减小，这是由于随着雷诺数的变大，边界层分离点更加靠近上游，导致前后压差变小。

而对于方形柱体，可以图10 不同雷诺数圆柱与方柱平均阻力系数看到上下游柱体阻力系数均随着雷诺数地增大而增大。

比较柱体不同形状的阻力系数。

方柱柱体上下游柱体所受到的阻力均比圆形柱体上下游柱体所受到的阻力大。

3.4升力系数本文给出了Re=3000时的圆柱绕流和方柱绕流上下游各圆柱体表面的脉图12 Re=3000圆柱绕流脉动升力系图11 Re=3000方柱绕流升力系数动升力系数变化曲线如图 12和图 11所示。

类似于阻力系数，升力系数由于漩涡的影响也将产生振荡。

圆柱绕流到达充分发展湍流要比方柱绕流需要的时间短。

图13 不同雷诺数圆柱与方柱平均升力系数由图13可以看到相对于圆柱绕流，方柱绕流两柱体所受到的升力要大一些。