圆柱绕流的数值模拟研究.docxVIP

圆柱绕流的数值模拟研究 摘要：选取直径为 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所示的计算域建立了 二维的计算模型。 小 ￡5D 1 圭o i 1 =3D J 甲 图1计算区域 1.2网格划分及边界条件设置 为提高模拟精度，计算区域采用分块网格划分与结构化网格划分相结合的技术。 计算区 域共分两块，尺寸见图 1所示。在圆柱区域采用 0型结构化网格（图2），尾流区域采用四边 形结构化网格分别划分（图 3）,使用GAMBIT对两块计算区域进行了网格划分，划分的结 果是网格总数为42946个。 对计算区域进行边界条件定义， 考虑到流入介质的为空气， 同时流速较低，就把介质假 定为不可压缩的流体。进而把左侧的入口定义为速度入口即： Velocity-inlet，右侧的出口假 定为充分发展的出流，即定义为： Outflow。其余的边界保持默认的壁面边界条件，同时定 义为绝热条件，即热流密度为 0。 图2圆柱区域0型结构化网格 图3整体网格划分情况 1.3计算方法 流动介质为空气，密度 p=1.225kg/m3，动力粘度 尸1.79 xi0-5pa?j，并保持为常数。计算 雷诺数（Re2000）均为小雷诺数。属层流范围，故采用 Laminar模型。采用基于压力基的分 离式求解器进行求解。计算中采用具有二阶隐式时间格式的非定常流动进行计算。 压力项与 速度项的耦合项计算采用 SIMPLE算法实现，压力项离散采用具有二阶精度的格式离散. 动 量方程采用二阶迎风格式离散。计算中压力、密度、体积力和动量项的欠松弛因子分别为 0.3、1、1 和 0.7。 2数值模拟 从0.03m/s开始逐渐增加入口初速度来对不同选用25C的空气为流经计算区域的流体, 速度下的圆柱绕流的流场进行模拟 从0.03m/s开始逐渐增加入口初速度来对不同 2.1.1涡量云图 0.03m/s速度下的涡量云图0.1m/s速度下的涡量云图0.2m/s速度下的涡量云图0.3m/s速度下的涡量云图0.5m/s速度下的涡量云图1m/s速度下的涡量云图 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时，整个计算域的涡量云图呈对称分布，在圆柱的后侧没 有明显的漩涡存在，空气流贴着圆柱流动，在贴近圆柱的上下侧涡量较其他区域大。速度 为0.1m/s时漩涡开始向圆柱的后侧推进，此时涡量较大的区域集中在圆柱的前半侧，在圆 柱右侧的中后部区域涡量图出现明显的远离现象，贴近圆柱的地方涡量很低，而外扩展的 话就又会变大，圆柱的后部区域出现比较明显的小漩涡。在入口速度达到 0.2m/s时圆柱后 侧的漩涡就比较的明显了，漩涡分布在圆柱的后侧离圆柱较近的区域，上下对称分布。 入口速度为0.3m/s时圆柱后侧的漩涡更加的明显，并且有细微的扩大趋势。当速度达 到0.5m/s时圆柱后侧的漩涡范围明显