算例5 辐射与自然对流模拟.pdfVIP

《数值计算与工程仿真》专刊—FLUENT HELP 算例精选中文版（一） 算例5 辐射与自然对流模拟 引言 在这个算例中，将会解决二维方箱中的辐射与自然对流相结合的问题，网 格采用四边形单元网格。 在这个算例中将会学到以下知识点： 1.应用 Fluent 中各种辐射模型 Rosseland, P-1, DTRM, 离散坐标模型(DO), 和 surface-to-surface (S2S)，并理解各个模型的适用范围。 2.使用 Boussinesq model 定义密度； 3.设定辐射与自然对流传热问题的边界条件； 4.将单一的墙划分为多个墙区域； 5.对已有的流体物性进行修改； 6.用隔离求解器求解； 7.显示速度矢量和流函数等值线，以及温度等值线； 前提条件 本指南是在你熟悉掌握 Fluent 的菜单操作，并且已经熟练进行指南 1 练习 的基础上继续的，有些步骤和操作过程将不会在以下内容中详细介绍。 问题描述 将被考虑的问题如图 5.1 所示，一个边长为 L 的正方形箱体，右墙温度为 2000K,左墙温度为 1000K,上下墙绝热，重力向下，由于热重引起密度梯度所以 发展为浮力流。箱体中的介质被认为是有吸收性和散射性的，因此墙壁间的辐 射交换因存在吸收被减弱，同时也因为介质的散射作用而增强了。所有墙壁被 认为是黑体，目的在于应用有效的辐射模型计算箱体中流场和温度场分布，以 及墙壁的热流量，并且对于不同光学深度 aL 比较 所表现出的特性。 — 99 — ─ 99 ─ www.myCFD.cn 《数值计算与工程仿真》专刊—FLUENT HELP 算例精选中文版（一） 工质普朗特数大约为 0.71，基于L 的雷诺数为 500000，这说明流动相当于 层状流动，应用 Boussinesq 假设来模拟浮力流动。普朗克数 为 0.02， 用于考虑传导与辐射的相对重要性，其中，T 0 = ( T h + T c)/2 。在这个算例中将 有三种 optical thickness 的情况会被考虑到，分别是 aL=0, aL=0.2, and aL=5 。 注意：物理属性和工作条件（重力加速度）都已经给定以适合于产生的想 要的普朗特数，雷诺数和普朗克数。如下图所示： 图5.1：问题描述 准备工作 1. 将网格文件 rad/rad.msh 拷至 fluent 的工作目录下（就像在指南 1 中描述 的一样）； 2. 起动fluent的二维单精度解算器。 — 100 — ─ 100 ─ www.myCFD.cn 《数值计算与工程仿真》专刊—FLUENT HELP 算例精选中文版（一） 第1 步： 网格 1. 读取网格文件rad.msh.。 File Read Case... 当网格读入的时候，在Fluent控制窗口会显示相应的信息，会报告网格 有2500个单元。 2. 检查网格质量。 Grid Check Fluent会对网格进行各种各样的检查，并会在控制窗口显示信息。特别 注意最小体积，确保它是正数值。 3. 显示出网格（如图5.2）。 Display Grid... — 101 — ─ 101 ─ www.myCFD.cn 《数值计算与工程仿真》专刊—