蒙特卡罗方法常用蒙特卡罗程序介绍
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
的第 9~80 列到后续卡片的第 1~80 列都可填写运行
信息。在标题卡之前用一个空行分隔符结束信息块。
整理ppt
9 08.03.2021
(2) 初始运行的输入文件
信息块 空行分隔符
标题卡
栅元卡 … 空行分隔符 曲面卡 … 空行分隔符 数据卡 … 空行分隔符
其它
选择项 仅一行,占用第 1~80 列。作为输出标题。 定义构成整个系统的各个基本介质单元以及 相应的物理信息。
KCODE,MPLOT,ZA,ZB,和 ZC)
其它
选择项
接续运行必须在运行行信息或信息块中给出C项选择, 即Cm,表示从RUNTPE文件中读出第m次转储的内容 接着运算,如果m未指定,则读最后一次转储的数据。
如果不需要改变内容,则不需要接续输入文件,仅需
运行RUNTPE以及在运行行加上C选择。
整理ppt
整理ppt
6 08.03.2021
3. MCNP程序输入的描述
MCNP 的 输 入 包 括 几 个 文 件 , 但 主 要 的 一个是由用户编写的INP文件,该文件包括描 述问题所必须的全部输入信息。文件采用卡片 结构,每行代表一张卡片,文件由一系列卡片 组成,对于任一特定的问题,只需用到INP全 部输入卡片的一小部分。
整理ppt
5 08.03.2021
3) MCNP程序
MCNP是美国Los Alamos国家实验室开发 的大型多功能通用蒙特卡罗程序,可以计算中
子、光子和电子的联合输运问题以及临界问题, 中子能量范围从10-11MeV至20MeV,光子和电 子的能量范围从1KeV至1000MeV 。程序采用 独特的曲面组合几何结构,使用点截面数据, 程序通用性较强,与其它程序相比,MCNP程 序中的减方差技巧是比较多而全的。
定义组成栅元的曲面信息。
其它数据,包括问题类型、源描述、材料描 述、计数描述,问题截断条件等。 选择项
整理ppt
10 08.03.2021
(3) 接续运行的输入文件
信息块 空行分隔符
选择项
CONTINUE 写在第 1~8 列
数据卡
… 空行分隔符
只允许部分数据卡。(FQ,DD,NPS,CTME,
IDUM,RDUM,PRDMP,LOST,DBCN,PRINT,
整理ppt
4 08.03.2021
2) EGS程序
EGS是Electron-Gamma Shower 的缩写, 它是一个用蒙特卡罗方法模拟在任意几何中, 能量从几个KeV到几个TeV的电子-光子簇射过 程 的 通 用 程 序 包 。 由 美 国 Stanford Linear Accelerator Center提供。EGS于1979年第一次 公开发表,提供使用。EGS4是1986年发表的 EGS程序的最新版本。
第八章 蒙特卡罗方法应用程序介绍
1. 蒙特卡罗方法应用软件的特点 2. 常用的通用蒙特卡罗程序简介 3. MCNP程序输入的描述
整理ppt
1 08.03.2021
第八章 蒙特卡罗方法应用程序介绍
建立完善的通用蒙特卡罗程序可以 避免大量的重复性工作,并且可以在程 序的基础上,开展对于蒙特卡罗方法技 巧的研究以及对于计算结果的改进和修 正的研究,而这些研究成果反过来又可 以进一步完善蒙特卡罗程序。
整理ppt
2 08.03.2021
1. 蒙特卡罗方法应用软件的特点
通用蒙特卡罗程序通常具有以下特点: 1) 具有灵活的几何处理能力 2) 参数通用化,使用方便 3) 元素和介质材料数据齐全 4) 能量范围广,功能强,输出量灵活全面 5) 含有简单可靠又能普遍适用的抽样技巧 6) 具有较强的绘图功能
整理ppt
3 08.03.2021
2. 常用的通用蒙特卡罗程序简介
1) MORSE程序 较早开发的通用蒙特卡罗程序,可以解决
中子、光子、中子-光子的联合输运问题。采 用组合几何结构,使用群截面数据,程序中包 括了几种重要抽样技巧,如俄国轮盘赌和分裂 技巧,指数变换技巧,统计估计技巧和能量偏 移抽样等。程序提供用户程序,用户可根据需 要编写源分布以及记录程序。
整理ppt
7 08.03.2021
MCNBiblioteka Baidu输入文件中物理量的单位
长度
厘米
能量
MeV
时间
10-8 秒
温度
MeV(kT)
原子密度
1024 个原子 / 厘米3
质量密度
克 / 厘米3
截面
10-24 厘米2
原子量
中子质量的1.008664967倍
阿伏加德罗常数 6.023×1023
整理ppt
8 08.03.2021
输入文件中,在标题卡之后及最后的空行 结束卡之前的任何地方都可插入注释卡。注释 卡必须是字母“C”写在 1~5 列中的任意位置, 且至少用一个空格隔开后面的注释内容。
整理ppt
12 08.03.2021
a) 行输入格式
栅元卡、曲面卡和数据卡的书写格式是相 同的。必须从 1~5 列开始填写这些卡片相应 的名字(或编号)和粒子标识符,后面填写用空 格分隔的数据项。如果 1~5 列为空,则表示 它是前一张卡片的继续卡。如果在一行的末尾 有一个用空格隔开的符号“&”,则表示下一 行是该行的继续卡,数据可填写在 1~80 列。 一个数据项必须在一张卡片上写完,不得跨到 下一张卡片上。完全空白的一行则为两组卡片 的分隔符。 对任何给定的带有粒子标识符的 类型卡只能有一张。需要整数的数据项必须填 写整数,其它数据可填写为整数或浮点数以及 MCNP能读的数据。
1) 输入文件的基本形式 (1) 信息块
(2)
信息块的卡片放在INP文件中标题卡之前。
信息块给出了MCNP的一些运行信息,信息块上各部
分的意思和运行行信息是一样的,当运行行信息与信
息块中所指定的信息相矛盾时,则忽略信息块中相应
的信息,而以运行行信息为准。
(3)
信息块是可选的,信息块的第一张卡片,
必须在第 1~8 列写上“MESSAGE:”,从第一张卡片
11 08.03.2021
(4) 卡片格式 INP输入文件的每一行(称之为一张卡片)
都限于使用第 1~80 列并构成卡片映象。大部 分输入卡片按行填写;然而,对数据卡允许按 列填写。 $ 符号为它所在那行数据的结束符, 在 $ 符号后面的内容作为注释,它可从 $ 符号 后面的任一列开始。
标题卡只占一行,整行都可填入用户需要 的信息,也可以是空行。但要注意在其它地方 使用空行是作为结束符或者分隔符。