实用MCNP3B教程.docVIP

大型通用中子-光子联合 输运蒙特卡罗模拟程序 郑华 编 大庆生产测井研究所 1998年10月 前言 本教程用于培训那些曾经、正在或将要从事核测井研究工作的人员，使他们掌握一种在计算机上模拟核测井过程的技术。 目前研究核测井问题的计算机模拟方法有两种，一是扩散理论方法，二是蒙特卡罗方法。扩散理论法对核测井过程建立玻尔兹曼扩散方程，然后求解。由于测井问题比较复杂，这个含有多个变量的微积分方程一般不能得到解析解，还要用数值方法进行求解，在求解时对粒子运动方向的角坐标可用级数展开或取离散角度处理，对空间和时间变量多采用差分处理。蒙特卡罗方法则对给定的问题建立相应的随机抽样模型，并用以一系列随机数跟踪大量粒子历程的方法完成对粒子输运的模拟。与确定论相比，蒙特卡罗方法能更好地适应复杂的几何条件。 蒙特卡罗方法模拟辐射输运的思想在40年代由美国Los Alamos实验室的科学家提出，76年开发出了通用程序MCNP。从80年代开始，美国一些核武器研究实验室和大学受雇于各大测井公司进行核测井数值模拟研究，编写出一批核测井专用蒙特卡罗模拟程序。专用程序的计算时效高，但针对性强，后来随着计算机速度的提高，一般人们还多使用通用程序进行计算。目前最常用的两个通用程序是MCNP和英国原子能署的MCBEND。 MCNP具有较强的通用性，可处理任意三维几何问题，提供了多种源分布和记数方式，使用精细的点截面数据。MCNP3版(83年)和3A版(85年)发行后，就成为用蒙特卡罗方法模拟核测井的最流行的通用程序，它解决了特征伽马谱线的问题，可以较好地模拟中子和光子联合输运，使用的主要核数据库是ENDF/B-4。3B是89年发行的版本，对3A版做了许多改进，并增加了几何重复定义功能。91年问世的4版可以联合模拟中子、光子、带电粒子(离子)的输运，可模拟探测器的测量结果。4版与以前版本相比在性能上又有较大改善，程序代码时效有显著提高，使用了更新的ENDF/B-6评价核数据库。 MCNP程序3B版传入中国已多年，可从中国原子能科学研究院计算机应用研究所得到拷贝。金文绵、李素梅等人将MCNP3B IBM 3033版移植到微机上，并翻译了原说明书第三章，对中国用户提供了很大帮助。 本教程主要内容来源于MCNP3B说明书的前三章，并加上部分本人应用心得。涉及数学方法和物理基础的部分仅作简要介绍，其中多数定义并不严格，目的是让读者易于接受。有关数学基础和蒙特卡罗方法在粒子输运问题中应用方面的详细知识可参阅裴鹿成和张孝泽的专著，中子、光子与物质相互作用的知识可从核物理教材中了解。 鉴于本教程有特定的应用范围，其中并不全面覆盖使用MCNP3B程序的所有知识。有关计算裂变核反应(如临界问题)、如何修改源程序和用户编写输入子程序(如定义用户源、用户记数)的内容本教程中一概不予介绍，还有一些我们不常用的MCNP功能(如几何描述重复定义、坐标变换等)也不介绍。对这些内容感兴趣的读者可以阅读MCNP3A和3B使用说明书。 由于本人才疏学浅，使用MCNP程序的经验又不多，编写过程也比较匆忙，本教程中不免存在错误和疏漏。如读者发现其中有差错之处或觉得部分叙述晦涩难懂，请与编者联系，电话是5970055(家)、5592052(办)。 郑华 1998/10/5 目录 页码 第一章 入门 1 §1.1 用蒙特卡罗方法模拟粒子输运 §1.2 MCNP简介 A. 程序的发展过程和应用领域 B. 程序的特点 §1.3 MCNP输入文件 2 A. 输入文件的基本形式 B. 一个简单的例子 C. 接续运行的输入文件 D. 卡片格式 3 E. 输入错误信息 §1.4 如何运行MCNP3B微机版 A. MCNP3B程序包的主要文件 B. 安装和运行MCNP3B的操作步骤 4 第二章 几何 §2.1 基础知识 概述 栅元定义中的一些概念 有关曲面的一些知识 §2.2 几何描述卡 5 栅元卡 曲面卡 §2.3 有效地构建几何 6 定义栅元的原则 检查几何错误 第三章 数理基础 §3.1 物理 粒子权重 粒子径迹 中子与物质作用 光子与物质作用 8 §3.2 记数 9 面流量记数 通量记数 栅元能量沉积记数 探测器记数 记数精度 10 §3.3 减小方差技巧 11 统筹考虑 能量截断 时间截断 几何分裂和轮盘赌 能量分裂和轮盘赌 暗含俘获和权重截断 强迫碰撞 12 DXTRAN 源变量偏倚 权重窗口 13 指数变换 相关抽样 点探测器 第四章 数据卡 §4.1 问题类型卡 14 §4.2 栅元参数和曲面参数卡 A. IMP 栅元重要性卡 B. VOL 栅元体积卡 C. AREA 曲面面积卡 D. PWT 光子产生权重卡 E. EXT 指数变换卡 F. VECT 矢量输入卡 G. FCL 强迫碰撞卡