蒙特卡洛法的研究随机表面光散射特性.pdfVIP

中国工程热物理学会 传热传质学 2004年学术会议 043076 蒙特卡洛法研究随机表面光散射特性 刘晓红阮立明 (哈尔滨工业大学能源科学与工程学院‘，黑龙江，哈尔滨，150001) Teh0451 E-Mail：ruanlm(础it．edu．cn 摘要 本文采用蒙特卡洛(Mc)数值模型生成高度起伏满足正态(Gauss)分布的随机粗糙表面． 在Ki。。hhoff近似的基础上，利用M—萌蟊；F算了一维随机粗糙表面的光教射，描述相关散射蒋翟：_采 用射线踪迹法有效地解决了遮挡效应和多次散射问题．并分析了表面轮廓统计参数对散射特性的影响． 关键词 随机表面，光散射特性，蒙特卡洛法 1．引言 粗糙面散射在光学、电磁学与声学等领域都有十分重要的研究与应用，比如光学界 面特性，电磁散射与波传播，海底声波探测等，这些都已成为现代通信、遥感、目标识 别、环境系统监测等不同学科共同感兴趣的课题，近年来取得了一系列重要的进展…。 研究粗糙表面散射的方法有Kirchhoff近似法，微扰法及全波法等，可以较好地解 释和预测随机表面的散射特征。微扰法只适用于弱粗糙表面，全波法的数学过程过于复 杂难以得到数值解。Kirchhoff近似的基本思路吲是：将粗糙表面视为斜率随机分布的平 面小面元的集合，将小面元近似看成光滑平面，在光与面元相互作用时只考虑镜反射， 不考虑其它方向的散射。在每个小面元上用入射场和反射场的和来近似实际的总场。它 的物理图像简单清晰，在数学上易于处理。更重要的是，它可以明确地将单次和多次散 射分离出来，有助于分析与多次散射直接相关的现象。这是微扰法和全波法难以实现的。 2．蒙特卡洛法模拟生成符合Gauss分布的粗糙表面 生成随机变量序列的基本思想pJ：由指定的分布函数出发，产生出服从该指定分布 特性的随机变量样本值序列。具有Gauss分布的随机变量的概率密度分布如下式所示I．1： 1螂 厶(x)=—杀一e驴 (1) q二孔。 令 Z=研+叮·Z (2) 则Z为N(O，1)分布的随机变量。由于不能直接求出积分分布函数的表达式，本文通过 【O，1】分布的随机变量转换得到Gauss分布的随机变量．具体过程参考文献【4】。 3．蒙特卡洛法模拟粗糙表面的散射特性 本文应用Kirchhoff近似计算粗糙表面的散射强度．但经典的散射场积分方程在计 算多次散射时，随机变量数目增加，积分收敛速度很慢．为简化计算而引入的近似不可 避免地降低了结果的精度。本文采用MC模拟来解决计算中的精度和效率问题15l。 另外，光线在系统内传播的过程中，遮蔽效应不容忽视．随机表面上遮蔽效应的分 析本身很复杂，早期的文献中大多采用一些经验公式嘲来描述遮挡效应，虽然比较简单， 774 但都只适用于一定的范围，不具有通用性。数值方法计算比较复杂，但它的出发点是粗 糙面的实际高度起伏，包含了粗糙面的所有特征。因此数值方法适用于任何高度起伏已 知的粗糙面，比经验公式更符合实际情况。本文采用射线追踪法【7l模拟一定入射角和散 射角下的遮挡过程，沿光线方向逐点考察阻挡情况，并在计算中去除那些被遮蔽的面元。 圈1一维随机粗糙表面一次和多次散射的几何关系 4．计算结果 形貌的粗糙表面计算了不同入射角度下散射强度在空间上的分布情况。通过计算散射方 向图得到了散射强度分布与表面的几何特征参数的关系。井对比和分析了同一粗糙表面 的一次散射和多次散射的结果。以下所有结果均达到与迭代次数无关的精度。 5．结果分析 不考虑波长、入射角度及入射光的偏振状态对散射