蒙特卡洛方法模拟中子吸收2014.5.12讲解.ppt

蒙特卡罗方法模拟中子吸收 为了求得圆周率π值，在十九世纪后期，有很多人作了这样的试验：将长为2l的一根针任意投到地面上，用针与一组相间距离为2a（ l＜a）的平行线相交的频率代替概率P，再利用准确的关系式： 求出π值 其中Ｎ为投计次数，n为针与平行线相交次数。这就是古典概率论中著名的蒲丰氏问题。 一些人进行了实验，其结果列于下表 ： 为了得到具有一定精确度的近似解，所需试验的次数是很多的，通过人工方法作大量的试验相当困难，甚至是不可能的。因此，蒙特卡罗方法的基本思想虽然早已被人们提出，却很少被使用。上世纪四十年代以来，由于电子计算机的出现，使得人们可以通过电子计算机来模拟随机试验过程，把巨大数目的随机试验交由计算机完成，使得蒙特卡罗方法得以广泛地应用，在现代化的科学技术中发挥应有的作用。 蒙特卡罗方法的优点 1、能够比较逼真地描述具有随机性质的事物的特点及物理实验过程。 2、受几何条件限制小。 3、程序结构简单，易于实现。 蒙特卡罗方法的主要应用范围 合作起源于Manhattan工程：利用ENIAC(Electronic Numerical Integrator and Computer)计算产额 目前国内外利用蒙特卡罗方法解决屏蔽问题很普遍。在粒子输运问题，统计物理，典型数学问题，真空技术，激光技术以及医学，生物，探矿等方面应用很广。 常用的通用蒙特卡罗程序简介 MORSE程序 较早开发的通用蒙特卡罗程序，可以解决中子、光子、中子－光子的联合输运问题。 EGS程序 用蒙特卡罗方法模拟在任意几何中，能量从几个KeV到几个TeV的电子-光子簇射过程的通用程序包。 MCNP程序 MCNP是美国Los Alamos国家实验室开发的大型多功能通用蒙特卡罗程序，可以计算中子、光子和电子的联合输运问题以及临界问题，中子能量范围从10-11MeV至20MeV，光子和电子的能量范围从1KeV至1000MeV 。程序采用独特的曲面组合几何结构，使用点截面数据，程序通用性较强，与其它程序相比，MCNP程序中的减方差技巧是比较多而全的。 MCNP程序的使用方法 要用MCNP程序来处理问题，需要按照MCNP规定的格式编写一个输入文件，该文件包括描述问题所必须的全部输入信息。文件采用卡片结构，每行代表一张卡片，文件由一系列卡片组成，对于任一特定的问题，只需用到INP全部输入卡片的一小部分。在这个文件中要根据要解决的问题给出发射粒子（即源的形状、位置，粒子的种类、能谱）、穿过的介质（包括化学成分、密度、几何形状）、要探测的计数（探测位置、能量区间）的信息，然后提交给MCNP程序进行处理即可得到结果，然后根据自己的需要从输出文件中撮自己感兴趣的信息。 输入文件的组成 整个输入文件由3大块组成，块之间由2个有且仅有1行的空行分隔。第1块用来定义Cell，每个Cell都有标号、组成材料、材料密度、几何结构；第2块用来定义那些定义Cell所需要的Surface，每个Surface也都有一个唯一的标号，而面类型、面参数用来定义一个面的具体位置；第3个块用来定义整个实验装置的一些全局特性：比如源的性质、所有的材料、要测量的信息等等。 举例应用 图1给出了我们要进行模拟的实验装置图，从“点中子源”发出的具有一定能谱的“中子”，在穿过“水介质”之后其能谱和运动方向都会发生一定的改变，通过放置在“水介质”背面的“探测器”可以得到“探测器”位置上的“中子”的能谱。 如何执行MCNP程序 输入文件已经写了，可以提交给MCNP程序进行计算了。这里需要你会进行简单的DOS操作。 首先进入MCNP程序的根目录，为自己的输入文件建立一个文件夹，比如L，把输入文件保存在里面，起一个名字，比如inp。这里要注意文件夹和文件的名字不要太长，加起来不要超过6个，否则MCNP处理的时候会出问题。 再回到MCNP根目录，执行 mcnp iprxz n=L\inp，就会进入MCNP程序啦。 先说一下mcnp执行时的几个参数的用途 我们的程序加入了参数p，所以直接按回车它首先进行几何结构绘图模式，如下 计算结束后，输入tally可以进入探测结果显示模式： 谢谢！ 通过察看如果发现几何结构和我们预想的是一样的，那么输入end就可以离开绘图模式，开始进行计算了，计算时间视结构的复杂性和计算的粒子数多少而定 至此，我们就通过MCNP解决了一开始提出的要求：一个点中子源，发出一定能谱的中子，经过水的散射和慢化后在后面的探测器上探测到的中子能谱。更多信息可以通过输出结果文件：输入文件名＋o得到。 蒙特卡罗方法简介 Monte Carlo方法： 蒙特卡罗方法，又称随机抽样或统计试验方法，属于计算数学的一个分支，它是在上世