激光粒度仪 MIE散射系数计算方法的研究济南微纳作者.docxVIP

散射系数计算方法的研究 王少清 任中京 张希明 何芳 江海鹰 山东建筑材料工业学院，济南 250022 摘要：介绍一种计算 Mie散射系数的新方法，给出计算实例。 关键词：Mie散射Mie系数 Mie计算 Study of calculati ng method of mie scatteri ng coefficie nt Wang shao qing Ren zho ngji ng Zhang ximing Hefang Jia nghai ying (Sha ngdong In stitute of buildi ng materiele. Jinan 25022) Abstract: A new method of calculat ing mie coefficie nt is in troduced,a nd calculat ing examples are give n Keywords :Mie scatteri ng ,Mie coefficie nt ,Mie calculati on 1引言 Mie理论是球形颗粒对单色光的散射场分布的严格解析解。目前在环保、动力、气象、 天文、两相流及粉体颗粒尺寸分布测量等工程技术领域中有广泛的应用。 利用单一颗粒或颗 粒群光散射场的测量数据，可以反推得散射颗粒或颗粒群的很多物理性质 ，如颗粒的尺寸和 颗粒的折射率等。但反推前必须事先计算出各种尺寸的颗粒在各种复折射率下的散射场分布 数据。由于散射解十分复杂，故虽然其解折表达式早在年速、 准确计算，以满足工程测量，尤 其是在线式工程测量方面的需要，一直是多年来没有得到彻底解决的一个难题。 1968年Dave最先报道了完整的 Mie散射计算方法，又针对部分计算提出新算法。 国内 也有人发表了他们自己的算法。但总的看来 ，这些算法均有各自的局限性。尤其是当颗粒尺 寸或折射率的虚部值较大时，往往计算速度过慢或产生溢出和不收敛等现象。 本文介绍散射 新算法。该算法的特点是不受颗粒尺寸及折射率的限制 ，不会产生溢出和不收敛现象，且具 有较快的计算速度。 2 Mie散射系数的计算公式 Mie散射计算的中心问题是计算 Mie散射系数an和bn,其表达式为 VF. 9 J. m a) 一 m V* wm) _ w 7^\ (ma) — W* =— ■(机小 2) 其中，a为颗粒的尺寸参数，定义为 m为颗粒在周围介质中的相对复D m为颗粒在周围介质中的相对复 m ~ nr L H- imi ± V 0) ⑷ (5) hO =— (m/2X— lr^J-.+! (jj) =-伍) 式中，Jn+1和Nn+1分别为半整数阶的第一类和第二类贝塞尔函数。 以上是计算 Mie散射系数的基本公式。与 Mie散射有关的物理量如散射场强度、颗粒 的消光系数、散射系数和吸收系数、 散射光中两偏振分量的相差、 颗粒群的比浊度、不对称 因子等)均可利用 Mie散射系数求得。 3计算产生溢出的原因 计算M ie散射系数须先计算和一般采用递推法。递推又分为向前递推 (即从n=0开始) 与向后递推(即从n=N开始至n=0为顶先设定值)。实验表明，向前递推总是快于向后递推。 妙疔=sin(in4 t叫心 —exp —exp )n j X COS tn yfjt)-卜 tCexp (— ) -卜 exp ((nt yU) —cos {iff 尹 i- ijn2a) —-^-Cexp (— + exp(irt2r7)] X Costm t xp (— m^a) —《叫d)闾门(叫G .M - Q 分析(8)( 9)两式可知，当■— 时，若颗粒尺寸 d很大或复折射率的虚部值 m2 很大，将使得乘积 m2d很大，可能使(8)(9)两式中的项的值超过计算机的数据限 ，从而产生 溢出。这是产生溢出的最要原因。另外在递推过程中 ，不恰当的算法也可能造成溢出。 4 an和bn的新算法 为解决上述问题，作者提出了新算法。将 an和bn公式变形如下，令 (ma) = + iL,f 其中賈屮 J和L*、匚柑分别衷示D 和h佈町的实部与虚部*于是叭和丹?可写为 = a-r + aBj f/r 一- m \ (fi} 叽=如+叽 __ 懾甲‘于二龙心}/；“( HIM、一 屮』即0》 (?) f** (fri a} t wi a) — d g、j和h”、.］分别卷示aH和叫的实部与虚 部°利用5)式可进一步推得’ bl + bl bl + bl (15) al(12) ~ (15)式应采用比值形式，即四个公式中的后半段形式，这样可避免递减过程中 与bl较大时乘法运算可能产生的溢出。在以上四式中 al TOC \o 1-5 \h \z 如=F.IQEm -桝岸1(亦 06) 巧=p?(