fluent的一个实例(波浪管道的内部流动模拟).

基于基于 FLUENT 的波浪管道热传递耦合模拟的波浪管道热传递耦合模拟 CFD 可以对热传递耦合的流体流动进行模拟 CFD 模拟可以观察到管道内 部的流动行为和热传递 这样可以改进波浪壁面复杂通道几何形状中的热传递 目的 目的 1 创建由足够数量的完整波浪组成的波浪管道 提供充分发展条件 2 应用周期性边界条件创建波浪通道的一部分 3 研究不同湍流模型以及壁面函数对求解的影响 4 采用固定表面温度以及固定表面热流量条件 确定雷诺数与热特性之间 的关系 问题的描述 问题的描述 通道由重复部分构成 每一部分由顶部的直面和底部的正弦曲面构成 如 图 图 1 管道模型 空气的流动特性如下 质量流量 m 0 816kg s 密度 1kg m3 动力粘度 0 0001kg m s 流动温度 Tb 300K 流体其他热特性选择默认项 流动初试条件 x 方向的速度 0 816m s 湍动能 1m2 s2 湍流耗散率 1 105m2 s3 所有湍流模型中均采用增强壁面处理 操作过程 操作过程 一 完整波浪管道模型的数值模拟 1 计算 Re uH v 0 816 1 0 0001 1 8160 Cf 2 0 0359Re 0 2 0 0359 8160 0 2 0 0059259 0628 0 0059259 0 816 0 2 0 f t C uu y uty v y 0 00159 2 创建网格 本例为波浪形管道 管道壁面为我们所感兴趣的地方所以要局部细化 入 口和出口处的边界网格设置如图 图 2 边网格 生成面网格 图 3 管道网格 3 运用 Fluent 进行计算 本例涉及热传递耦合 所以在 fluent 中启动能量方程 如图 图 4 能量方程 设定条件 湍流模型选择标准 k e 模型 近壁面处理选择增强壁面处理 图 5 湍流模型 设定材料 密度为 1 动力粘度改为 0 0001 如图 图 6 材料设定 设定边界条件 入口速度为 0 816 湍动能为 1 湍流耗散率为 100000 出 口为自由出口 壁面温度为固定温度分别为 300k 500k 图 7 边界条件 初始化 并计算 图 8 残差 残差中的 e 和 k 并没有减小 没有达到 10 3一下 并且由于网格很大 计 算时间很长 图 9 压力分布图 随着流体流动 管道中压力分布趋于平稳 波浪管道中波谷的压力最低 在入口处的压力较高 图 10 速度分布图 从图中可以看出 在管道 7 11 个波浪处 流动已经充分发展 贴近上壁处 速度最大 在波谷出的速度最小 甚至接近于零 图 11 温度分布图 贴近波浪壁面出的温度较高 流动充分发展后 由于换热作用 管道后部 流体温度逐渐升高 在波峰与波谷之间的流体温度最高 如图 图 12 速度矢量图 可以看出 波浪壁面出流体出现反流 在波谷出反流的流体最多 速度在 波谷出最小 接近于 0 出现滞留区 若要观察波峰 波谷处流体流动速度 需要在波峰 波谷处创建两条直线 观察直线上的速度 因为管道 7 11 节处流动充分发展 所以在第十一节波峰 波谷处建立两条直线 如图 图 13 波峰 波谷 图 14 波峰 波谷的速度 波峰贴近两侧壁面出的速度梯度很大 在管道中间速度随高度增加而增加 在 0 7m 左右达到最大 波谷处靠近上壁面的速度梯度很大 但是由于有波谷存 在 波谷处的速度梯度不大 在谷中速度先增大再减小 在 0 处左右达到最下 随后逐渐增加 在 0 7m 高度左右速度最大 高度在 0 5m 处以上波峰波谷处的 速度基本相等 二 周期性波浪管道模型的数值模拟 图 15 周期性网格 网格密度与完全管道网格相同 在 fluent 中输入以下指令 创建周期性网格 grid modify zones grid modify zones make periodic Periodic zone inlet Shadow zone outlet Rotational periodic if no translational yes no Create periodic zones yes yes Auto detect translation vector yes yes computed translation deltas 1 000000 0 000000 all 100 faces matched for zones 6 and 5 zone 5 deleted created periodic zones grid modify zones 边界条件中可以看到已经没有 outlet inlet 也变成了 periodic 周期性的 这 里要设置周期性边界条件 质量流量为 0 816 其他设置与之前相同 图 16 边界条件 图 17 残差 Energy 并不收敛 反而随计算而发散 其他参数都收敛 计算量很小 计 算速度明显提高 图 18 压力图 图 18 压力图 速度图 温度图 速度矢量图 管道中间压力最大并比较均匀 在波浪波谷出压力最低 温度分布不太正 确 图 19 波峰波谷速度图 波峰贴近两侧壁面出的速度梯度很大 与完全模型相似 波峰在管道中间 速度随高度增加先增加在减小 但是变化不大 其他分布趋势与完全模型相似 图 20 周期模型与完全模型波峰速度图比较 周期性模型 管道中速度基本相等 没有呈现完全模型的变化趋势 图 21 周期模型与完全模型波谷速度图比较 两模型的波谷速度的变化趋势大致相同 在小于 0 3m 处 周期性模型速度 大于完全模型 0 3m 之后 周期性模型速度小于完全管道模型 三 周期性波浪管道不同湍流模型的数值模拟 应用周期模型 再选分别应用 RNG 与 Realizable 湍流模型计算 实验数据 如图 图 21 RNG 模型波峰波谷速度图 图 22 Realizable 模型波峰波谷速度图 将之前的数据放在一起进行比较 图 23 不同湍流模型波峰速度图 图 24 不同湍流模型波谷速度图 从图中可以看出 realizable 和标准 k e 模型波峰的速度分布大致相同 RNG 模型与完全模型更接近 变化更大 三种湍流模型波谷的速度分布与波峰分布 特点相似 四 整体模型壁面采用固定热流率的数值模拟 重新选用整体模型 只需改变波浪面即 wall2 的边界条件即可 改为热流 率 1000w m2 其他条件不变 图 25 不同湍流模型波谷速度图 图 26 固定壁面温度与固定热流率波峰速度图比较 图 27 固定壁面温度与固定热流率波谷速度图比较 从图中可以看出固定壁面温度与固定热流率波峰与波谷的速度分布还是比 较相似的 不同的是在 0 6m 处固定热流率波峰波谷速度比固定壁面温度的速度 稍高 在 0