三维圆管流动状况的数值模拟分析.docxVIP

三维圆管流动状况的数值模拟分析 三维圆管流动状况的数值模拟分析 三维圆管流动状况的数值模拟分析 三维圆管流动状况的数值模拟分析 在工程和生活中，圆管内的流动是最常见也是最简单的一种流动，圆管流动有层流和紊流两种流动状 况。层流，即液体质点作有序的线状运动，彼此互不混掺的流动；紊流，即液体质点流动的轨迹极为紊乱， 质点相互掺混、碰撞的流动。雷诺数是判别流体流动状态的准则数。本研究用CFD软件来模拟研究三维 圆管的层流和紊流流动状况，主要对流速分布和压强分布作出分析。 物理模型 三维圆管长l2000mm，直径d100mm。 流体介质：水，其运动粘度系数 1 106m2/s。 Inlet：流速入口，1 /s，2 0.1m/s Outlet：压强出口 Wall：光滑壁面，无滑移 2在ICEMCFD中建立模型 2.1首先建立三维圆管的几何模型Geometry 2.2做Blocking 因为截面为圆形，故需做“O”型网格。 2.3划分网格mesh 注意检查网格质量。 在未加密的情况下，网格质量不是很好，如下图 因管流存在边界层，故需对边界进行加密，网格质量有所提升，如下图 2.4生成非结构化网格，输出fluent.msh等相关文件 数值模拟原理 3.1层流流动 1 /s，从Inlet υd 当水流以流速 方向流入圆管，可计算出雷诺数 Re 500 ，故圆管内 ν 流动为层流。 假设水的粘性为常数（运动粘度系数 1 106m2/s）、不可压流体，圆管光滑，则流动的控制方程 如下： ①质量守恒方程： ( u) ( v) (w) 0 (1-1) t x y z ②动量守恒方程： ( u) ( uu) ( uv) ( uw) ( u) ( u) ( u) p (1-2) t x y z x x y y z z x ( v) ( vu) ( vv) ( vw) ( v) ( v) ( v) p (1-3) t x y z x x y y z z y (w) ( wu) ( wv) ( ww) ( w ( w ( w p t x y z ) ) ) (1-4) x x y y z z z 式中， 为密度，u、 、w是流速矢量在 x、y和z方向的分量，p为流体微元体上的压强。 方程求解：对于细长管流， FLUENT建议选用双精度求解器，流场计算采用 SIMPLE算法，属于压强 修正法的一种。 3.2紊流流动 υd 10000 当以水流以流速 2 0.1m/s，从Inlet方向流入圆管，可计算出雷诺数Re ，故圆管内 ν 流动为紊流。 假设水的粘性为常数（运动粘度系数 1 106m2/s）、不可压流体，圆管光滑，则流动的控制方程 如下： ①质量守恒方程： ( u) ( v) (w) 0 (1-5) t x y z ②动量守恒方程： (u) ( uu) (uv) ( uw) ( u) ( u) ( u) t x y z x x y y z z (1-6) (u2) (uv) (uw)] p [ x y z x ( v) ( vu) (vv) ( vw) ( v) ( v) ( v) t x y z x x y y z z (1-7) (uv) (v2) (vw)] p [ x y z y ( w) ( wu) ( wv) ( ww) x ( w) y ( w) ( w) t x y z x y z z (uw) (vw) (w2)] p (1-8) [ x y z z ③湍动能方程： ( k) ( ku) ( kv) ( kw) [( t ) k)] [( t ) k)] t x y z x k x y k y k)] (1-9) [( t ) Gk z k z ④湍能耗散率方程： ( ) ( u) ( v) ( w) [( t ) )] [( t ) )] t x y z x k x y k y C1 2 (1-10) t [( ) )] Gk C2 z k z k k 式中， 为密度，u、 、w是流速矢量在 x、y和z方向的分量，p为流体微元体上的压强。 方程求解：采用双精度求解器，定常流动，标准 k 模型，SIMPLEC算法。 4在FLUENT中求解计算层流流动 4.1导入并检查网格 注意调整Scale大小。因在ICEM中作网格时，已采用的是以“米”为单位的长度，故不需更换单位。 网格显示流动沿X方向，共存在283575hexahedralcells，范围DomainExtents: x-coordinate:min(m)=0.000000e+000,max(m)=2.000000e+000 4.2设置求解器 本模型基于压强计算，可采取绝对流速计算，Solver求解器可采取默认设置。 雷诺数Re υd 500，故圆管内流动为层流， Viscous设置