版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
...wd......wd......wd...米散射〔Miescattering〕;又稱“粗粒散射〞。粒子尺度接近或大于入射光波長的粒子散射現(xiàn)象。德國物理學家米(GustavMie,1868—1957)指出,其散射光強在各方向是不對稱的,順入射方向上的前向散射最強。粒子愈大,前向散射愈強。米散射當球形粒子的尺度與波長可比較時,必須考慮散射粒子體內電荷的三維分布。此散射情況下,散射粒子應考慮為由許多聚集在一起的復雜分子構成,它們在入射電磁場的作用下,形成振蕩的多極子,多極子輻射的電磁波相疊加,就構成散射波。又因為粒子尺度可與波長相比較,所以入射波的相位在粒子上是不均勻的,造成了各子波在空間和時間上的相位差。在子波組合產生散射波的地方,將出現(xiàn)相位差造成的干預。這些干預取決于入射光的波長、粒子的大小、折射率及散射角。當粒子增大時,造成散射強度變化的干預也增大。因此,散射光強與這些參數的關系,不象瑞利散射那樣簡單,而用復雜的級數表達,該級數的收斂相當緩慢。這個關系首先由德國科學家G.米得出,故稱這類散射為米散射。它具有如下特點:①散射強度比瑞利散射大得多,散射強度隨波長的變化不如瑞利散射那樣劇烈。隨著尺度參數增大,散射的總能量很快增加,并最后以振動的形式趨于一定值。②散射光強隨角度變化出現(xiàn)許多極大值和極小值,當尺度參數增大時,極值的個數也增加。③當尺度參數增大時,前向散射與后向散射之比增大,使粒子前半球散射增大。當尺度參數很小時,米散射結果可以簡化為瑞利散射;當尺度參數很大時,它的結果又與幾何光學結果一致;而在尺度參數比較適中的范圍內,只有用米散射才能得到唯一正確的結果。所以米散射計算模式能廣泛地描述任何尺度參數均勻球狀粒子的散射特點。19世紀末,英國科學家瑞利首先解釋了天空的藍色:在清潔大氣中,起主要散射作用的是大氣氣體分子的密度漲落。分子散射的光強度和入射波長四次方成反比,因此在發(fā)生大氣分子散射的日光中,紫、藍和青色彩光比綠、黃、橙和紅色彩光為強,最后綜合效果使天穹呈現(xiàn)藍色。從而建設了瑞利散射理論。20世紀初,德國科學家米從電磁理論出發(fā),又稱粗進一步解決了均勻球形粒子的散射問題,建設了米散射理論,粒散射理論。質點半徑與波長接近時的散射,特點:粗粒散射與波長無關,對各波長的散射能力一樣,大氣較混濁時,大氣中懸浮較多的的塵粒與水滴時,天空呈灰白色。米散射理論是由麥克斯韋方程組推導出來的均質球形粒子在電磁場中對平面波散射的準確解。一般把粒子直徑與入射光波長相當的微粒子所造成的散射稱為米散射。米散射適合于任何粒子尺度,只是當粒子直徑相對于波長而言很小時利用瑞利散射、很大時利用夫瑯和費衍射理論就可以很方便的近似解決問題。米散射理論最早是由G1Mie在研究膠體金屬粒子的散射時建設的。1908年,米氏通過電磁波的麥克斯韋方程,解出了一個關于光散射的嚴格解,得出了任意直徑、任意成分的均勻粒子的散射規(guī)律,這就是著名的米氏理論[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實現(xiàn)假設散射體為均勻球體,如圖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);而對大于微米至毫米的大粒子那么近似為意義明確的夫朗和費衍射規(guī)律了。Mie散射理論給出了球型粒子在遠場條件下的散射場振幅an、bn以及粒子內部電磁場振幅cn、dn的計算表達式,通常稱為Mie散射系數式中m表示微粒子外部介質的相對折射率,x=κa,a為球的半徑,κ=2π/λ稱為波數,μ為相對磁導率,即球的磁導率與介質磁導率的比值,jn(x)和h(1)n(x)分別為第一類虛宗量球Bessel函數和Hankell函數。散射系數,消光系數及偏振狀態(tài)下散射相位函數:散射截面σ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時即瑞利散射情況,顆粒的內部平均場強為常數,其值為:ImprovedMiescatteringalgorithmsW.J.WiscombeMie計算存在的問題就是若何最有效地構造Mie計算,同時保證準確性和防止數值的不穩(wěn)定性和病態(tài)。Mie計算以耗時著稱,首先無窮項級數N的求和,例如:100的水滴在0.5的可見光散射情況下,大約需1260項求和。其次,典型的計算都希望能對一系列半徑〔如對尺寸分布求積分〕、一系列波長〔如對太陽光譜求積分〕及一系列折射率求和〔如通過散射參量反推折射率〕。當折射率虛部mIm很大時,用向后循環(huán)法求An很不穩(wěn)定。而向前遞推總是穩(wěn)定的〔但向后遞推安全時,總是優(yōu)先選擇,因為其計算速度很快〕。得出允許向后遞推的經歷標準:用正確的向前地推與相對應的向后地推做比較,當發(fā)現(xiàn)對和g的相對誤差超過10-6時,認為計算失敗。對于一對確定的(x,mRe),我們采用向后遞推尋找第一個循環(huán)失敗的研究說明:對于確定的,,的值隨著x的增加很快趨向于一個確定值。對如果在任意角度下、的實部和虛部的相對誤差超過時,認為對和的向后遞推失敗?!捕藭r,并不受影響,因為當,的相對誤差到達時,的相對誤差總維持在以下?!硨蛯ι⑸鋸姸群推冗B分式算法總結: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序列
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
- 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年度個人房屋抵押民間借貸借款協(xié)議文本3篇
- 2025年度平板車租賃運輸合同風險評估與控制
- 個人租房合同標準版可打印
- 二零二五年度處理男子外遇妻子懷孕離婚撫養(yǎng)權糾紛調解合同
- 二零二五年度智慧家居項目報建代理與智能家居合同2篇
- 2025年度鐵礦石進出口關稅及稅費繳納合同規(guī)范
- 二零二五版礦產資源租賃居間代理合同3篇
- 二零二五年度企事業(yè)單位廉潔從業(yè)監(jiān)督員聘任合同
- 2025阿里巴巴智能物流機器人研發(fā)及銷售合同3篇
- 重慶文化藝術職業(yè)學院《應用文寫作》2023-2024學年第一學期期末試卷
- 《大學生職業(yè)發(fā)展與就業(yè)指導》課程標準
- 第23課《出師表》課件(共56張)
- GB/T 3953-2024電工圓銅線
- 發(fā)電機停電故障應急預案
- 接電的施工方案
- 常用藥物作用及副作用課件
- 幼兒阿拉伯數字描紅(0-100)打印版
- 社會組織等級評估報告模板
- GB/T 12173-2008礦用一般型電氣設備
- 新媒體研究方法教學ppt課件(完整版)
- 2020新版?zhèn)€人征信報告模板
評論
0/150
提交評論