偏振辐射传输的蒙特卡洛方法求解.docVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：123456 偏振辐射传输的蒙特卡洛方法求解基金项目：国家自然科学基金 (No. 基金项目：国家自然科学基金 (No.51176040)，新世纪优秀人才支持计划 (NCET-09-0067) 贲勋，易红亮，谈和平 （哈尔滨工业大学，哈尔滨，150001） (Tel: 0451 Email: yihongliang@hit.edu.cn) 摘 要 辐射的偏振特性有重要应用，与标量辐射传递相比，其求解更加困难。本文论述了偏振辐射传输的基本原理，给出了采用蒙特卡洛方法（MCM）求解偏振辐射传输的基本思路，着重介绍了如何采用拒绝法（rejection method）确定Stokes矢量的散射方向；提出了矢量辐射传递系数（vector radiation transfer coefficient）的概念。结合一维算例，验证了该方法应用于Rayleigh散射辐射传递模型的正确性。 关键词 偏振辐射；Stokes矢量；蒙特卡洛法；拒绝法 0 前言 可见光、红外辐射等本质上都是电磁波，其横波的特性决定了光波具有偏振性 REF _Ref325185095 \r \h [1]。光的偏振特性在国民经济、科学研究以及国防等领域有着广泛的应用。对于偏振信号的捕捉识别与模拟，正在成为辐射传输领域研究的重要方向。偏振辐射传输在数值模拟中已有多年的探索，离散坐标法 REF _Ref326000522 \r \h [2]、谱元法 REF _Ref326000593 \r \h [3]、有限体积法 REF _Ref326000663 \r \h [4]等方法被应用其中。由于偏振辐射传输方程的复杂性，使得基于方程离散的方法进行的研究表现出一定的局限性，尤其针对非均匀介质更是如此。蒙特卡洛法（MCM）在辐射传输研究中有着广泛的应用，对于传统的标量辐射传输方程(SRTE)求解，该方法已经发展得相当成熟 REF _Ref326000715 \r \h [5]，并表现出精度高、易于适应复杂几何形状的优越性。目前为止，将MCM方法应用在求解矢量辐射传输方程（VRTE）中的研究还较少，且体系不够完整、成熟。将MCM高精度、适应性强的优势应用在VRTE的求解中有重大的研究价值的和意义。 1 蒙特卡洛法求解偏振辐射传输基本原理 1.1 托克斯矢量的定义 描述光的偏振性的方法很多，例如琼斯(Jones)矢量法、三角函数表示法、邦加(Poincare)球图示法等，本文采用斯托克斯矢量法来描述光的偏振特性。这种方法可以准确的描述出光的强度及偏振状态 ，对于完全偏振光、部分偏振光和非偏振光都适用。Stokes矢量具有可加性，一束由许多不同偏振状态的独立光线组成的光的Stokes矢量，是每个不同偏振状态的独立光线Stokes矢量的累加。斯托克斯矢量定义如下： (1) 其中表示强度，表示偏振光中垂直或水平方向偏振的比例，表示偏振光中与垂直或水平方向呈方向偏振的比例，表示偏振光中圆偏振的比例。各分量的定义及关系为： (2) 上式中的表示在频率相关的间隔内取时间平均，上标*表示共轭复数。、为垂直于辐射传播方向的电矢量的两个正交分量。 1.2 Stokes矢量的矩阵化运算规则 Stokes矢量是一个的矩阵，为了清晰准确的描述Stokes在传输过程中发生的旋转、散射等过程，通常采用矩阵运算来表示上述过程。 1.2.1 旋转矩阵 Stokes矢量相对于某一参考面定义的，在传输过程中，为了计算方便，经常要改变参考面的选择，这样就需要将Stokes矢量由一个参考面旋转到另一个参考面，其变换规则如下： (3) 其中代表旋转到新的平面上的Stokes矢量，代表沿着光传播的方向看过去，顺时针旋转的的角度 REF _Ref326000793 \r \h [6]（逆时针旋转为）。为Stokes矢量的旋转矩阵。 1.2.2 散射相矩阵 辐射在传输过程中遇到粒子会发生散射，Stokes矢量会发生改变，在数学上表示是通过与散射相矩阵相乘实现的。 散射相矩阵是一个的矩阵，通常来说散射相矩阵由散射角和粒子物性决定的。对于常见的Rayleigh和Mie散射来说，散射相矩阵可以分别化简为式(4)、(5)