《蒙特卡罗方法及应用》6.MCNP在粒子输运中的应用课件.ppt

上传人(卖家):三亚风情 文档编号:3366459 上传时间:2022-08-24 格式:PPT 页数:63 大小:495.50KB
下载 相关 举报
《蒙特卡罗方法及应用》6.MCNP在粒子输运中的应用课件.ppt_第1页
第1页 / 共63页
《蒙特卡罗方法及应用》6.MCNP在粒子输运中的应用课件.ppt_第2页
第2页 / 共63页
《蒙特卡罗方法及应用》6.MCNP在粒子输运中的应用课件.ppt_第3页
第3页 / 共63页
《蒙特卡罗方法及应用》6.MCNP在粒子输运中的应用课件.ppt_第4页
第4页 / 共63页
《蒙特卡罗方法及应用》6.MCNP在粒子输运中的应用课件.ppt_第5页
第5页 / 共63页
点击查看更多>>
资源描述

1、1.MCNP应用概述应用概述2.MCNP一般应用步骤一般应用步骤3.MCNP输入文件描述输入文件描述4.例子例子 MCNP是现成的蒙特卡罗应用软件之一。是现成的蒙特卡罗应用软件之一。一般情况下,使用者只是应用,无需再一般情况下,使用者只是应用,无需再去按照蒙特卡罗基本思想编制软件。去按照蒙特卡罗基本思想编制软件。在这种情况下,使用者要做的工作有:在这种情况下,使用者要做的工作有:准备准备MCNP软件(复制、安装软件)软件(复制、安装软件)建模、编写输入文件,建模、编写输入文件,运行、计算,运行、计算,分析结果。分析结果。按照如下步骤可快速入门按照如下步骤可快速入门 (课堂讲解、演示操作)(课堂

2、讲解、演示操作)(1)准备软件,)准备软件,MCNP4C(2)根据任务建模)根据任务建模(3)编写输入文件,)编写输入文件,t.inp(4)运行、计算)运行、计算,MCNP i=t.inp o=t.o(5)提取目标信息、作图与分析,)提取目标信息、作图与分析,t.o MCNP4C易学易用、不需要安装,易学易用、不需要安装,复制到硬盘中即可使用;复制到硬盘中即可使用;在文件夹在文件夹MCNP4C中运行。中运行。MCNP5版本需要安装。版本需要安装。最新版本最新版本MCNP6,还没见国内用户公开使用。,还没见国内用户公开使用。根据工作任务建根据工作任务建立模型,根据模型建立模型,根据模型建立输入文

3、件,模型与立输入文件,模型与输入文件一一对应。输入文件一一对应。根据工作任务建根据工作任务建立模型,根据模型建立模型,根据模型建立输入文件,模型与立输入文件,模型与输入文件一一对应。输入文件一一对应。根据工作任务建立模型,根据模型建立根据工作任务建立模型,根据模型建立输入文件,模型与输入文件一一对应。输入文件,模型与输入文件一一对应。MCNP的输入包括几个文件,但主要的一的输入包括几个文件,但主要的一个是由用户编写的个是由用户编写的INP文件,是文本文件。文件,是文本文件。该文件包括描述问题所必须的全部输入信该文件包括描述问题所必须的全部输入信息。文件采用卡片结构,每行代表一张卡片,息。文件采

4、用卡片结构,每行代表一张卡片,文件由一系列卡片组成。文件由一系列卡片组成。输入文件,就是按照规定的格式和规则,输入文件,就是按照规定的格式和规则,填写模型的几何参数、物理参数和输出参数。填写模型的几何参数、物理参数和输出参数。占用每行前占用每行前80列。列。文本文件,采用文本编辑器,不要采用文本文件,采用文本编辑器,不要采用word等字处理软件。文件名,如:等字处理软件。文件名,如:lr.inp(1)inp文件的组成主要有三块,即三个部分:文件的组成主要有三块,即三个部分:最上一块最上一块栅元卡栅元卡 中间一块中间一块曲面卡曲面卡 最下一块最下一块数据卡数据卡 且相邻的两块之间必须空一行。且相

5、邻的两块之间必须空一行。(1)输入文件构成)输入文件构成信息块选择项标题卡仅一行,占用第 180 列。作为输出标题。栅元卡栅元卡定义构成整个系统的各个基本介质定义构成整个系统的各个基本介质单元以及相应的物理信息。单元以及相应的物理信息。曲面卡曲面卡定义组成栅元的曲面信息。定义组成栅元的曲面信息。数据卡数据卡其它数据,包括问题类型、源描述、其它数据,包括问题类型、源描述、材料描述、计数描述,问题截断条材料描述、计数描述,问题截断条件等。件等。其它选择项表表1,MCNP输入文件中物理量的单位输入文件中物理量的单位长度长度厘米厘米能量能量MeV时间时间10-8 秒秒温度温度MeV(kT)原子密度原子

6、密度1024 个原子个原子/厘米厘米3质量密度质量密度克克/厘米厘米3截面截面10-24 厘米厘米2原子量原子量中子质量的中子质量的1.008664967倍倍阿伏加德罗常数阿伏加德罗常数6.0231023 栅元卡部分栅元卡部分 第一列,栅元号第一列,栅元号 第二列,材料号第二列,材料号 随后,描述该行栅元的密度、所有曲面,随后,描述该行栅元的密度、所有曲面,重要性,重要性,$后是注释后是注释栅元卡部分栅元卡部分 曲面卡部分曲面卡部分 数据卡部分数据卡部分 问题类型(问题类型(MODE)卡)卡 如果不给出MODE卡,则缺省形式是MODE N,即缺省值是中子输运问题。格式:MODEx1 xi xi

7、 N,中子输运。P,光子输运。E,电子输运。ERG 能量(能量(Mev)POS 抽样位置的参考点抽样位置的参考点 PAR 特殊类型源特殊类型源 1 MODE N或或N P或或N P E为中子为中子 2 MODE P或或P E为光子为光子 3 MODE E为电子为电子Fn 计数卡计数卡 F8:P或或F8:E或或F8:P,E 探测器中脉冲数的能量分布探测器中脉冲数的能量分布En 记数能量卡记数能量卡例:例:E11 0.1 1 20 把把F11流量计数分成流量计数分成4个能量箱:个能量箱:(1)从截断能量到)从截断能量到0.1MeV,(2)从)从0.1到到1.0MeV,(3)从)从1.0到到20.0

8、MeV,(4)整个能量范围,即从截断能到)整个能量范围,即从截断能到20.0MeV。例例2,栅元卡部分,栅元卡部分 反推几何模型反推几何模型 晶体的几何描述几种方法?晶体的几何描述几种方法?栅元用各相关曲面的布尔运算表示,布栅元用各相关曲面的布尔运算表示,布尔算符包括:尔算符包括:交交(用空格表示(用空格表示)、)、并并(用冒号(用冒号:表示)、表示)、非非(用(用#表示)。表示)。缺省的运算顺序是先缺省的运算顺序是先非非,其次是,其次是交交,最,最后是后是并并,使用括号可控制布尔运算的次序。,使用括号可控制布尔运算的次序。非运算有两种形式:(1)#n,n是某个栅元号,#n表示一个由不在栅元n

9、内的点组成的空间区域。(2)#(-),括号内是对某一个栅元进行描述的曲面栅元关系组,这一形式定义的几何区域由不属于括号内描述区域的点组成的空间。例如:3 0 -1 2 -4$定义栅元3#3$与下行相同#(-1 2 -4)例例2,曲面卡部分,曲面卡部分 曲面编号可以不连续曲面编号可以不连续 表表3 MCNP曲面卡曲面卡助记符 类型说 明方 程卡片上的数据项PPXPYPZ平面一般垂直X轴垂直Y轴垂直Z轴SOSSXSYSZ球面球心在原点一般球心在X轴球心在Y轴球心在Z轴C/XC/YC/ZCXCYCZ圆柱面平行于X轴平行于Y轴平行于Z轴轴心在X轴轴心在Y轴轴心在Z轴RzRyRxRzyxRRzzyxRz

10、yyxRzyxxRzzyyxxRzyx0)(0)(0)(0)()()(022222222222222222222DDDDCBADzDyDxDCzByAx0000RRRRyxRzxRzyRyxRzxRzyRyyxxRzzxxRzzyy0000)()(0)()(0)()(222222222222222222表表3.1 MCNP曲面卡(续一)曲面卡(续一)助记符 类型说 明方 程卡片上的数据项K/XK/YK/ZKXKYKZ圆锥面平行于X轴平行于Y轴平行于Z轴轴心在X轴轴心在Y轴轴心在Z轴1只用于单叶锥面SQ椭球面双曲面抛物面轴平行于X、Y或Z轴GQ圆柱面圆锥面椭球面双曲面抛物面轴不平行于X、Y或Z轴

11、1111110)(0)(0)(0)()()(0)()()(0)()()(222222222222222222tztytxtzyxtzyxtzyxzztyxyytzxxxtzyzztyyxxyytzzxxxxtzzyyzyxGFEDCBAGzzFyyExxDzzCyyBxxA0)(2)(2)(2)()()(222KJHGFEDCBAKJzHyGxFzxEyzDxyCzByAx0222表表3.1 MCNP曲面卡(续二)曲面卡(续二)助记符 类型方 程卡片上的数据项TXTYTZ椭圆或圆形的圆环面对称轴平行于X、Y或Z轴X Y Z P由点定义的曲面CBAzyxCBAzyxCBAzyxCAyyxxBzz

12、CAzzxxByyCAzzyyBxx01)()()(01)()()(01)()()(222222222222222222 曲面的指向在曲面的指向在MCNP中具有重要意义。中具有重要意义。如果点(如果点(x,y,z)在一个曲面的方程计算)在一个曲面的方程计算值为正,则称该点对于这个曲面是正向的;值为正,则称该点对于这个曲面是正向的;反之则为负向的。反之则为负向的。对于球、柱、锥及环,曲面外部是正向对于球、柱、锥及环,曲面外部是正向的。对于垂直于坐标轴的平面(的。对于垂直于坐标轴的平面(PX、PY或或PZ),大于相应平面截距的点是正向的。),大于相应平面截距的点是正向的。在栅元的几何说明中,关于曲

13、面的指向在栅元的几何说明中,关于曲面的指向是很重要。是很重要。假定曲面假定曲面 s 的曲面方程为的曲面方程为 f(x,y,z)0,则,则对于对于f(x,y,z)0的区域对于曲面的区域对于曲面 s 具有正的指具有正的指向;向;对于对于f(x,y,z)0的区域对于曲面的区域对于曲面 s具有负具有负的指向。的指向。正指向的区域用正指向的区域用+s表示,表示,“+”号可不写;号可不写;负指向的区域用负指向的区域用-s表示。表示。例例2,数据卡部分,数据卡部分 为书写方便,可以使用如下书写功能:为书写方便,可以使用如下书写功能:I.nR功能,表示将它前面的数据重复功能,表示将它前面的数据重复n次。次。例

14、如:例如:2 4R 等同于等同于 2 2 2 2 2II.xM功能,它表示的数值为前面的数据乘上x。例如:1 1 2M 2M 4M 2M 1 1 2 4 16 32III.nI功能,表示在与其前后相邻的两个数之间,功能,表示在与其前后相邻的两个数之间,插入插入n个线性插值点。对于个线性插值点。对于 X nI Y 的结构,的结构,如果如果X和和Y是整数,且是整数,且XY刚好是刚好是n+1的整倍的整倍数,则产生标准的整数插值,否则产生实数数,则产生标准的整数插值,否则产生实数插值,但插值,但Y值直接存储。值直接存储。例如:例如:1.5 2I 3.0 1.5 2.0 2.5 3 2.0可能不精确可能

15、不精确而而 1 4I 6 1 2 3 4 5 6 都是精确定整数都是精确定整数必须满足以下规则:i.nR前面必须放一个数或者放由R或M产生的数据项。ii.nI 前面必须放一个数或者放由R或M产生的数据项,而它的后面还必须有一个常数。iii.xM 前面必须放一个数或者放由R或M产生的数据项。例如:1 3M 2R1 3 3 31 3M I 41 3 3.5 41 3M 3M1 3 91 2R 2I 2.51 1 1 1.5 2.0 2.51 R 2M1 1 21 R R1 1 11 2I 4 3M 1 2 3 4 121 2I 4 2I 101 2 3 4 6 8 101 4I 3M错误!错误!I

16、V.NPS 历史数截断卡历史数截断卡格式格式:NPS N N 粒子的历史数。粒子的历史数。该卡指定要计算的粒子历史数,当计算该卡指定要计算的粒子历史数,当计算的粒子历史数达到指定值时,的粒子历史数达到指定值时,MCNP将中止将中止计算。计算。举例,举例,nps 100000 nps 10000000 两者运行结果的区别?两者运行结果的区别?V.CTIME 计算时间截断卡计算时间截断卡格式:格式:CTIME x x 计算时间(分钟)。计算时间(分钟)。该卡指定该问题在计算机上花费的时间该卡指定该问题在计算机上花费的时间限制,超过该时间,限制,超过该时间,MCNP将中止计算。将中止计算。所有中止计

17、算的条件只要达到一个就中所有中止计算的条件只要达到一个就中止计算。止计算。举例举例 在栅元卡上可定义栅元参数以代替在输入文件中数据卡部分定义的栅元参数。格式为:关键词值。这里允许的关键词有:带有粒子标识符的IMP、VOL、PWT、EXT、FCL、WWN、DXC、NONU、PD和TMP,以及关于重复结构的4个栅元参数卡:U卡、TRCL卡、LAT卡和FILL卡。例如:10 16 -4.2 1 -2 3 IMP:N=4 IMP:P=8解释之?例如例如:10 16 -4.2 1 -2 3 IMP:N=4 IMP:P=8 表示栅元表示栅元10由曲面由曲面1的正面、曲面的正面、曲面2的负面的负面和曲面和曲

18、面3的正面的交集组成,填充质量密度为的正面的交集组成,填充质量密度为 4.2 克克/厘米厘米3 的的16号材料。该栅元的中子重要性号材料。该栅元的中子重要性为为4,光子重要性为,光子重要性为8。在简写格式LIKE n BUT中,分别表示栅元的介质号和密度。例如:2 3 -3.7 -1 IMP:N=2 IMP:P=4 3 LIKE 2 BUT TRCL=1 IMP:N=10表示栅元3除了有不同的中子重要性和位置以外,其它方面与栅元2完全一样。即栅元3的定义及其材料和密度与栅元2一样,它们的光子重要性也一样。(3)减方差减方差MCNP运用以下卡片来减小方差:助记符卡片类型IMP栅元重要性ESPLT

19、能量分裂和俄国轮盘赌PWT次级光子权重EXT指数变换VECT方向矢量定义FCL强迫碰撞助记符卡片类型WWE权重窗的能量或时间间隔WWN权重窗的边界WWP权重窗的参数WWG权重窗生成器WWGE权重窗生成器的能量或时间间隔MESH分层重要性网格权窗生成器PD探测器贡献DXCDXTRAN贡献BBREM韧致辐射偏倚因子IMP或WWN卡之一是必须的,其它卡片是可选的。I.IMP 栅元重要性卡格式:IMP:nx1 x2 xi xI n 中子为 N,光子为 P,电子为 E。N,P、P,E 或 N,P,E 也是允许的,如果它们的值相同。xi 栅元 i 的重要性,i 1,2,II 问题中的栅元总数。缺省:在一个

20、MODE N P 问题中,若省略了IMP:P卡,则所有栅元的光子重要性都置 1,除非其中子重要性为零,这时其光子重要性也为零。栅元重要性卡用于输入各个栅元的重要性。零重要性可以用来终止粒子的历史。IMP 卡是必须的,除非用了WWN卡。例:IMP:N1 2 2M 0 1 20R 表示栅元1的中子重要性为1,栅元2为2,栅元3为4,栅元4为0,栅元5至25为1。此时如果一个粒子从栅元2进入栅元3,则粒子分裂为两个粒子,每个粒子权重为原来的一半。若粒子从栅元3进入栅元2,则以50的概率进行俄国轮盘赌终止粒子的历史,如果未能终止,则粒子的权重加倍。如果粒子进入“0”重要性栅元,则粒子被杀死。如果粒子进

21、入真空栅元,则即不分裂也不终止,然而,当粒子离开真空进入一个非真空栅元时,则要根据其重要性比值的增减进行分裂或俄国轮盘赌。源变量源变量说明(续二)变量变量说说 明明缺缺 省省POS抽样位置的参考点。0,0,0RAD 抽样位置离开POS或AXS的径向距离。0EXT栅元源:沿AXS方向离开POS的距离。曲面源:离开AXS的角度余弦。0AXS用于EXT和RAD的参考矢量。没有方向X抽样位置的X坐标。无XY抽样位置的Y坐标。无YZ抽样位置的Z坐标。无Z源变量源变量说明(续三)变量变量说说 明明缺缺 省省CCCCookie-cutter栅元。无ARA曲面的面积(当平面源对点探测器要求直接贡献时才需要此卡

22、)。无WGT源粒子权重。1EFF位置选舍抽样效率判据。0.01PAR源发射的粒子的类型。1中子2光子3电子1:包括中子的问题。2:无中子但有光子问题。3:仅有电子的问题。II.SIn 源信息卡格式:SIn选项选项 I1 Ik n 分布号(n 1999)。选项选项 Ii 的说明,允许的值是:省略或H:直方图分布的区间边界,仅用于标量。L:离散的源变量值。A:定义概率分布密度的点。S:分布号。I1 Ik 源变量值或分布号。缺省:SIn H I1 Ik III.SPn 源概率卡格式:SPn选项选项 P1 Pk 或:SPn f a b n 分布号(n 1999)。选项选项 Pi 的说明,允许的值是:省

23、略:对H或L分布与D相同,对A分布为概率密度。D:H或L分布的各个区间的概率值(不用归一)。C:H或L分布的各个区间的累积概率值。V:仅用于栅元分布。概率与栅元体积成 正比(若有Pi,要乘以Pi)。P1 Pk 源变量概率。f 内部函数的标识符。a,b 内部函数的输入参数。缺省:SPn D P1 Pk IV.SBn 源偏移卡格式:SBn选项选项 B1 Bk 或:SBn f a b n,选项选项,f,a,b 均和SPn卡一样,只有一点不同,即内部函数 f 只允许21和31。B1 Bk 源变量偏移概率。缺省:SBn D B1 Bk SPn卡的第一种格式的第一项是正数或非数值,表明该卡及它的SIn卡定

24、义了一个概率分布函数。SPn卡的第二种格式的第一项是负数,表示使用一个内部解析函数产生源变量的一个连续的概率密度函数。VII.一般源的例子例1:SDEF ERG=D1 POS=x y z WGT=wSI1H E1E2 Ek SP1D 0 P2 PkSB1D 0 B2 Bk这是个位于(x,y,z),权重为w的各向同性点源,能量采用偏移分布抽样。例2:SDEF SUR=m AXS=i j k EXT=D6 SB6-31 1.5$SI6 -1 1;SP6 -31 0这是曲面m上的源。发射方向用余弦分布抽样确定。通过分布6对与方向(i,j,k)的夹角余弦按指数偏移抽样确定源粒子在曲面上的位置。粒子权重

25、的最大和最小值是e1.54.48和e-1.50.223。例3:SDEF SUR=m NRM=-1 DIR=D1 WGT=wSB1-21 2$SI6 0 1;SP6 -21 1这是球面m上向内发射的源。如果wr2(r为球面m的半径),则该源和VOID卡、VOL卡以及计数类型2和4一起可以估计曲面的面积和栅元的体积。SB1卡的方向偏移使更多粒子朝着最感兴趣栅元所在的球m的中心发射。同时,该偏移顺便提供了球m体积的零方差估计。例4:SDEF SUR=m POS=x y z RAD=D1 CCC=n SI1r$SI1 0 r;SP1 -21 1这是沿正法线方向从曲面m上发射的单向源。源位置在以(x,y

26、,z)为中心,r 为半径的曲面上,用户必须确保点(x,y,z)在曲面m上。如果所抽取的点不在Cookie-cutter栅元n内,则舍弃该点并重新抽样。I.Fna 计数类型卡助记符类型说明Fn单位*Fn单位F1:(N、P、E)面流粒子MeVF2:(N、P、E)面通量粒子/cm2MeV/cm2F4:(N、P、E)体通量粒子/cm2MeV/cm2F5a:(N、P)点或环探测器通量 粒子/cm2MeV/cm2F6:(N、P、N,P)平均沉积能量MeV/克109 J/克F7:N平均裂变沉积能量 MeV/克109 J/克F8:(P、E、P,E)探测器探测的能谱 脉冲MeV+F8:E沉积电荷电荷无i.曲面和

27、栅元计数(类型1、2、4、6和7)简单格式:Fn:pl S1 Sk 一般格式:Fn:pl S1 (S2 S3)(S4 S5)S6 S7 n 计数号(n 1999)。pl N 或 P 或N,P 或 E。Si 用于计数的问题曲面号或栅元号或 T。这里只允许用栅元卡上列出的栅元以及界定这些栅元是曲面。上面的简单格式建立了该计数的 k 个曲面或栅元区间,对每个曲面或栅元单独列出结果。对于一般格式,区间按单个曲面或栅元以及若干个曲面或栅元的集合来划分。括号表示内部所有项的集合,括号内的项还可以用在其它集合或单独使用。T 是卡上列出的其它所有项的合集的简写。例1:F2:N 1 3 6 T 该卡指定了4个中

28、子通量计数,分别为穿过曲面1、3、6的平均通量以及穿过所有这三个曲面的平均通量。例2:F1:P (1 2)(3 4 5)6 该卡指定了3个光子流量计数,这三个计数分别是对曲面 1 和 2 的合集;曲面 3、4 和 5 的合集以及曲面6的计数。例3:F371:N (1 2 3)(1 4)T 该卡指定了3个中子流量计数,分别是对曲面 1、2和 3 的合集;曲面 1和 4 的合集以及曲面 1、2、3 和 4 的合集的计数。注意,这里的T 对重复出现的曲面1只使用了一次。iii.脉冲幅度计数(类型8)简单格式:Fn:pl S1 Sk 一般格式:Fn:pl S1 (S2 S3)(S4 S5)S6 S7

29、n 计数号(n 1999)。pl P或 E或P,E。Si 用于计数的问题栅元号或 T。脉冲幅度计数可以记录在一个探测器中由辐射产生的脉冲的能量分布。该卡后面列出的是栅元号或栅元集合,与F4卡一样。集合计数是累计计数而非平均计数。只允许有栅元区间以及能量区间。不论 pl 取何值,该计数类型将记录光子和电子,即粒子类型P、E或P,E对计数类型8是等价的。III.En 计数能量卡格式:EnE1 Ek n 计数号。Ei 计数 n 的第 i 个能量间隔的上界。缺省:如果没有此卡,则建立一个整个能量范围的区间,除非用E0卡更改了缺省值。可以使用E0卡对所有计数建立一个缺省的能量间隔结构。MCNP会自动给出

30、覆盖所有能量区间的总计数,但如果在该卡的末尾填写字符“NT”,则不记录该总计数。如果在该卡的末尾填写字符“C”,则记录的是计数累计值,最后一个是覆盖所有能量的总计数。例子:E11 .1 1 20该例子为F11计数划分了四个能量间隔。1)截止能量0.1 MeV2)0.11 MeV3)120 MeV4)整个能量范围特殊处理的描述以及参数:高斯能量展宽GEB a b c参数指定了物理辐射探测器能量展宽的半高宽。,E为能量。2cEEba半高宽II.ELPT 栅元能量截断卡格式:ELPT:n x1 x2 xi xIn 中子为 N,光子为 P,电子为 E。xi 栅元 i 的截断能量(MeV)。I 问题中的

31、栅元总数。该卡为每个栅元指定能量截断的下限。当粒子能量低于栅元指定能量下限或低于全局能量下限(由CUT卡定义)时,粒子的历史被终止。5)计算结果评价计算结果评价 MCNP在计算结果的同时,还会计算该结果的相对误差R,R估计值的方差/估计值。根据R的范围可以知道计算结果的可信程度:R的范围计数的质量0.51垃圾0.20.5不可信(只有少许的贡献记录)0.10.2可疑的0.10除点探测器外,一般可信0.05对点探测器一般可信6)接续运行接续运行这个问题可以分为三类:这个问题可以分为三类:第一类:第一类:第一次运行没有正常结束,中途中断,第一次运行没有正常结束,中途中断,后来在原有基础上接续运行。后

32、来在原有基础上接续运行。mcnp r=runtpe c 即可。即可。运行结束,另外生成了一个系统默运行结束,另外生成了一个系统默认的输出文件认的输出文件outp,最终的输出结果就,最终的输出结果就在这个在这个outp文件中。文件中。6)接续运行接续运行第二类:第二类:第一次运行已经正常结束,但还需要增第一次运行已经正常结束,但还需要增加一定的历史数。加一定的历史数。建立建立LB0.inp,LB0.bat两个文件,运行两个文件,运行LB0.bat文件。生成文件。生成LB0.inpo和和runtp*(*代代表一个字母表一个字母),等一段时间,等一段时间(如如nps 800000),运行正常结束。运

33、行正常结束。如果想再追加运行如果想再追加运行 100000个粒子,那么个粒子,那么设为设为nps 900000 6)接续运行接续运行第二类:第二类:另外写一个程序另外写一个程序inp文件,加入名字为文件,加入名字为11,内容为:内容为:continue nps*其中,其中,*为运行性的为运行性的nps数,即原来的数,即原来的nps+现在需要增加的现在需要增加的nps。接着再在接着再在dos下输入下输入mcnp i=11.inp r=runtpe c即可,运行完毕后,另外生成了即可,运行完毕后,另外生成了一个系统默认的输出文件,例如一个系统默认的输出文件,例如outp,最终,最终的输出结果就在这个的输出结果就在这个outp文件中。文件中。6)接续运行接续运行第三类:提前取数据的问题第三类:提前取数据的问题 历史数设定很大,计算时间长,计算途历史数设定很大,计算时间长,计算途中取数。中取数。如果使用的是如果使用的是mcnp4c的话可以通过:的话可以通过:ctrl+c先暂停运行先暂停运行 q 退出退出,这个时候就可以输出计数结果。,这个时候就可以输出计数结果。s会出现状态窗口,会出现状态窗口,我的问题是如何从状态窗口恢复到计算状态?我的问题是如何从状态窗口恢复到计算状态?直接回车即可。直接回车即可。

展开阅读全文
相关资源
猜你喜欢
相关搜索
资源标签

当前位置:首页 > 办公、行业 > 各类PPT课件(模板)
版权提示 | 免责声明

1,本文(《蒙特卡罗方法及应用》6.MCNP在粒子输运中的应用课件.ppt)为本站会员(三亚风情)主动上传,163文库仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。
2,用户下载本文档,所消耗的文币(积分)将全额增加到上传者的账号。
3, 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(发送邮件至3464097650@qq.com或直接QQ联系客服),我们立即给予删除!


侵权处理QQ:3464097650--上传资料QQ:3464097650

【声明】本站为“文档C2C交易模式”,即用户上传的文档直接卖给(下载)用户,本站只是网络空间服务平台,本站所有原创文档下载所得归上传人所有,如您发现上传作品侵犯了您的版权,请立刻联系我们并提供证据,我们将在3个工作日内予以改正。


163文库-Www.163Wenku.Com |网站地图|