Monte Carlo方法及相关软件在实验核物理中的应用教材课程.ppt

Monte Carlo方法及应用软件在实验核物理中的应用;内容 ;将一根长度为l的针，随机的投放在两条间距为d(d＞l)的两条平行线中间，计算针与两条平行线相交的概率。 通过解析法，容易得到 针与平行线相交的概率 由此可以得到圆周率 ;;;;Monte Carlo;蒙特卡罗方法又称统计模拟方法，是利用随机数进行数值模拟的方法 可以将蒙特卡罗法看作利用随机试验的方法计算积分，所计算的积分可看做服从某种分布的密度函数为f(r)的随机变量g(r)的数学期望 而通过某种试验得到N个观察值r1，r2...rN等，将N个随机变量的值g(r1), g(r2),…, g(rN)的算术平均值作为积分估计值 ;Monte Carlo;优点 对于具有随机性质的事件或物理过程其物理意义逼真 受几何限制小 收敛速度与问题的维数无关 误差容易确定 程序结构简单，易于实现 缺点 收敛速度慢 误差具有概率性 在粒子输运问题中，计算结果与系统大小有关 ;蒙特卡罗算法主要组成部分 概率密度函数(pdf)— 必须给出描述一个物理系统的一组概率密度函数; 随机数产生器—能够产生在区间[0,1]上均匀分布的随机数 抽样规则—如何从在区间[0,1]上均匀分布的随机数出发,随机抽取服从给定的pdf的随机变量; 模拟结果记录—记录一些感兴趣的量的模拟结果 误差估计—必须确定统计误差（或方差）随模拟次数以及其它一些量的变化； 减少方差的技术—利用该技术可减少模拟过程中计算的次数； 并行和矢量化—可以在先进的并行计算机上运行的有效算法 ;蒙特卡罗算法模拟的主要步骤 ;Monte Carlo;实际;Monte Carlo模拟;在核物理中的应用 而其中粒子输运问题是蒙特卡罗模拟的基本核物理过程 ;蒙卡模拟粒子输运问题的主要步骤 ;蒙卡模拟粒子输运问题的主要步骤 多粒子情况下常用的有字典编辑分支法 ;蒙卡模拟粒子输运问题的主要步骤 ;蒙卡模拟粒子输运问题的主要步骤 以容易实现，所用计算时间少为标准;在通过对粒子输运问题的进行模拟的基础上，考虑实际情况进一步对其结果进行分析，可以针对不同应用得出不同的物理结果，如屏蔽问题中的光通量，核辐射探测器中的探测效率、??量沉积谱、能量分辨率等。;Monte Carlo应用软件;全名Multigroup?Oak Ridge Stochastic Experiment，是美国橡树岭国家实验室从60年代开始研制的大型、多功能、多群中子-光子耦合输运程序 特点 具有三维几何能力 具有多种功能，可解决中子、光子、中子-光子耦合输运问题 使用多群截面数据 模块结构 包括几种有效地减小方差技巧和计算技巧 程序具有很大的灵活性 ;全名Monte Carlo Neutron andPhoton Transport Code，由美国洛斯阿拉莫斯国家实验室编制的具有目前最高水平的中子-光子输运程序，被称为“超级蒙特卡罗程序”。 特点 程序中的几何可以是三维任意态 使用精细的点截面数据 可用于中子、光子、中子-光子耦合输运问题，及临界系统特征值问题 减小方差技巧比较全面 程序的通用性很强 使用Fortran语言;全名为electron-gamma shower，由美国Stanford Linear Accelerator Center提供，用以模拟在任意几何中能量从几keV到几TeV的电子-光子簇射过程，目前最新版本为EGS-5 特点 元素和材料介质齐全 带电粒子和光子的输运均采用随机游动模式进行 光子的能量可以从1keV到几千GeV 反应类型非常齐全 使用Mortran语言 ;全名FLUktuierende Kaskade，是INFN（意大利核物理研究所）和CERN（欧洲粒子物理研究中心）联合开发的一种蒙特卡罗模拟计算程序，质子、电子加速器屏蔽设计，量热计，活化，计量学，探测器设计，宇宙射线，中微子物理及放射治疗等 特点 采用微观模型，物理过程严谨 约60种不同的粒子可以使用FLUKA进行模拟，能量范围宽 使用改进的CG（Combinatorial Geometry）软件包，可处理复杂几何结构 针对于大多数的应用，FLUKA使用者不需要编写程序 完全免费;全名GEometry ANd Tracking，由上世纪90年代CERN（欧洲核子研究组织）和KEK（日本高能加速器研究中心）组织了来自欧洲、俄罗斯、加拿大、日本和美国的几十个实验室、高校和研究机构的超过100名科学家和工程师进行合作编译出来的具有开放结构的蒙卡程序 特点： 采用面向对象的类的程序结构设计 采用C++语言编译程序，用户可自定义的开放代码 重新编写物理过程，不借助于外部程序 采用商业化的通用技术 完全免费 ;对于简单探测器相应问题的编程思路;谢谢