球谐函数法求解各向异性散射辐射传输.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093261 球谐函数法求解各向异性散射辐射传输 王雁鸣 谈和平 易红亮 （哈尔滨工业大学，能源科学与工程学院456 信箱，哈尔滨 150001） (Tel: 0451 Email: hit_wangym@) 摘要：本文针对各向异性散射介质内的红外辐射传输开展数值方法研究。应用球谐函数展开推导了一 维吸收、发射、散射性灰介质的辐射传输近似方程组，采用差分方法结合三对角矩阵解法建立了任意 阶辐射传输近似方程数值解法。验证计算结果表明，高阶球谐函数数值解与理论解吻合的很好，该方 法可以求解任意非线性各向异性散射辐射传输问题。 关键词：辐射传输 球谐函数法 各向异性散射 0 前言 目前，已发展了很多求解吸收散射性介质内红外辐射传输的数值方法，如热流法、 [1] 区域法、蒙特卡洛法、离散传递法、离散坐标法、有限容积法等 。球谐函数法(Spherical Harmonics Method) ，又称 PN 近似法。该方法首先由Jeans[2]在研究星际辐射传输时提出， 后来经Murray[3]，Davison[4]和 Kourganoff[5]进一步完善并应用在中子传输问题中，近几 年来，Modest[6,7]将其应用于三维辐射传热中。国内目前对于该方法在各向异性散射参与 [8] 性介质内红外辐射传输求解问题中的应用尚不多见 。 本文对一维吸收、发射、散射性灰介质的辐射传输方程进行球谐函数形式展开，推 导了任意阶PN 近似方程，并采用差分方法及三对角矩阵法求解各向异性散射辐射传输问 题。 1 各向异性散射介质内辐射传输球谐函数方法 1.1 PN 近似方程 一维吸收、发射、散射性灰介质的辐射传输方程（RTE ）为： dI κ Ω =− κ + κ I + κ I + s I Φ Ω,Ω d Ω (1a) dL ( a s ) Ω a b,Ω 4π ∫Ωi 4π Ωi ( i ) i 其中κ 、κ 为吸收、散射系数，方程整理可得： a s dI τ ,µ 1 ( ) ω ′ ′ ′ µ =−I τ ,µ + 1−ω I τ + I τ ,µ Φ µ ,µ d µ (1b) dτ ( ) ( ) b ( ) 2 ∫− 1 ( ) ( ) 其中µ cosθ ，光学厚度τ κ =+ κ L ，反照率ω κ / κ =+ κ 。 ( a s ) s ( a s ) 将辐射强度按球谐函数展开： 资助项目：国家自然科学基金重大国际合作研究项目(No.50620120442)资助 N I (τ ,µ ) ∑ I (τ ) P (µ ) (2) n n n