咨询热线:15190210002

三维仿真       

基于蒙特卡罗的光散射传输三维动态仿真

摘  要: 以典型偏振蒙特卡罗模型为基础,采用斯托克斯-穆勒(Stokes-Mueller)形式描述偏振光散射传输过程,利用VC和OpenGL工具设计实现了光散射传输过程的三维动态仿真系统,并对系统仿真的数据结果进行了分析。结果表明,此仿真系统可以直观地在三维空间中查看光子在介质中的散射传输过程和接收器的累积过程;且由仿真系统所得的数据结果分析,偏振信息具有一定的对称性。
关键词:蒙特卡罗; 三维仿真; 偏振特性; 斯托克斯

天空中存在许多不同偏振度、不同偏振化方向的散射光,便形成了特定的包含大量偏振信息的偏振分布模式[1],这种偏振分布模式蕴含丰富的方位信息和大气参数信息。这种偏振特性而主要是由于光在传输过程中与介质中粒子散射碰撞而产生的,研究光在散射介质中的传输和散射过程对偏振导航、探测都有重要的意义[2-3]。
    光在散射介质中的传输是一个非常复杂的过程。典型的传输理论为辐射传输理论,该理论一般难以得到解析解。对散射介质中的光传输问题也可以用统计方法进行描述。目前解决光线在大气中的辐射传输问题比较有效的方法是蒙特卡罗方法,其思路清晰,能方便地应用到实际计算中。WILSON和ADAM首次将蒙特卡罗思想引入激光同组织相互作用方面[4],随后,蒙特卡罗模拟方法就被广泛地应用到光在组织中的传输模拟。SEBASTIAN B应用斯托克斯-穆勒(Stokes-Mueller)形式模拟了偏振光在半无限介质中的传输,将蒙特卡罗模拟方法的应用扩展到了偏振领域[5]。JESSICA总结并给出了半无限介质中的偏振光传输的三种蒙特卡罗模拟方法[6],蒙特卡罗方法得到的结果依赖于对大量光子的统计,由于模拟得到的结果与实验结果比较符合,因此蒙特卡罗方法被越来越多的科研工作者所采用。
    本文采用蒙特卡罗统计方法模拟光束在散射介质中的传输和散射过程。同时,利用开源图形函数库OpenGL工具在VC++环境下设计实现了三维动态仿真,提供了三维空间中光散射传输的可视化系统,并对光散射传输后的偏振特性进行了分析。
1 偏振蒙特卡罗模型
    英国物理学家Stokes于1852年提出利用一种四维向量(即Stokes矢量),来描述光的偏振状态,由于其具有可加性,因此非常适用于对光线的分析工作。用Stokes矢量来描述光的偏振状态:S=[I Q U V]T,其中,I表示光术的总强度,±Q反映线偏振光水平或垂直的程度;±U反映线偏振光与Q成 ±45° 的程度;±V反映光束中圆偏振光的左旋或右旋的程度。其偏振度可以表示为:


    散射介质层可以看作是一个无限大但厚度有限的面板,在本文中厚度取4/us(us为散射系数),介质层中为同一介质。如图1所示,沿着Z轴正方向从Z=0(即X-Y面)到Z=4/us面为介质层,可以认为在发射点O有一个光子发射器沿着Z轴正方向发射光子,在Z=4/us面上放置一个正方形接收器(图中灰色区域所示),本文中接收器的边长取14/us。为了较真实地反应一束光穿过介质层之后偏振态的变化情况,在本文中,如果没有特殊说明,取1e6个光子连续发射。接收器取100×100个方格点采样接收。
    蒙特卡罗方法的算法流程图如图2所示,其中灰色步骤只有在偏振蒙特卡罗模型中需要。

下面对蒙特卡罗方法的关键步骤做简要解释(可参考JESSICA总结的偏振蒙特卡罗子午面模型[6]):
    (1)初始参考平面为X-Z平面;
    (2)启动时刻光的E矢量是在X-Z平面内定义的;
    (3)根据已知的散射物质和拒绝方法的相位函数,随机确定角度?琢(散射角度)和角度?茁(入射平面和散射平面的夹角)。
    (4)根据散射参数确定穆勒矩阵计算出射光;
    (5)在光子散射出介质层后,若打在接收器上,进行累计至所有光子散射完。
    另外,当光子其权重低于阈值时,光子运动终止。
2 光散射传输的三维动态仿真系统设计实现
    按照上述的偏振蒙特卡罗模型,利用OpenGL图形库和VC++工具设计实现光在介质中散射传输的三维动态仿真系统,图3所示为仿真系统的功能框图。



苏ICP备000214454号