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

Monte Carlo方法及应用软件在实验核物理中的应用 内容 Monte Carlo简介 核技术特点及Monte Carlo方法的应用 Monte Carlo应用软件 Geant4具体应用实例 将一根长度为l的针，随机的投放在两条间距为d(d＞l)的两条平行线中间，计算针与两条平行线相交的概率。 通过解析法，容易得到 针与平行线相交的概率 由此可以得到圆周率 Buffon投针问题 Buffon投针问题 Buffon投针问题 Buffon投针问题 Buffon投针问题 试验者 时间(年) 针长 投针次数 相交次数 π的估计值 Wolf 1850 0.80 5000 2532 3.15956 Smith 1855 0.60 3204 1218 3.15665 Fox 1884 0.75 1030 489 3.15951 Lazzarini 1925 0.83 3408 1808 3Buffon投针问题 Monte Carlo 蒙特卡罗方法又称统计模拟方法，是利用随机数进行数值模拟的方法 可以将蒙特卡罗法看作利用随机试验的方法计算积分，所计算的积分可看做服从某种分布的密度函数为f(r)的随机变量g(r)的数学期望 而通过某种试验得到N个观察值r1，r2...rN等，将N个随机变量的值g(r1), g(r2),…, g(rN)的算术平均值作为积分估计值 Monte Carlo Monte Carlo 根据概率论知识，可知 优点 对于具有随机性质的事件或物理过程其物理意义逼真 受几何限制小 收敛速度与问题的维数无关 误差容易确定 程序结构简单，易于实现 缺点 收敛速度慢 误差具有概率性 在粒子输运问题中，计算结果与系统大小有关 Monte Carlo 蒙特卡罗算法主要组成部分 概率密度函数(pdf)— 必须给出描述一个物理系统的一组概率密度函数; 随机数产生器—能够产生在区间[0,1]上均匀分布的随机数 抽样规则—如何从在区间[0,1]上均匀分布的随机数出发,随机抽取服从给定的pdf的随机变量; 模拟结果记录—记录一些感兴趣的量的模拟结果 误差估计—必须确定统计误差（或方差）随模拟次数以及其它一些量的变化； 减少方差的技术—利用该技术可减少模拟过程中计算的次数； 并行和矢量化—可以在先进的并行计算机上运行的有效算法 Monte Carlo 蒙特卡罗算法模拟的主要步骤 Monte Carlo Monte Carlo 数值 解法 Monte Carlo 实际 模型/ 原理 理论结果 模拟结果 测试结果 测试 理论分析 模拟 对比 对比 理解 模型错误 原理正确 原理错误 一致 不一致 一致 不一致 Monte Carlo Monte Carlo模拟 + 统计特性 高样本量 非常有利于Monte Carlo方法的应用 在核物理中的应用 而其中粒子输运问题是蒙特卡罗模拟的基本核物理过程 Monte Carlo模拟 蒙卡模拟粒子输运问题的主要步骤 Monte Carlo模拟 确定粒子输运过程中涉及的物理过程 蒙卡模拟粒子输运问题的主要步骤 多粒子情况下常用的有字典编辑分支法 Monte Carlo模拟 确定所用的蒙卡技巧（单粒子） 蒙卡模拟粒子输运问题的主要步骤 Monte Carlo模拟 确定粒子状态参数及状态序列 蒙卡模拟粒子输运问题的主要步骤 以容易实现，所用计算时间少为标准 Monte Carlo模拟 确定粒子输运过程中状态分布的抽样方法 在通过对粒子输运问题的进行模拟的基础上，考虑实际情况进一步对其结果进行分析，可以针对不同应用得出不同的物理结果，如屏蔽问题中的光通量，核辐射探测器中的探测效率、能量沉积谱、能量分辨率等。 Monte Carlo模拟 Monte Carlo应用软件 全名Multigroup?Oak Ridge Stochastic Experiment，是美国橡树岭国家实验室从60年代开始研制的大型、多功能、多群中子-光子耦合输运程序 特点 具有三维几何能力 具有多种功能，可解决中子、光子、中子-光子耦合输运问题 使用多群截面数据 模块结构 包括几种有效地减小方差技巧和计算技巧 程序具有很大的灵活性 Monte Carlo应用软件——MORSE 全名Monte Carlo Neutron andPhoton Transport Code，由美国洛斯阿拉莫斯国家实验室编制的具有目前最高水平的中子-光子输运程序，被称为“超级蒙特卡罗程序”。 特点 程序中的几何可以是三维任意态 使用精细的点截面数据 可用于中子、光子、中子-光子耦合输运问题，及临