第五讲 二维室内气流数值计算实例（学生用）.pdfVIP

第七章 二维室内气流数值计算实例 §7.1 二维室内气流标准算例描述 As proposed in IEA Annex 20, a two-dimensional case, for which detailed experimental data are available, has been specified by Nielsen (1990) to test different CFD codes. This report presents the results computed by the PHOENICS- 84 code (Rosten and Spalding, 1987) and the comparison between the computed and measured results. 一、等温状态计算 The computations are carried out for two test cases: (1) an isothermal case and (2) a non-isothermal case as shown in Figure 1. The geometry of the test case is L/H= 3.0, h/H= 0.056, and t/H = 0.16, where H = 3.0 m. The inlet conditions are specified by a Reynolds number： where kinematic viscosity v for air is 15.3x106 m2/s at 20°C. The inlet velocity can then be determined as u0 = 0.455 m/s The inlet conditions for turbulent kinetic energy k and dissipation rate e are given by： 耗散率（dissipation rate ）： Case 1 边界条件： 1 / 17 U=0.455m/s K0=1.5(0.04u0)^2=4.969e-4 E0=k0^1.5(0.04u0)^2=1.07e-3 二、非等温状态： For the non-isothermal case, the aim is to predict flow with strong buoyancy effect. A constant heat flux is added along the floor so that the cold jet may deflect from the ceiling before it reaches the end wall. The buoyancy effect is weighted by the Archimedes number defined as ： where ρ, g and ∆T0 are volume expansion coefficient, gravitational acceleration and temperature difference between return and supply, respectively. 三、建模与计算 1、利用GAMBIT 建立混合器计算模型 （1）利用坐标网格创建节点； （2）在两个节点之间创建直线 （边）； (3) 由边创建面； （4）对各条边定义网格节点的分布； (5) 在面上创建网格； （6） 定义边界类型； （7） 为FLUENT5/6 输出网格文件。 2、利用FLUENT—2D 求解器进行求解 （1）读入网格文件； (2) 确定长度单位为m; （3）确定流体材料及其物理属性; (4) 确定边界条件； 2 / 17 （5）计算初始化并设置监视器； （6）使用非耦合、隐式求解器求解； （7）利用图形显示方法察看流场与温度场； 三、结果分析： The grid of 45X29, sta