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

《顶盖驱动流数值模拟分析.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、《数值传热学》作业：顶盖驱动流数值模拟分析西安科技大学能源学院安全技术及工程申敬杰顶盖驱动流数值模拟分析顶盖驱动流作为经典的数值计算模型，常常用来考核源程序和计算思想的正确性。这种流动边界条件简单，而且不涉及模型的影响，便于直接评价差分格式的性能。1.引言数值传热学，又称计算传热学，是指对描写流动与传热问题的控制方程采用数值方法，通过计算机求解的一门传热学与数值方法相结合的交叉学科。数值传热学的基本思想是把原来在空间与时间坐标中连续的物理量的场（如速度场，温度场，浓度场等），用一系列有限个离散点上的值的集合来代替，通过一定的原则建立起这

2、些离散点变量值之间关系的代数方程（称为离散方程）。求解所建立起来的代数方程已获得求解变量的近似值。由于实验方法或分析方法在处理复杂的流动与换热问题时，受到较大的限制，例如问题的复杂性，即无法做分析解，也因为费用的昂贵而无力进行实验测定，而数值计算的方法正具有成本较低和能模拟复杂或较理想的过程等优点，数值传热学得到了飞速的发展。特别是近年来，计算机硬件工业的发展更为数值传热学提供了坚实的物质基础，使数值模拟对流动与传热过程的研究发挥了重要的作用。目前，比较著名的数值模拟分析应用软件有FLUENT、CFX、STAR-CD、和PHOENICS

3、等，而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/sL=0.12

4、m图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：设置边界类型操作：→l

5、在Name栏输入边界名称wall-1，将Type栏选为Wall，在Entity栏选取Edges，并选中方腔顶部边线。l在Name栏输入边界名称wall-2，将Type栏选为Wall，在Entity栏选取Edges，并选中方腔其它三条边线。Step5：输出网格文件操作：Film→export→mesh打开对话框如图5所示，选中Export2-Dmesh前面的复选框，输出网格文件。图5网格文件输出对话框4.求解计算求解计算的操作步骤如下：Step1：启动Fluent选择2d单精度求解器，点击Run，如图6所示。Step2：设置Fluent运

6、算模式操作：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．

7、显示网格操作：Display→Grid...如下图所示：图9显示生成网格Step4：选择计算模型1．设置求解器操作：Define→Models→solve...保留默认设置。点击OK确认。Step5：定义流体材料性质操作：Define→Materials...打开对话框如图10所示，在该对话框中可定义材料的物理属性，也可以从材料数据库中选择其它材料，或者创建新的材料。在本例中，只要选择空气就可以了，在定义空气的物理性质进行一些调整。（1）将Density的值改为1；（2）将Viscosity中的值改为2e-05；（3）点击Change/

8、Create，然后关闭对话框。图10材料设置对话框Step5：设置边界条件操作：Define→BoundaryConditions...打开“BoundaryConditions”设置对话框如图11所示。图