FLUENT算例 三维圆管紊流流动状况的数值模拟分析.doc

三维圆管紊流流动状况的数值模拟分析 在工程和生活中，圆管内的流动是最常见也是最简单的一种流动，圆管流动有层流和紊流两种流动状况。层流，即液体质点作有序的线状运动，彼此互不混掺的流动；紊流，即液体质点流动的轨迹极为紊乱，质点相互掺混、碰撞的流动。雷诺数是判别流体流动状态的准则数。本研究用CFD软件来模拟研究三维圆管的紊流流动状况，主要对流速分布和压强分布作出分析。 1 物理模型 三维圆管长，直径。 流体介质：水，其运动粘度系数。 Inlet：流速入口，， Outlet：压强出口 Wall：光滑壁面，无滑移 2 在ICEM CFD中建立模型 2.1 首先建立三维圆管的几何模型Geometry 2.2 做Blocking 因为截面为圆形，故需做“O”型网格。 2.3 划分网格mesh 注意检查网格质量。 在未加密的情况下，网格质量不是很好，如下图 因管流存在边界层，故需对边界进行加密，网格质量有所提升，如下图 2.4 生成非结构化网格，输出fluent.msh等相关文件 3 数值模拟原理 紊流流动 当以水流以流速，从Inlet方向流入圆管，可计算出雷诺数，故圆管内流动为紊流。 假设水的粘性为常数（运动粘度系数）、不可压流体，圆管光滑，则流动的控制方程如下： ①质量守恒方程： ②动量守恒方程： ③湍动能方程： ④湍能耗散率方程： 式中，为密度，、、是流速矢量在x、y和z方向的分量，p为流体微元体上的压强。 方程求解：采用双精度求解器，定常流动，标准模型，SIMPLEC算法。 4 在FLUENT中求解计算紊流流动 4.1 FLUENT设置 除以下设置为紊流所必须设置的外，其余选项和层流相同，不再详述。 ①Viscous设置 雷诺数，故圆管内流动为紊流，Viscous设置为Realizable K-epsilon模型，其余默认。 ②Boundary设置 Inlet设置为速度入口，为，Turbulence设置为Intensity and Hydraulic Diameter方法，即 Outlet设置为自由出口Outflow，如设置成压力出口，则之后计算会存在问题（已验证）。 ③Solution设置 采用双精度求解器，定常流动，Realizable 模型，SIMPLEC算法。 4.2 开始迭代 设置迭代次数为300，实际比这个少，迭代收敛时会自动停止。 5 紊流计算结果及分析 计算293步后，已收敛，自动停止运算。 残差监视窗口为 5.1 显示流速等值线图 5.1.1 入口和出口截面的流速分布图 分布在Surface里选择inlet及outlet （1）Velocity of inlet 可见，入口处流速分布不明显，基乎都等于入口流速，只是外层靠近壁面处流速几乎为零。 （2）Velocity of outlet 可见，出口截面流速分布较为明显，和层流一样，显同心圆分布，内层流速偏大，外层靠近壁面处流速几乎为零，分层更为严重，边界层很薄。 5.1.2 Y轴和Z轴方向流速截面 圆管内各个截面的流速分布均不相同，可以认为紊流还没达到稳定状态，在此不再分析各个截面的流速分布，仅对整个圆管的流速作出分析。 截面沿圆管长度X方向截取，可看到对称的效果。 （1）Velocity of y-0 整根圆管： （2）Velocity of z-0 整根圆管： 以上两个截面流速分布图的效果是一样的，可以看出圆管水流紊流入口段及之后的流速发展趋势，而且显示流速变化的规律更为明显。 （3）入口段 与层流入口段的流速分布相比，可以明显的看出紊流入口段的流速分布不太明显，且基本没有分层，符合紊流流动的基本规律。流速分布也不像层流流速那样显明显抛物线分布，而是更加平滑，越超后发展发展越平滑，到底是什么曲面，之后再加分析。紊流过流断面的流速对数分布比层流的抛物面分布均匀的多，符合的规律，即 （4）出口段 出口段的流层分布很明显，切趋于均匀，但仔细观察圆管轴心的速度，其实速度分布并未达到均匀，可见紊流并未达到充分发展的状况。 5.1.3轴向流速的变化 执行Plot→XY Plot，选择Y Axis Function里的Velocity和Velocity Magnitude，选择Surfaces里圆管的对称轴line-x，可得到轴向流速分布散点图。 由上图可以看出，在圆管的轴上，进口段流速分布变化较大，从进口流速急剧上升到最大流速。之后又下降。但实际经验表明，紊流应该在进口段后达到稳定状态，轴向流速应该趋于恒定，可见此模拟实验设置长度不够，使流动并未达到充分紊流。 紊流入口段长度有经验公式可以算的，即 由此可见，紊流的边界层厚度的增长比层流边界层要快，因此紊流的进口段要