圆柱绕流数值模拟.docx

摘要：选取直径为D=10mm的圆柱及6Dx3D的计算区域，利用GAMBIT进行模型的 创建模型，对计算区域采用分块网格划分与结构化网格划分相结合的技术进行网格划分对 0.03m/s〜1.0m/s的低流速情况下的圆柱绕流进行模拟研究，结果发现在速度达到0.1m/s前圆 柱后侧没有出现明显的漩涡，在速度大于0.1m/s后漩涡开始出现，当速度达到0.5m/s时漩 涡的范围最大最后利用FLUENT的网格自适应技术对入口速度为0.5m/s的情况进行了网 格加密，发现网格自动加密可以改进网格分布情况，但对计算结果的影响程度有限关键词：网格划分；圆柱绕流；涡量；网格自适应钝体绕流中尤其以圆柱体的绕流问题最为经典和引起人们的注意圆柱绕流属于非定常 分离流动问题，在工业工程中的应用非常广泛圆柱绕流同时也是一个经典的流体力学问题， 流体绕圆柱体流动时，过流断面收缩，流速沿程增加，压强沿程减小，由于黏性力的存在， 就会在柱体周围形成附面层的分离，形成圆柱绕流而由于圆柱的存在，会在圆柱迎水面产 生壅水现象，同时也增加了圆柱的受力，使得圆柱绕流问题变得十分复杂研究圆柱绕流问题在工程实际中也具有很重要的意义如在水流对桥梁、海洋钻井平台 支柱、海底输运管线、桩基码头等的作用中，风对塔建筑、化工塔设备、高空电缆等的作用 中，都有重要的工程应用背景。

因此，对圆柱绕流进行深入研究，了解其流动机理和水动力 学规律，不仅具有理论意义，还具有明显的社会经济效益1数学模型与计算方法1.1几何模型结合本文研究目标，取圆柱直径D=10mm,计算区域为6DX3D的矩形区域，如图1所 示上游尺寸1.5D，下游尺寸4.5D使用GAMBIT建模软件按照图1所示的计算域建立了 二维的计算模型4-5D -■11壬O i ；图1计算区域1.2网格划分及边界条件设置为提高模拟精度，计算区域采用分块网格划分与结构化网格划分相结合的技术计算区 域共分两块，尺寸见图1所示在圆柱区域采用O型结构化网格（图2）,尾流区域采用四边 形结构化网格分别划分（图3），使用GAMBIT对两块计算区域进行了网格划分，划分的结 果是网格总数为42946个对计算区域进行边界条件定义，考虑到流入介质的为空气，同时流速较低，就把介质假 定为不可压缩的流体进而把左侧的入口定义为速度入口即：Velocity-inlet,右侧的出口假 定为充分发展的出流，即定义为：Outflow其余的边界保持默认的壁面边界条件，同时定 义为绝热条件，即热流密度为0图2圆柱区域O型结构化网格 图3整体网格划分情况1.3计算方法流动介质为空气，密度p=1.225kg/m3，动力粘度卩=1.79xlO-5Pa・s，并保持为常数。

计算 雷诺数(Rev2000)均为小雷诺数属层流范围，故采用Laminar模型采用基于压力基的分 离式求解器进行求解计算中采用具有二阶隐式时间格式的非定常流动进行计算压力项与 速度项的耦合项计算采用SIMPLE算法实现，压力项离散采用具有二阶精度的格式离散.动 量方程采用二阶迎风格式离散计算中压力、密度、体积力和动量项的欠松弛因子分别为0.3、1、1 和 0.72数值模拟选用25°C的空气为流经计算区域的流体，从0.03m/s开始逐渐增加入口初速度来对不同 速度下的圆柱绕流的流场进行模拟2.1结果与分析2.1・1涡量云图0.03m/s速度下的涡量云图0.1m/s速度下的涡量云图0.2m/s速度下的涡量云图0.3m/s速度下的涡量云图0.5m/s速度下的涡量云图1m/s速度下的涡量云图图4 0.03m/s~1m/s下的涡量图由上图可知在速度为0.03m/s时，整个计算域的涡量云图呈对称分布，在圆柱的后侧没 有明显的漩涡存在，空气流贴着圆柱流动，在贴近圆柱的上下侧涡量较其他区域大速度为O.lm/s时漩涡开始向圆柱的后侧推进，此时涡量较大的区域集中在圆柱的前半侧，在圆 柱右侧的中后部区域涡量图出现明显的远离现象，贴近圆柱的地方涡量很低，而外扩展的 话就又会变大，圆柱的后部区域出现比较明显的小漩涡。

在入口速度达到0.2m/s时圆柱后 侧的漩涡就比较的明显了，漩涡分布在圆柱的后侧离圆柱较近的区域，上下对称分布入口速度为0.3m/s时圆柱后侧的漩涡更加的明显，并且有细微的扩大趋势当速度达 到0.5m/s时圆柱后侧的漩涡范围明显的扩大，扩展到圆柱后侧很大区域，区域长度和圆柱 直径相当，并且在圆柱后侧紧近圆柱的部分也出现较小的漩涡，但是漩涡的整体强度却没 有明显的增大从涡量云图上看在速度达到lm/s时漩涡的范围有所减少，但是强度确明显 的得到了增强，在圆柱的后侧有两个非常明显的漩涡存在，漩涡的影响范围涉及到了紧近 圆柱后侧的区域再则是由于模拟区域的上下边界都定义为了壁面，从上面6幅涡量云图 可以看见入口区域，贴近壁面的范围也出现了轻微的漩涡，并且在速度为0.5m/s时漩涡最 明显速度云图0.03m/s初速下的速度云图 0.1 m/s初速下的速度云图0.2 m/s初速下的速度云图0.3 m/s初速下的速度云图0.5 m/s初速下的速度云图 1 m/s初速下的速度云图图5 0.03m/s〜lm/s下的速度云图从上图可以看到速度云图总体呈对称分布，在入口速度为0.03m/s时高速区域集中在圆 柱的上下两侧，在圆柱的前后是速度的低速区，并且有两处速度很小的区域。

当速度达到 0.1m/s时圆柱后侧的低速区域得到扩大，低速区域的范围呈现箭头状，越往外速度的值越来 越大在入口速度为0.2m/s的条件下，可以从图上圆柱后也有明显的低速区域，但是速度 明显大于0.1m/s时的速度同时可以看到圆柱后侧有三个点状的低速区，围绕着一个速度 较高的椭圆区域从0.3m/s的云图可以看到圆柱后方有三个明显的低速区，在这三个低速区中后面的面 积较大，而前面的两个低速区较小，并且三个低速区域有连在一起的趋势到了 0.5m/s时 我们可以清楚的看到圆柱后侧的低速区域连在了一起，整体呈现出狭长的v字形状，同时包 裹着一个近椭圆的区域，在椭圆中速度从中间往外逐渐变小，可以看到有一块明显的黄色区域同时由于入口速度的提高圆柱两侧的高速区域较前面有了很大的扩展，两侧呈狭长的矩 形分布，并且在贴近计算域边界的壁面处也出现了明显的条状的低速区但是当速度达到 1m/s时圆柱上下两侧的高速区域宽度有所增加而长度却明显减少，同时圆柱后侧的低速区 域也减少很多，三个低速区也明显的分离开来3网格自适应技术的应用对初速度为0.5m/s的情况应用FLUENT的网格自适应技术对速度梯度大于1.0xl0-6n/s 的地方进行了网格加密，网格前后效果对比如下图：未加密的网格 自动加密后的网格图6网格效果对比图上图仅仅截取了圆柱周围的部分网格图，对比两个网格可以清楚的看到在后图里在贴近圆柱的区域里网格明显的密于前图。

同时在圆柱的后部区域也有三个条状的加密区域加密前的速度云图 加密后的速度云图对比加密前后的速度云图，可以发现圆柱后部上侧区域的低速区域有点不同，同时后图 的圆柱后的中轴线上的速度较前图有所降低，更能反映实际的情况4结论使用FLUENT软件对0.03m/s~1m/s间的入口速度下的定常流动的圆柱绕流进行了数值 模拟的对比研究结果发现：（1） 在入口速度小于O.lm/s时，圆柱后侧没有出现漩涡，在速度为0.1m/s时圆柱后侧 开始有不明显的漩涡出现当速度大于0.1m/s后就会出现明显的漩涡，同时在0.5m/s的条 件下漩涡的范围最大，在1m/s时漩涡强度大于0.5m/s的情况，但范围却小于前者2） 使用FLUENT的网格自适应技术对入口速度为0.5m/s的情况依据速度梯度大 小对网格进行了加密，结果发现虽然可以对高梯度的区域进行网格加密，但是再次计算 的结果却没有达到预想的效果，不过总的来讲网格自适应技术也加密了网格增进了模拟 的精度3） 由于仅仅对速度小于1m/s的情况进行了对比研究，对高速度高雷诺数的情况 未有涉及，结果具有局限性参考文献：[1] 夏雪渝，等.工程分离流动力学[M].北京：北京航空航天大学出版社，1991.[2] 孟元元.圆柱绕流的数值模拟研究[D]•甘肃农业大学，2010.[3] 张立.小雷诺数下圆柱绕流的数值模拟[J].力学季刊,2010,04:543-547.[4] 杨纪伟，付晓丽.圆柱绕流研究进展[J].中国水运（下半月），2008，05:156-158.[5] 张立.小雷诺数下圆柱绕流的涡态演化[J].机械科学与技术,2012,04:679-684.⑹李国亚.有限水深横流中近壁水平圆柱绕流的实验研究[D].武汉大学,2004.[7] 段志强•低雷诺数下尾部隔板影响圆柱绕流的数值研究[D] •重庆大学,2012.[8] 何鸿涛.圆柱绕流及其控制的数值模拟研究[D].北京交通大学,2009.。