采取蒙特卡罗法分析离散坐标法的假散射.pdfVIP

( ) 第 30 卷 第 6 期 　 　　　华 　中 　科 　技 　大 　学 　学 　报 自然科学版 　　　　　Vol. 30 　No. 6 2002 年 　6 月 　　　J . Huazhong Univ. of Sci. Tech. (Nature Science Edition) 　　　　J un. 　2002 采用蒙特卡罗法分析离散坐标法的假散射 李宏顺 Gilles Flamant 陆继东 华中科技大学 　 IMP/ CNRS , 66125 　 华中科技大学 能源与动力工程学院 FontRomue , France 煤燃烧国家重点实验室 摘要 : 采用蒙特卡罗法获得了离散坐标方程的不含任何假散射的高精度解 ,其基本思路是 : 当采用一个离散 坐标格式来计算一个漫射表面的热辐射时 ,就意味着用有限个离散方向去“代表”2π立体空间上的无穷多个 方向. 因此 ,不妨假定存在着这样一个虚拟表面 :该表面的确只沿着该离散坐标格式的离散方向上发射热辐 射 ,并且在每个方向上的热流也完全遵循离散坐标方程. 然后将蒙特卡罗法应用于此虚拟表面 , 由此所获得的 解即为此虚拟表面的高精度解 ,不言而喻它也等效于离散坐标方程的高精度解. 在此基础上 ,分析了离散坐标 法的假散射对计算结果的影响. 关 　键 　词 : 辐射传热 ; 离散坐标法 ; 假散射 中图分类号: T K124 　　文献标识码 : A 　　文章编号 : (2002) ( ) 　　近 20 年来 ,离散坐标法 DOM 在辐射传热 的真实解 ,显然前者与后者可能不相等. 的数值计算中得到了越来越广泛的应用 , 目前它 迄今为止尚未见有任何离散坐标方程的精确 已被成功地应用到各种辐射传热问题中,而其本 解的文献报导. 本文的基本思想是采用蒙特卡罗 身也在被不断改进 ,新的求解方法、新的离散坐标 法来获取不含任何假散射的离散坐标方程的高精 格式都在不断地被提出. 然而 , 目前已公认离散坐 度解. 为便于阐述 ,先考虑如图 1 (a) 所示的一灰 标法有着两个缺点 : 一是射线影响 , 二是假散 体表面 ,此表面各向同性地向各个方向发射热辐 射[1 ] . 为了克服前一缺点 , 一些学者提出了各种 办法 , 例如 Ramankutty 等[2 ] 、Liou 等[3 ] 、Baek [4 ] 等 采用了改进的微分近似法的基本思想 ,也就 是将辐射强度分为壁面强度和介质强度两个部 分 ;此外 ,Li 等[5 ] 提出了三种方法来降低射线追 图 1 　原理示意图 [6 ] ( 踪法中的射线影响 ,而 Cumber 提出了两种措施 射. 现在 ,假定用一 L SO S 4 离散坐标 该系列可 ) ( ) 从许多文献中得到 来近似地描述这一表面, 由于 来降低离散传输法 D TM 的射线影响 ,这些方 法均可借鉴于离散坐标法. 对于假散射这一问题 , S 4 在每个卦限中仅仅只有 3 个离散方向, 因而此 一些学者试图通过改进空间差分格式来