FLUENT算例三维圆管紊流流动状况的数值模拟解析总结计划.docxVIP

FLUENT算例三维圆管紊流流动状况的数值模拟解析总结计划 FLUENT算例三维圆管紊流流动状况的数值模拟解析总结计划 PAGE PAGE16 FLUENT算例三维圆管紊流流动状况的数值模拟解析总结计划 PAGE 三维圆管紊流流动情况的数值模拟解析 在工程和生活中，圆管内的流动是最常有也是最简单的一种流动，圆管流动有层流和紊流两种流动状 况。层流，即液体质点作有序的线状运动， 相互互不混掺的流动； 紊流，即液体质点流动的轨迹极为凌乱， 质点相互掺混、碰撞的流动。雷诺数是鉴识流体流动状态的准则数。本研究用 CFD软件来模拟研究三维 圆管的紊流流动情况，主要对流速分布和压强分布作出解析。 物理模型 三维圆管长l 2000mm，直径d 100mm。 流体介质：水，其运动粘度系数 1 106m2/s。 Inlet：流速入口，1 /s，2 /s Outlet：压强出口 Wall：圆滑壁面，无滑移 2在ICEMCFD 中建立模型 第一建立三维圆管的几何模型 Geometry 做Blocking 由于截面为圆形，故需做“ O”型网格。 划分网格 mesh 注意检查网格质量。 在未加密的情况下，网格质量不是很好，以以下图 因管流存在界线层，故需对界线进行加密，网格质量有所提升，以以下图 生成非结构化网格，输出 等相关文件 数值模拟原理 紊流流动 当以水流以流速 2 /s，从Inlet 方向流入圆管，可计算出雷诺数 Re υd 10000，故圆管内 ν 流动为紊流。 假设水的粘性为常数（运动粘度系数 1 106m2/s）、不可以压流体，圆管圆滑，则流动的控制方程 以下： ①质量守恒方程： ( u) ( v) ( w) 0 (0-1) t x y z ②动量守恒方程： ( u) ( uu) ( uv) ( uw) ( u) ( u) ( u) t x y z x x y y z z (u2) (uv) (uw) p (0-2) [ x y z ] x ( v) ( vu) ( vv) ( vw) ( v) ( v) ( v) t x y z x x y y z z (uv)(v2) (vw) p (0-3) [ x y z ] y ( w) ( wu) ( wv) ( ww) ( w) ( w) ( w) t x y z x x y y z z (uw) (vw) (w2) p (0-4) [ x y z ] z ③湍动能方程： ( k) ( ku) (kv) (kw) x [( t ) k)] [( t ) k)] t x y z k x y k y k)] (0-5) [( t) Gk z k z ④湍能耗散率方程： ( ) ( u) ( v) ( w) [( t ) )] [( t ) )] t x y z x k x y k y C1 2 (0-6) t) [( )] G k C z z k 2 k k 式中， 为密度，u、 、w是流速矢量在 x、y和z方向的重量，p为流体微元体上的压强。 方程求解：采用双精度求解器，定常流动，标准 k 模型，SIMPLEC算法。 4在FLUENT 中求解计算紊流流动 FLUENT 设置 除以下设置为紊流所必定设置的外，其他选项和层流相同，不再详述。 Viscous设置 雷诺数Re υd 10000，故圆管内流动为紊流，Viscous设置为RealizableK-epsilon模型，其他默认。 ν Boundary设置 Inlet设置为速度入口，为 2 /s，Turbulence设置为IntensityandHydraulicDiameter 方法，即 Outlet设置为自由出口 Outflow，如设置成压力出口，则此后计算会存在问题（已考据） 。 Solution设置 采用双精度求解器，定常流动， Realizable k 模型，SIMPLEC算法。 开始迭代 设置迭代次数为 300，本质比这个少，迭代收敛时会自动停止。 紊流计算结果及解析 计算293步后，已收敛，自动停止运算。 残差监察窗口为 显示流速等值线图 入口和出口截面的流速分布图 分布在Surface里选择inlet及outlet 1）Velocityofinlet 可见，入口处流速分布不明显，基乎都等于入口流速 2 /s，可是外层凑近壁面处流速几乎为 零。 2）Velocityofoutlet 可见，出口截面流速分布较为明显，和层流相同，显同心圆分布，内层流速偏大，外层凑近壁面处流速几乎为零，分层更加严重，界线层很薄。 5.1.2 Y轴和Z轴方向流速截面 圆管内各个截面的流速分布均不相同，可以认为紊流还没达到牢固状态，在此不再解析各个截面的流速分布，仅对整个圆管的流速作出解析。 截面沿圆管长度 X方向截取，可看到对称的收效。 1）Veloc