文档介绍：圆柱绕流的数值模拟研究摘要: 选取直径为 D=10mm 的圆柱及 6D × 3D 的计算区域,利用 GAMBIT 进行模型的创建模型, 对计算区域采用分块网格划分与结构化网格划分相结合的技术进行网格划分。对 ~ 的低流速情况下的圆柱绕流进行模拟研究, 结果发现在速度达到前圆柱后侧没有出现明显的漩涡,在速度大于后漩涡开始出现,当速度达到时漩涡的范围最大。最后利用 FLUENT 的网格自适应技术对入口速度为的情况进行了网格加密,发现网格自动加密可以改进网格分布情况,但对计算结果的影响程度有限。关键词: 网格划分;圆柱绕流;涡量;网格自适应钝体绕流中尤其以圆柱体的绕流问题最为经典和引起人们的注意。圆柱绕流属于非定常分离流动问题, 在工业工程中的应用非常广泛。圆柱绕流同时也是一个经典的流体力学问题, 流体绕圆柱体流动时,过流断面收缩,流速沿程增加,压强沿程减小,由于黏性力的存在, 就会在柱体周围形成附面层的分离, 形成圆柱绕流。而由于圆柱的存在, 会在圆柱迎水面产生壅水现象,同时也增加了圆柱的受力,使得圆柱绕流问题变得十分复杂。研究圆柱绕流问题在工程实际中也具有很重要的意义。如在水流对桥梁、海洋钻井平台支柱、海底输运管线、桩基码头等的作用中, 风对塔建筑、化工塔设备、高空电缆等的作用中, 都有重要的工程应用背景。因此, 对圆柱绕流进行深入研究, 了解其流动机理和水动力学规律,不仅具有理论意义,还具有明显的社会经济效益。 1 数学模型与计算方法几何模型结合本文研究目标, 取圆柱直径 D=10mm , 计算区域为 6D × 3D 的矩形区域, 如图 1所示。上游尺寸 , 下游尺寸。使用 GAMBIT 建模软件按照图 1 所示的计算域建立了二维的计算模型。图1 计算区域网格划分及边界条件设置为提高模拟精度, 计算区域采用分块网格划分与结构化网格划分相结合的技术。计算区域共分两块, 尺寸见图 1 所示。在圆柱区域采用 O 型结构化网格(图 2), 尾流区域采用四边形结构化网格分别划分(图3), 使用 GAMBIT 对两块计算区域进行了网格划分, 划分的结果是网格总数为 42946 个。对计算区域进行边界条件定义, 考虑到流入介质的为空气, 同时流速较低, 就把介质假定为不可压缩的流体。进而把左侧的入口定义为速度入口即: Velocity-inlet ,右侧的出口假定为充分发展的出流,即定义为: Outflow 。其余的边界保持默认的壁面边界条件,同时定义为绝热条件,即热流密度为 0。图2 圆柱区域 O 型结构化网格图3 整体网格划分情况计算方法流动介质为空气,密度ρ= 3 ,动力粘度μ= × 10 -5 Pa?s ,并保持为常数。计算雷诺数(Re < 2000) 均为小雷诺数。属层流范围,故采用 Laminar 模型。采用基于压力基的分离式求解器进行求解。计算中采用具有二阶隐式时间格式的非定常流动进行计算。压力项与速度项的耦合项计算采用 SIMPLE 算法实现,压力项离散采用具有二阶精度的格式离散. 动量方程采用二阶迎风格式离散。计算中压力、密度、体积力和动量项的欠松弛因子分别为、1、1和