圆柱绕流.docxVIP

二维圆柱绕流 摘要:采用有限体积法求解二维N-S方程，对雷诺数1,10,100的二维圆柱非定常流场进行了数值模拟，对比各雷诺数下其流动情况发现，在Re=1时，圆柱上下游的流线前后对称，此Re数范围的绕流称为斯托克斯区；随着Re的增大，圆柱上下游的流线逐渐失去对称性。当Re=10时，沿圆柱表面流动的流体在到达圆柱顶点（90度）附近就离开了壁面，分离后的流体在圆柱下游形成一对固定不动的对称漩涡（附着涡），涡内流体自成封闭回路而成为“死水区”；随着Re的增大，死水区逐渐拉长圆柱前后流场的非对称性逐渐明显，此Re数范围称为对称尾流区。圆柱下游流场不再是定常的，圆柱后缘上下两侧有涡周期性地轮流脱落，形成规则排列的涡阵，这种涡阵称为卡门涡街。 圆柱绕流研究 圆柱绕流是一个经典的流体力学问题，流体绕圆柱体流动时，过流断面收缩，流速沿程增加，压强沿程减小，由于粘性力的存在，，就会在柱体周围发生边界层的分离，形成圆柱绕流。圆柱绕流现象比较复杂，因此，对圆柱绕流研究具有重要的基础理论意义。 研究圆柱绕流问题在工程实际中也具有非常重要的意义。如水流对桥梁、海上钻井平台支柱、海上输运管线、桩基码头等的作用中，风对塔设备、化工塔设备，高空电缆等的作用中，都有着重要的工程应用背景。因此，对圆柱绕流进行深入研究，对其流动机理进行分析，不仅具有理论意义，还有明显的社会经济效益。 数值方法 因为本文主要求解雷诺数Re=1,10,100时的圆柱绕流情况，需要求解二维非定常不可压的N—S方程组：本文采用有限体积法对上述微分方程进行离散，然后用SIMPLE算法对离散方程进行求解，计算中时间推进采用二阶隐式格式，空间离散采用三阶精度的QUICK格式。控制方程如下： (1) (2) 网格划分及模拟工况3.1计算网格 计算的区域大小为上游边界距圆柱圆心为2.5D，下游边界距圆柱圆心10D，顶部和底部边界距圆柱圆心2.5D，如图1所示。流场的网格采用结构网格，其中圆柱壁面网格采用O型网格如图2所示。 Outlet 10D 2.5D D 2.5D Inlet 图1. 计算模型 图2. 网格划分 数值计算边界条件：人口为速度入口，出口是压力出口，圆柱及上下表面均为壁面无滑移条件。3.2模拟工况 雷洛数Re反映了惯性力与粘性力的比值： Re=pulμ=ulν （3）Re，特征尺寸D，运动粘度ν如表1所示： 表1. 进口流速 Re特征直径D(mm)运动粘度ν（m2/s）流速v(m/s)11010-610-51010-410010-3结果分析雷洛数分别为1、10、100时的流线图与涡量图，如下所示： 图a. 速度矢量图 图b. 涡量云图 图3. Re=1时流线图和涡量图 图a. 速度矢量图 图b. 涡量云图 图4. Re=10时流线图和涡量图 图a. 速度矢量图 图b. 涡量云图 图5. Re=100时流线图和涡量图 5. 结论 6. 参考文献