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

三维圆管流动状况的数值模拟分析三维圆管流动状况的数值模拟分析 在工程和生活中 圆管内的流动是最常见也是最简单的一种流动 圆管流动有层流和紊流两种流动 状况 层流 即液体质点作有序的线状运动 彼此互不混掺的流动 紊流 即液体质点流动的轨迹极为 紊乱 质点相互掺混 碰撞的流动 雷诺数是判别流体流动状态的准则数 本研究用 CFD 软件来模拟研 究三维圆管的层流和紊流流动状况 主要对流速分布和压强分布作出分析 1 物理模型物理模型 三维圆管长 直径 2000mml 100mmd 流体介质 水 其运动粘度系数 62 1 10 m s Inlet 流速入口 1 0 005m s 2 0 1m s Outlet 压强出口 Wall 光滑壁面 无滑移 2 在在 ICEM CFD 中建立模型中建立模型 2 1 首先建立三维圆管的几何模型 Geometry 2 2 做 Blocking 因为截面为圆形 故需做 O 型网格 2 3 划分网格 mesh 注意检查网格质量 在未加密的情况下 网格质量不是很好 如下图 因管流存在边界层 故需对边界进行加密 网格质量有所提升 如下图 2 4 生成非结构化网格 输出 fluent msh 等相关文件 3 数值模拟原理数值模拟原理 3 1 层流流动 当水流以流速 从 Inlet 方向流入圆管 可计算出雷诺数 故圆管内 1 0 005m s 500 d Re 流动为层流 假设水的粘性为常数 运动粘度系数 不可压流体 圆管光滑 则流动的控制方 62 1 10 m s 程如下 质量守恒方程 MERGEFORMAT 1 1 0 uvw txyz 动量守恒方程 MERGEFORMAT 1 2 uuuuvuwuuup txyzxxyyzzx MERGEFORMAT 1 3 vvuvvvwvvvp txyzxxyyzzy MERGEFORMAT 1 4 wwuwvwwwwwp txyzxxyyzzz 式中 为密度 是流速矢量在 x y 和 z 方向的分量 p 为流体微元体上的压强 u w 方程求解 对于细长管流 FLUENT 建议选用双精度求解器 流场计算采用 SIMPLE 算法 属于压 强修正法的一种 3 2 紊流流动 当以水流以流速 从 Inlet 方向流入圆管 可计算出雷诺数 故圆管 2 0 1m s 10000 d Re 内流动为紊流 假设水的粘性为常数 运动粘度系数 不可压流体 圆管光滑 则流动的控制方 62 1 10 m s 程如下 质量守恒方程 MERGEFORMAT 1 5 0 uvw txyz 动量守恒方程 MERGEFORMAT 1 6 2 uuuuvuwuuu txyzxxyyzz uu vu wp xyzx MERGEFORMAT 1 7 2 vvuvvvwvvv txyzxxyyzz u vvv wp xyzy MERGEFORMAT 1 8 2 wwuwvwwwww txyzxxyyzz u wv wwp xyzz 湍动能方程 MERGEFORMAT 1 9 tt kk t k k kkukvkwkk txyzxxyy k G zz 湍能耗散率方程 MERGEFORMAT 1 10 2 1 2 tt kk t k k uvw txyzxxyy C GC zzkk 式中 为密度 是流速矢量在 x y 和 z 方向的分量 p 为流体微元体上的压强 u w 方程求解 采用双精度求解器 定常流动 标准模型 SIMPLEC 算法 k 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 500 d Re 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 监视切面 4 7 1 首先切取所需面 以网格 Grid 为单位 在 X 方向 在 0 到 2m 之间 每隔 0 2m 切一平面 以来监视流速和压强的变化 在 Y 方向 取 Y 0 的位置切面 相当于横剖圆柱截面 在 Z 方向 取 Z 0 的位置切面 相当于沿 X 轴方向竖剖圆柱截面 注意标清切面名称 以供查找 4 7 2 设置监视窗口 因不需监视太多所需切面 故建立 4 个监视窗口即可 需将 Plot 和 Write 选取 设为时间步长 再 Define 内容 例监视 1 监视 Inlet 切面的流速 可设置为 4 8 开始迭代 设置迭代次数为 200 实际比这个更少 迭代收敛时会自动停止 5 层流计算结果及分析层流计算结果及分析 计算 120 步后 已收敛 自动停止运算 残差监视窗口为 5 1 显示流速等值线图 打开 Display Contours 选择 Velocity 和 Velocity magnitude 5 1 1 入口和出口截面的流速分布图 分布在 Surface 里选择 inlet 及 outlet 1 Velocity of inlet 可见 入口处流速分布不明显 基乎都等于入口流速 只是外层靠近壁面处流速几 1 0 005m s 乎为零 符合圆管层流流动规律 也符合边界层理论 2 Velocity of outlet 出口截面流速分布较为明显 显同心圆分布 内层流速偏大 外层靠近壁面处流速几乎为零 边界 层很薄 分层更为严重 层流显现的更为明显 且趋于稳定状态 5 1 2 圆管内不同截面的流速分布图 下述截面均为距 inlet 从 0 2m 到 1 8m 的截面 1 Velocity of inlet 0 2 2 Velocity of inlet 0 4 3 Velocity of inlet 0 6 4 Velocity of inlet 0 8 5 Velocity of inlet 1 6 Velocity of inlet 1 2 7 Velocity of inlet 1 4 8 Velocity of inlet 1 6 9 Velocity of inlet 1 8 上述图像为圆管内部 X 轴方向不同截面的流速分布 可看出流速在截面上从入口到出口的变化 水 流在圆管内部的流速分层很明显 靠近壁面处流速接近于零 有一很薄的边界层 流速在边界层内很快 上升 到最大流速 在圆管中央的一大片圆形区域内 流速基本一致 达到最大 且中心流速最大 为 流速在截面的变化规律可以看出 在 0 到 1 2m 之间 每个截面的流速分布都不同 当离 Inlet max u 1 2m 远之后 流速在截面的分布基本一致 说明层流达到了稳定状态 这符合圆管流动进口段及流中层 流分布规律 以上图像因只能看到沿 X 轴截面的流速分布 故下面讨论从 Y 轴和 Z 轴方向看圆管的整体 流速分布 5 1 3 Y 轴和 Z 轴方向流速截面 截面若均沿圆管长度 X 方向截取 可看到对称的效果 1 Velocity of y 0 整根圆管 入口段 出口段 2 Velocity of z 0 整根圆管 以上两个截面流速分布图的效果是一样的 可以看出圆管水流入口段及之后的流速发展趋势 而且 显示流速变化的规律更为明显 由数值模拟实验设置了入口均匀流速 可以认为在进口处的流速分布是 均匀的 进入管内后 靠近壁面的流动受到阻滞 流速降低 形成边界层 且边界层的厚度逐渐加大 以致尚未受管壁影响的中心部分的流速加快 进口段的流动是流速分布不断变化的非均匀流动 且边界 层的厚度在进口段逐渐增加 之后的流动是各个截面流速分布均相同的均匀流动 由于为层流流动 故 流速分层现象很明显 但平均流速为多少 最大流速为多少 进口段长度为多少 等等问题需要再进行讨论 max uL 5 2 轴向流速的变化 沿 X 轴截取轴线 执行 Plot XY Plot 选择 Y Axis Function 里的 Velocity 和 Velocity Magnitude 选择 Surfaces 里圆管 的对称轴 line x 可得到轴向流速分布散点图 由上图可以看出 在圆管的轴上 进口段流速分布变化较大 从进口流速急剧上升 1 0 005m s 到最大流速 max 0 00848m su 层流入口段长度有经验公式可以算的 即 MERGEFORMAT 1 11 0 058LdRe 可算得入口段长度约为 1 18m 由上图显示效果可以看出 流速在离入口 1 1m 到 1 2m 之间 即入口段长 度约为 1 1 1 2m 符合理论计算结果 5 3 截面流速分布散点图 取流动充分发展后 离 Inlet 1 6m 远的截面 x coordinate 1 6 其流速分布如下图 注意 Plot direction 的选取 可以看处流速沿半径 Y 方向成抛物线分布 与理论公式抛物面公式相符 即 MERGEFORMAT 1 12 22 0 4 gJ urr 取沿 Y 方向中心轴线的流速分布 即 5 4 显示压强分布图 在 Contours 里选取 Pressure 和 Static pressure 在 Surfaces 里选择 int solid 即管道内部流体整体 以两种方式显示 Pressure of int solid top Pressure of int solid isometric 由以上两图可以看出圆管内部压强分布从管口处向延伸方向逐渐减小 可知流速相应增大 符合流 速大 压强小的流动定律 也符合圆管流动压降的原理 另外从入口处的压强分布可以看出 在圆管任 何截面上 其压强分布也不是均匀的 也有分层现象 5 5 轴向压强的变化 执行 Plot XY Plot 选择 Y Axis Function 里的 Pressure 和 Pressure Magnitude 选择 Surfaces 里圆管 的对称轴 line x 可得到轴向压强分布散点图 圆管层流中的压降 理论上存在下述公式 MERGEFORMAT 1 13 4 128 V Lq p d 即压降与流体的粘度 管道长度 流体的流量成正比 在本模拟实验中 由于流体的粘度 流体的流量 不变 可认为压降与长度成正比 即与成正比 由上图可以看出 除了入口段压强分布因流速急剧p L 上升而下降过快外 其余部分均可看做是一条直线 即随的增加而降低 是正比关系 p L 5 6 总结报告 5 6 1 系统总流量 Mass Flow Rate kg s inlet 0 039138829 int solid 7 8335983 outlet 0 039138853 wall 0 Net 2 4280432e 08 5 6 2 入口出口流速积分 Integral Velocity Magnitude m s m2 inlet 3 8362443e 05 outlet 3 9154264e 05 Net 7 7516706e 05 5 6 3 入口出口压强积分 Integral Static Pressure pascal m2 inlet 0 00038809504 outlet 0 Net 0 00038809504 6 在在 FLUENT 中求解计算紊流流动中求解计算紊流流动 6 1 FLUENT 设置 除以下设置为紊流所必须设置的外 其余选项和层流相同 不再详述 Viscous 设置 雷诺数 故圆管内流动为紊流 Viscous 设置为 Realizable K epsilon 模型 其余默10000 d Re 认 Boundary 设置 Inlet 设置为速度入口 为 Turbulence 设置为 Intensity and Hydraulic Diameter 方法 即 2 0 1m s Outlet 设置为自由出口 Outflow 如设置成压力出口 则之后计算会存在问题 已验证 Solution 设置 采用双精度求解器 定常流动 Realizable 模型 SIMPLEC 算法 k 6 2 开始迭代 设置迭代次数为 300 实际比这个少 迭代收敛时会自动停止 7 7 紊流计算结果及分析紊流计算结果及分析 计算 293 步后 已收敛 自动停止运算 残差监视窗口为 7 1 显示流速等值线图 7 1 1 入口和出口截面的流速分布图 分布在 Surface 里选择 inlet 及 outlet 1 Velocity of inlet 可见 入口处流速分布不明显 基乎都等于入口流速 只是外层靠近壁面处流速几乎 2 0 1m s 为零 2 Velocity of outlet 可见 出口截面流速分布较为明显 和层流一样 显同心圆分布 内层流速偏大 外层靠近壁面处 流速几乎为零 分层更为严重 边界层很薄 7 1 2 Y 轴和 Z 轴方向流速截面 圆管内各个截面的流速分布均不相同 可以认为紊流还没达到稳定状态 在此不再分析各个截面的 流速分布 仅对整个圆管的流速作出分析 截面沿圆管长度 X 方向截取 可看到对称的效果 1 Velocity of y 0 整根圆管 2 Velocity of z 0 整根圆管 以上两个截面流速分布图的效果是一样的 可以看出圆管水流紊流入口段及之后的流速发展趋势 而且显示流速变化的规律更为明显 3 入口段 与层流入口段的流速分布相比 可以明显的看出紊流入口段的流速分布不太明显 且基本没有分层 符合紊流流动的基本规律 流速分布也不像层流流速那样显明显抛物线分布 而是更加平滑 越超后发 展发展越平滑 到底是什么曲面 之后再加分析 紊流过流断面的流速对数分布比层流的抛物面分布均 匀的多 符合的规律 即 1 ln u yC uK 4 出口段 出口段的流层分布很明显 切趋于均匀 但仔细观察圆管轴心的速度 其实速度分布并未达到均匀 可见紊流并未达到充分发展的状况 7 1 3 轴向流速的变化 执行 Plot XY Plot 选择 Y Axis Function 里的 Velocity 和 Velocity Magnitude 选择 Surfaces 里圆管 的对称轴 line x 可得到轴向流速分布散点图 由上图可以看出 在圆管的轴上 进口段流速分布变化较大 从进口流速急剧上升到 2 0 1m s 最大流速 之后又下降 但实际经验表明 紊流应该在进口段后达到稳定状态 轴向 max 0 1369m su 流速应该趋于恒定 可见此模拟实验设置长度不够 使流动并未达到充分紊流 紊流入口段长度有经验公式可以算的 即 MERGEFORMAT 1 11 2540Ld 由此可见 紊流的边界层厚度的增长比层流边界层要快 因此紊流的进口段要短些 而且长度主要受来 流扰动的程度有关 与雷诺数无关 扰动越大 进口段越短 可算得入口段长度约为 3m 由上图显示效 果可以看出 轴向流速一直在变化 并未达到最大且稳定的速度 故紊流未发展充分 改进实验应加大 圆管长度 7 1 4 出口截面的流速分布散点图 因紊流并未充分 故选取出口截面来进行分析 注意 Plot direction 的选取 可见截面流速分布已很平滑 与层流出口截面的流速分布截然不同 若紊流充分发展 则截面流速 散点图最高处几乎为一条直线 说明圆管内大多数流体流速趋于稳定 几乎没有分层 取沿 Y 方向中心轴线的流速分布 即 7 2 显示压强分布图 在 Contours 里选取 Pressure 和 Static pressure 在 Surfaces 里选择 int solid 即管道内部流体整体 Pressure of int solid top 和层流圆管内压强分布一样 进口压强大 出口压强小 即存在压降 另外在圆管任何截面上 其 压强分布是均匀的 没有分层现象 这点和层流截面压强分布很不同 7 3 轴向压强的变化 执行 Plot XY Plot 选择 Y Axis Function 里的 Pressure 和 Pressure Magnitude 选择 Surfaces 里圆管 的对称轴 line x 可得到轴向压强分布散点图 圆管紊流中的压降 虽然不存在理论上的经验公式 但从上图可以看出 紊流的压降和层流类似 除了入口段压强分布因流速急剧上升而下降稍快外 其余部分均可看做是一条直线 即随的增加而p L 降低 是正比关系 7 4 总结报告 5