mcnp剂量率模拟笔记_兵器核科学_工程科技_专业资料.docx

.LOG采用美国洛斯阿拉莫斯(oLSAlomos)编制的蒙特卡罗软件MCNP,对核工业航测遥感中心二级放射性计量站的地面圆柱体放射性模型(Fi220cmx60cm,底面直径*高度，以下简称圆柱体模型) 表面及上方不同高度处的空气吸收剂量率进行了模拟计算。在空气吸收剂量率模拟计算中,首先按U-238，Ra-226，Th-232和K-40，分别模拟计算得到各衰变系列或核素的剂量率贡献,然后加和得到总吸收剂量率。质能吸收系数：伽马射线在物质中，穿过单位厚度以后，其能量被物质吸收的份额读到page3，2.MC模拟计算结果F4 通量计数，再定义计数乘子卡可转换成剂量Q1：如果想计算一定源强下（比如1000个/秒），某栅元内的剂量率，该怎么写？方案一：sdefpos = 0 0 0 wgt = 1000 f4：n 1 de4 df4 方案二：sdefpos = 0 0 0 f4：n 1 fm4 1000 de4 df4 是用方案一呢还是用方案二呢，还是其他方法？此外，方案一的结果与方案二的结果一样吗？//补充格式: DEn A E1 ... EkDFn A F1 ... Fkn: 记数号Ei: 第i个能量 (MeV) Fi: 第i个能量的函数值A: 取LOG或LIN插值方法(缺省为LOG) 用户可以使用DEn和DFn卡填写与能量相关的函数来修改常规记数，主要用于F6型记数的剂量计算。这两个卡必须有对应的数据项，对能量点之间的函数用对数或线性插值给出，如果粒子能量超出DEn卡指定的范围，则使用此卡最低或最高能量。可以使用DE0和DF0卡作为缺省卡。de为能量值，df为该能量点对应的注量-剂量转换系数。因为不可能给出任意一点的de,df，这样mcnp(x)会根据已知的de,df数据，进行插值，计算任意一点的df值。lin：线性插值；log：对数插值。一般用F4配合DE和DF卡DE 1MeV 2MeV DF a b （自己定义）若a=有效能量转移系数*1MeV*单位转换系数 a=有效能量转移系数*2MeV*单位转换系数吸收剂量=a*（F4在0-1MeV）至于吸收剂量的单位，看单位转换系数怎么计算了单位质量的能量沉积应该就是吸收剂量，MCNP给*F6出的是单个抽样粒子吸收剂量，考虑半衰期/衰变常数和分支比，就可以换算为吸收剂量率A1：方案一是指MCNP计算一个粒子来代表实际的1000个粒子。即认为这1000个粒子的踪迹是一模一样的啊！！！！方案二是MCNP计算一个粒子代表实际的一个粒子。最后归一，再乘以源活度。所以方案二是实际情况模拟，方案一是存在偏差的。但为什么还有方案一呢？那是为了省点时间，基本上误差范围内也允许。所以你的时间允许情况下，选方案二吧，呵呵。另外，计算出的结果是否一致或有差别，你可以自己实验一下。不要用真空环境或空气材料，使用水、人体模型、土壤或几种混合材料验证一下。nps设在1E+08。粒子能量最好大于2MeV。我赶脚MCNP很多问题要靠自己多编程计算、思考。就理解深了。Q2：最近用MCNP做剂量计算时，遇到下面一个问题：我利用DE、DF配合F4计数卡和Fm乘子卡算一点处的吸收剂量值，并且利用算出的剂量率反推源活度时发现比实际的活度低很多，比如我在Fm处设置的源活度为4.3Ci，但由此算出的值反推回去的源活度只有3Ci左右，是什么造成这样的源活度“缺失”呢？对于这点，我做了实际实验，实验佐证了我的算法没问题，但MCNP的模拟却造成这样的结果，是我模拟时哪一点没设置好吗，没考虑到吗？请大家帮着分析一下，尤其是做过剂量计算的坛友们，谢谢！物理模型很简单：一个点源，计算距点源30厘米处的剂量率，并且反推点源的活度，我的思路是：已知点源活度，将其设置为一个Fm卡乘子，这样得出的30厘米处剂量率，并且根据此剂量率反推出源的活度，推出的活度和我设置的活度（乘子）作比较，我的预先设想是这个反推出的活度值应该我设置的源的活度值（已知值）一样大小的，但实际模拟计算出的情况却差别很大，具体如我顶楼所说的，我搞不清原因，因为我在用MCNP算剂量这方面没什么经验，所以发贴问一下大家，这是怎么回事？这是我在算剂量的哪个细节没考虑到造成的？或者哪个概念我理解的不对呢？希望能及时得到大家的指点，谢谢。A：点源？什么样的点源？四面八方均匀发射的吗？计算一点的剂量，用错了卡，不是F4，要用F5. 还有，你描述的模型似乎不是很清楚，呵呵。你de，df卡里面的参数是从哪里来的？对于源点和剂量计算点的相对位置，源的体积大小与剂