1、 第 卷 第 期 年 月：多 圆 柱 绕 流 的 数 值 模 拟 及 力 学 分 析赵 红 晓（同济大学 航空航天与力学学院，上海）摘 要：基于计算流体力学理论，采用 对 种不同雷诺数下多圆柱绕流模型进行了仿真研究，得到了速度云图和升阻力系数随时间变化的曲线。仿真结果表明：尾涡区的长度和紊动强度随来流的雷诺数增大而增强。圆柱的阻力系数随雷诺数的增加而降低，流动由层流变为紊流，且湍动更强烈；第 列小圆柱升力系数随雷诺数的增大而逐渐减小。从第 列至第 列升力系数振荡幅度逐列增加，第列圆柱的升力系数振荡幅度较大，最大值为其他列圆柱的 倍以上；时的速度云图和实验装置的流动图像接近，与数值模拟实验装置的
2、流动规律相似。关键词：多圆柱绕流；数值模拟；力学分析中图分类号：.文献标志码：文章编号：（）（，）：，：；收稿日期：基金项目：同济大学实验教改项目（）作者简介：赵红晓（），女，河南郏县人，博士，高级工程师，研究方向为固体力学、实验力学。：；：引 言 圆柱的绕流现象在工程中广泛存在，如管排、桥墩、电缆和高耸建筑等结构周围的流动问题。年 月 日，美国华盛顿州塔克马海峡大桥因风振致毁，事故原因正是空气绕流形成的卡门涡街引起的桥梁共振。作为典型钝体的一种，近代学者对圆柱绕流的研究一直没有停息过。年，冯卡门观察到在一定条件下的定常来流绕过某些物体时，物体两侧会周期性地脱落出旋转方向相反、并排列成有规则的
3、双列线涡，并对涡街的诱导速度、稳定性和阻力等做了分析。第 期赵红晓：多圆柱绕流的数值模拟及力学分析等揭示了圆柱后尾流形态和雷诺数之间的对应关系。研究了双圆柱串列在流体流动速度很小时，不同的间距比 （圆柱间距与圆柱直径之比）对两柱间隙中的速度分布和下游柱尾流场中的速度分布具有明显影响。等发现双圆柱绕流在间距比 时尾流产生的涡街迅速消失，间距比为.时两条尾流的形状随时间随机发生变化，出现涡脱落方向相同、相反等 种典型的尾流形状。等研究了双柱并列（非等直径），间距比为.之间时流体形态发生交替偏斜，间距比大于.时，两柱后分别形成各自的涡脱落。张爱社等利用有限元法对等边三角形布置的等直径 圆柱绕流问题进
4、行了研究，表明小圆柱的位置、大小、个数以及雷诺数变化等对主圆柱有影响。等研究了 个并列小圆柱置入尾流区，可有效地抑制涡街的产生，并使主圆柱的升力系数降低了 ，压力、阻力减小也很明显。张洪军等研究了附属小圆柱对并列双圆柱尾流流场的影响。等研究了放置于主圆柱后的旋转小圆柱可以有效地抑制涡街的生成。由于流动的复杂性，目前以多圆柱群为对象的研究相对较少，工程情况中圆柱以圆柱群的形式出现的案例却十分常见，本文选取不等直径多圆柱群实验装置进行仿真研究。数值模拟.控制方程本文采用 软件进行数值模拟，选用 算法（半隐式方法 ）进行计算。本文研究中，低雷诺数下水的流动，可近似为不可压缩流体，则二维连续性方程可简
5、化为：（）动量方程为动量守恒定律在流体力学中的表示形式：()（）()（）式中：、为 个速度分量；为压强；为流体密度；为流体动力黏度。升力系数 和阻力系数 是流体力学中的无因次量，用来表示物体在流体中的升力和阻力大小（）（）式中：、分别为物体受到的升力和阻力；为无限远场的来流速度；为物体在来流方向上的投影面积。雷诺数 （）式中，、分别为流场的特征速度和特征长度。雷诺数物理上表示惯性力和黏性力量级的比。本文均设置左侧边界为速度入口，右侧边界为压力出口，远场壁面为对称边界条件，圆柱壁面采用无滑移壁面。流体介质为水，密度 ，动力黏度 .（）。.计算模型多圆柱绕流研究实验装置如图 所示，以还原实验装置中
6、的流动规律，划分模型进行多组仿真计算。图 流态演示装置断面示意图 为避免过于庞大的模型造成不必要的计算时间浪费，需对模型进行适当的简化和划分。如图 所示，模型保留了实验装置中近壁面绕流和多圆柱群绕流的流域，并在入口和出口处做了适当的加长。与实验室现象进行对比，观察能否正确模拟实际流动，保留了实验装置中的栅格，多圆柱排列方式、尺寸和编号如图 所示。（）模型形状及尺寸（）（）栅格形状及尺寸（）图 实验装置物理模型及尺寸.网格划分及计算方法验证采用结构化网格的方法进行网格划分，建立 与圆柱面以及栅格前缘、尾缘的圆角关联，将其第 卷（）多圆柱排列方式及尺寸（）（）多圆柱编号图 多圆柱群尺寸及编号示意余
7、边与壁面和栅格长边关联，并在湍流情况中，运用式（）估算第 层网格到柱面的距离 进行相应加密。最终得到的网格具有较好的质量和长宽比，网格划分如图 所示。图 实验装置网格划分 网格划分与层流验证模型时相同。为保证圆柱面上第 层网格满足，估算确定第 层网格到柱面的距离 .（）代入特征长度 .，得到.，加密后网格如图 所示。（）局部圆柱网格（）圆柱周围网格加密图 圆柱周围及局部网格加密图 层流模型和紊流模型验证参考文献，使用相同的物理模型进行层流验证。左侧为速度入口，右侧为压力出口，远场壁面为对称边界条件，圆柱壁面采用无滑移壁面。计算区域，采用结构化网格，计算参数见表，计算结果见表。表 层流验证模型仿
8、真参数雷诺数计算模型圆柱直径 来流速度（）时间步长 残差层流 表 层流模型验证结果与文献数据对比数据来源文献本文 由表 可知，计算结果与文献基本一致，可使用层流模型对小雷诺数的情况进行模拟。参考文献中使用相同的物理模型进行紊流验证。左侧为速度入口，右侧为压力出口，远场壁面为对称边界条件，圆柱壁面采用无滑移壁面。计算区域，靠近圆柱加密网格，采用结构化网格，计算参数见表，计算结果见表。表 湍流验证模型仿真参数雷诺数计算模型圆柱直径 来流速度（）时间步长 残差 分离涡 表 湍流模型验证结果与文献数据对比数据来源平均阻力系数，斯特劳哈尔数，本文文献 计算结果与文献基本一致，可使用分离涡模型进行模拟。仿
9、真计算分别采用以上验证计算模型和划分网格的方法分别对多圆柱绕流模型（见图（）进行数值模拟，表 中雷诺数均为小圆柱雷诺数，参数设置见表。表 实验装置中多圆柱绕流仿真参数雷诺数计算模型大圆柱直径 小圆柱直径 来流速度（）时间步长 残差层流 分离涡 分离涡 由速度云图图（）可见，在 时，以大圆柱的前缘为中心，水流在驻点受阻后，被迫向两边流动，在靠近圆柱尾部的边界上，主流开始与圆柱体分离，边界层分离后，在分离区的下游形成回流区，即为尾涡区。后流经栅格区域，在流经栅格的各狭缝中，由于截面积的减少，流动皆有一定的加速。流体流经多列小圆柱时，列与列之间的尾涡出现掺混而紊动更强烈，涡的长度由于受到限制变短。第
10、 列圆柱则直接在后方产交替脱落涡。由图（）速度云图可见，当 时，大圆柱尾涡出现明显的两列涡，且脱落形成长度明显变短，涡街的宽度变窄，尾流呈复杂的湍流特征涡。小圆柱的第 列圆柱限制了后面第、列圆柱分离剪切层的卷起，使其不能产生强烈的旋涡脱落，第、列圆柱的涡在后方混合，形成更为复杂的尾流。各圆柱后方仍有一定范围的回流泡，旋涡脱落现象在下游第 列圆柱后尾涡脱落成偏斜的尾流模式，随着雷诺数的增大，第 期赵红晓：多圆柱绕流的数值模拟及力学分析（）（）（）图 不同雷诺数下实验装置中多圆柱绕流速度云图第 列后的紊动就越强烈。当 时，从图（）的速度云图可观察到其流动与 时十分相似，但紊动更强烈，尾流模式和圆柱
11、受力不具有周期性，是由于尾涡区的长度和紊动强度与来流的雷诺数有关，雷诺数越大，紊动越强烈。当 时，由图（）、（）可见，大圆柱受到最大阻力系数最大，小圆柱中第 列（从左向右看）圆柱阻力次之，第 列圆柱所受阻力略小于第 列圆柱，位于中心线上的圆柱、仍相比同列受到更大的阻力。第 和第 列圆柱所受阻力比第 列略小，第 列圆柱所受阻力略大于第 列。第 列圆柱升阻力系数的振幅比前几列大。当 时，由图（）、（）和（）可见，大圆柱的最大阻力系数仍是最大，但比 时小。圆柱 阻力系数次之，第 列圆柱的阻力系数比第 列大，位于中心线上的圆柱、仍相比同列其他圆柱受到更大的阻力，各圆柱的阻力系数均有一定幅度的振荡，后方
12、的圆柱较前方有更大的振荡幅度。（）时大圆柱和小圆柱 ，和 的阻力系数随时间变化曲线（）时圆柱 和 的阻力系数随时间变化曲线（）时大圆柱和小圆柱 ，阻力系数随时间变化曲线（）时圆柱 的阻力系数随时间变化曲线（）时圆柱 和 的阻力系数随时间变化曲线（）时大圆柱和小圆柱 ，的阻力系数随时间变化曲线（）时圆柱 的阻力系数随时间变化曲线（）时圆柱 和 的阻力系数随时间变化曲线注：图中：为大圆柱的阻力系数；分别为小圆柱 的阻力系数；分别为小圆柱 的阻力系数；分别为小圆柱 的阻力系数；分别为小圆柱 的阻力系数。图 多圆柱绕流阻力系数随时间变化曲线第 卷 当 时，由图（）、（）和（）可见，各圆柱的阻力系数均有
13、一定幅度的下降和振荡，后方的圆柱较前方有更大的振荡幅度。同样位于中心线上的圆柱、仍相比同列受到更大的阻力。由图 中不同雷诺数下的阻力系数可见，大圆柱的升力系随着雷诺数的增加而降低，小圆柱的绕流阻力系数总体上也随着雷诺数的增加而降低，随着雷诺数的增加，流动由层流变为紊流，湍动更强烈，各阻力系数均随时间产生微小振荡，振幅逐列增加。当 时，由图（）（）可见，第 列小圆柱升力系数大约在 .这一范围内振荡，其次是圆柱、，圆柱、。圆柱、升力系数在零上下振荡，对称轴上的大圆柱和小两圆柱、所受升力为零。第、列和第、列同侧圆柱所受升力仍具反向性。当 时，由图（）（）可见，升力系数的振荡幅度也逐列增加，第 列圆柱
14、的升力系数振幅较大，最大值为其他圆柱的 倍以上，第 列升力系数次之，第 列最小。不过前 列圆柱的振荡较小，除去中心线上的圆柱、，其他圆柱受力方向也没发生（）时大圆柱和小圆柱 的升力系数随时间变化曲线（）时小圆柱 的升力系数随时间变化曲线（）时小圆柱 的升力系数随时间变化曲线（）时小圆柱 的升力系数随时间变化曲线（）时大圆柱和小圆柱 的升力系数随时间变化曲线（）时小圆柱 的升力系数随时间变化曲线（）时小圆柱 的升力系数随时间变化曲线（）时小圆柱 的升力系数随时间变化曲线（）时大圆柱和小圆柱 的升力系数随时间变化曲线（）时小圆柱 的升力系数随时间变化曲线（）时小圆柱 的升力系数随时间变化曲线（）时
15、小圆柱 的升力系数随时间变化曲线注：图中 为大圆柱的升力系数；分别为小圆柱 的升力系数；分别为小圆柱 的升力系数；分别为小圆柱 的升力系数；分别为小圆柱 的升力系数。图 多圆柱绕流升力系数随时间变化曲线 第 期赵红晓：多圆柱绕流的数值模拟及力学分析改变，而大圆柱和小圆柱第 列圆柱随着尾流的复杂化，其升力系数开始较大幅度的振荡，且方向会随时间改变，说明紊动强度增加，能量损失更大。当 时，由图（）（）可见，规律和雷诺数 时相似，升力系数的振荡幅度较之雷诺数 时更大，流体紊动强度继续增加。对比 时的速度云图和实验装置的流动图像（见图），可发现两者较好符合，认为仿真结果与实际流动是很接近的，可用数值模
16、拟的方法替代实验装置来演示紊流的流态。（）（）图 实验装置流动图像对比 结 语本文对实验装置简化模型在 种雷诺数下进行了数值计算和流态、阻力系数和升力系数分析，可知，多圆柱绕流在雷诺数 时的速度云图和实验装置的流动图像符合，认为仿真结果与实际流动是很接近的，可采用模拟的方法去演示流体的紊动机理，为流体力学虚拟仿真开发提供有益的探索。参考文献（）：，（）：陈薛浩，朱志夏 基于 的海底管线附近流场分析 上海交通大学学报，（）：何国建，方红卫，府仁寿 桥墩群对河道水流影响的三维数值分析 水动力学研究与进展：辑，（）：，（）：彭宇清，应旭永 高层建筑典型断面绕流的非定常数值模拟工程建设与设计，（）：武
17、玉涛，任华堂，夏建新 圆柱绕流研究进展及展望 水运工程，（）：，：，（）：，（）：，（）：，（）：张爱社，张 陵 等边布置三圆柱绕流的数值分析 应用力学学报，（）：，（）：张洪军，苏中地，尹招琴 并列双圆柱尾迹被动控制中流场特性的实验研究 水动力学研究与进展：辑，（）：，：，（）：史里希廷 边界层理论 徐燕侯等译 北京：科学出版社，顾 罡 二维单圆柱、双圆柱绕流问题和三维垂荡板运动的数值模拟 上海：上海交通大学，（上接第 页）陈聚祥，黎映静 基于数据包络分析的福建三级甲等医院运营效率评价 中国卫生统计，（）：吴锦涛，董 刚，栗保明 基于多线程并行的动态自适应建表加速算法在气相爆轰模拟中的应用
18、兵工学报，（）：王宇宙，徐忠伟，程嗣怡，等 基于包络信息的 信号与侦收信道带宽定量关系分析 现代雷达，（）：白广争，李 亮，白广良 基于通信的列车控制系统数据采集与解析软件的设计与实现 城市轨道交通研究，（）：周佳佳，李 涛，黄小康 多核同时多线程处理器的线程调度器设计 电子技术应用，（）：丁小磊，宋 俐，沈文琪 基于数据包络分析法的江苏省基本公共卫生服务实施效率评价 中国全科医学，（）：贾 康，郑 帅，荆冲，等 基于离散曲线包络的刮削加工仿真数学模型 计算机集成制造系统，（）：黄晓霞，吴 燕，王君燕，等 基于数据包络分析法的上海市闵行区社区卫生服务机构公共卫生服务效率评价 中国全科医学，（）：田 佳，胡 威，李 琳，等 基于粒子群优化算法的多核处理器任务 调度 研 究 计 算 机 应 用 研 究，（）：陈 俊，周宽久，贾敏 多线程并行程序数据竞争静态检测方法 计算机工程与设计，（）：王 静，王春梅，智 佳，等 面向有效载荷高速数据流的数据处理方法 计算机工程与设计，（）：
