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

蒙特卡罗方法模拟中子吸收 蒙特卡罗方法简介 Monte Carlo方法： 蒙特卡罗方法，又称随机抽样或统计试验方法，属于计算数学的一个分支，它是在上世纪四十年代中期为了适应当时原子能事业的发展而发展起来的。 Monte Carlo名字的由来： Monte Carlo是摩纳哥（monaco)的首都，该城以赌博闻名 第二次世界大战期间，冯·诺依曼和乌拉姆两人把从事的研究与原子弹有关的秘密工作，以赌城名蒙特卡罗作为秘密代号来称呼。他们的具体工作是对裂变物质的中子随机扩散进行直接模拟。用赌城名比喻随机模拟，既风趣又贴切，很快得到了人们的普遍赞许和接受。 为了求得圆周率π值，在十九世纪后期，有很多人作了这样的试验：将长为2l的一根针任意投到地面上，用针与一组相间距离为2a（ l＜a）的平行线相交的频率代替概率P，再利用准确的关系式： 求出π值 其中Ｎ为投计次数，n为针与平行线相交次数。这就是古典概率论中著名的蒲丰氏问题。 Monte Carlo方法举例：Buffon投针实验 一些人进行了实验，其结果列于下表 ： 为了得到具有一定精确度的近似解，所需试验的次数是很多的，通过人工方法作大量的试验相当困难，甚至是不可能的。因此，蒙特卡罗方法的基本思想虽然早已被人们提出，却很少被使用。上世纪四十年代以来，由于电子计算机的出现，使得人们可以通过电子计算机来模拟随机试验过程，把巨大数目的随机试验交由计算机完成，使得蒙特卡罗方法得以广泛地应用，在现代化的科学技术中发挥应有的作用。 蒙特卡罗方法的优点 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给出了我们要进行模拟的实验装置图，从“点中子源”发出的具有一定能谱的“中子”，在穿过“水介质”之后其能谱和运动方向都会发生一定的改变，通过放置在“水介质”背面的“探测器”可以得到“探测器”位置上的“中子”的能谱。 标题卡 栅元卡 空行分隔符 曲面卡 空行分隔符 数据卡 仅一行，占用第 1～80 列。作为输出标题。 定义构成整个系统的各个基本介质单元以及相应的物理信息。 其它数据，包括问题类型、源描述、材料描述、计数描述，问题截断条件等。 第1行：A sample 这一行是一个注释性质的语句，里面的内容没有任何实际作用，仅仅是用来简单的向他人说明这个输入文件要解决的问题。 第2行－5行：定义了4个Cell。 Cell是MCNP输入文件中的基本单元，它代表实验装置中某一特定的物体，Cell的定义中包含了该