《蒙特卡罗方法及应用实验讲义2016》由会员分享,可在线阅读,更多相关《蒙特卡罗方法及应用实验讲义2016(18页珍藏版)》请在金锄头文库上搜索。
1、规苹*烟$评滋国迪济累陶H蕭H洲2016.8实验一 蒙特卡罗方法基本思想一、实验目的1、了解蒙特卡罗方法方法的基本思想;2、掌握蒙特卡罗方法计算面积、体积的方法;3、掌握由已知分布的随机抽样方法。二、实验原理Monte Carlo 方法,又称统计模拟方法或计算机随机模拟方法,是一种基于 “随机数”进行数值模拟的方法,一种采用统计抽样理论近似求解物理或数学问 题的方法。如待求量可以表述成某些特征量的期望值、某些事件出现的概率或两者的函 数形式,那么可采用蒙特卡罗方法求解。在求解某些特征量的期望值或某些事件 出现的概率时,必须构建合符实际的数学模型。例如采用蒙特卡罗方法计算某函 数所围面积时,构建
2、的数学模型是构造一已知面积的可均匀抽样区域,在该区域 投点,由伯努利定理大数定理可知,进入待求区域投点的频率依概率 1 收敛于该 事件出现的概率(面积之比)。由已知分布的随机抽样方法指的是由已知分布的总体中抽取简单子样。具体 方法很多,详见教材第三章。三、实验内容1、安装所需计算工具(MATLAB、fortran、C+等);2、学习使用 rand(m,n)、unifrnd(a,b,m,n)函数3、求解下列问题:3.0、蒲丰氏投针求圆周率。3.1、 给定曲线y =2 - x2和曲线y3 = x2,曲线的交点为:P1( - 1, 1 )、P2( 1, 1 )。曲线围成平面有限区域,用蒙特卡罗方法计
3、算区域面积;z Jx2 + y23.2、 计算所围体积z 1 + 41 x 2 y 2其中 0 = (x, y, z) | 1 x 1,1 y 1,0 z 106)与抽样时间。实验二 由已知分布的随机抽样方法一、实验目的1、掌握由已知分布的随机抽样方法。2、用编程语言实现某具体随机抽样方法。二、实验原理由已知分布的随机抽样方法指的是由已知分布的总体中抽取简单子样。具体 方法很多,本实验综合直接抽样方法、挑选抽样方法和替换抽样方法,以散射方 位角余弦分布的抽样为例。实验原理详见教材对应章节。1. 连续型分布的直接抽样方法对于连续型分布,如果分布函数F(x)的反函数FT(x)存在,则直接抽样方X
4、二 F-1( )F2. 挑选抽样方法为了实现从己知分布密度函数f(x)抽样,选取与f(x)取值范围相同的分布密 度函数h(x),如果f ( x)M = sup gh( x)-g x g则挑选抽样方法为:A /(兀)?M-h(Xh)3. 替换法抽样方法为了实现某个复杂的随机变量 y 的抽样,将其表示成若干个简单的随机变 量xi,x2,,xn的函数y= g(x1,x2, ,xn)得到 x1, x2, , xn 的抽样后,即可确定 y 的抽样,这种方法叫作替换法12n抽样。三、实验内容蒲丰氏问题的算法如何产生任意的(x,e)?x在0, a上任意取值,表示x在0, a上是均匀分布的,其分布密度函数为:
5、1/ a, 0 x afi( x) = o,其他类似地,0的分布密度函数为:f仆/兀,00兀f2( ) 0,其他因此,产生任意的仗,0)的过程就变成了由f/x)抽样x及由厶他抽样0的过程了。由此得到:x ag10 兀g2其中0刍均为(0,1)上均匀分布的随机变量。每次投针试验,实际上变成在计算机上从两个均匀分布的随机变量中抽样得到(x,0),然后定义描述针与平行线相交状况的随机变量s(x,0),为s( X, 0) = 0,x l - sin 0,其他如果投针N次,则s = s(x,0 )N N i ii1是针与平行线相交概率P的估计值。事实上,P = JJ s(x,0) f (x) f (0
6、)dxd012d0J i sin 0 dx2L0 n 0 ana于是有2l2ln =沁aP asN四、实验报告1、给出源程序程序并解释语句的含义;2、作出抽样框图、试验累积频率与理论累积频率关系图,并给出抽样次数 (106)与抽样时间。实验三 MCNP 方法在实验核物理中的应用一、实验目的1、了解 MCNP 程序运行流程;2、掌握MCNP输入文件编写规范;3、理解模拟内容、并能编写输入文件、运行,并获得计算结果;二、实验原理MCNP是一种常见的粒子输运模拟软件,软件的安装、运行和输入文件编写 方法详见相关参考资料。MCNP输入文件编写完成后,先确认输入模型是否正确,在DOS环境下进 行,打开运
7、行 DOS 环境,进行以下操作:DOS命令操作命令含义Mcnp i=name.inp o=name.o打开画图框PX vx输出模型在 x=vx 面上的切面PY vy输出模型在 y=vy 面上的切面PZ vz输出模型在 z=vz 面上的切面FACTOR m将输出图放大1/m倍Extent a b切面沿两坐标轴方向分别放大ORIGIN X Y Z定义画图中心位置(X,Y,Z)三、实验内容1、学习 MCNP 程序常见各种运行方法;2、编写以下问题的输入文件;2.1 对课堂讲解的实例,模拟溴化镧探测器对点源的能谱,实验做一遍。2.2 有一 HPGe 探测器,结构如图 1 所示。分别给出位于探测器轴心、
8、距离 探测器晶体中心25cm处的i37Cs源、60Co、1311源对应特征y射线的探测效率(计 算时相应特征射线的源粒子至少为107个),并给出三者混合源(活度比为1:1:3) 的能谱图(源发射总粒子数大于3x108个)。外死层锂厚度0.09cmAl 壳005cm冷指Cu的尺寸056cmx175cm图 1 HPGe 探测器结构图四、实验报告1、给出 2.1和 2.2 的 MCNP 输入文件并解释每一行的含义2、分别运行实例,给出实验结果,并对结果进行分析。实验四MCNP模拟计算y射线造成的剂量一、实验目的1、掌握应用软件MCNP、应用范围以及在辐射剂量计算和防护中的作用;2、进一步掌握MCNP
9、程序基本用法;3、利用MCNP解决一个简单的求解Y射线在空气、组织等效材料(肌肉)中 造成的剂量沉积的计算问题,并进行结果分析,得出结论;4、利用MCNP程序解决实际工作中碰到的实际问题;二、实验内容1、学习MCNP程序的基本组成、操作方法以及问题描述文件的写法;2、利用MCNP程序计算简单的Y射线源在空气、肌肉模型中的剂量沉积分 布,并对计算结果进行分析并绘图,得出结论,调整数据重新计算,并与理论计 算结果进行比较;三、内容简介1、MCNP程序的计算流程如下图1所示:图1 MCNP程序的计算流程2、MCNP输入文件inp的基本结构通过这个文件描述并建立一个蒙特卡罗计算问题,对问题的几何结构、
10、材料 记数要求等给以描述,如果需要,便可直接运行。该文件的格式如下:信息块I 空行分隔符(optional可选)标题卡栅元卡1栅元卡2栅元卡 n空行分隔符曲面卡1曲面卡 2。曲面卡 n 空行分隔符 数据卡1 数据卡2 数据卡 n空行分隔符(optional)其它选择项 (optional)其中栅元卡用来描述由不同的封闭曲面分割的立体空间区域,并用独有的数 字 ID 号加以标示,同时在各个栅元卡中说明包围该区域的曲面类型 (曲面卡)、 填充该区域的材料类型(材料卡)以及对应的材料密度等;曲面卡是用来描述不同类型曲面的,并用独有的数字 ID 号加以标示,最终 曲面卡被应用在栅元卡中,并利用交(与)
11、、联(或)、补(非)这些逻辑运算符号联合 不同曲面组成所需要的复杂的栅元。在mcnp中支持的常见曲面类型见参考文献3,4。数据卡类型很多,主要有粒子类型标识卡mode、重要性卡imp、通用源卡 sdef、粒子计数器卡Fn、材料描述卡Mn以及粒子截断卡(nps或ctme)等,数据 卡的类型涉及到了方方面面,类型很多,具体请见参考文献3,4。F面利用一个简单的例子来配合说明mcnp中的输入卡(inp)的编写格式。3、一个简单的说明例子z为说明如何填写INP文件,这里例举一个简 单问题。如图3所示,在一个边长10cm的石墨立 方体3中有两个半径0.5cm的球形空间,球1中充 满氧气,球2是铁球。在球
12、1中置一 14Me V各向 同性中子点源,计算球2外表面与能量相关的中 子通量。建立的INP文件如下:SAMPLE PROBLEM INPUT DECK11-0.0014 -72 2 -7.86-8图3 例子的几何示意图3 3 -1.601 -2-3 4-5 6 7 84 0-1:2:3:-4:5:-6空行PZ -52PZ53PY54PY-55PX56PX-57S0 -4 2.5 .58S0 44 .5J空行MODE pIMP:p1 1 1 0SDEFPOS=0 -4 -2.5 ERG=14F2:n8E01E-5 1E-4 1E-3 .01 .1 1 2 3 4 5 6 7 8 9 10 11
13、 12 13M18016 1M226000 1M36000 1NPS100000J空行本例中没有信息块,第一行是标题卡,之后至空格前为栅元块。栅元卡上依 次填写栅元号、材料号、密度和构成栅元界面的曲面号(带正负号),这里定义了 4个栅元:栅元1由球面7围成,里面填充材料1(1602气体),密度是0.0014g/cm3; 栅元2由球面8围成,填充材料2(铁),密度7.86g/cm3;栅元3由平面1、2、3、4、 5、6围成,不包括球面7、8以内的空间,填充材料3(石墨),密度1.6g/cm3;栅元 4是栅元3以外的空间,为真空。曲面卡上需要填写曲面号、曲面类型和曲面参数,本例中定义了 8个曲面,