流体力学fluent报告——圆柱绕流.pdfVIP

亚临界雷诺数下串列双圆柱与方柱绕流的数值模拟 摘要 ：本文运用Fluent 软件中的 RNG k-ε 模型对亚临界雷诺数下 二维串列圆柱和方柱绕流 问题进行了数值研究，通过结 对比，分析 了雷诺数、柱体形状对柱体绕流阻力、升力以及涡脱 率的影响。一 般而言，Re 数越大，方柱的阻力越大，圆柱体则不然；而 Re 越大，两种 柱体的升力均越大。相对于圆柱，同种条件下，方柱受到的阻力要大； 相反地，方柱涡脱落 率要小。Re 越大，串列柱体的 Sr 数越接近于单 圆柱体的 Sr 数。 关键字 ：圆柱绕流、升力系数、阻力系数、斯特劳哈尔数 在工程实践中，如航空、航天、航海、体育运动、风工程及地面交 通等广泛的实际领域中，绕流研究在工程实际中具有重大的意义。当 流体流过圆柱时 , 由 于漩涡脱落，在圆柱体上产生交变作用力。这 种作用力引起柱体的振动及材料的疲劳，损坏结构，后 严重。因此， 近些年来，众多专家和学者对于圆柱绕流 问题进行过细致的研究，特 别是圆柱所受阻力、升力和涡脱落以及涡致振动问题。 沈立龙等[1]基于 RNG k⁃ε 模型，采用有限体积法研究了亚临界雷 诺数下二维圆柱和方柱绕流数值模拟，得到了圆柱和方柱绕流阻力系 数 Cd 与 Strouhal 数随雷诺数的变化规律。姚熊亮等[2]采用计算流体 软件 CFX 中 LES 模型计算了二维不可压 缩均匀流中孤立圆柱及串 列双圆柱的水动力特性。使用非结构化网格六面体单元和有限体积法 对二维 N- S 方程进行求解。他们着重研究了高雷诺数时串列双圆柱 在不同间距比时的压力分布、阻力、升力及 Sr 数随 Re 数的变化趋势。 费宝玲等[3]用 FLUENT 软件对串列圆柱绕流进行了二维模拟，他们选 取间距比 L/D(L 为两圆柱中心间的距离，D 为圆柱直径)2、3、4 共 3 个 间距进行了数值分析。计算均在 Re = 200 的非定常条件下进行。计 算了圆柱的升阻力系数、尾涡脱落 率等描述绕流 问题的主要参量， 分析了不同间距对圆柱间相互作用和尾流特征的影响。 圆柱绕流的一个重要特征是流动形态取决于雷诺数。Lienhard[4] 总结了大量的实验研究结果并给出了圆柱体尾流形态随雷诺数变化 的规律。当 Re5 时，圆柱上下游的流线呈对称分布,流体并不脱离圆 柱体，没有旋涡产生。此时与理想流体相似，若改变流向，上下游流形 仍相同。当 5Re40 时,边界层发生分离,分离剪切层在圆柱体背后形 成一对稳定的“附着涡”。当 40Re150 时，流动保持层流状态并且流 体旋涡交替地从圆柱后部作周期性的脱落并在尾流中形成两列交叉 排列的涡，即卡门涡街。从 150Re300 开始，旋涡内部开始由层流向 湍流转捩，直至增加至 3x105 左右，此时圆柱体表面附近的边界层仍 为层流，整个涡街逐渐转变为湍流，及 e3xl05 称为亚临界区域。当 5 6 3xl0 Re3.5x10 时,边界层的流动也逐渐趋于湍流状态，尾流中没有 明显的涡街结构，称为临界状态。[5] 圆柱绕流的另一个显著特征是斯特劳哈尔数是雷诺数的函数。早 在 1878 年，捷克科学家 Strouhal[6]就对风吹过金属丝时发出鸣叫声作 过研究，发现金属丝的风鸣音调与风速成正比,同时与弦线之粗细成 反比，并提出计算涡脱落 率 f 的经验公式: f D 19.7 Sr   0.198(1 ) v Re 式中即斯特劳哈尔数 Sr 由 Re 所唯一确定。 本文运用 Fluent 软件中的 RNG k-ε 模型对亚临界雷诺数下二维串 列圆柱和方柱绕流 问题进行了数值研究，通过结 对比，分析了雷诺 数、柱体形状对柱体绕流阻力、升力以及涡脱 率的