非结构化网格中辐射传热的数值计算.pdfVIP

2w 非结构化网格中辐射传热的数值计算 1 2 1 张敏 ，John C. Chai ，许彬 1 南京理工大学动力工程学院，南京（210094 ） 2School of Mechanical and Production Eng., Nanyang Tech. University, Singapore （639798 ） 摘 要：用基元有限体积法和非结构化网格求解吸收/散射介质空间的辐射传热问题。空间 离散采用三角形非结构化网格，方向角离散采用四边形结构化网格，并且采用二阶迎风格式 计算。计算结果分别与精确解和结构化网格的数值解相比较，并得到了满意的一致性。虽然 这两种网格都能得到满意的结果，但在矩形域中用结构化网格比用非结构化三角形网格更精 确些。 关键词：非结构网格，有限体积法，辐射传热 引言 尽管数值辐射传热研究的历史不长，但在近十几年有了飞速的发展。利用有限体积法 （FVM ）在非结构化网格上进行数值计算，是近几年才开始的研究领域。Raithby 和 Chai 是最先采用FVM 进行数值辐射传热研究的学者之一。1990 年，Raithby 和Chui 给出了在笛 卡尔直角坐标系中的数值辐射传热的计算方法[1] ，接着又进行了在圆柱坐标系(Chui et al.,1992)[2]和非正交坐标系(Chui and Raithby,1993)[3] 中数值辐射传热的研究。1994 年，Chai 把上面的过程推广到可穿透的吸收、辐射和散射介质之中，以及平射光束和不规则的几何形 [4—7] 状中 。 上世纪末，Murthy 、Mathur(1998)[8]、Vaidya(1999)[9]、Raithby(1999)[10]和Liu et al.(1999)[11] 给出了在非结构化网格中 FVM 的辐射传热计算。近年来不少商用软件也采用了这些结果 (Liu et al.,1999 ；Murthy and Mathur,1998 ；Vaidya,1998[12]) 。这些都预示着这些研究在未来的 工业界有着广阔的应用前景。对此感兴趣的读者可查阅2002 年Chai 和Patankar 的综述报告 [13] 。 1. 控制方程和边界条件 1.1 辐射传热方程 [13] 穿过灰体介质的辐射输送能量方程（Radiative Transfer Equation ，RTE ）可写成 ： r ˆ dI (r ,s) r r ˆ r ˆ −β (r )I (r ,s) + SR (r , s) (1.1) ds 其中，衰减（extinction ）系数 S β 和源项函数 分别为： R r r r β (r ) κ (r ) + σ s (r ) (1.2) r r ˆ r r σ s (r ) r ˆ′ ′ S R (r , s) κ (r )I b (r ) + 4π