《张根发-蒙特卡罗方法课件-mcnp程序专题讲座.ppt》由会员分享，可在线阅读，更多相关《张根发-蒙特卡罗方法课件-mcnp程序专题讲座.ppt（166页珍藏版）》请在金锄头文库上搜索。
1、入门 Introduction 输入文件基础 InputFileBasics 几何描述 GeometryDefinition 源的描述 SourceDefinition 计数描述 TallyDefinition 材料 物理成分和数据 Materials Physics Data 高级计数描述 AdvancedTallies 高级几何描述 AdvancedGeometries 各种简化 VarianceReduction 临界问题 Criticality MCNP程序专题讲座 MCNP内容提要 入门 Introduction 输入文件基础 InputFileBasics 几何描述 Geometr
2、yDefinition 源的描述 SourceDefinition 计数描述 TallyDefinition 材料 物理成分和数据 Materials Physics Data 高级计数描述 AdvancedTallies 高级几何描述 AdvancedGeometries 各种简化 VarianceReduction 临界问题 Criticality 返回首页 MCNP入门 熟悉Unix或Dos操作系统用Unix或Dos命令行运行MCNP绘制MCNP几何图象MCNP输入文件结构 UNIX UNIX操作系统是一个强大的多用户 多任务操作系统 支持多种处理器架构 按照操作系统的分类 属于分时操作
3、系统 最早由KenThompson DennisRitchie和DouglasMcIlroy于1969年在AT T的贝尔实验室开发 目前它的商标权由国际开放标准组织所拥有 只有符合单一UNIX规范的UNIX系统才能使用UNIX这个名称 否则只能称为类UNIX UNIX like UNIX的系统结构可分为两部分 操作系统内核 由文件子系统和进程控制子系统构成 最贴近硬件 系统的外壳 贴近用户 外壳由Shell解释程序 支持程序设计的各种语言 编译程序和解释程序 实用程序和系统调用接口等组成 UNIX系统采用树状目录结构 具有良好的安全性 保密性和可维护性 Linux Linux是一套免费使用和自
4、由传播的类Unix操作系统 是一个基于POSIX和UNIX的多用户 多任务 支持多线程和多CPU的操作系统 它能运行主要的UNIX工具软件 应用程序和网络协议 它支持32位和64位硬件 Linux继承了Unix以网络为核心的设计思想 是一个性能稳定的多用户网络操作系统 Linux的基本思想有两点 第一 一切都是文件 第二 每个软件都有确定的用途 其中第一条详细来讲就是系统中的所有都归结为一个文件 包括命令 硬件和软件设备 操作系统 进程等等对于操作系统内核而言 都被视为拥有各自特性或类型的文件 至于说Linux是基于Unix的 很大程度上也是因为这两者的基本思想十分相近 DOS DOS是Dis
5、kOperationSystem 磁盘操作系统 的简称 是个人计算机 PC 上的一类操作系统 DOS是命令模式下的人机交互界面 人通过这个界面使用一些接近于自然语言或其缩写的命令来运行和控制计算机 就好像两个人相互沟通 早期的DOS系统是由微软公司为IBM的个人计算机开发的 称为MS DOS 从1980年到1995年的15年间 DOS在IBM PC兼容机市场中占有举足轻重的地位 自1998之后Windows完全脱离了DOS 虽然DOS过时了 就像当年的软驱一样 命令行的批处理BAT程序却存活下来 而且活得很好 因为它高效 简单 实用 WindowsXP WindowsVista和Windows
6、7在 附件 中有一个 命令提示符 CMD 其模拟了一个DOS环境 运行MCNP mcnpi o options Options 选项 i处理输入文件默认值p绘图x处理截面默认值r粒子传输默认值z标绘计数结果标绘截面 运行MCNP mcnpi o 选项 默认文件名inp输入文件outpASCII输出文件runtpe二进制重启文件通过命令行改变默认值 mcnpinp exloupt exlorun exlr mcnpname exl 练习1a 运行一次 mcnpi demo1 什么文件被创建 再运行一次 mcnpi demo1 这时候什么文件被创建 不要使用默认文件名 始终清楚地定义文件名 或者i
7、 inNameo outNamer runNamen baseNamei innamen baseOutName这将会防止你覆盖先前地计算结果这将会帮助你知道哪一个结果是正确的 练习1b 运行一次 mcnpn demo1 什么文件被创建 再运行一次 mcnpi demo1n demo1 这时候什么文件被创建 绘制几何图象 计算机上的二维几何图象显示能够用来检查几何问题的很多方面 栅元和表面序号材料密度材料位置几何错误用红色虚线显示经常绘图检查几何结构 练习2a 绘制图象 mcnpi demo1n demo1 ip 概念图象放大全景显示改变方向验证材料 栅元 表面 密度等 MCNP内容提要 入门
8、 Introduction 输入文件基础 InputFileBasics 几何描述 GeometryDefinition 源的描述 SourceDefinition 计数描述 TallyDefinition 材料 物理成分和数据 Materials Physics Data 高级计数描述 AdvancedTallies 高级几何描述 AdvancedGeometries 各种简化 VarianceReduction 临界问题 Criticality 返回首页 输入文件基础 了解MCNP中的物理单位了解MCNP输入文件三个主要部分了解MCNP输入文件的格式规定了解MCNP输入文件的简写特点 MC
9、NP中的物理单位 MCNP输入文件 标题卡栅元卡要求空行分隔曲面卡要求空行分隔数据卡推荐空行作为结束 MCNP输入文件 Demoproblem 1forNEEP602MCCourse Spring2002 CCellcardsCcell Mat DensSurfacecombinations11 1 122 5 23012 3403CSurfacecards1sph30022sph 20013sph00015CDatacardsIMP n1110m110012 080161 0m230071 0sdefnps10000 MCNP输入文件格式 每行最多80个字符不含控制字符 比如 Tab注释行
10、标题卡之后的任何位置都可插入 第一列是字母 C 且随后四个空格 从输入数据之后的 符号后开始以上三种情况可以单独或同时存在 输入简写 nR 表示将它前面的一个数据重复n遍例如 24R 22222nI 表示在与它前后相邻的两个数之间插入n个线性插值点 例如 15I7 1234567xM 表示它前面的数据与x之积例如 54M 520nJ 表示从它所在位置跳过n项不指定的数据而使用缺省值 输入简写规则 如果n R I J 中的n缺省 则假设n 1 如果xM中的x缺省 则致命错误 nR前面必须放有一个数或由R或M产生的数据项 nI前面必须放有一个数或由R或M产生的数据项 后面还要跟有一个常数 xM前面
11、必须放有一个数或由R或M产生的数据项 nJ前面可以放除了I以外的任何内容 MCNP内容提要 入门 Introduction 输入文件基础 InputFileBasics 几何描述 GeometryDefinition 源的描述 SourceDefinition 计数描述 TallyDefinition 材料 物理成分和数据 Materials Physics Data 高级计数描述 AdvancedTallies 高级几何描述 AdvancedGeometries 各种简化 VarianceReduction 临界问题 Criticality 返回首页 MCNP几何描述 几何基础快速开始曲面组
12、合曲面Macrobodies栅元特性例子 几何基础 懂得四种定义曲面的方法懂得怎样由曲面创建栅元了解Macrobodies的定义细节懂得进行曲面变换懂得何时使用特殊曲面 几何基础 universe 根据材料和特性被分成不同的区域整个无穷的universe必须包括在几何模型之内几何的基本单位是栅元所有的栅元都由闭合曲面定义所有的曲面都能将universe分成两部分 曲面 由方程定义曲面由方程及参数确定例如 一个球心在原点半径为R的球jsoR平行于y轴半径为R的圆柱jc yxzR垂直于z轴的平面jpzz 常用面的定义 面定义格式如下 面序号面助记符 Mnemonic 入口参数 entries 如图
13、所示的面1 9等面 可描述如下 1py 1 5 平面 序号为1 平面方程为y 1 5cm2py13py24py2 55py3 56cy1 5 圆柱面 7cy28cy2 59cy1 栅元中的复合曲面 栅元中的点和曲面的关系通过栅元对曲面的坐向联系起来 和 曲面将universe分为两个半区布尔算符 将不同的半区与创建的栅元联系起来交 Intersection 联 Union 余 Complement 坐向 栅元中所有的点都通过坐向与定义栅元的曲面联系起来 坐向说明了栅元中的点在曲面的那一边 正的坐向对于开放曲面 平面 点在曲面的坐标轴正方向 对于封闭的曲面 球 圆柱等 点在曲面以外 负的坐向对于
14、开放曲面 平面 点在曲面的坐标轴负方向 对于封闭曲面 球 圆柱等 点在曲面以内 栅元的复合曲面 交 Intersection 同时满足两个坐向的空间算符输入 在两个曲面号中用空格2 1只表示同时满足坐向 2和坐向 1的空间区域 栅元的复合曲面 联 Union 任意满足两个坐向之一的空间算符输入 在两个曲面号中用冒号 2 1表示任意满足坐向 2和坐向 1之一的空间区域 栅元的复合曲面 余 Complement 表示栅元之外的空间算符输入 在曲面号前用 5表示栅元5之外的空间 余以后的区域可以和其他区域进行交和并的运算 2 5代表曲面2之内且在曲面5之外的区域 栅元 栅元输入卡包括三个部分栅元号
15、1 99999栅元号 写在第1 5列上栅元内容材料号材料密度 0 表示原子密度 0 表示质量密度复合曲面 Ccell Mat DensSurfacecombinations11 1 122 5 23012 3403 几何区域 cell 定义 区域号物质号物质密度区域定义 若干面围成 例111 7 91 23 45 6 区域120 7 1 区域232 8 97 8 区域3 说明 图中所示的区域1为一个长方体 由平面1 6围成 该区域为密度7 9g cm3的物质1 区域2位于球面7以内但不包括区域1 该区域为真空 区域3位于球面7以外且在球面8以内 由密度为8 9g cm3的物质2填充 例2 10
16、 1 220 12 21 例3 12 18 3 1 2 6 4 5 9 21 8 92 3 734 11 33 4 8 四类MCNP曲面 方程定义曲面平面 球 圆柱 圆锥 圆环 任意的二次曲面Macrobodies基于闭合图元的复杂曲面用点定义对称曲面平面 线性曲面或二次曲面由三个点定义一般平面一般平面 MCNP曲面 方程定义 基本格式jnalistj 曲面号 1 99999 写在第1 5列上n 缺省值为0 表示不进行任何坐标变换 0 用TRn卡对曲面坐标变换 0 曲面j是伴随曲面n的周期边界a 方程助记名list 方程描述的数据项 MCNP曲面 圆锥 圆锥的等式定义了两个 叶 参数中额外的条目是用来区分 正叶 和 负叶 的只有在圆锥平行于轴的时候才有效 MCNP曲面 Macrobodies 有限的 模块 构成的曲面BOX任意指向的正交框RPP直角平行六面体所有的表面垂直于各自的轴SPH球与方程表示的球是一样的RCC直圆柱体轴与底面垂直 但是方向任意RHP HEX 直六面棱柱与RCC相似但底为任意的六边形 MCNP曲面 Macrobodies 坐向 与其他封闭曲面相类似 正的坐向 点在
