辐射传输介质散射相函数的蒙特卡洛算法.pdf

1910 光谱学与光谱分析 第31卷 LO 1．O 0．8 O．8 O．6 0．6 口oI苗日届o∞量山 O．4 0-4 口。召u_县Q甓‘山 O．2 O．2 O 0 0 加80 120 160 0 40 80 120 160 Scatteringangle Scatteringangle volume fuctions诵mdifferentradiuses 飚1Controlscatteringphase partical Ca)：r=3．5脚；(b)：r=5／zm；(c)：r=8．2p“ (d)：r=lO．6um；(e)：r=12．5tan；(f)：r=15tan 影响加剧，计算中必须考虑其影响；当数密度较小时，MCM 性的算法。基于该算法得出的辐射介质散射相函数能更准确 计算的相函数与Mie相函数很接近，说明此时可按独立散射 的描述多重散射带来的影响： 计算。 (1)当控制体内粒子数少于1×109时，复分散系可以当 当人射波长确定，改变复分散系粒子半径时，会导致相 作单次散射计算。控制体的相函数随着粒子数密度的增大而 函数的变化。由于粒子尺度参数：[一2矿A，实际改变了x，趋于平缓，当数密度大于l×1011时需考虑多重散射的影响。 而不对称因子g与消光因子Q。都是z的函数，因此使得相 当水雾粒子半径较小时，控制体内需要更多的粒子数才能体 函数变化。根据图1(a)一(f)的变化规律，当粒子半径较小时现出所重散射的影响。若将水雾用于红外隐身，通过提高数 须有更的粒子数密度才能体现出多重散射的影响，而粒子半 密度加强多重散射，可有效降低透过率。 径较大时，较小的粒子数密度便可体现处多重散射的影响。 (2)只要能确定复分散系控制体的曲面方程，采用MCM 根据计算，粒径3～5“m时粒子数N1×1011时最好按多可以计算任意形状的控制体相函数。 (3)该算法的不足之处在于当提高样本容量时会大大增 重散射计算；当粒子半径8～10．6肛m时，粒子数N5×109 需按多重散射计算；当粒子半径12～15 加计算量，准确性与计算时间成为一对矛盾。因此选取合适 gxn时需按多重散射 计算。