圆柱坐标系下辐射换热离散坐标法的影响因素分析 discrete ordinates method for radiative transfer equation in cylindrical coordinate.pdfVIP

第9卷第3期 热 科 学 与 技 术 V01．9NO．3 2010年9月 Scienceand Sep．2010 ofThermal Journal Technology 文章编号：1671—8097(2010)03—0218—07 圆柱坐标系下辐射换热离散坐标法的影响因素分析 薛然然， 刘玉英。， 金劲睿， 张荣春 (北京航空航天大学能源与动力工程学院，北京100191) 摘要：以辐射平衡条件下的一维无限长同心圆环之间的发射、吸收、非散射介质的辐射换热问题为例，探讨 了辐射传递方程的离散坐标求解方法，并分析了角积分方案、空间差分格式及角度差分格式对圆柱坐标系下 离散坐标法的影响。研究结果表明：角积分方案近似阶数越高，计算结果的精度越高，但随着积分方案阶教的 增加，所需计算时间也相应增大；对于同一近似阶数的角积分方案．F。角积分方案的计算精度高于SⅣ角积分 方案；空间差分格式对计算精度的影响较小，其中阶梯格式精度最高．正格式所需计算时间最长，菱形格式所 需时间最短；角度差分格式对计算精度的影响较大。其中菱形格式计算精度最高，阶梯格式和正格式的计算精 度相当，同时三种角度差分格式所需计算时阃相当。因此，角积分方案一定的情况下，建议空间与角度差分格 式分别选用阶梯格式和菱形格式。 关键词：圆柱坐标系；辐射传递方程；离散坐标法；角分布项；差分格式 中图分类号：V235．11文献标识码：A 系中辐射传递方程的离散坐标法与直角坐标系中 0 引 言 的解法有很大的不同‘7|。本文以无限长同心圆筒 辐射换热是航空发动机燃气轮机、火箭喷射 之间吸收、发射、无散射介质的辐射换热问题为 器、工业炉、锅炉等高温燃烧设备中重要的传热方 例，介绍了一维圆柱坐标系中辐射传递方程的离 式之一[卜3。。辐射换热的计算以求解辐射传递方程 散坐标法，编制了DOM计算程序，通过具体算例 (radiativetransfer equation，RTE)为基础，由于 验证了程序的可行性，并对影响离散坐标法计算 RTE是一个复杂的积分一微分方程，只在某些简 精度的角积分方案及差分格式的选取进行了讨论 单或极限情况下才可能得到精确解或近似解，因 分析。 此实际的工程应用中多采用数值解法，如区域法、 1 圆柱坐标系下的离散坐标法 热流法、离散传递法、有限容积法及离散坐标法 ordinates (discrete method，DOM)等。DOM因易 1．1 圆柱坐标系下的辐射传递方程 与流动及燃烧计算耦合求解、对计算机的性能要 以图1所示的一维柱坐标系为例，发射、吸 求较低、可以很方便地处理各向异性散射项等优 收、非散射灰体介质在位置r(r，么，z)处、沿辐射 点，得到广泛的研究和发展[4-8]。 传输方向；的辐射传递方程可表示为Ⅲ 与直角坐标系相比，在诸如圆柱、球等曲线坐 等导(订(r，；))一了1磊(矿(r，；))； 标系中，即便光子传递方向不变，沿直线传播，在 rdr rd曲。 (1)