圆柱坐标系下辐射换热离散坐标法的影响因素分析.pdfVIP

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