黑曲霉孢子凝聚粒子红外透过率计算讲稿.docVIP

黑曲霉孢子凝聚粒子红外透过率计算 目录 1 课题背景及研究意义 3 1.1.1 课题背景 3 1.1.2 课题研究意义 5 1.1.3本文主要工作 5 2、课题研究方法 5 2.1.1黑曲霉孢子凝聚粒子红外透过率计算的理论基础 5 （1）单次散射计算模型 5 （2）蒙特卡洛仿真模型 7 2.1.2 计算程序结构图 8 2.1.3计算程序介绍 8 （1）CCA（团簇-团簇凝聚）模拟 8 （2）DDA（离散偶极子近似法）模拟 10 （3）MC（蒙特卡洛）计算模型 11 （4）MFC界面程序介绍 14 3、黑曲霉孢子红外波段透过率的计算 17 3.1.1用CCA构建凝聚粒子空间结构的具体步骤 18 （1）为笔记本安装JAVA虚拟机 18 （2）运行CCA模型程序 20 （3）调整运行结果 21 （4）matlab计算孔隙率 22 3.1.2用DDA计算凝聚粒子单次散射参量的具体步骤 23 （1）算出凝聚粒子的等效半径 23 （2）ddscet.exe程序的参数调整 24 （3）ddscat.exe程序的运行 26 (4)用matlab计算消光系数 27 3.1.3用蒙特卡洛方法仿真烟幕透过率的具体步骤 27 （1）调整程序运行参数 27 （2）运行MCML程序 28 3.1.4计算结果及小结 29 （1）计算结果 29 （2）小结 30 1 课题背景及研究意义 1.1.1 课题背景 在自然界中，存在着许多由微小粒子组成的形状不规则的复杂结构凝聚粒子群，如大气中的气溶胶粒子、星际尘埃粒子和烟幕粒子等。在大气气溶胶中，细小的微生物颗粒由于静电、碰撞粘附而形成的微生物凝聚粒子群是其中一个重要组成部分。在环境监测、大气遥感、目标探测以及烟幕消光机理的认识等方面，研究具有随机取向的微生物凝聚粒子群所构成的复杂介质环境对激光传输的影响，具有非常重要的意义。【1】 目前，针对复杂结构凝聚粒子群消光性能的研究工作已经展开。前人利用扩散限制凝聚模型（DLA）【2】模拟了纳米石墨凝聚粒子的可能结构，基于离散偶极子近似方法（DDA）计算其红外消光截面。【3】；采用团簇-团簇凝聚模型（CCA）给出了随机分布烟尘团簇粒子的空间结构，并用蒙特卡罗方法研究了激光信号在随机分布烟尘团簇粒子中的传输衰减特性【4】。但是上述研究工作均针对纳米量级的无机气溶胶凝聚粒子，在微米量级的微生物粒子方面，虽有文献用有限时域差分法（FDTD）【5】或DDA计算【6】其单次散射特性，但是对于微生物凝聚粒子群对激光衰减特性方面的研究还有待展开。 我在前期学习了相关知识，大量阅读了相关文献材料，对其计算步骤及原理有了基本的了解。首先采用团簇-团簇凝聚模型模拟了微生物凝聚粒子的空间结构，利用离散偶极子近似方法得出微生物凝聚粒子的单次散射参量，通过蒙特卡洛方法仿真了10.6μm激光在微生物凝聚粒子群中的透过率值，最后，通过C++MFC平台设计一个计算界面，简化繁琐的计算操作步骤，实现时间上的提高。 1.1.2 课题研究意义 1.1.3本文主要工作 2、课题研究方法 2.1.1黑曲霉孢子凝聚粒子红外透过率计算的理论基础 （1）单次散射计算模型 用团簇-团簇凝聚（CCA）模型模拟了黑曲霉孢子凝聚粒子的空间结构（每个凝聚粒子包含40个半径为1.5μm的原始微粒），如图2.1所示，图2.2为放大5000倍的黑曲霉孢子扫描电镜图。 图2.1 图2.2 从图2.1可以看出，微生物凝聚粒子的空间结构是不规则的，局部与整体结构相似，这种自相关性说明其具有分型特征，结构有孔隙。其原始微粒数目M和回旋半径Rg的关系为 , (1) 式中,分别为第i个和第j个原始微粒在空间的坐标。凝聚粒子的孔隙率为P： , (2) 其中，为原始微粒半径。求出凝聚粒子的孔隙率之后就可以给出具有孔隙结构的凝聚粒子的等效折射率： , (3) 式中m为该种微生物粒子的复折射率，在计算中用等效折射率代替m，所得结果更能反映P对微生物凝聚粒子群透过率的影响。 微生物凝聚粒子的空间结构确定后，其单次散射参量用离散偶极子近似（DDA）计算【11】。Bruce T.Draine指出，DDA算法的使用条件是，m是目标的复折射率【12】。生物细胞粒子的折射率m一般在之间【13】，符合DDA使用条件。 （2） 2.1.3计算程序介绍 （1）CCA（团簇-团簇凝聚）模拟 模拟气溶胶凝聚粒子空间结构最常用的是采用CCA模型来模拟。CCA模型的基本思想是：让所有原始气溶胶单体粒子同时进入一个区域并做随机运动，一旦粒子与粒子相遇就