DRESOR法求解一维各向异性散射性介质中的辐射传递方程.pdfVIP

摘 要 辐射传热是一种重要的传热方式，而辐射传递方程（Radiative Transfer Equation ） 求解是目前辐射换热研究领域里的热门话题。由于 RTE 的积分/微分性质，同时参与 辐射换热介质的发射、吸收、散射等特性使得辐射传递方程的求解更加复杂。经过近 百年的发展，辐射换热数值计算领域涌现出了多种辐射换热计算方法，其中常用的四 种基本方法为：离散坐标法（DOM ，Discrete Ordinate Method ）、球形谐波法（Spheric Harmonics ，PN ）、区域法（Zonal ）、Monte Carlo 方法（MCM ）。MCM 可以综合考虑 辐射传热模拟过程中所有重要的影响而不需要近似处理，同时，它可以用重复多次的 简单关系来替代一定程度上的数学复杂问题，这样，复杂的物理问题可以很简单的编 制计算程序。因此，Monte Carlo 法对辐射换热研究具有重要的意义。 为充分利用Monte Carlo 法在辐射传热过程模拟方面的优越性，同时能够给出具 有较高方向分辨率的辐射强度空间分布，我们在 READ 法的基础上提出了一种求解 RTE 的新方法－DRESOR 法。作为DRESOR 法求解RTE 系统研究的一部分，本文主 要运用该方法对具有真空边界的一维灰性平行平板介质中的RTE 进行求解，其中，介 质具有发射、吸收、各向异性散射特性。为对各向异性散射进行处理，同时加快能束 追踪的计算过程，本文提出了散射分布因子的概念，将以入射方向和散射方向间夹角 为变量的相函数转化为可以直接利用的以系统坐标系下的角度为变量的系数矩阵。对 不同光学厚度、不同散射率及不同相函数下的RTE 进行求解，求解结果代入RTE 积 分形式精确解中进行验证计算。验证计算结果和 DRESOR 法求解结果的对比证明了 该方法的正确性。给出了具有较高方向分辨率的辐射强度空间分布、介质内热源、计 算节点热流密度及投射辐射分布。从计算结果可以看出，光学厚度、散射相函数和散 射率均对辐射强度分布有影响，从而对热流密度、投射辐射及介质内热源均有影响。 对相同的计算节点，在相同的角度下辐射强度在不同相函数作用下表现出不同的分布 规律，辐射强度分布在小光学厚度下受相函数影响较小对辐射强度分布影响较小，且 在散射率发生改变时不同相函数对辐射强度分布的影响相对来说也比较小。 关键词：辐射强度 DRESOR 法 各向异性散射 散射分布因子 平行平板介质 I Abstract The radiative heat transfer is an important heat transfer mode ,and solving RTE is a fashional aspect in the studies of radiative heat transfer. The RTE governing the propagation of radiative intensity in participating media is an integro-differential equation, and the media is always with absorption, emitting and anisotropic scattering characters, so solving RTE is complicated. Hundred years later, many methods for calculating radiative transfer come forth. The four primary methods are Discrete Ordinate Method, Spheric Harmonics, Zonal and Monte Carlo Method. MCM can consider all the major influences in the procedure of simulating radiative heat tra