米散射理论基础.docx
- 文档编号:9523591
- 上传时间:2023-02-05
- 格式:DOCX
- 页数:13
- 大小:21.80KB
米散射理论基础.docx
《米散射理论基础.docx》由会员分享,可在线阅读,更多相关《米散射理论基础.docx(13页珍藏版)》请在冰豆网上搜索。
米散射理论基础
米散射(Miescattering);又称“粗粒散射”。
粒子尺度靠近或大于入射光波长的粒子散射现象。
德国物理学家米(GustavMie,1868—1957)指出,其散射光强在各方向是不对称的,顺入射方向上的前向散射最强。
粒子愈大,前向散射愈强。
米散射
当球形粒子的尺度与波长可比较时,一定考虑散射粒子体内电荷的三维分
布。
此散射状况下,散射粒子应试虑为由很多齐集在一同的复杂分子构成,它们在入射电磁场的作用下,形成振荡的多极子,多极子辐射的电磁波相叠加,就构成散射波。
又因为粒子尺度可与波长对比较,所以入射波的相位在粒子上是不均匀的,造成了各子波在空间和时间上的相位差。
在子波组合产生散射波的地方,将出现相位差造成的干预。
这些干预取决于入射光的波长、粒子的大小、折射率及散射角。
当粒子增大时,造成散射强度变化的干预也增大。
所以,散射光强与这
些参数的关系,不象瑞利散射那样简单,而用复杂的级数表达,该级数的收敛相当迟缓。
这个关系第一由德国科学家G.米得出,故称这种散射为米散射。
它拥有以下特色:
①散射强度比瑞利散射大得多,散射强度随波长的变化不如瑞利散射那样激烈。
跟着尺度参数增大,散射的总能量很快增添,并最后以振动的形式趋于必定值。
②散射光强随角度变化出现很多极大值和极小值,当尺度参数增大时,极值的个数也增添。
③当尺度参数增大时,前向散射与后向散射之比增大,使粒子前半球散射增大。
当尺度参数很小时,米散射结果能够简化为瑞利散射;当尺度参数很大时,它的结果又与几何光学结果一致;而在尺度参数比较适中的范围内,只实用米散射才能获得独一正确的结果。
所以米散射计算模式能宽泛地描绘任何尺度参数均匀球状粒子的散射特色。
世纪末,英国科学家瑞利第一解说了天空的蓝色:
在洁净大气中,起主要散射作用的是大气气体分子的密度涨落。
分子散射的光强度和入射波长四次方成反比,所以在发生大气分子散射的日光中,紫、蓝和青色彩光比绿、黄、橙和红色彩光为强,最后综合成效使天穹体现蓝色。
进而成立了瑞利散射理论。
世纪初,德国科学家米从电磁理论出发,进一步解决了均匀球形粒子的
散射问题,成立了米散射理论,又称粗粒散射理论。
质点半径与波长靠近时的散射,特色:
粗粒散射与波长没关,对各波长的散射能力相同,大气较浑浊时,大气中悬浮许多的的尘粒与水滴时,天空呈灰白色。
米散射理论是由麦克斯韦方程组推导出来的均质球形粒子在电磁场中对平面
波散射的精准解。
一般把粒子直径与入射光波长相当的微粒子所造成的散射称为
米散射。
米散射合适于任何粒子尺度,不过当粒子直径相对于波长而言很小时利
用瑞利散射、很大时利用夫琅和费衍射理论就能够很方便的近似解决问题。
米散
射理论最早是由G1Mie在研究胶体金属粒子的散射时成立的。
1908年,米氏经过电磁波的麦克斯韦方程,解出了一个对于光散射的严格解,得出了任意直径、任意成分的均匀粒子的散射规律,这就是有名的米氏理论[4-
6]。
依据米散射理论,当入射光强为I0,粒子四周介质中波长为λ的自然光平行入射到直径为D的各向同性真球形粒子上时,在散射角为θ,距离粒子r处的散
射光和散射系数分别为:
从上式中能够看到,因为是各向同性的粒子,散射光强的散布和φ角没关。
同时,
1
上式中:
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为散射光强;为散射角;为偏振
2
光的偏振角。
式中:
S1()和S2()是振幅函数;an和bn是与贝塞尔函数和汉克尔函数相关的函数;n和n是连带勒让得函数的函数,仅与散射角相关。
此中
式中:
n()和n()分别是贝塞尔函数和第一类汉克尔函数;n()和
n()是n()和n()的导数;为无因次直径,
D,D为颗粒的
实质直径;是入射光的波长;m是散射颗粒相对于四周介质的折射率,它是一
个复数,虚部是颗粒对光的汲取的量化。
由以上公式可见,Mie散射计算的重点
是振幅函数S1()和S2(),它们是一个无量乞降的过程,理论上没法计算。
求解振幅函数的重点是计算an和bn,所以Mie散射的计算难点是求解an和bn。
Mie散射理论的数值计算
经过以上剖析可知,Mie散射计算的核心是求解an和bn,我们编制程序也
是环绕它进行编写。
在an和bn的表达式中n(),n(),n()和n()
知足以下递推关系:
这些函数的初始值为;
3
与散射角相关的n()和n()知足以下递推公式:
有了这些递推公式能够很方便地经过计算机程序求解。
可是对于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。
则散射光垂直于散射面和平行于散射面的两
个重量的振幅函数为:
4
以上式中:
Jn+1/2(z)和Yn+1/2(z)分别为半整数阶的第一类,第二类贝塞尔函数。
(1)
让德函数。
在数值模拟过程中选用初始下:
5
微粒子对光的散射和汲取是电磁波与微粒子互相作用的重要特色,而微粒对电磁辐射的汲取与散射与粒子的线度有亲密关系,对于不一样线度的粒子一定应用不一样的散射理论。
Mie散射理论主要用于从亚微米至微米的尺寸段;在微米以下至纳米的光散射则近似为形式更清晰简单的瑞利散射定律,散射光激烈依靠于光波长λ(I~λ-4);而对大于微米至毫米的大粒子则近似为意义明确的夫朗和费衍射规律了。
Mie散射理论给出了球型粒子在远场条件下的散射场振幅an、bn以及粒子内部电磁场振幅cn、dn的计算表达式,往常称为Mie散射系数
式中m表示微粒子外面介质的相对折射率,x=κa,a为球的半径,κ=2π/λ
称为波数,μ为相对磁导率,即球的磁导率与介质磁导率的比值,jn(x)和h
(1)n分
(x)
别为第一类虚宗量球Bessel函数和Hankell函数。
散射系数,消光系数及偏振状态下散射相位函数:
6
散射截面σsca(散射率Qsca)、汲取截面σabs(汲取率Qabs)、消光截面σext(消光
率Qext)、后向散射截面σb(后向散射率Qb)以及辐射压力σpr(辐射压力效率Qpr)。
其表达式以下:
此中i为sca、abs、ext、pr分别表示散射、汲取、消光、辐射压力。
依据能量守恒定律有:
Qpr(辐射压力效率的计算公式):
Qb(后向散射系数):
这些都是无量级数乞降,在实质计算过程中一定取有限项,Bohren和Huffman给出了级数项最大值弃取的标准:
对于单位振幅入射波经微粒散射后,其散射场振幅的大小与散射角相关,在球坐标系下,远场散射振幅的大小为:
此中S1和S2为散射辐射电场在垂直及平行于散射面的两个偏振重量。
微球内部场振幅计算公式
颗粒内部电场强度为:
7
此中M
(1)o1n和N
(1)e1n为矢量波球谐函数,在球坐标系中定义以下:
汲取截面Qabs
拥有消耗介质颗粒的汲取截面为:
此中ε″是粒子相对介电常数的虚部,经整理可得:
式中mn、nn为:
实质上由Mie散射理论可知,上式中的积分项为电场强度的平方对角度θ、φ全空间积分的均匀值,即:
于是汲取效率为:
8
式中x′=rk=z/m。
当xn1时即瑞利散射状况,颗粒的内部均匀场强为常数,其值
为:
ImprovedMiescatteringalgorithms
Mie计算存在的问题就是如何最有效地结构Mie计算,同时保证正确性和防止数值的不稳固性和病态。
Mie计算以耗时著称,第一无量项级数N的乞降,比如:
m的水滴在0.5m的可见光散射状况下,大概需1260项乞降。
其次,典型
的计算都希望能对一系列半径(如对尺寸散布求积分)、一系列波长(如对太阳光谱求积分)及一系列折射率乞降(如经过散射参量反推折射率)。
当折射率虚部mIm很大时,用向后循环法求An很不稳固。
而向前递推老是稳固的(但向后递推安全时,老是优先选择,因为其计算速度很快)。
得出同意向后递推的经验标准:
用正确的向前地推与相对应的向后地推做比较,当发现对
和g的相对偏差超出10-6
时,以为计算失败。
对于一对确立的
(x,mRe),我们采纳向后递推找寻第一个循环失败的
研究表示:
对于确立的
,,
的值跟着x的增添很快趋势于一个确立值。
对
9
假如在任意角度下S1、S2的实部和虚部的相对偏差超出105时,以为对S1和
S2的向后递推失败。
(而此时,QscaQext其实不受影响,因为当S1,S2的相对误
差达到105时,QscaQext的相对偏差总保持在1010以下。
)
对S1和S2
对散射强度和偏正度
连分式算法总结:
10
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(α),依据贝塞耳差积公式:
11
由以上二式整理得:
上式中Cn的计算是采纳近似于Dn的连分式的形式,计算中可调用同一函数计算。
若已知初值:
这样便可计算出各级Jn和Yn。
WilliamJ.Lentz对于连分式的文章:
此中。
以为基础,采纳贝塞尔函数比值的连分式表示法:
,利
用此法可产生全部的,只管耗时,但能减少储存需求。
同时可经过计算
高阶值,使用下边的递推公式,从后往前算出其余值。
不像一般的函数,贝塞尔函数的比值一旦超出可控制的界限,就不再增添,初始
的高阶值决定了全部低阶值的正确性,所以,采纳新方法计算正确的初始比值是必需的。
12
处于分母地点的+号表示分母上加上一个特别的连分式。
近似于上式中的表示
形式。
定义一种新的符号:
Lentz给出了n阶部分收敛值为:
比如:
实变量,虚数计算过程:
米散射学习当前所碰到的困难:
究竟如何的计算结果才算正确,如何能找到一个米散射计算结果正确又有效的数据库,来考证自己算法及程序的正确性。
倒退式算法的总结:
13
Dn的计算采纳Dave的倒推式:
因为Dn函数有很强的收敛性,对于Dn的倒推计算的初值的选用有很强的任意性。
因为当n→∞时Dn(mα)→0,所以能够取0作为初值。
倒推起点选用大一
些,能够保证Dn函数的收敛完整,可是同时却增添了计算时间。
所以一定选用一个最正确的选择标准。
经过试算,作者以为最正确的上限为
这里m1是复折射率的实部.
相同,对于贝塞耳函数Jn的计算也能够用倒推的方法计算产生:
上式是一个一般的Jn的递推式,知道了Jn和Jn-1,能够顺利地计算出全部的J
n序列值。
为了防止计算Jn的繁琐而又能发挥递推式的迅速的长处
采纳下边的
方法:
假定N→∞时,取某一个递推初始值为:
JN*()0,JN1*(
),
此中ε是一个很小的数,如可取10-6。
将初值代入上式,就能够算出全部的J*。
察看同一自变量的J*和J序列,发现它们对应项之间有固定的倍数关系。
如定义这个倍数为β,那么
因为J1(α)的计算是特别便利的(J1=sinα/α2-cosα/α),所以β=J1/J1*,计算出Jn*(α)能够算出Jn(α)。
和Dn的计算相同,Jn的倒推开端点的公式为:
对于贝塞尔函数的倒退过程在另一文件中的描绘:
14
利用初始值
Jn
sin
Jn
J0
sin
0
Jn
cos
sin
0
J1
15
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 散射 理论基础