Mie散射系数的新算法

王少渄�/span>任中京　张希明　何芳江海�?/span>

(山东建筑材料工业学院，济南　250022)

提要９�/span>介绍了一种计箖�/span>Mie散射系数的新方法，给出了计算实例、�/span>

关键词：米氏散射＋�/span>Mie系数＋�/span>Mie计算

New algorithm of Mie scattering parameters

Wang Shaoqing Ren Zhongjing　Zhang Ximing 　HeFang Jiang Haiying

(Shandong Institute o f Building Materials＋�/span>Jinan　250022)

Abstract９�/span>A new algorithm of calculating Mie Parameters was introduced. Some calculation results done with this algorithm were given.

Key words:Mie scattering＋�/span>Mie parameters＋�/span>Mie calculation

1　引言

Mie理论是球形颗粒对单色光的散射场分布的严格解析觢�/span>^[1]，目前在环保、动力、气象、天文、两相流及粉体颗粒尺寸分布测量等工程技术领域中有广泛的应用。利用单一颗粒或颗粒群光散射场的测量数据。可以反推得散射颗粒或颗粒群的很多物理性质，如颗粒的尺寸、颗粒的折射率等^[2]。但反推必须事先计算出各种尺寸的颗粒在各种复折射率下的散射场分布数据、�/span>1968平�/span>Dave^[3]最先发表了完整皃�/span>Mie散射计算方法，以名�/span>Lentz^[4]咋�/span>Wiscomb^[5]又针对部分计算提出了新算法。国内也有人发表了他们自己的算法^{[6�?/span>7�?/span>8]}。但总的看来，这些算法均有各自的局限性。尤其是当颗粒尺寸或折射率的虚部值较大时，往往计算速度过慢或产生溢出和不收敛等现象。本文介绍作者发展的Mie散射新算法。该算法的特点是不受颗粒尺寸及折射率的限制，不会产生溢出和不收敛的现象，且具有较快的计算速度、�/span>

2　Mie散射系数的计算公弎�/span>

Mie散射计算的中心问题是计算Mie散射系数a_n咋�/span>b_n，其表达式为^[9]

其中为颗粒的尺寸参数，定义为α=πd/λ＋�/span>d为颗粒直径，λ为入射光在颗粒周围介质中的波长，耋�/span>m为颗粒在周围介质中的相对复折射率，即m= m¹+ im²( m²< 0) ,式中i为虚数单位。而�?/span>_n( Z)和ν�/span>_n( Z) ( Z既表礹�/span>α，又表示mα)的表达式丹�/span>

Ψ_n( Z)=＇�/span>πZ/2）J_n+1(Z)

ξ_n( Z)=Ψ_n( Z)+iΧ_n(Z)

Χ_n(Z)= ̶ (πZ/2)( ̶ 1)_n-1J_-(n-1)(Z)

=̶ (πZ/2)N_n+1(Z)

3 计算中产生溢出的原因

计算Mie散射系数须先计算Ψ_n咋�/span>Χ_n。一般采用递推的方法。递推又分为向前递推(即从n= 0开姊�/span>)与向后递推(即从n= N开始至n=0, N为预先设定倻�/span>)。实验表明，向前递推总是快于向后递推。�?/span>_n咋�/span>Χ_n的初值为

分析以上两式可知，当m₂≟�/span>0旵�/span>,若颗粒尺寷�/span>d很大,或复折射率的虚部倻�/span>m₂很大，将使得乘积m₂d很大，可使两式中的项exp( - m₂α) =exp ( - πm₂d/λ)的值超过计算机的数据限，从而产生溢出，这是产生溢出的重要原因。另外，在递推过程中，不恰当的算法也可能造成溢出、�/span>

4　a_n咋�/span>b_n的新算法

为解决上述问题，作者提出了新的算法。将a_n咋�/span>b_n的公式变形如上�/span>:代�/span>

其中Lnr�?/span>Lnj分别表示Ln(m)的实部与虚部。将＇�/span>2）式代入＇�/span>1）式，并�?/span>a_nr�?/span>a_nj咋�/span>b_nr�?/span>b_nj分别表示的实部与虚部。如此可推得

在上述四个公式中采用比值的形式是非常重要的，这样可避免递推过程中当a_i不�/span>b_i较大时乘法运算可能产生的溢出。这是本算法的一个重要特点。在以上四式�?/span>

由于

均为实变量函数，计算不产生溢出。关键是 咋�/span>

的算法如何处理，才能保证计算中不产生溢出。在Lentz的算法中是采用连分式计算Ln的值，其精度的保证是由在大量计算基础上得出一个截断项�?/span>N与参�?/span>a不�/span>m的经验公式而实现的。这样的经验公式，一则有实用上的局限性，再则也会带来截断误差。文献（6）对此经验公式做了改进，但仍陷于a=1~100＋�/span>m1=1~2＋�/span>m2=0~1的范围。下面介绍本文作者发展的关于Ln的新算法。该算法的特点是不受a叉�/span>m值的限制，不会产生溢出或不收敛等病态现象，且具有较快的计算速度。令

上面导出皃�/span>(3) - (20)诸式，构成了Mie系数an咋�/span>bn的完整算法。由亍�/span>an咋�/span>bn是从n= 1开始计算，利用初值公弎�/span>(16) - (20)即可算得任意级数皃�/span>an咋�/span>bn的值，故没有舍入误差的问题。从(16)式可见，因为y= m₂ɑ≣�/span>0,故无讹�/span>m₂咋�/span>ɑ取何值均不会产生溢出，再加上(3)的各式中采用了比的形式，又避免了计算过程中的溢出，这就从根本上解决了溢出的问题、�/span>

5　计算实例

利用以上算法编制了计算颗粒散射场强度和消光系数的计算机程序。当波长为λ的单位振幅的平面自然光入射颗粒时，颗粒的散射光强为^[9]

Ξ为计算机在双精度下的最小数据限、�/span>

国�/span>1示出散射光强的一组计算实例。其中取m₁=1.33＋�/span>m₂=-0.4，κ�/span>=0.6328三图分别对应于颗粒的直径丹�/span>d=0.001,1.0咋�/span>30μm。耋�/span>d图为颗粒直径d=100μm时散射花样的局部放大。可见，随着颗粒尺寸的增大，前向散射迅速加强，并且出现了复杂的旁瓣、�/span>

实部(a)与虚�?/span>(b)的变化情况。可见随着m₁咋�/span>m₂的增大，虽然颗粒的尺寸保持不变，但散射也加强，且后向散射随着m₁咋�/span>m₂的增大而加强、�/span>

国�/span>3给出有关消光系数的计算结果，其中a)咋�/span>b)分别表示消光系数随折射率实部与虚部的变化情况。可见随着颗粒直径的增大，消光系数趋近亍�/span>2;折射率的增大，尤其是折射率虚部的增大，使这一趋近变得更快和更明显。另外当折射率的虚部m₂= 0时，消光系数随颗粒直径的增大而振荡；但当m2≟�/span>0时，振荡迅速消失、�/span>

光强最大值处所对应皃�/span>F_M(Z)的值即由上式确定。由上式也可见：此极限情况下的焦移大小主要由S₀/f叉�/span>N_a所决定、�/span>

参考文献：

[1]M. Bor n and E. Wolf, Principles of optics ( sixt hedition ) , Pergaman Press ( OXFORD, NEW YORK, PARIS) , pp611(1980)

[2]Peter Chylek, V. Ramasw amy , A. Ashkin and J.M . Dziedzic,‛�/span>Simultaneous determination o f refractiv endex and size o f spherical dielectric particles from light scattering data“�/span>,Appl. Opt . ,Vo l. 22, No . 15, p2302- 2307(1983)

[3]J. V. Dave, Repo rt 320 - 3237 ( IBM Scientific center , 1968)

[4]W. J. Lentz, Appl. Opt. , 15, 668( 1976)W. J. Wiscome, Appl. Opt, Vol. 19, No. 9, pp1505( 1980)

[6]顾冠亮等,有关光散射物理量的数值计箖�/span>,上海机械学院学报, 1984平�/span>4朞�/span>pp21

[7]余其铮等, Mie散射算法的改迚�/span>,哈尔滨工业大学学�?/span>, 1987平�/span>3朞�/span>, pp21

[8]郑刚筈�/span>, M ie散射的数值计箖�/span>,应用激先�/span>, Vol.12, No . 5, pp220( 1992)

[9]H. C. v an de Hulst , Lig ht scatter ing by small particles,Do ver Publica tion lnc. , New Yor k, chp. 13( 1981)