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

三维圆管流动状况的数值模拟分析 在工程和生活中，圆管内的流动是最常见也是最简单的一种流动，圆管流动有层流和紊流两种流动状况。层流，即液体质点作有序的线状运动，彼此互不混掺的流动；紊流，即液体质点流动的轨迹极为紊乱，质点相互掺混、碰撞的流动。雷诺数是判别流体流动状态的准则数。本研究用CFD软件来模拟研究三维圆管的层流和紊流流动状况，主要对流速分布和压强分布作出分析。 1 物理模型 三维圆管长，直径。 流体介质：水，其运动粘度系数。 Inlet：流速入口，， Outlet：压强出口 Wall：光滑壁面，无滑移 2 在ICEM CFD中建立模型 2.1 首先建立三维圆管的几何模型Geometry 2.2 做Blocking 因为截面为圆形，故需做“O”型网格。 2.3 划分网格mesh 注意检查网格质量。 在未加密的情况下，网格质量不是很好，如下图 因管流存在边界层，故需对边界进行加密，网格质量有所提升，如下图 2.4 生成非结构化网格，输出fluent.msh等相关文件 3 数值模拟原理 3.1 层流流动 当水流以流速，从Inlet方向流入圆管，可计算出雷诺数，故圆管内流动为层流。 假设水的粘性为常数（运动粘度系数）、不可压流体，圆管光滑，则流动的控制方程如下： ①质量守恒方程： ②动量守恒方程： 式中，为密度，、、是流速矢量在x、y和z方向的分量，p为流体微元体上的压强。 方程求解：对于细长管流，FLUENT建议选用双精度求解器，流场计算采用SIMPLE算法，属于压强修正法的一种。 3.2 紊流流动 当以水流以流速，从Inlet方向流入圆管，可计算出雷诺数，故圆管内流动为紊流。 假设水的粘性为常数（运动粘度系数）、不可压流体，圆管光滑，则流动的控制方程如下： ①质量守恒方程： ②动量守恒方程： ③湍动能方程： ④湍能耗散率方程： 式中，为密度，、、是流速矢量在x、y和z方向的分量，p为流体微元体上的压强。 方程求解：采用双精度求解器，定常流动，标准模型，SIMPLEC算法。 4 在FLUENT中求解计算层流流动 4.1 导入并检查网格 注意调整Scale大小。因在ICEM中作网格时，已采用的是以“米”为单位的长度，故不需更换单位。 网格显示流动沿X方向，共存在283575 hexahedral cells，范围Domain Extents: x-coordinate: min (m) = 0.000000e+000, max (m) = 2.000000e+000 y-coordinate: min (m) = -4.995393e-002, max (m) = 4.995393e-002 z-coordinate: min (m) = -4.995393e-002, max (m) = 4.995393e-002 4.2 设置求解器 本模型基于压强计算，可采取绝对流速计算，Solver求解器可采取默认设置。 雷诺数，故圆管内流动为层流，Viscous设置为Laminar。 4.3 定义材料 因本研究采用水流动，故需使Material type定义为Fluent，设置成水。 4.4 设置边界条件 4.4.1 将Solid定义为Fluid，并设置成水 4.4.2 定义进口Inlet Inlet定义为流速入口Velocity-inlet，并设置入口流速为0.005m/s。 4.4.3 定义出口Outlet 出口为压强出口Pressure Outlet，默认设置。 4.4.4 定义壁面Wall。 设置为默认。 4.5 设置操作条件 因为圆管截面较小，故可不考虑重力选项。压强选项默认为一个大气压。 4.6 求解方法的设置与控制 4.6.1 求解参数的设置、 在Solution controls中，将Momentum设置为Second order upwind，其他保持默认。 4.6.2 设置监视残差 注意点选Plot。 4.6.3 流场初始化 Compute from设置为Inlet。 监视切面 4.7.1 首先切取所需面 以网格Grid为单位， 在X方向，在0到2m之间，每隔0.2m切一平面，以来监视流速和压强的变化； 在Y方向，取Y=0的位置切面，相当于横剖圆柱截面； 在Z方向，取Z=0的位置切面，相当于沿X轴方向竖剖圆柱截面。 注意标清切面名称，以供查找。 4.7.2 设置监视窗口 因不需监视太多所需切面，故建立4个监视窗口即可，需将Plot和Write选取，设为时间步长，再Define内容。