顶盖驱动流数值模拟分析.docVIP

《数值传热学》作业：

顶

盖

驱

动

流

数

值

模

拟

分

析

西安科技大学能源学院平安技术及工程

申敬杰201112612

顶盖驱动流数值模拟分析

顶盖驱动流作为经典的数值计算模型，常常用来考核源程序和计算思想的正确性。这种流动边界条件简单，而且不涉及模型的影响，便于直接评价差分格式的性能。

1.引言

数值传热学，又称计算传热学，是指对描写流动与传热问题的控制方程采用数值方法，通过计算机求解的一门传热学与数值方法相结合的交叉学科。数值传热学的根本思想是把原来在空间与时间坐标中连续的物理量的场〔如速度场，温度场，浓度场等〕，用一系列有限个离散点上的值的集合来代替，通过一定的原那么建立起这些离散点变量值之间关系的代数方程〔称为离散方程〕。求解所建立起来的代数方程已获得求解变量的近似值。

由于实验方法或分析方法在处理复杂的流动与换热问题时，受到较大的限制，例如问题的复杂性，即无法做分析解，也因为费用的昂贵而无力进行实验测定，而数值计算的方法正具有本钱较低和能模拟复杂或较理想的过程等优点，数值传热学得到了飞速的开展。特别是近年来，计算机硬件工业的开展更为数值传热学提供了坚实的物质根底，使数值模拟对流动与传热过程的研究发挥了重要的作用。

目前，比拟著名的数值模拟分析应用软件有FLUENT、CFX、STAR-CD、和PHOENICS等，而FLUENT是国内外比拟流行的商用CFD软件包，该软件以其市场占有率高、计算准确、界面友好、使用简单、应用领域广、物理模型多而获得较高的市场占有率和用户的肯定。

2.物理模型

在一个正方形的二维空腔中充满等密度的空气，方腔每边长为0.12m，取雷诺数为Re=12000，由Re=vd/υ，方腔的当量直径，计算知d=0.12m，又υ=15.7×10﹣6m2/s，那么顶盖驱动流的速度v=1.57m/s，即其顶板以1.57m/s的速度向右移动，同时带动方腔内流体的流动，流场内的流体为紊流。计算区域示意图如图1所示。

v=1.57m/s

L=0.12m

图1计算区域示意图

3.在Gambit中建立模型

在Gambit中建立模型的操作步骤如下：

Step1：启动Gambit并选择求解器为Fluent5/6。

Step2：创立面

操作：→→

翻开对话框如图2所示。输入长度和宽度12，在Direction中选择XYCentered。

图2创立面设置对话框

Step3：划分面网格

操作：→→

翻开对话框如图3所示，Shift+鼠标左键选择正方形面，Internalsize=0.5，其它保存默认，点击Apply确认。划分后的网格如图4所示。

图3网格划分设置对话框图4计算区域网格图

Step4：设置边界类型

操作：→

在Name栏输入边界名称wall-1，将Type栏选为Wall，在Entity栏选取Edges，并选中方腔顶部边线。

在Name栏输入边界名称wall-2，将Type栏选为Wall，在Entity栏选取Edges，并选中方腔其它三条边线。

Step5：输出网格文件

操作：Film→export→mesh

翻开对话框如图5所示，选中Export2-Dmesh前面的复选框，输出网格文件。

图5网格文件输出对话框

4.求解计算

求解计算的操作步骤如下：

Step1：启动Fluent

选择2d单精度求解器，点击Run，如图6所示。

Step2：设置Fluent运算模式

操作：define→models→viscousmodels...设置如图7所示：

图6启动求解器图7设置运算模式

图8网格尺寸设置对话框

Step3：导入并检查网格

1．读入网格文件

操作：File→Read→Case...

找到文件后，单击OK按键确认。

2．检查网格

操作：Grid→Check

网格读入后，一定要进行网格检查，注意最小体积不能为负值。

3．网格比例设置

操作：Grid→Scale...

在Gambit中，生成网格使用的单位是cm，在GridWasCreatedIn下拉菜单中，选取cm，如图8所示，然后单击Scale，关闭对话框。

4．显示网格

操作：Display→Grid...如下列图所示：

图9显示生成网格

Step4：选择计算模型

1．设置求解器

操作：Define→Models→solve...

保存默认设置。点击OK确认。

Step5：定义流体材料性质

操作：Define→Materials...

翻开对话框如图10所示，在该对话框中可定义材料的物理属性，也可以从材