amie米散射理论基础Word文件下载.docx

amie米散射理论基础Word文件下载.docx

《amie米散射理论基础Word文件下载.docx》由会员分享，可在线阅读，更多相关《amie米散射理论基础Word文件下载.docx（10页珍藏版）》请在冰豆网上搜索。

③当尺度参数增大时，前向散射与后向散射之比增大，使粒子前半球散射增大。

当尺度参数很小时，米散射结果可以简化为瑞利散射；

当尺度参数很大时，它的结果又与几何光学结果一致；

而在尺度参数比较适中的范围内，只有用米散射才能得到唯一正确的结果。

所以米散射计算模式能广泛地描述任何尺度参数均匀球状粒子的散射特点。

19世纪末，英国科学家瑞利首先解释了天空的蓝色：

在清洁大气中，起主要散射作用的是大气气体分子的密度涨落。

分子散射的光强度和入射波长四次方成反比，因此在发生大气分子散射的日光中，紫、蓝和青色彩光比绿、黄、橙和红色彩光为强，最后综合效果使天穹呈现蓝色。

从而建立了瑞利散射理论。

20世纪初，德国科学家米从电磁理论出发，又称粗进一步解决了均匀球形粒子的散射问题，建立了米散射理论，粒散射理论。

质点半径与波长接近时的散射，特点：

粗粒散射与波长无关，对各波长的散射能力相同，大气较混浊时，大气中悬浮较多的的尘粒与水滴时，天空呈灰白色。

米散射理论是由麦克斯韦方程组推导出来的均质球形粒子在电磁场中对平面波散射的精确解。

一般把粒子直径与入射光波长相当的微粒子所造成的散射称为米散射。

米散射适合于任何粒子尺度,只是当粒子直径相对于波长而言很小时利用瑞利散射、很大时利用夫琅和费衍射理论就可以很方便的近似解决问题。

米散射理论最早是由G1Mie在研究胶体金属粒子的散射时建立的。

1908年,米氏通过电磁波的麦克斯韦方程,解出了一个关于光散射的严格解,得出了任意直径、任意成分的均匀粒子的散射规律,这就是着名的米氏理论[4-6]。

根据米散射理论,当入射光强为I0,粒子周围介质中波长为λ的自然光平行入射到直径为D的各向同性真球形粒子上时,在散射角为θ,距离粒子r处的散射光和散射系数分别为:

从上式中可以看到,因为是各向同性的粒子,散射光强的分布和φ角无关。

同时,上式中:

i1、i2为散射光的强度函数;

s1、s2称为散射光的振幅函数;

a为粒子的尺寸参数（a=πD/λ）;

m=m1+im2为粒子相对周围介质的折射率,当虚部不为零时,表示粒子有吸收。

对于散射光的振幅函数,有:

式中an、bn为米散射系数,其表达式为:

其中:

是半奇阶的第一类贝塞尔函数;

是第二类汉克尔函数;

Pn（cosθ）是第一类勒让德函数;

P

（1）n（cosθ）是第一类缔合勒让德函数。

Mie散射理论

Mie散射理论是麦克斯韦方程对处在均匀介质中的均匀颗粒在平面单色波照射下的严格数学解。

由Mie散射知道,距离散射体r处p点的散射光强为

式中:

为光波波长;

I0为入射光强;

Isca为散射光强;

为散射角;

为偏振光的偏振角。

和

是振幅函数;

an和bn是与贝塞尔函数和汉克尔函数有关的函数;

和

是连带勒让得函数的函数,仅与散射角

有关。

其中

分别是贝塞尔函数和第一类汉克尔函数;

是

的导数;

为无因次直径,

D为颗粒的实际直径;

是入射光的波长;

m是散射颗粒相对于周围介质的折射率,它是一个复数,虚部是颗粒对光的吸收的量化。

由以上公式可见,Mie散射计算的关键是振幅函数

它们是一个无穷求和的过程,理论上无法计算。

求解振幅函数的关键是计算an和bn,所以Mie散射的计算难点是求解an和bn。

Mie散射理论的数值计算

通过以上分析可知,Mie散射计算的核心是求解an和bn,我们编制程序也是围绕它进行编写。

在an和bn的表达式中

满足下列递推关系:

这些函数的初始值为;

与散射角有关的

满足下列递推公式:

有了这些递推公式可以很方便地通过计算机程序求解。

但是对于n的大小,因为计算机不可能计算无穷个数据,所以n在计算之前就要被确定。

　散射理论基础与Matlab实现

若散射体为均匀球体,如图1所示,照射光为线偏振平面波,振幅为E,光强I0,沿z轴传播,其电场矢量沿x轴振动。

散射体位于坐标原点O,P为观测点。

散射光方向（OP方向）与照射光方向（z轴）所组成的平面称为散射面,照射光方向至散射光方向之间的夹角θ称为散射角,而x轴至OP在xy平面上投影线（OP′）之间的夹角φ称为极化角。

观测点与散射体相距r。

根据经典的Mie散射理论,散射粒子的尺度参数为α=2πa/λ,其中a为球形粒子的半径,散射粒子相对周围介质的折射率为m=m1+i*m2。

则散射光垂直于散射面和平行于散射面的两个分量的振幅函数为:

以上式中：

Jn+1/2（z）和Yn+1/2（z）分别为半整数阶的第一类,第二类贝塞尔函数。

P

（1）n（cosθ）为一阶n次第一类缔合勒让德函数;

Pn（cosθ）为第一类勒

让德函数。

在数值模拟过程中选取初始下:

微粒子对光的散射和吸收是电磁波与微粒子相互作用的重要特征,而微粒对电磁辐射的吸收与散射与粒子的线度有密切关系,对于不同线度的粒子必须应用不同的散射理论。

Mie散射理论主要用于从亚微米至微米的尺寸段;

在微米以下至纳米的光散射则近似为形式更明晰简单的瑞利散射定律,散射光强烈依赖于光波长λ（I～λ-4）;

而对大于微米至毫米的大粒子则近似为意义明确的夫朗和费衍射规律了。

Mie散射理论给出了球型粒子在远场条件下的散射场振幅an、bn以及粒子内部电磁场振幅cn、dn的计算表达式,通常称为Mie散射系数

式中m表示微粒子外部介质的相对折射率,x=κa,a为球的半径,κ=2π/λ称为波数,μ为相对磁导率,即球的磁导率与介质磁导率的比值,jn（x）和h

（1）n（x）分别为第一类虚宗量球Bessel函数和Hankell函数。

散射系数,消光系数及偏振状态下散射相位函数：

散射截面σsca（散射率Qsca）、吸收截面σabs（吸收率Qabs）、消光截面σext（消光率Qext）、后向散射截面σb（后向散射率Qb）以及辐射压力σpr（辐射压力效率Qpr）。

其表达式如下:

其中i为sca、abs、ext、pr分别表示散射、吸收、消光、辐射压力。

按照能量守恒定律有:

Qpr（辐射压力效率的计算公式）:

Qb（后向散射系数）：

这些都是无穷级数求和,在实际计算过程中必须取有限项,Bohren和Huffman给出了级数项最大值取舍的标准:

对于单位振幅入射波经微粒散射后,其散射场振幅的大小与散射角有关,在球坐标系下,远场散射振幅的大小为:

其中S1和S2为散射辐射电场在垂直及平行于散射面的两个偏振分量。

微球内部场振幅计算公式

颗粒内部电场强度为:

其中M

（1）o1n和N

（1）e1n为矢量波球谐函数,在球坐标系中定义如下:

吸收截面Qabs

具有损耗介质颗粒的吸收截面为：

其中ε″是粒子相对介电常数的虚部,经整理可得:

式中mn、nn为:

实际上由Mie散射理论可知,上式中的积分项为电场强度的平方对角度θ、φ全空间积分的平均值,即:

于是吸收效率为:

式中x′=rk=z/m。

当xn1时即瑞利散射情况,颗粒的内部平均场强为常数,其值为:

Mie计算存在的问题就是如何最有效地构造Mie计算，同时保证准确性和避免数值的不稳定性和病态。

Mie计算以耗时着称，首先无穷项级数N的求和，例如：

100

的水滴在0.5

的可见光散射情况下，大约需1260项求和。

其次，典型的计算都希望能对一系列半径（如对尺寸分布求积分）、一系列波长（如对太阳光谱求积分）及一系列折射率求和（如通过散射参量反推折射率）。

当折射率虚部mIm很大时，用向后循环法求An很不稳定。

而向前递推总是稳定的（但向后递推安全时，总是优先选择，因为其计算速度很快）。

得出允许向后递推的经验标准：

用正确的向前地推与相对应的向后地推做比较，当发现对

和g的相对误差超过10-6时，认为计算失败。

对于一对确定的（x,mRe），我们采用向后递推寻找第一个循环失败的

研究表明：

对于确定的

，，

的值随着x的增加很快趋向于一个确定值。

对

如果在任意角度下

、

的实部和虚部的相对误差超过

时，认为对

的向后递推失败。

（而此时，

并不受影响，因为当

，

的相对误差达到

时，

的相对误差总维持在

以下。

）

对散射强度

和偏正度

连分式算法总结：

Mie散射计算的核心是计算an和bn

其中ψn（α）=αJn（α）,ξn（α）=αJn（α）+iαYn（α）,Jn和Yn分别是第一和二类贝塞耳函数,α

称为当量直径,α=2πr/λ,r是球形颗粒的真实半径,λ是入射光的波长,m为折射率

式中ρ为函数任一自变量。

贝塞耳函数递推关系式：

Mie散射计算中Jn、Yn、Dn的计算是关键和难点。

对于Dn,我们采用的是Lentz的连分式的算法:

Lentz证明有如下关系:

其中，

。

我们注意到当

所以可以利用上式累积相乘直到满足精度要求。

（可根据精度要求例如10-7来确定所要达到的k值）

对于Jn、Yn的生成本文也采用连分式的算法。

具体方案如下:

令Cn=Jn-1（α）/Jn（α）,根据贝塞耳差积公式：

由以上二式整理得：

上式中Cn的计算是采用类似于Dn的连分式的形式,计算中可调用同一函数计

算。

若已知初值：

这样就可计算出各级Jn和Yn。

WilliamJ.Lentz关于连分式的文章：

其中

以

为基础，采用贝塞尔函数比值的连分式表示法：

，利

用此法可产生所有的

，尽管耗时，但能减少存储需求。

同时可通过计算

高阶

值，使用下面的递推公式，从后往前算出其他值。

不像一般的函数，贝塞尔函数的比值一旦超过可控制的边界，就不再增长，初始

的高阶

值决定了所有低阶值的准确性，因此，采用新方法计算准确的初始比值是必要的。

处于分母位置的+号表示分母上加上一个特殊的连分式。

类似于上式

中的表示形式。

定义一种新的符号：

Lentz给出了n阶部分收敛值为：

例如：

实变量

，虚数

计算过程：

米散射学习目前所遇到的困难：

到底怎样的计算结果才算正确，如何能找到一个米散射计算结果准确又有效的数据库，来验证自己算法及程序的正确性。

倒退式算法的总结：

Dn的计算采用Dave的倒推式:

由于Dn函数有很强的收敛性,对于Dn的倒推计算的初值的选取有很强的随意性。

因为当n→∞时Dn（mα）→0,所以可以取0作为初值。

倒推起点选取大一些,可以保证Dn函数的收敛完全,但是同时却增加了计算时间。

所以必须选取一个最佳的选择标准。

通过试算,作者认为最佳的上限为

这里m1是复折射率的实部.

同样,对于贝塞耳函数Jn的计算也可以用倒推的方法计算产生:

上式是一个普通的Jn的递推式,知道了Jn和Jn-1,可以顺利地计算出所有的Jn序列值。

为了避免计算Jn的繁琐而又能发挥递推式的快速的优点,采用下面的办法:

假设N→∞时,取某一个递推初始值为:

其中ε是一个很小的数,如可取10-6。

将初值代入上式,就可以算出所有的J*。

观察同一自变量的J*和J序列,发现它们对应项之间有固定的倍数关系。

如定义这个倍数为β,那么

由于J1（α）的计算是非常便利的（J1=sinα/α2-cosα/α）,所以β=J1/J1*,计算出Jn*（α）可以算出Jn（α）。

和Dn的计算一样,Jn的倒推起始点的公式为:

关于贝塞尔函数的倒退过程在另一文献中的描述:

利用初始值