求解高方向分辨率辐射强度两种数值方法的比较分析.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083197 求解高方向分辨率辐射强度两种 数值方法的比较分析 黄志锋 周怀春 程 强 （华中科技大学煤燃烧国家重点实验室，武汉，430074） Tel:027 Email: hczhou@mail.hust.edu.cn 摘 要 针对一维平行平板介质的辐射传递过程求解算例，分别用DRESOR 法和DOS+ISW 法进行求 解，对两种数值方法得到的高方向分辨率的辐射强度代入到辐射传递方程进行检验。结果表明，在介 质温度均匀时，两种方法计算得到的辐射强度都具有很好的准确性，在介质温度不均匀时，DOS+ISW 法计算得到的辐射强度存在明显的误差，而DRESOR 法依然能够得到精度较高的辐射强度。同时，对 两种方法计算得到的辐射热流也进行了比较，DOS+ISW 法只有在选取较多的初始离散方向时才能得 到较好的结果。 关键词 辐射强度 DRESOR 法 DOS+ISW 法 高方向分辨率 0 引言 辐射传热问题的求解一直以来是国内外众多学者所关注的问题。到目前已经发展了 多种数值计算方法，如离散坐标法、离散传递法、区域法、球形谐波法、蒙特卡洛法等。 这些方法都各有优势，在实际应用过程中需根据实际问题选取合适的方法。获取高方向 分辨率的辐射强度在许多理论研究和实际应用领域都有重要的作用，例如辐射传热逆问 [1] 题 。文献[2]提出了一种基于蒙特卡洛法的辐射传递问题新方法，它继承了蒙特卡洛法 [3] 的优点，同时还可以求解得到高方向分辨率的辐射强度。DOS+ISW 法 是通过向原有 的离散坐标系列中加入一些权值无穷小的新方向，通过对新的离散坐标系列进行求解也 可以获得任意方向的辐射强度分布。 本文利用DRESOR 法和DOS+ISW 法同时研究了一维平行平板介质内的辐射传递 过程，对计算得到的辐射强度和辐射热流进行了比较分析，以期对两种方法在求解辐射 传热问题过程中的精确性作出评价分析，便于在研究或实际应用过程中有针对性地选择 合适的方法。 1 数值方法介绍 1.1 DRESOR 法求解高方向分辨率辐射强度简介 介质具有吸收、发射、各向异性散射特性的辐射传递方程为：  ˆ ˆ ˆ ˆ s ˆ ˆˆ s I (r, s)  I (r, s)  I (r, s)  I (r, s ) (s , s)d  (1)   b   4 4  i  i i 资助项目：国家自然科学基金 （No No）和 “111”引智计划 （B06019 ）资助项目 ˆ 它描述了辐射强度沿着路径方向 上的变化. s [4] 其积分形式为 ：  s  s  s  I (r , s) I (r , s)exp   ds  S (r , s)exp   ds  ds (2)  w  0   0   0    这里, s r rw , s r r ,积分方向转为沿s (从点r 指向壁面).