实用MCNP3B教程

大型通用中子-光子联合 输运蒙特卡罗模拟程序

郑华 编

大庆生产测井研究所

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 基础知识

A.概述

B.栅元定义中的一些概念

C.有关曲面的一些知识

§2.2 几何描述卡 5

A.栅元卡

B.曲面卡

§2.3 有效地构建几何 6

A.定义栅元的原则

B.检查几何错误

第三章数理基础

§3.1 物理

A.粒子权重

B.粒子径迹

C.中子与物质作用

D.光子与物质作用8 §3.2 记数9

A.面流量记数

B.通量记数

C.栅元能量沉积记数

D.探测器记数

E.记数精度10 §3.3 减小方差技巧11

A.统筹考虑

B.能量截断

C.时间截断

D.几何分裂和轮盘赌

E.能量分裂和轮盘赌

F.暗含俘获和权重截断

G.强迫碰撞12

H.DXTRAN

I.源变量偏倚

J.权重窗口13 K.指数变换

L.相关抽样

M.点探测器

第四章数据卡

§4.1 问题类型卡14 §4.2 栅元参数和曲面参数卡

A. IMP 栅元重要性卡

B. VOL 栅元体积卡

C. AREA 曲面面积卡

D. PWT 光子产生权重卡

E. EXT 指数变换卡

F. VECT 矢量输入卡

G. FCL 强迫碰撞卡15

H. WWE 权窗能量卡

I. WWN 权窗边界卡

J. WWP 权窗参数卡

K. WWG 权窗产生器卡

L. WWGE 权窗产生器能量卡

M. PDn 探测器贡献卡

N. DXC DXTRAN贡献卡16 §4.3 源的描述

A. SDEF 通用源卡

B. SIn 源信息卡17

C. SPn 源概率卡

D. SBn 源偏倚卡

E. DSn 相关源分布卡18

F. SCn 源注释卡

§4.4 记数方式的指定19

A. Fna 记数类型卡

B. FCn 记数注释卡20

C. En 记数能量卡

D. Tn 记数时间卡

E. Cn 记数余弦卡

F. EMn 记数能量乘子卡

G. TMn 记数时间乘子卡

H. CMn 记数余弦乘子卡

I. DEn/DFn 剂量能量/剂量函数卡

J. CFn 记数栅元标志卡

K. SFn 记数界面标志卡

L. FSn 记数片段划分卡

M. SDn 记数片段的体积/面积卡

N. FQn 记数打印层次卡21 O. TFn 记数涨落打印卡

P. DD 探测器和DXTRAN诊断指定卡

Q. DXT DXTRAN的参量卡

R. FTn 记数特殊处理卡

§4.5 材料的指定22

A. Mm 材料成份卡

B. DRXS 离散反应截面卡23

C. AWTAB 原子量卡

D. VOID 否定材料卡

§4.6 能量和热处理方式的指定

A. PHYS 能量物理截断卡

B. ESPLT 能量分裂和轮盘赌卡

C. TMP 自由气体热处理卡

D. THTME 热时间卡

E. MTm S(α,β)材料卡24 §4.7 问题截断卡

A. CUT 截断卡

B. NPS 历史截断卡

C. CTME 计算时间截断卡

§4.8 外围卡

A. PRDMP 打印及转储周期卡

B. LOST 丢失粒子卡

C. DBCN 调试信息卡

D. PRINT 打印控制卡25 §4.9 MCNP输入文件综述

A. 输入卡

B. 存储限制第五章经验26§5.1 一般应用步骤

§5.2 需注意的问题

附录

连续能量中子截面库ENDL851数据目录

中子热截面库BMCC1数据目录

离散中子截面库D91数据目录27 光子截面库MCPLIB1数据目录

特殊材料S(α,β)热截面库TMCC1数据目录

参考文献

MCNP 3B 实用教程郑华编

第一章入门

MCNP3B版是一套模拟中子和光子联合输运的通用蒙特卡罗程序，它具有连续能量、三维几何和与时间相关的处理能力。本章仅简要介绍该程序的使用方法，详细内容将在后续各章中叙述。

§1.1 用蒙特卡罗方法模拟粒子输运

中子和光子在物质中输运的宏观表现是大量粒子与原子核微观作用的平均结果，蒙特卡罗方法通过逐一模拟和记录单个粒子的历程来求解输运问题。要得到比较合理的平均结果需要跟踪大量的粒子，至于单个粒子在其生命中的某一阶段如何度过，可以在已知统计分布规律的前提下通过抽取随机数来决定。这就象掷骰子赌博一样，因而得名蒙特卡罗方法*。

图1.1显示了模拟

中一个中子射入物质后

的随机历程。首先根据

中子与物质作用的物理

规律(分布函数)，选取

一个随机数决定中子在

何处与原子核碰撞，本

例中在1点碰撞；然后再

用抽取随机数的方法决

定中子与原子核发生了哪种反应，这里抽出的是非弹性散射反应；散射中子的能量和向哪个方向飞行也是用抽取随机数的方法从已知分布函数中决定的；碰撞过程中是否产生光子以及光子的能量、飞行方向等参数还是要通过抽取随机数从已知分布中决定，这里产生了一个光子。跟踪光子，确定它在7点与原子核碰撞并被吸收。散射后的中子在2点与原子核发生(n,2n)反应，其中一个出射中子射向探测器，另一个中子在3点被吸收。在2点的碰撞还产生了一个光子，它在5点又与原子核发生了一次散射反应，并离开物质。这一入射中子的历史过程结束了，有一个中子到达了探测器，感兴趣的结果被记录下来。跟踪越来越多的入射粒子历程后，平均结果就能反映出宏观效果。通过以上描述，读者不难领略蒙特卡罗方法如何通过跟踪粒子历程的方法计算问题，也了解了随机数在蒙特卡罗计算中的独特作用。

§1.2 MCNP简介

A. 程序的发展过程和应用领域

40年代美国Los Alamos实验室的Fermi、von Neumann 和Ulam等人提出用蒙特卡罗方法模拟辐射输运的思想。47年Fermi发明了第一台用蒙特卡罗方法计算中子链式反应的机器。从50年代开始，von Neumann领导一个小组研究输运问题的蒙特卡罗处理方法，编写出模拟中子输运的程序MCS。63年蒙特卡罗方法描述语言标准化。65年完成的中子输运程序MCN有了很大改进，使用了标准的截面库，并且具有复杂几何描述功能。后来，Los Alamos实验室又开发了模拟光子输运的程序MCG(高能)和MCP(能量低至1keV)。73年MCN和MCG合并成MCNG，为MCNP的雏形。MCNP于76年开发成功，77年6月发行。3B版之前还发行过1A、1B、2、2A、2B、2C、2D、3和3A版。

* Monte Carlo，又译蒙特卡洛，摩纳哥的世界著名赌城。

MCNP软件包(a general Monte Carlo code for Neutron and Particle transport)是一套通用的、模拟三维空间中连续能量的中子、粒子联合输运的程序，其名字早先来源于the analog Monte Carlo method for Neutron and Proton's transport 的缩写。MCNP3版(83年)和3A版(85年)发行后，这一软件在核测井领域逐渐成为最流行的通用程序，程序模拟结果和模型井实验结果较好地吻合，此时程序使用的主要核数据库是ENDF/B-4。88年发行的3B版程序增加了几何重构功能。91年4版问世，加入了模拟带电粒子(离子)输运部分，可以模拟探测器的测量结果，使用了新的ENDF/B-6评价核数据库。

MCNP程序的应用范围十分广阔，主要包括：反应堆设计、核临界安全、辐射屏蔽和核防护、探测器的设计与分析、核测井、个人剂量与物理保健、加速器靶的设计、医学物理与放射性治疗、国家防御、废物处理、射线探伤等。

B. 程序的特点

MCNP3B版程序能处理的中子能量范围是10-11MeV 至20MeV，光子能量范围1keV至100MeV。程序中物理量的单位规定为：

物理量单位

长度cm

能量MeV

时间shake, =10-8sec

温度MeV (kT)

原子密度atoms/barn-cm

质量密度g/cm3

截面barns (10-24cm2)

加热量MeV/collision

此外，原子质量按中子质量为1.0计算，这种单位下阿佛伽德罗常数是0.59703109；程序运行时间以分钟为单位。

MCNP程序的源代码是用FORTRAN语言编写的，3B 版主程序有6万多行。程序包中携带了大量的核反应数据库文件，我们得到的库文件约占80Mb磁盘空间。

MCNP3B版具有很强的通用性，主要体现在：

1.可以处理任意三维几何结构的问题。在输入文件INP中，空间被曲面(surface)分割成相互邻接的区域，称为栅元(cell)，可以给栅元填充各种物质。栅元的界面可以是平面、二阶曲面或某些四阶曲面(如椭圆环状面)。有关MCNP几何描述的知识将在第二章介绍。

2.可以模拟中子输运、光子输运和二者联合输运。

3.用户可以非常方便地在任何位置指定体源、面源、线源或点源，设置源粒子位置、能量、时间、飞行方向等参数的分布。源的指定方法将在第四章中介绍。

4.程序提供多种记录模拟结果方法，包括通过某一界面的粒子流量或通量、进入某一栅元的通量、沉积能量和点通量。模拟结果在MCNP中称为记数(tally)，可以按位置、能量、时间、粒子来向和粒子种类记数。有关记数的指定方法将在第四章中介绍。

5.程序包携带了大量核反应数据库文件，包括连续和离散的中子截面库、光子点截面库、热中子点截面库等，几乎可对所有天然物质进行计算。程序能比较精细地模拟中子和光子输运过程，并对一些特定的物理过程允许用户选择使用哪种方式进行处理，如对热中子处理可选用自由气体模型或S(α,β)模型，对低能光子处理可以考虑或忽略相干散射等。

6.为了提高计算时效，给用户提供了许多可选用的减小方差(variance)技巧，主要包括：重要抽样、权重截断和轮盘赌、时间和能量截断、模拟俘获、指数变换、强迫碰撞、能量分裂和轮盘赌、源的偏倚、点探测器记数、确定

图1.1 单个中子随机历程示意

论输运、权窗等。有关减小方差技巧的理论基础在第三章中阐述，应用方法在第四章中介绍。

7.用户可通过设置源粒子数或运行时间来通知程序何时终止运行，还可以在原有计算结果的基础上接续运行程序。一些结果不会因计算的意外中断而丢失。

8.在输出文件OUTP中给用户提供丰富的信息，包括输入列表、使用的截面表、粒子生成和丢失表、栅元中的粒子活动情况、中子诱发光子表、记数和记数涨落表等，还可以根据用户要求给出其它信息。

9.提供了简单的问题调试工具。

§1.3 MCNP输入文件

若用户不修改源程序，MCNP的输入文件包括截面数据库文件、XSDIR文件、INP文件等。INP文件是用户要填写的主要输入文件，一般把该文件特指为输入文件。OUTP 是MCNP的主要输出文件，其它输出文件还有转储文件RUNTPE、运行信息文件OUTPUT等。

A. 输入文件的基本形式

INP文件有初始运行和接续运行两种形式，用户须在INP中描述问题的几何、材料、源、记数和其它要求。INP 由一些被空行分隔的输入块组成，主要的输入块是信息块、标题和栅元块、曲面块和数据块等。输入块又由一些被称为卡的输入行组成。初始运行输入文件的格式如下：

信息卡

} 可选

空行分隔符

标题卡

栅元卡

.

.

.

空行分隔符

曲面卡

.

.

.

空行分隔符

数据卡

.

.

.

空行结束符

其它

信息卡的1~8列应填写MESSAGE:，后面跟着用空格分隔的参数项。可用A=B参数项更改输出文件名，如OUTP = MYOUT。信息块是可选的。

在信息块之后的第一行是问题的标题卡，它仅限于一行，占用1~80列，可以是任何信息，将作为OUTP文件中各个输出表的标题被复制。

用户在栅元块和曲面块中描述问题的几何。栅元由栅元卡描述。空间必须由彼此相邻的栅元填满，栅元之间不能重叠，也不能出现无栅元的空区，否则会出现错误。构建栅元的曲面由曲面卡定义，曲面卡在曲面块中给出。曲面卡和栅元卡的填写方法见第二章。

曲面块之后是数据块，在数据块中用户描述源、记数方式、材料等。数据卡在第四章详细介绍。

无论数据块后是否有空行结束符，MCNP都能运行。用户可以把希望保留的一些附加信息写在数据块的空行之后，MCNP会将它们复制到OUTP文件末尾。

B. 一个简单的例子

为说明如何填写INP文件，这里例举一个简单问题。如图1.2所示，在一个边长10cm的石墨立方体3中有两个半径0.5cm的球形空间，球1中充满氧气，球2是铁球。在球1中置一14MeV各向同性中子点源，计算球2外表面与能量相关的中子通量。建立的INP文件如下：

SAMPLE PROBLEM INPUT DECK

1 1 –0.0014 -7

2 2 –7.86 -8

3 3 –1.60 1 –2 –3

4 –

5

6

7 8

40 -1:2:3:-4:5:-6

←空行

1 PZ -5

2 PZ 5

3 PY 5

4 PY -5

5 PX 5

6 PX -5

7 S 0 –4 –2.5 .5

8 S 0 4 4 .5

←空行

MODE N

IMP:N 1 1 1 0

SDEF POS=0 –4 –2.5 ERG=14

F2 8

E0 1E-5 1E-4 1E-3 .01 .1 1 2 3 4 5 6 7 8 9 10 11 12 13

M1 8016 1

M2 26000 1

M3 6000 1

NPS 10000

←空行

本例中没有信息块，

第一行是标题卡，之后至

空格前为栅元块。栅元卡

上依次填写栅元号、材料

号、密度和构成栅元界面

的曲面号(带正负号)，这里

定义了4个栅元：栅元1由

球面7围成，里面填充材料1(16O2气体)，密度是0.0014g/cm3；栅元2由球面8围成，填充材料2(铁)，密度7.86g/cm3；栅元3由平面1、2、3、4、5、6围成，不包括球面7、8以内的空间，填充材料3(石墨)，密度1.6g/cm3；栅元4是栅元3以外的空间，为真空。

曲面卡上需要填写曲面号、曲面类型和曲面参数，本例中定义了8个曲面，前6个为与原点距离5cm垂直于各坐标轴的平面，后两个是半径0.5cm的球面，球心分别在(0,-4,-2.5)和(0,4,4)。

数据块中指定了问题类型、源、记数方式、材料和运行粒子数，各卡数据项的意义如下：

MODE卡问题类型是中子输运

IMP卡4个栅元的中子重要性分别是1 1 1 0

SDEF卡位于(0,-4,-2.5)、能量14MeV的各向同性点源

F2卡在曲面8上做中子通量记数

E0卡对记数能量分区，1~14MeV之间间隔为1MeV，

1MeV~10-5MeV之间间隔为一个数量级

M1卡材料1是16O核素

M2卡材料2是Fe元素

M3卡材料3是C元素

NPS卡运行源粒子数10000

以上例子仅用于说明INP文件格式，有关各输入卡的详细内容，读者在此不必深究。

C. 接续运行的输入文件

为了继续一个早先被终止的计算，或想对早先终止的计算进行重新编辑打印(指向文件输出，即PRINT语句，不是向打印机输出)，可以接续运行MCNP。

接续运行INP文件中必须有信息块，信息卡中必须有C 或C m参数，其中m表示从第m次转储开始接续运行，m缺省表示从最后一次开始。接续运行输入文件的格式为：

4

3

2

1

8

7

43

2

1

图1.2 例子的几何示意图

MESSAGE: C m

空行分隔符

CONTINUE

数据卡

.

.

空行结束符

其它

此时INP中不能有标题卡、栅元块和曲面块。数据块第一行必须是CONTINUE，允许使用的数据卡只是NPS、CTME、PRINT、FQ、DD、PRDMP、LOST和DBCN卡。NPS卡参数为累积运行的总源粒子数，CTME卡参数相对于本次运行的开始时间。如果给NPS置负参数，表示仅从转储文件中提取数据重新打印。

接续运行的另一个必要条件是上次运行的转储文件RUNTPE还没被删除。因每次运行都生成新的转储文件，若想对一次运行结果多次重新打印，则应把RUNTPE文件做备份。各次转储时运行的源粒子数和时间可从OUTPUT 文件中查到。

D. 卡片格式

INP文件的每一行都限于使用1~80列。卡片都可以按行填写，数据卡也允许按列填写。

a) 行输入格式

通常卡片的1~5列用于填写栅元号、曲面号或数据卡的助记名，6~72列填写卡片参数，73~80列为注释，$符号之后也为注释。序号或卡片助记名可以写在1~5列的任何地方。带有粒子标识符助记名可能需要5列以上，但冒号必须写在6列以内。如果1~5列空白，表示本行为上一行的接续行。6列之后可以写数据项，多个数据项之间用空格分隔。一个数据项必须在一行上写完，不得接续到下一行。相同编号的卡片只能有一张。

填写卡片参数时可以使用以下书写功能：

1.nR功能，表示它前面的数据再重复n遍。

2.nI功能，表示与其前后相临两个数之间插进n个线性

插值点。

3.xM功能，表示数值等于它前面数据的x倍。

4.nJ功能，表示从本项开始的n个数据使用缺省值。

b) 列输入格式

列输入只能用于数据块中，对栅元参数和源的描述比较有用。按行输入的栅元重要性、体积、权窗等数据项可读性较差，而且增加或删除栅元时要在行输入卡上仔细寻找相应项。列输入的可读性有很大提高，删除或增加与某一栅元相对应的数据项时也比较方便。

列输入格式的第一行以#开始，#可以放在1~5列的任意位置，卡片助记名逐个放在该行6列以后，在这些助记名之下按列给出数据项。同一个列输入格式块中的卡片必须是同一类卡片，比如都是栅元参数卡、都是曲面参数卡或都是源参数卡等，在#号下面的1~5列放置栅元号、曲面号或源分布号。

c) 粒子标识符

粒子标识符是卡片助记名的一部分，:N表示中子，:P 表示光子，有时也能遇到表示中子-光子联合输运的:N P。下面一些数据卡需要粒子标识符：IMP、EXT、FCL、WWN、WWE、WWP、WWCE、DXT、DXC、Fna、PHYS、ESPLT 和CUT。

d) 缺省值

MCNP许多输入卡的参数项有缺省值，用户不必每次都给出这些参数，如果卡片输入项有固定顺序，可以使用nJ功能跳过n个输入项。如果卡片上所有数据项都想缺省，只给出卡片助记名即可。有些卡片不给出也有缺省值，如MODE N卡就可以省略。

E. 输入错误信息

MCNP对输入文件作广泛的检查，如果用户违反输入规定，将在输出文件OUTP中打印致命错误信息(FATAL ERROR)，程序停止运行。第一个致命错误是真实的，后面的致命错误不一定是真实的。

MCNP还给出许多警告性的提示(WARNING)，对这些提示用户也不可忽视。

如果遇到任何代码损坏，如零作除数，MCNP将给出BAD TROUBLE信息，并终止运行。

§1.4 如何运行MCNP3B微机版

通常编辑INP文件后，键入RUN386 MCNP并回车，就可运行MCNP程序。3B微机版源程序用NDP FORTRAN 386编译器编译，要求386以上的微机和DOS3.3之后的操作系统，不能在Windows系统或保护方式下运行。

A. MCNP3B程序包的主要文件

1.MCNP.DA T：MCNP3B源程序MCNP.F的缩写型文件。文

件中所有公用说明语句只列出一次，便于修改。

2.PRPR.F：辅助程序，把MCNP.DAT展成MCNP.F。

3.MAKXSF.F：该程序用于对NCMP所需截面数据加工处

理，从原有的截面数据文件中挑选若干常用的截面数据，生成新的较小文件，可以节省磁盘空间和缩短程序运行时的数据读取速度。

4.SPECS：是MAKXS.F程序的输入文件，用来指定加工截

面数据文件的任务清单，包括加工格式、输出文件名和核素。

例如，某SPECS文件的内容是：

XSDIR1 XSDIR2

MCPLIB2 2 2048 512

1000.01P

2000.01P

1600.01P

BMCCS2 2 2048 512

1003.03C

2004.03C

92235.11C

表示原截面库目录文件名是XSDIR1，新截面库目录文件是XSDIR2(为二进制的直接存取文件)，MCPLIB2中包含

H、He和S的光子截面数据，BMCCS2中包含3H、4He和

235U的连续能量中子截面数据，它们必须在XSDIR1目录规定的数据库中存在。

5.XSDIR1：它也是MAKXS.F程序的输入文件，是现有全

部截面数据库目录。

6.BMCCS1：连续能量中子截面数据。

7.D91：离散能量中子截面数据。

8.MCPLIB1：光子截面数据。

9.TMCCS1：热处理中子的数据。

10.ENDL851：连续能量中子截面数据。

11.531DOS1：中子剂量数据。

12.532DOS1：中子剂量数据。

13.LLLDOS1：中子剂量数据。

14.PROB1~PROB18：例题。

15.OUTP1~OUTP18：例题结果。

B. 安装和运行MCNP3B的操作步骤

1.拷贝原盘上的文件PRPR.F、MAKXSF.F、MCNP.DA T、

XSDIR1及所需截面数据库文件。

2.安装NDP FORTRAN 386编译器。

3.建立MCNP输入文件INP。

4.建立MAKXSF的输入文件SPECS。

5.编译PRPR.F和MAKXSF.F：

F77 PRPR.F

F77 MAKXSF.F

6.运行MAKXSF.EXP：

RUN386 MAKXSF

生成新的截面库目录文件XSDIR2和新截面数据库文件。7.展开MCNP.DA T：

RUN386 PRPR

生成MCNP.F文件。

8.编译运行MCNP.F：

COPY XSDIR2 XSDIR

F77 MCNP.F

RUN386 MCNP

运行结束生成结果文件OUTP、转储文件RUNTPE、运行信息文件OUTPUT等。

第二章几何

本章介绍MCNP几何构建方法，INP文件中栅元卡和曲面卡格式，以及用户定义栅元时应注意的事项。

§2.1 基础知识

A. 概述

用户在三维正交坐标系(见图2.1)

下定义曲面，用曲面围成栅元，在栅

元中填充物质。栅元可由一些曲面包

围区域的交集和并集构成。曲面可通

过设置解析方程参数或指定曲面上一

些已知点定义。

B. 栅元定义中的一些概念

a)坐向(sense)

定义栅元时，一个重要的概念是坐标点相对于曲面的坐向。假设s=f(x,y,z)=0是一曲面S的方程，对于任意一点(x’,y’,z’)，代入s=f(x,y,z)后，若s=0则该点在曲面S上。若该点不在曲面S上，有两种情况：s>0时称该点相对于曲面S 具有正的坐向，反之称该点相对于S具有负的坐向。例如，球面以外的点相对于球面具有正坐向，沿坐标轴正向离开垂直于坐标轴平面的点相对该平面具有正坐向。

任何曲面都把空间分成两部分，一部分相对于该曲面具有正坐向，另一部分具有负坐向。填写栅元卡时，在曲面号之前用正负号表示栅元中的点相对于该曲面的坐向，正号可省略。

b)交和并(intersection & union)

交集是两个集合的公共部分，如图

2.2a所示。并集是两个集合的合集，见图

2.2b。MCNP中交算符是隐含的，在栅元

卡上两个曲面号之间只有空格就表示交

运算。并算符用冒号: 表示。

若想仅用交算符定义栅元，栅元内所有点对特定的界面必须有相同的坐向。如图2.3所示，1、2、3、4、5面围成一个实体，用密度为1.0g/cm3的

材料1填充，外部是真空。实体内

P1和P2点相对于4面有负坐向，P3

点相对于4面有正坐向；P1点相对

于3面有正坐向，P2和P3点相对于3

面有负坐向。这时最好借助另一平

面6把实体分成两部分，从而顺利

定义栅元：

1 1 -1.0 1 –

2 –

3 6

2 1 -1.0 1 –6 –4 5

然而，只用交算符还不便定义实体之外的区域，这时用并算符比较方便，实体外可由3、4面以上区域的交集和其它几个面以外区域的并集构成一个栅元：

3 0 -5 : -1 : 2 : 3 4

当然实体之内的区域也可以借助并算符定义，把1、2栅元合并(这样定义会降低计算速度)。由于交运算优先，所以要用括号把并运算括起来：

1 1 -1.0 1 -

2 (-

3 : -4) 5

c) 补(complement)

补算符#的用法有两种：

1. #n 表示当前定义的栅元是栅元n其余的部分。

2. #(...) 表示括号内曲面号定义区域之外的部分。

引入坐向概念后，补运算就不是交运算和并运算之外的新概念，而是交和并的另一种表示方法。例如，对上面第一种情况的处理过程是，去掉#号，把n用括号括上，n 中的交算符变成并算符，n中的并算符变成交算符而且加上背靠背的括号) (，曲面之前的坐向取相反符号。数学中常用?表示交运算、用+表示并运算、ˉ表示补运算，有b

a

b

a?

=

+，b

a

b

a+

=

?，此处a和a则用坐向指定。

使用补算符时要格外小心，用户常因对补算符理解不透而弄出几何逻辑错误，因此建议不使用补算符。另外，#(n)是不合规定的。

C. 曲面定义中的一些概念

a) 模糊面(ambiguity surface)

有时需借助栅元界面之外的曲

面来描述栅元。如图2.4所示，假设图

中实体绕y轴旋转对称，栅元1和2相对

于球面都有正坐向、相对于平行于y

轴的柱面都有负坐向。这时，可用垂

直于y轴的平面a(y=0)来区分栅元1和2，平面a被称为模糊面。模糊面是为了使栅元描述不模糊而设置的曲面，它不是本栅元界面，但通常是其它栅元的界面，有时需要借助不止一个模糊面。

图2.1 迪卡尔坐标系

图2.2 交和并

图2.3 例子

图2.4 模糊面的作用

b) 叶(sheet)

锥面和椭环面分多叶。通常用方程定义的锥面含有两个叶，相对于锥面旋转对称轴所平行的坐标轴而言，一个具有正的展开方向，另一个具有负向斜率。一个使用二叶锥面描述的栅元常要借助模糊面来区分栅

元在锥面的哪一个叶。为了减少使用模糊

面，可以在曲面卡上只定义锥面的一个叶，

用+1表示正斜率的叶，-1表示负斜率的叶。

这种描述方法仅适用于旋转对称轴平行于

坐标轴的锥面。

D. 思考题

1. MCNP用什么方法描述几何。

2. 如何用曲面描述栅元几何。

3. 模糊面有什么作用。

§2.2 几何描述卡

在输入文件INP中，栅元卡用于栅元块，曲面卡用于曲面块。材料卡出现在数据块中，将在第四章中介绍它。

A. 栅元卡(cell card)

格式: j m d geom param

j: 在1~5列上填写的栅元标号，可以不连续。MCNP

按照读入的顺序对栅元另行编号，称为栅元的程

序编号，从1开始单调上升。为避免程序编号和

栅元标号混淆，建议按顺序给出栅元标号m: 该栅元的材料号，是材料卡Mm中相应的序号。

真空栅元m=0

d: 栅元材料的密度。填正值表示原子密度，单位atoms/barn-cm。负值表示质量密度，g/cm3。对

于真空栅元，该项不填写

geom: 栅元几何说明，列出界定该栅元的曲面的号(有符号，表坐向)，和描述这些曲面包围区域之间

关系的布尔算符。布尔算符包括交算符(缺省)、

并算符( : )和补算符(#)。缺省运算顺序是，先补，

再交，最后并。可用括号改变运算次序param: 可选的栅元参数说明，形式为：关键字=某一值

在栅元卡上可以定义栅元参数，以替代在输入块中定义栅元参数。这里允许的关键字是：带有粒子标识符的IMP、VOL、PWT、EXT、ECL、WWN、DXC和TMP。栅元卡参数项上的等号可用空格代替。若在栅元卡上指定了某种栅元参数，则在数据块中不能重复指定。

B. 曲面卡(surface card)

a) 由方程定义曲面

格式: j n list

j: 在1~5列上填写的曲面标号，可以不连续

n: 方程助记名

list: 按规定次序给出的方程数据项

下表列出了MCNP识别的曲面的类型、方程、助记名和数据项次序：

b) 用点定义对称曲面

格式: j n list

j: 在1~5列上填写的曲面标号，可以不连续

n: =X、Y或Z，

list: 一至三个点的坐标，每个点用一对值表示

类型为X、Y或Z的曲面卡用坐标点而不是用方程来分

别描述关于X、Y或Z轴对称旋转的曲面。如果这个曲面由

一叶以上组成，指定的点必须都在一叶上。每对坐标值定

义曲面上的一个点，例如，在Y卡上可以给出

j Y y1 r1y2 r2

其中(y i, r i=(x i2+z i2)1/2)是i点的坐标。给出点数不同，曲面类

型也不同：i) 一对坐标定义一个平面(PX、PY、PZ)；ii) 两

对坐标定义一个线性曲面(PX、PY、PZ、CX、CY、CZ、

KX、KY、KZ)；iii) 三对坐标定义一个二次曲面(除以上类

型还有SO、SX、SY、SZ和SQ)。

c) 用点定义一般平面

格式: j P x1 y1 z1 x2 y2 z2 x3 y3 z3

j: 在1~5列上填写的曲面标号，可以不连续

x1 y1 z1: 定义该平面的坐标点

MCNP对用户指定的P型曲面，将检查数据项个数，

若是4项，则作一般平面方程Ax+By+Cz-D=0的系数理解；

若是9项，则作为平面上三个点的坐标值理解，所产生的平

面方程系数遵循以下原则：1) 应使坐标原点相对于该平面

具有负坐向；2) 当平面通过原点(D=0)，应使(0,0,∞)相对

该平面具有正坐向；3) 若以上两项无法满足(D=C=0)，应

使(0,∞,0)相对该平面具有正坐向；4) 若以上三项无法满足

(D=B=C=0)，应使(∞,0,0)相对该平面具有正坐向；5)若以

上四项都无法满足，那么用户指定的三个点一定处于一条

直线上，则打印警告信息。

C. 思考题

1. 栅元卡的格式如何。为什么最好按顺序指定栅元标号。卡上材料号

对应什么，密度数据项填负值表示什么。

图2.5 二叶锥面

2. 如何用曲面卡定义曲面。

§2.3 有效地构建几何

即使是十分有经验且细心的用户，在构建比较复杂的几何时，也经常会犯一些错误。有些错误程序可以检查出来，有些则不能，以致造成错误结果。几何构建还直接影响程序运行效率，不合适的几何构建将浪费机时。

A.定义栅元的原则

原则1：栅元界面要少

尽管可用大量的并算符

定义一个几何上比较复杂的

栅元，但这样做并不明智。

问题在于，每次计算栅元内

碰撞之间的径迹时，都计算

径迹和栅元所有界面是否相

交以判断粒子是否进入其它栅元，太多的栅元界面会浪费机时。如，图2.6a中的几何只是一些平行的柱面，较容易定义，然而两个大柱面之间的栅元有14个界面(包括前后界面)；图2.6b则是比较有效的栅元定义，复杂界面的栅元被分隔成一些小栅元，每个栅元的界面有所减少。

原则2：栅元大小要合适

构建物理空间时，除按不同材料区域分割栅元外，还必须考虑效率因素。为了节省运算时间，栅元不能定义得太大或太小。在比较重要的区域，建议调整栅元大小，使粒子进入相邻栅元后数量约减少一半，以利于使用几何分裂和轮盘赌技巧。

原则3：尽量用简单曲面定义栅元

如前所述，MCNP频繁计算碰撞点到界面距离。点到二阶以上曲面的距离一般不能用解析法求解，用牛顿下山法求解时，多次叠代很浪费时间。如果问题条件允许，尽量使用二阶以下曲面分段替代二阶以上曲面。

原则4：避免曲面和曲面相切

即使使用双精度变量，也不能计算出点到曲面的确切距离，因而在曲面和曲面相切处有时会因计算精度不够导致栅元重叠或空缺，造成粒子丢失。虽然这种粒子丢失几率很小，一般不会影响记数结果，但计算大量粒子问题时可能会因为丢失粒子数量超过限制，使程序停止运行。

B.检查几何错误

MCNP检查INP文件的数据时，不能检查出几何逻辑错误。而运行中粒子丢失时，才能检查出栅元的重叠或栅元之间的空缺。因此在正式运算前，应采取以下办法检查几何逻辑错误：

1)增加一个VOID卡。这个卡将废弃这个问题的材料

说明，将全部栅元设为真空，并把加热记数转化为

通量记数。

2)对这个问题增加一个大球面，这个球面包围需计算

的系统，球面把系统之外的栅元分成两个栅元。球

内面是现在问题的边界，内部所有栅元的重要性都

置为1。

3)源描述改为：SDEF SUR=m NRM=-1，其中m是刚

才定义大球面的标号。

由于上述改动使粒子在计算系统中没有碰撞，一个短时间运行就会使大量粒子的轨迹穿过系统，可能在出现几何错误的地方引起粒子丢失。

粒子丢失时，OUTP文件上将产生丢失粒子打印表。表中列出丢失粒子跨越的所有曲面，并告诉你粒子在什么位置上向什么方向运行后丢失。你可以从这些信息中推断出粒子丢失的原因。

C. 习题

1. 请描述如图

2.7所示井

况模型，为高1m、半径

70cm的圆柱体。地层半径

为10cm至70cm，材料号为

1，密度2.65g/cm3。水泥

环厚30mm，材料号为2，密度1.95g/cm3。套管外径140mm，厚7.72mm，材料号为3，密度7.86g/cm3。井内介质材料号为4。为缩短计算时间，地层栅元的大小取高5cm、环距7.5cm，其它栅元高度也取5cm。本题中共有220个栅元。

第三章数理基础

本章介绍粒子在物质中输运的处理，记数的设置，以及减小方差的技巧。

§3.1 物理

A. 粒子权重(particle weight)

若严格模拟实际过程，则一个MCNP粒子仅代表一个实际粒子，且具有为1的权重。然而为了提高计算效率，MCNP允许用户使用许多并非实际过程的偏倚技巧。偏倚就是使一些感兴趣区域(包括位置、时间、能量等方面的区域)中粒子数增加，另一些区域中粒子数减少。为使计算结果不出现偏差，人为地改变粒子数后，相应地要调整每个粒子的权重。所谓粒子权重，就是每个模拟粒子携带的一个数w，它代表该粒子对最终结果有w份贡献，或者说该粒子相当于w个实际粒子。w不一定是整数，有时小于1。

不论模拟多少个源粒子，MCNP都给出相对于一个源粒子的记数。如果用户希望得到一定源强条件下的结果，可在SDEF卡上改变源粒子权重，指定WGT=w，这时记数结果是相对于发射w个源粒子。

B. 粒子径迹(particle tracks)

径迹是源粒子和诱发粒子在其历程中的运行轨迹，粒子在两次与原子核碰撞之间的运行线路是直线，粒子自身之间的作用忽略不计。在任何给定材料的栅元中，MCNP 计算沿径迹方向到下次碰撞点的距离而不是计算单位体积内碰撞几率。假设微观总截面是∑t，要抽样的径迹长度是l，可以证明，l和抽取出的[0,1)区间随机数ξ之间的关系是l=-lnξ/∑t(3.1)

C. 中子与物质作用

中子不带电，它与物质中原子的相互作用主要表现为与原子核作用，与核外电子的作用可忽略。不同能量的中子与原子核的作用有不同特点：1)能量在0~1keV之间的慢中子与原子核主要发生弹性散射和俘获反应。热中子能量kT=0.025eV(k是玻尔兹曼常数，T是绝对温度，取293K)，它是和原子热运动达到平衡时的能量，其速度分布接近麦克斯韦分布。能量≥0.5eV的中子称为超热中子。能量在1eV 至1keV之间的中子称为共振中子，这种中子被原子核强烈地共振吸收，吸收截面很大。2)能量在1keV~0.5MeV之间的中能中子与原子核的主要作用形式是弹性散射。3) 快中子的能量在0.5MeV~10MeV之间，与原子核主要发射弹性

b a

图2.6 栅元界面示意

图2.7 井况模型

散射和非弹性散射反应。4)能量高于10MeV 的特快中子与原子核发生弹性散射、非弹散射反应等。MCNP3B 模拟中子的能量范围是10-11~20MeV 。

当能量很高的中子与大块物质作用时，经过多次与原子核碰撞，能量不断减小，最后或者被俘获，或者变成热中子。低能中子和原子核碰撞时受热运动的影响，也受碰撞原子核周围原子的影响。MCNP 一般用自由气体模型处理原子核热运动，这种模型下弹性散射截面在零温度截面基础上作修正；其它中子反应的截面不随温度变化。对于某些常用原子，MCNP 还提供一些温度条件下的考虑了晶格和化学键影响的精细S(α,β)表。如果中子能量足够低，使用S(α,β)表后，总截面是常规截面表中俘获截面以及S(α,β)表中弹性截面、非弹性截面的总和。否则，总截面从常规截面表中得出，随温度变化作调节。

MCNP 模拟中子输运时，要跟踪径迹反复考虑如下事项：a) 碰撞前中子的能量、飞行方向和径迹长度；b)与哪种核碰撞，被碰撞核的速度与方向(除非使用S(α,β)处理)；c) 碰撞过程是哪种反应，如弹性散射、(n ,xn )、(n ,p )、(n ,α)、俘获等；c) 碰撞过程选择地产生光子，模拟每个光子的初始能量、运动方向、权重；d) 碰撞后中子(或几个中子)的能量、飞行方向、权重，下一次在什么位置和哪种核碰撞。如此反复，直至中子能量足够低或权重足够小。

1) 选择碰撞核

设介质由n 种核组成，ξ是一个[0,1)上的随机数，当

∑∑∑-===∑<∑≤∑

1

1

1

1

k i n i k

i ti

ti ti

ξ

(3.2)

成立时，第k 种核被选择为碰撞核，∑ti 是栅元中第i 种核的微观总截面。

2) 靶核速度抽样

如果不使用S(α,β)处理，需要考虑靶核速度。对轻核，若中子能量高于400kT (10eV)，对重核高于40kT (1eV)，靶核的速度对碰撞截面的影响可忽略不计。其它情况下，要按麦克斯韦分布抽取靶核速度，

)exp(4)(22232/1V V V P β-βπ=-

(3.3) 其中β=(AM n /2kT )1/2，A 是靶核质量，M n 是以MeV-sh 2

/cm 2为单位的中子质量。靶核运动方向从各向同性分布中抽取。然后得到动能为E 的中子的等效散射截面

??μ

σ=

σ2)()(1)(d dV V P v v v E rel rel s n eff s (3.4) 其中v rel 是中子和靶核间的相对速度，σs (v rel )为相对速度的散

射截面(从截面数据库中查出，未作温度修正)，μ是中子入射和出射方向夹角的余弦。相对速度的计算公式是，

2/122

)2(T n n rel V v V v v μ-+=

(3.5)

其中v n 是中子速度，V 是靶核速度，μT 是二者速度夹角的余弦。

3) 弹性散射和非弹性散射

如果碰撞时中子未被俘获，它就参与弹性散射或非弹散射反应，广义的非弹散射可以是(n ,n’)、(n ,np )、(n ,2n )等。弹性散射以)/(in el el σ+σσ的概率被选择，非弹散射以余下的概率被选择，其中σel 是弹性散射截面, σin 是非弹散射截面。

出射粒子方向：对弹性和非弹散射，抽取出射粒子飞行方向的方法是一样的，入射和出射方向夹角的余弦μ从截面文件的角分布表中抽取。对每种特定的中子入射能量，角分布表都有32个分立的μ箱。如果入射能量是E ，介于有角分布表的能量E n 和E n+1之间，则以(E -E n )/(E n+1-E n )的几率从n 表中抽取角分布，以(E n+1-E )/(E n+1-E n )的几率从表n +1中抽取角分布。抽取随机数ξ，找到第i 个余弦箱，使得i -1≤32ξ

μ=μi +(32ξ-i )(μi+1-μi ) (3.6) 出射和入射方向之间的方位角则从各向同性分布抽取，

出射粒子能量：对于弹性散射和分立能级非弹散射，出射粒子的能量由两体动力学规律确定；对其它的非弹散射，出射粒子能量的确定要依靠截面表。

弹性散射出射粒子能量等于

2

2)1(21A A A E E cm

in

out +μ++= (3.7) 其中，E in 是入射中子能量，μcm 是质心系下入射和出射夹角余弦。

非弹散射情况比较复杂，要按截面比例分别考虑每个出射粒子的权重，抽取出射能量的规则有24条，常用的一些规则是：

规则1 象抽取入射和出射方向夹角余弦一样，出射能量从

一些等几率的能量箱中抽取。

规则2 分立光子能量。假设能级差是EG ，则瞬发光子能

量是EG ，俘获辐射光子能量是EG +[A /(A +1)]E n ，其中A 是原子质量，E n 是入射中子能量。

规则3 (n ,n’)反应。用’表示质心系，无’表示实验室系，

])

1([)1(2A

A Q E A A E in out +++='

])1())(1(2[2

2/1''

+'+μ++=A E E A E E E out in in out out

(3.8) 规则4 表分布。对每种特定的入射能量，表包括三部分，

即次级能量E (I )、概率密度函数PDF (I )和累积密度函数CDF (I )。常用内插方法得到出射能量，抽取随机数ξ，若CDF (K )≤ξ

?-+-+++=')()1()

()1(2)({[)(2K E K E K PDF K PDF K PDF K E E

))

()1()

()1(/()}())]((2/1K E K E K PDF K PDF K PDF K CDF -+-+--ξ?

只有PDF (K+1)=PDF (K )时才使用直方分布。

规则5 一般蒸发谱。函数g (x )随x 变化关系用表给出，能量

随入射粒子能量E in 的关系也用表给出，规则是

))

(()(in out out in E T E

g E E f =→

(3.10) 密度函数由以下方式抽取 )()(in out E T E ξι=

(3.11)

其中T (E in )是入射能量的表函数，ι(ξ)是随机数的表函数。

4) S(α,β)处理

若中子能量足够低，且用户指定使用S(α,β)表，则碰撞通过S(α,β)处理被模拟，而不经过第3)步。中子与原子核碰撞受到原子热运动及附近其它原子的影响，S(α,β)同时考虑了二者的影响。

a)对于非弹散射，对应入射能量表的特定能量，出射能量分布由等概率的确定能量表示，覆盖4eV 到10-5eV 的能区，伴随着一系列角分布数据。若入射能量为E ，出射能量E ′从分布

()[]∑=++?ρ--ρ-'δ=<<'N i j i j i i i E E E N E E E E P 1

,1,11 1)| ( (3.12) 抽取，其中E i 和E i +1是在入射能量表中的相邻能量，N 是等概率出射能量的数目，E i,j 是对应于入射能量E i 的第j 个离散出射能量，ρ是选中入射能量表E i 的概率

i

i i E E E

E --=ρ++11

(3.13) 对应入射能量E i 和出射能量E i,j ，有一套等几率分立余弦箱{μi,j,k (k =1...ν)}，以1/ν的概率抽取其中的第k 个，有

k j i k j i ,,1,,)1(+μρ-+ρμ=μ (3.14) b)对于弹性散射，上述过程也适用于数据源于非连续

近似的情况。对于连续近似，数据常被存储为一系列参数D K ，

σel =D K /E ，当E Bk ≤E

E 是入射能量，E Bk 是由晶格参数得出得Bragg 能量。如果E Bk ≤E ≤E Bk +1，P i =D i /D k (i=1...k )表示一积累概率分布，从中抽得i ，认为散射发生在第i 个Bragg 晶格的边缘。散射余弦从关系式μ=1-2E Bi /E 中得到。

5) 俘获

有两种处理俘获的方法：模拟和暗含。

若用模拟处理，粒子以σa /σt 的概率被杀死，σa 和σt 是在入射中子能量下碰撞核的吸收截面和总截面。吸收截面是所有的(n ,x )截面总和，x 是除中子以外的任何粒子。

若用暗含处理，中子权重被减小到原来的(1-σa /σt )倍。如果新的权重低于问题的截断权重(在CUT 卡上给出)，使用轮盘赌减少粒子，且每一个存活的粒子都具有更大的权重。缺省情况下，MCNP 使用暗含俘获处理。

对于模拟俘获，被杀死粒子的能量和权重沉积在碰撞栅元中。对于暗含俘获，σa /σt 倍的入射能量和权重沉积下来。二者宏观效果一致。

也可以在确定中子飞行距离时作暗含俘获处理，这时抽样的飞行距离是到散射点的距离而不是到碰撞点的距离，用下式抽取l ，

l =-ln ξ/∑s (3.16) 在散射点中子权重减小到原来的)exp(l s ∑-倍。这种处理是把指数变换技巧和碰撞暗含俘获相结合。

6) 光子的产生

在中子-光子联合输运问题中，若碰撞核光子产生截面不为零，就会产生光子，光子权重为

WGT p =WGT n σγ/σT (3.17) 其中，WGT p 是光子权重，WGT n 是中子权重，σγ是光子产生截面，σT 是总中子截面。

假设中子从栅元S 进入栅元I ，二栅元的中子重要性分别是IN S 和IN I ，PWT 卡上I 栅元光子权重下限置为WM I ，若WGT P >WM I ?IN S /IN I ，就可能导致光子分裂，产生光子数N P =WGT P ?IN I /(5WM I ?IN S )+1，N P ≤10，光子权重为WGT P /N P ；反之，用轮盘赌杀死部分光子，光子存活几率WGT P ?IN I /(WM I ?IN S )，权重增大为WM I ?IN S /IN I 。

假设入射中子能量下第i 种反应的光子产生截面是σi ，可以使用抽取随机数ξ的方法用

∑∑∑-===σ<σξ≤σ11

1

1

k i n i k

i i

i i

(3.18)

选择反应k ，从而得到光子能量。诱发光子的方向从与能量相关的角分布表中抽取，抽样方法和前述次级中子飞行方向抽样方法类似(见式3.6)。

D. 光子与物质作用

MCNP 模拟光子的能量范围是1keV~100MeV 。有两种处理光子输运的方法，可通过PHYS:P 卡选择。

1) 简单的物理处理过程

这种处理对高Z 核和深穿透问题不适用，整个物理处理过程只有光电效应、电子对效应和康普顿散射，总截面σt =σpe +σpp +σs 。

光电效应被认为是纯吸收过程，由光子权重相应减少的暗含俘获来实现。在每次碰撞中，权重WGT σpe /σt 和能量E ?WGT σpe /σt 在相应的俘获箱中记录下来，而没有被俘获的权重WGT (1-σpe /σt )被迫经历电子对效应或康普顿散射：

以几率σpp /(σt -σpe )选择电子对效应，认为在碰撞点产生正负电子对后，动能WGT (E -1.022MeV)转化为热能沉积下来，正电子立即湮灭，各向同性地发射一个携带2WGT 权重的0.511MeV 光子。

以几率σs /(σt -σpe )抽样与自由电子碰撞的康普顿散射，确定出射光子能量E ’和方向余弦μ，能量WGT (E ’-E )沉积下来。微分截面由Klein-Nishina 公式给出：

K (α,μ)d μ=πr 0(α'/α)2[α'/α+α/α'+μ2-1]d μ (3.19)

其中r 0=2.817938?10-13

cm 为经典电子半径，α、α'为入射和散射光子的折合能量，α=E /mc 2，α'=α/[1+α(1-μ)]。 2) 详细的物理处理过程

这一过程在MCNP 中是缺省的，包括光电效应、电子对效应、相干散射和非相干散射，记录产生的荧光光子，并且考虑电子束缚的影响，从而对Thomson 和Klein-Nishina 微分截面进行修正。总截面σt =σpe +σpp +σscoh +σs 。，其中σpe 是光电截面，σpp 是对产生截面，σscoh 是相干(Thomson)散射截面，σs 是非相干(Compton)散射截面。

a) 平均自由程

对每种元素，MCPLIB 截面库中有一些与对数能量L i (Z )=ln[E i (Z )]相对应的截面数据L j,i (Z )=ln[σj,i (Z )]，其中j =1、2、3和4，分别表示非相干散射、相干散射、光电效应和电子对效应。

能量为E 的光子在特定物质中的平均自由程是

∑∑=σ=λZ

j j Z Z N 4

1

)

()(1

(3.20) Z 对物质中所有元素求和，N (Z )是原子密度，σj (Z )是截面。抽取随机数ξ∈[0,1)，确定到下次碰撞点距离l =-λln ξ。发生碰撞后，再用两个随机数ξ1和ξ2确定元素Z 和过程j ，前者由ξ1/λ与N (Z )σt (Z )的部分求和相比较得到，后者由ξ2σt (Z )与σj (Z )部分求和比较得到。

b) 非相干(Compton)散射

设非相干散射微分截面σ1(Z ,α,μ)d μ=l (Z ,v )K (α,μ)d μ，其中l (Z ,v )是修正公式(3.19)的因子，以表的形式给出，v 也是α和μ的函数。和简单物理处理一样，从修正Klein-Nishina 公式出发抽取出射光子能量E ’和方向余弦μ。

c) 相干(Thomson)散射

相干散射能量不变，只需抽样μ。微分截面可表示为σ2(Z ,α,μ)d μ=C 2(Z ,v )T (μ)d μ，T (μ)d μ=πr 02(1+μ2)d μ，C 2(Z ,v )以表的形式给出。

对于高能光子，相干散射主要是前向散射，μ?1。使用点探测器时，由于探测器一般并不位于散射方向，所以详细处理比简单处理结果的方差更大。

d) 光电效应

与简单处理相比，详细处理的光电效应还包括荧光光子的发射。光子敲出壳层电子后，高层电子向下跃迁后发

射荧光光子，仅考虑有1~2个能量高于1keV 的荧光光子出射，剩余能量沉积在碰撞点。Z <12时光电效应按热吸收过程处理，因为荧光光子能量低于1keV ；31>Z ≥12时可能发射一个荧光光子，Z >31时可能发射两个荧光光子，能量可能是K α1、K α2、K β1，Z >37时还可能是K β2，等几率抽样。荧光光子方向从各向同性分布中抽取。

e) 对产生

与简单处理相同。 f) 其它

MCNP 不考虑韧致辐射。认为能量低于0.001MeV 的截面全是纯吸收截面。

建议不同时使用点探测器和详细的光子物理处理，使用简单物理处理能够明显提高点探测器的计算效率。对能谱计算，当荧光辐射比较重要，或问题几何对能量低于1MeV 的光子有几个平均自由程厚，则推荐使用详细物理处理。两种处理的比较见下表。

§3.2 记数

记数是对计算结果的记录，MCNP 提供一些标准记数

方式，以助记符F1、F2、F4、F5a 、F6等区分，用户可以根据自己的意愿加以选择。

记数都被归一化成对应一个源粒子。

A. 面流量记数

F1型记数记录如下物理量

????μμμ=A t E dA dEdtd t E r J ),,,(1F (3.21)

????μμμ?=A t E dA dEdtd t E r J E ),,,(1F *

(3.22)

此记数是流过界面的粒子数(*F1是能量)，流量与通量的关系是J (r ,E ,t ,μ)=|μ|φ(r ,E ,t )A ，其中r 、E 、t 和μ分别表示粒子通过曲面时的位置、能量、时间和方向余弦，φ是通量函数，A 是面积。用户通过FSn 、En 、Tn 和Cn 卡指定面积、能量、时间和角度的积分区间。

B. 通量记数

F2和F4型记数可表示为 ??φ=t E dEdt t E r ),,(4,2F (3.23)

??φ?=t E dEdt t E r E ),,(4,2F *

(3.24)

用户可在En 和Tn 卡上指定积分区间。

1) 栅元通量的径迹长度估计(F4)

粒子通量的定义是φ(r ,E ,t )=vN (r ,E ,t )，其中v 是粒子速度，N 是粒子密度(等于粒子权重/栅元体积)。对时间积分后通量可表示为

?==φV WT V Wvt dt t E r l //),,( (3.25) 其中W 是粒子权重。由于径迹长度T l 在分子上，这类记数就叫做栅元径迹长度记数。

2) 界面通量(F2)

界面通量可以理解为厚度δ极薄情况下的栅元通量 |)|/()/(|)|/(/lim 2F 0

μ=δμδ==→δA W A W V WT l (3.26)

C. 栅元能量沉积记数

栅元加热和能量沉积记数也与径迹长度有关

??φρρ=t E g a dEdt t E r E H ),,()(/6F (3.27) 其中ρa 是原子密度(atoms/barn-cm)，ρg 是质量密度(g/cm 3)，

H (E )是加热函数，

)()(E H E H avg T σ= (3.28) σT 是中子或光子的总截面，中子和光子H avg (E )的形式不同。

对于中子，

∑γ-+-=i

i i i out i avg E E Q E E E p E E H )]()()[()(,, (3.29)

其中)(E p i 是反应i 的几率，)(,E E i out 是平均中子出射能量，Q i 是反应能，)(,E E i γ平均光子出射能量。

对于光子，

∑=-=3

1

))(()(i out i avg E E E p E H

(3.30)

其中i =1、2和3，分别表示非相干散射、对产生和光电效应，

)(E E out 是平均光子出射能量。

D. 探测器记数

在MCNP 中，某点处的通量可用点探测器或环探测器记数来估计。

1) 点探测器

点探测器用确定论方法计算空间中某点处的通量，记录源粒子和碰撞产生粒子对该点的贡献。设p (μ,?)d Ω是粒子射向、被散射到或生成后飞向(μ,?)方向d Ω立体角的几率，R 是源点或碰撞点到探测器的距离，则p (μ,?)d Ω?e -λ是沿(μ,?)方向d Ω立体角无碰撞地到达探测器的几率，其中?∑=λR

t ds s 0)(是平均自由程的个数，])(exp[0?∑-R

t ds s 是一

束单能粒子通过介质的衰减几率，s 是朝向探测器方向的位

移，∑t (s )是在s 点的总截面。由于d Ω=dA /R 2，p (μ,?)= p (μ)/2π，则在探测器处单位面积dA 上的上述来源的有权重粒子数(即通量分布函数)等于

)2/()(),,,(2R e Wp t E r πμ=μφλ-

(3.31) 其中W 是粒子权重。

点探测器也被称为“下次事件估计器”，因为估计点通量记数时，下一个事件是一个无碰撞直接到达点探测器的事件。这种过程可以理解为伪粒子的输运。

注意(3.31)式中分母上的R 2造成R →0时φ→∞。解决此问题的方法是设置一个包围探测点的邻域球，半径R 0可以以cm 或平均自由程为单位给出。如果以cm 为单位设置R 0，在R 0内的通量用球内平均通量替代，

?∑π∑--μ=φ=<φt

t R R Wp V dV R R 3

0002))

exp(1)((3/)( (3.32) 当∑t =0，探测器在真空中，球内无碰撞，则

3

02)(3)0,(R R Wp R R t πμ==∑<φ (3.33) 如果以平均自由程为单位设置R 0，则λ0=R 0=∑t R ，有

3

2

002))exp(1)((3)(πλ∑λ--μ=

<φt Wp R R (3.34) 对于探测器在真空中或稀薄介质(如空气)中，可将R 0

设置为0。典型情况下，R 0设置为平均自由程或它的几分之一较合适。如果以cm 为单位设置R 0，应考虑球内平均粒子能量的平均自由程。一般邻域球内不应出现多种物质。

2) 环探测器

环探测器和点探测器类似，但环探测器在一个环上记数。F5X 、F5Y 、F5Z 探测器分别表示探测器环相对于x 、y 、z 轴旋转对称。在环状对称系统中，绕对称轴环上各点的通量相等，用户自然可以在环上选择一点记录抽样结果来计算环上平均通量。然而环探测器的优点表现在，MCNP 通过位置偏倚技巧选择环上探测点，使靠近环的碰撞得到优化处理，在探测器总记数不变的前提下，经常抽样贡献较大的事件，减小相对误差。所以环探测器一般比点探测器有更高的效率。

3) 伪粒子和探测器可靠性

探测器是使用了确定论捷径的混合蒙特卡罗方法，从源点或碰撞点到探测位置的输运被确定估计所替代。可以理解为存在着伪粒子，伪粒子不再参与碰撞，也不影响产生它们的粒子的权重和其它效应，只用于估计对探测器的贡献。伪粒子估计通量时需要一个衰减因子，由对源点或碰撞点到探测器之间径迹的求和得到。

探测器记数虽然效率较高，但也经常误导用户。例如，考虑一个各向同性单能点源和探测器附近有一个散射区域的情况(见图3.1)，如果粒子未射向散射区，每次到达探测器的伪粒子提供相同的贡献，记数统计涨落为零。第一个粒子射向散射区后，探测器记数及其误差突然跳变，随着后来的源伪粒子到达探测器，统计误差又逐渐减小，第二个射向散射区的事件发生后，探测器记数的误差又突然变大。这种周期性的误差变大变小是探测器的一个特点。要较好地解决这一问题，需对散射区充分抽样。建议用户仔细研究探测器诊断打印表。探测器记数相对误差大于5%的结果是值得怀疑的。

4) 探测器记数的减小方差技巧

一些适用于模拟粒子的减小方差技巧不适用于伪粒子，例如，虽然伪粒子在零重要性区域被杀死，但不根据栅元重要性或权重窗口分裂。MCNP 提供两个专门适用于伪粒子的轮盘赌技巧。

PD 卡可用于指定每个栅元的伪粒子产生几率。对栅元i 的输入项p i 应介于0~1之间，伪粒子产生几率为p i ，权重为1/p i 。缺省时p i =1，即每个源或碰撞事件产生一个伪粒子。设置p i =0可以使栅元i 中不产生伪粒子，从而去掉该栅元对探测器的直接贡献。

DD 卡同时控制探测器诊断打印和向探测器输运伪粒子的轮盘赌。DD 卡上的参数项k 用于指定启动轮盘赌的权重下限w c : 若k ≤0，w c =-k ；若k >0，从第201个伪粒子开始，w c 等于以前伪粒子的平均权重的k 倍。缺省k =0.1。

5) 直接贡献和总贡献

MCNP 的输出文件中将自动打印对探测器的总贡献和直接贡献。直接贡献是来源于源事件的伪粒子记数，总贡献则还包括碰撞事件伪粒子的贡献。在N P 联合输运方式

下，中子碰撞产生光子是光子源事件。

E. 记数精度

1) 平均值，方差和标准偏差

记数结果是抽样许多历史过程贡献的平均值，假定P (x )是选择一个随机过程的概率密度函数，x 是该过程产生的要估计的值，x i 是从P (x )中抽取的第i 个历史的贡献，共抽样N 个粒子，则x 的近似期望值为

∑==N i i x N

x 1

1

(3.35) x 就是MCNP 输出文件中给出的记数平均值。

x 值分布的方差σ2是其离散性的量度，此方差的估计值S 2等于

21

2

1

2

2)(1~

1)(x x N N x x S N i i N

i i ---=

∑∑==

(3.36) σ称为标准偏差，是方差的平方根，其估计值S 是实际抽样值x i 标准偏差的总体估计。

x 的方差由下式估计，

N S S x /22=

(3.37)

由于x S 与1/N 成正比，减小一个数量级x S 要计算100倍原来的粒子数目。也可以固定N 通过S 的减小来减小x S ，即通过减小方差的技巧来实现。

在MCNP 的输出文件中，对应每个记数平均值x ，还给出其相对误差R ，

2

1

21122

/122111??????

????????-??? ??=???????????? ??-==∑∑==N x x x x N x S R N i i N i i x (3.38) MCNP 建议R <0.05的结果才是可信的。

2) FOM

FOM(figure of merit)反应了MCNP 的计算效率，定义为FOM=1/R 2T ，R 是记数平均值的相对误差估计，T 是计算时间。因R 2∝1/N ，T ∝N ，故理想情况下FOM 应为常数。

在MCNP 输出文件的末尾给出对应一定运行粒子数的FOM 。FOM 是一个很重要的统计参数，它还是运行结果可靠性的指示器。对于一个好的运行结果，运行初期FOM 有一定统计涨落，运行后期FOM 趋于某一固定值。只要FOM 还随运行粒子数的增加不断跳变，结果就不可靠。

影响效率的因素有三种：历史贡献的效率、非零记录的离散度和每个历史所用的计算时间。其中前两个影响R ，后一个影响T 。R 可以分成两部分，即非零记录的效率部分

2eff R 和本身具有的离散度2int R ，R 2=2eff R +2int R 。定义f 为产生

非零记录的那部分历史， 有

)/()1(2

fN f R eff -=

(3.39)

fN x i x R N

i N

i i

112

12

int -

∑∑=

?

?

? ??==

(3.40)

所以可以通过一种或多种减小方差的技巧来抽样那些产生

非零记录的粒子，增加f 值，减小2

eff R 和R 。同时，粒子权

重也作相应的调整，使结果保持不变。而较小的f 值则需模拟大量的粒子来补偿，以取得精确的计算结果。

3) 精度和准确性

图3.1 探测器记数

精度和准确性是两个不同的概念。精度反应结果绕平均值涨落的情况，而准确性则反应平均值与期望值的偏差。模拟结果的统计误差小并不一定意味着结果准确。

影响MCNP 计算准确性的因素有两方面：程序本身和用户。程序本身的不精确是由于物理处理方法不完善、截面数据存在误差等原因造成的，有待于在今后新版本中不断完善。后一个原因是把实际问题模型化时带进偏差，主要包括源的能量及角分布指定得有问题，或不恰当的几何描述和材料物质描述等。用户也可能选用不合适的记数方式，或，也是经常地，滥用减小方差的技巧。

影响MCNP 记数结果精度的因素有三方面：记数类型、减小方差技巧和计算源粒子数。

§3.3 减小方差技巧

A. 统筹考虑

减小方差的技巧用于缩短记数达到足够精度的计算时间。影响MCNP 计算时效的主要因素是记数类型和粒子随机行走的抽样方法，二者都影响着每个粒子历程对记数结果平均贡献的大小。MCNP 记数结果可以表示为

????>=<),,(),,(t v r T t v r N dt v d r d T

(3.41) 其中),,(t v r N 是粒子密度函数，),,(t v r T

是记数函数。记

数函数),,(t v r T

是δ函数，粒子达到记数区域时T =1，否则T =0。如果只简单地依物理规律抽样事件，大量粒子可能根本就不到达记数区域(包括位置、能量、时间和方向等)，使得计算效率很低。实际上，可以多抽样一些T =1的事件，少抽样T =0的事件，为使保持不变，相应地调整粒子权重，这就是减小方差的技巧。

对用户而言，合理选择和适当应用MCNP 减小方差的技巧并非易事，因为每种减小方差方法有各自的优缺点和适用性。然而，还是有一些供用户遵循的一般性的原则。首先，也是最重要的一点，用户应该非常熟悉要计算的物理问题和可能使用的各种减小方差方法。如果简单抽样不能有效地计算，用户应该分析问题出在哪里，以便相应地解决。其次，选择了减小方差技巧后，考虑如何赋以适当的参数，这对用户来说是更难的事。最初设置的参数可能来源于猜测或计算同类问题的经验，一般应持较保守的态度，即使用较小的偏倚而不是较大的偏倚。然后，进行一个短计算，仔细研究输出文件提供的信息，包括：使用的技巧是否改进了对记数结果的抽样；使用的几种技巧之间是否相互干扰；FOM 是否变大并且随计算粒子数变化正常；有没有其它看起来明显不对的地方等等。经过反复调整和短计算后，用户还可以试着使用权窗产生器让程序自动调整参数。

无论多么熟练的用户都应当对使用减小方差的技巧持保守和慎重的态度，过分偏倚往往是用户的通病。应当牢记减小方差技巧就是减少了对某些情况的抽样，如果这种减少达到了某种不恰当的程度，会降低结果的可靠性。用户应当保证对问题的所有区域都有与其重要程度相当的足够的抽样。

下面逐一介绍MCNP 减小方差的技巧。

B. 能量截断

用户设置一截断能量，低于此能量的粒子将被杀死。能量截断截掉了粒子径迹，所以能节省时间。它的优点是简单，并不直接和方差相关，但存在以下几点问题：

1) 低能粒子也可能产生高能粒子。尽管用户可能对记录低能粒子不感兴趣，但低能粒子有时也是重要的。 2) 对整个问题区域使用同一截断能量。一些区域中低能粒子可能不重要，另一些区域也许比较重要。

3) 也许会因能量截断出现结果系统偏小。如果计算大量粒子，记数统计涨落趋于零，容易误导用户以为计算结果以足够准确。因此，能量截断不如能量窗口好。

C. 时间截断

用户设置一时限，粒子生存时间超限后，终止其径迹。时间截断应只用于与时相关问题。

D. 几何分裂和轮盘赌

几何分裂和俄罗斯轮盘赌是最早的、使用最广泛的技巧。若使用得当，可节省大量机时。它的原理是，当粒子向重要方向迁移时，粒子数增加，反之粒子数减少。伴随着粒子数的增减，权重相应调整，因而结果是无偏的。

在输入文件中，用户通过IMP 卡给每个栅元设置一个重要性。如当前粒子具有权重WGT ，从重要性为I 的栅元要进入重要性为I’的栅元。当I>I’，粒子有I’/I 的几率存活，存活粒子权重变为WGT ?I/I’。当I

栅元重要性可以设置成任何非负实数，但相邻栅元的重要性不能相差太多(除重要性为0外，一般不大于4倍)，否则粒子数过分分裂也会降低计算实效。建议在穿透方向上栅元厚度设置为粒子数减小一半的距离，栅元重要性设为相差2倍，总体上使穿透方向上粒子数基本不变。

其它需要注意的事项是：粒子在真空栅元中不分裂也不赌；粒子进入零重要性栅元后就被杀死。

E. 能量分裂和轮盘赌

能量分裂和轮盘赌经常同时使用，但用户可以根据需要只使用其中一种。能量分裂/轮盘赌与空间位置无关，如果问题具有位置-能量依赖性，使用与位置-能量相关的权窗是更好的选择。

1) 分裂

有时一些能量区域比其它能区更重要，可以在重要能区使粒子数分裂，伴随相应的权重调整。粒子可以在能量降低或升高时分裂，可以设置5级能量分裂界限。

能量分裂可能导致方差减小或增大。目前MCNP 权重截断不考虑粒子数是否依能量分裂过。

2) 轮盘赌

多数情况下，尤其是中子进入热能区时，随能量降低粒子数增多，这些低能粒子要经历大量耗时的碰撞。当低能粒子在问题中比较重要时，不能使用能量截断技巧，但可以对粒子数用轮盘赌技巧加以控制。

采用轮盘赌技巧后，当粒子能量低于设置能限时，以设定的几率存活。存活粒子的权重相应被调大，以保证整体权重无偏。

F. 暗含俘获和权重截断

1) 暗含俘获

暗含俘获、存活偏倚和权重减小的吸收是同义词。选

中被碰撞原子核i 后，如果σti 是总微观截面，σai 是微观吸收截面，并不以σai /σti 的几率吸收粒子，粒子总是存活，但权重减小为WGT (1-σai /σti )。

暗含俘获有两个优点：1) 粒子不在到达记数区域之前被杀死；2) 不对存活粒子权重抽样，节省时间。暗含俘获也存在两个缺点：1) 引入了粒子权重涨落；2) 跟踪每个粒子的时间增加，当粒子权重较小时可能浪费机时。

2) 权重截断

当粒子权重低于用户设置的权重下限时，启动轮盘赌，粒子或者被杀死，或者权重增大到用户设置值：在栅元j 中的权重截断不仅与用户在CUT 卡上设置的WC1和WC2参数有关，而且和该栅元重要性有关。当粒子权重低于R j ?WC2后权重截断就起作用，其中R j 是源栅元和当前栅元重要性之比。粒子存活几率是WGT /(R j ?WC1)，存活粒子权重调整为R j ?WC1。和能量截断、时间截断不一样，权重截断使用轮盘赌，因而结果是无偏的。

权重截断是为配合几何分裂/轮盘赌和暗含俘获而设计的。为了理解权重截断的必要性，可以从反面考虑如果没有权重截断会发生什么情况，设想粒子在很大的介质中输运，不使用能量截断或时间截断，粒子经过多次碰撞后损失了大量权重，跟踪大量小贡献粒子自然会浪费机时。权重截断设计成与栅元重要性之比有关，是为了兼顾在重要性大的栅元中粒子分裂后权重减小。MCNP 的许多技巧导致权重改变，由于权重截断最初是为几何分裂和暗含俘获设计的，同时使用其它技巧时要格外小心。

设置好权重截断的界限不是件易事，需要思考和试验。用户需要猜测多大的权重是值得保留的，然后在它附近做调整。

G. 强迫碰撞

强迫碰撞是一种在指定栅元增加碰撞抽样的技巧。由于探测器伪粒子和DXTRAN 粒子只在碰撞或源抽样过程中产生，因此增加碰撞次数可使这两种粒子增多。也可以应用强迫碰撞技巧在较薄的栅元中充分抽样。

强迫碰撞把进入某栅元的一个粒子分裂成碰撞部分和不碰撞部分。不碰部分还是一个粒子，直接离开本栅元，权重变为W 0?exp(-DLS ?QPL )，其中W 0是分裂前的权重，DLS 是沿粒子运动方向的栅元长度，QPL 是栅元内物质的微观截面。碰撞部分具有剩余权重W 0?(1-exp(-DLS ?QPL ))，并使用轮盘赌控制这种粒子的数量，存活几率f 由用户设置，存活粒子权重调整为W 0?(1-exp(-DLS ?QPL ))/f 。当强迫碰撞栅元相邻或碰撞次数过多时，控制碰撞部分的粒子数是有必要的。粒子在飞行方向上碰撞的位置PMF 由抽取随机数ξ确定

)]]exp(1[1ln[1

QPL DLS QPL

PMF ?--ξ--= (3.42)

强迫碰撞常产生很多小权重粒子，需注意权重截断和权窗对这种粒子的影响：

1) 如果在强迫碰撞卡上设置负输入项。强迫碰撞后次级粒子正常抽样。关闭权重截断。权窗在探测器粒子和DXTRAN 粒子产生前不起作用。为使伪粒子和DXTRAN 粒子能达到记数区，如果使用权窗，其下限应设置为碰撞粒子权重或零。

2) 如果在强迫碰撞卡上设置正输入项。碰撞后产生探测器粒子和DXTRAN 粒子，权重截断和权窗都起作用。存活的碰撞粒子可以继续参与强迫碰撞。如果使用权窗，它将限制进入栅元粒子的权重。

H. DXTRAN

DXTRAN 是确定性输运(deterministic transports)的缩写，它是类似探测器记数的混合蒙特卡罗方法，一般用于不能很好抽样的小区域，因为粒子射向小区域的几率很小。这时，用户可设置一包围感兴趣区域的DXTRAN 球。在球外每次碰撞中，产生一特殊的DXTRAN 粒子，按确定几率散射向并达到DXTRAN 球表面，不再参与碰撞，权重按确定性输运计算。碰撞产生的正常粒子按蒙特卡罗方法继续处理，权重也不做特殊调整，只是正常粒子试图进入DXTRAN 球时被杀死，从而可保证整体结果无偏。在球的内部，DXTRAN 粒子按普通粒子处理。因此，DXTRAN 是分裂、轮盘赌和从非模拟几率密度函数抽样的综合。

DXTRAN 的主要缺点是，额外地跟踪DXTRAN 粒子花费机时。有时DXTRAN 粒子在球内多次碰撞，导致权重很低。DXTRAN 球内粒子不被常规权重截断和权重窗口所控制，但有两种控制低权DXTRAN 粒子数量的方法：

1) DXC 卡与栅元有关，用户可以指定栅元i 中碰撞或源抽样产生DXTRAN 粒子的几率p i ，产生的DXTRAN 粒子的权重乘以p i -1，由此可以产生权重较大的DXTRAN 粒子。用户可通过设置各栅元的p i 使进入DXTRAN 球粒子的权重大致相等，但选择合适的p i 较困难，一般在权重涨落主要来源于按物质吸收系数指数衰减时应用此方法。

2) 权重涨落的另一个、也经常是主要的来源是粒子向DXTRAN 球所在立体角度散射的几率不同，因此采用DD 卡控制DXTRAN 粒子数更方便。前面在§3.2 D 4)中介绍的控制探测器伪粒子的DD 卡也同样适用于DXTRAN 粒子。DD 卡上的参数项k 用于指定启动轮盘赌的权重下限w c ：若k ≤0，w c =-k ；若k >0，从第201个DXTRAN 粒子开始，w c 等于以前DXTRAN 粒子的平均权重的k 倍(缺省k =0.1)。

使用DXTRAN 技巧需注意的一些事项有：1) 在光密介质中不宜使用DXTRAN ，不要指望用DXTRAN 解决深穿透问题。2) 当一个问题中设置多个DXTRAN 球时，往往产生相当多的权重极低的DXTRAN 粒子。3) 因为常规粒子不能进入DXTRAN 球，不要在DXC 卡上对DXTRAN 粒子设置0存活几率。4) 最好使产生DXTRAN 粒子和杀死常规粒子的权重基本平衡。5) DXTRAN 技巧和探测器记数都是向一小区域输运非常规粒子，因此使用探测器记数需注意的事项也使用于DXTRAN ，见§3.2 D 3)。

I. 源变量偏倚

MCNP 允许对任何或全部通用源(SDEF 卡上)的参数项做偏倚，这就是说，允许在每个源变量的重要区域多抽样一些粒子，相应地减小这种粒子的权重。例如，在一个屏蔽问题中，对于一个有确定能谱的各向同性点源，可以在记数方向上多抽样一些高能源粒子，并相应地调整每个源粒子的权重。

1) 通过指定外部抽样频率来偏倚

可以在SB 输入卡上指定不同变量区间的源粒子产生几率，从而对某一源变量进行源粒子偏倚。如果这种硬性指定的几率与应在该变量区间产生粒子的实际几率不相符，相应源粒子的权重调整为实际几率(在SP 卡上定义)除以硬性指定的几率(在SB 卡上指定)，从而使该区间中粒子总权重保持不变。

2) 通过指定内部函数来偏倚

象指定源变量分布几率一样，也可以通过指定内部函数来做源粒子抽样偏倚。这时SB 卡的输入类似于SP 卡的输入(详细指定方法将在第四章讲述)，第一个参数项是用负数指定的函数号，例如：a) 方向偏倚。对源粒子出发方向的

抽样可以按指数函数偏倚。b) 沿圆柱伸展方向偏倚。对源粒子出发位置的抽样，即出发位置沿圆柱体对称轴方向离开中心的距离，可以按内部函数偏倚。c) 沿球半径方向偏倚。对源粒子出发位置的抽样，即出发位置距球心的距离，可以按内部函数偏倚。等等。

J. 权重窗口

权重窗口是与位置和能量相关的分裂/轮盘赌技巧。用户可以指定一些能量区间，然后对每个栅元的每个能量区间指定权重上下限，其中上限是下限的固定倍数。如果在特定位置-能量区域中，粒子权重低于设置下限，则启动轮盘赌；若高于设置上限，则分裂。可以在粒子碰撞或穿越界面时检查其权重。

虽然权重窗口可以单独使用，且实践证明它是一种很有效的减小方差技巧，但设计它的初衷是为了弥补其它减小方差技巧造成权重变化过大的缺陷。尤其是，在某些情况下，由轮盘赌导致粒子权重过大，这种生存几率很小且权重相对较大的粒子对记数稳定性起负面作用；相反，分裂过程产生一些很小权重的粒子，对记数贡献也很小，跟踪它们会造成机时浪费。因此保证抽样过程中粒子权重的基本平衡是十分重要的。虽然不能强求每个粒子历史对记数的贡献完全相等，但通过设置与重要性成反比的权重窗口，可以使每个径迹的贡献大致相当。与重要性成反比的系数和窗口的宽度是需要用户调节的主要参数。实践中发现，当窗口上限是下限的5倍时比较好，并且计算效率对窗口的宽度并不十分敏感。与重要性成反比的系数一般选择成使一些参考栅元中粒子数比较合理，最方便的选择方法是使源栅元中粒子权重处于窗口中。

对某些位置的某些能量区域，可以通过把下限设置为零来关掉权窗。需要注意的是，如果碰撞中权窗关闭，权重截断工作，有时可能杀死过多的粒子。

1) 权窗和几何分裂的区别

虽然二者都使用分裂和轮盘赌技巧，但存在一些显著区别：a) 权窗与位置和能量有关，几何分裂只与位置有关；b) 权窗依粒子权重决定如何处理，几何分裂不管粒子携带的权重有多大；c) 权窗技巧的判据是权重上下限，几何分裂的判据是重要性之比；d) 权窗在栅元界面和碰撞点起作用，几何分裂只在界面起作用；e) 权窗可控制其它技巧造成的权重涨落，几何分裂保持这种涨落。

2) 权窗产生器

用户使用权窗技巧时需要设置随位置和能量变化的权窗上下限，如果问题中栅元较多，或想在较多能量区间中操纵粒子权重，就得输入很多参数。用户可能需要凭借猜测、直觉、经验和不断摸索才能将这一技巧有效地利用起来。由于权窗上限都是下限的固定倍数，下限又设置成反比于区域重要性，可以把下限随栅元-能量区间的变化看成二维离散函数，称为权窗重要性函数。MCNP 提供的权窗产生器能自动计算权窗重要性函数，从而简化了用户的工作。

如果想让粒子径迹对记数的贡献相等，应使相空间中某点的重要性正比于进入该区间的粒子对记数-的贡献。模拟一定数量的粒子后，可按下式计算重要性，

进入该栅元的总权重

的贡献

进入该栅元粒子对记数重要性=

(3.43)

然后按重要性倒数的分布指定权窗，并在输出文件中提供输入文件需要的权窗卡。

虽然实践证明权窗产生器很有用，但它决不是解决所

有重要抽样问题的万金油。用户有时要站在自己的角度思考问题，权窗产生器计算的权窗就不太适用；多数情况下，用户还要判断计算出的权窗是否合理。用户需牢记，产生的参数是统计结果，有可能带很大误差，除非在设置的相空间中已充分抽样。为此，用户应在使用权窗产生器时预先给相空间粗略设置比较合理的重要性，经过计算-设置-再计算的多次过程才能得到较好的权窗；要对计算结果加以判断，对那些值得怀疑的参数重新设置。另外，若未充分划分相空间，也会导致权窗产生器失效。

K. 指数变换

指数变换从非模拟概率密度函数中抽样碰撞距离，在特定方向上人为地减小微观截面并且在相反方向上增大微观截面，截面按下式修改

)1(*μ-∑=∑p t t

(3.44)

其中*t ∑是变换后的截面，∑t 是实际总截面，p 是控制偏倚程度的指数变换参数，μ是和参考方向夹角的余弦。|p |<1，可取常数，在暗含俘获中取p =∑a /∑t 。

为保证结果无偏，粒子碰撞或穿越界面时需修改其权重，乘以因子w c 。碰撞时，

)1/()exp(μ-μ∑-=p s p w t c (3.45) 穿越界面时，

)exp(s p w t c μ∑-=

(3.46)

其中s 是碰撞距离或到界面距离。

对于高吸收性介质中的深穿透问题，方差一般随p 从0变到某一p ’逐渐减小，再随p 从p ’变到1逐渐增大。习惯把p

p ’为过偏倚。用户需要凭经验和不断试验寻找合适的p 。以下是一些应用经验：

1) 若不把指数变换和权窗同时使用，MCNP 将提出警告。大量实践表明，单独使用指数变换经常得到不可靠的记数结果和误差估计，因此需要合适的权窗加以控制。

2) 指数变换适用于高吸收性介质；不适用于高散射性介质。

3) 对于中子在水泥或土壤中输运，变换参数p 取0.7比较合适。对于光子在高Z 介质中输运，更高的偏倚，如p 取0.9，也许是合适的。

L. 相关抽样

非相关抽样往往使问题微小变化的结果差异被统计误差所掩盖，相关抽样技巧使得对这种小量变化的估计成为可能。MCNP 总是从固定的序列中挑选伪随机数。假设第i -1个粒子使用的第一个随机数是序列中的第n 号数，不管第i -1个粒子使用了多少个随机数，第i 个粒子所使用的第一个随机数都是序列中的第n +4297号数，然后该粒子的后续随机数又在序列中此后固定步长的位置挑选。

如果计算微绕问题，对原始和微绕后情况而言，相对应的每一个新历史的起始随机数都相同，而后模拟粒子径迹的随机数序列也相同，除非微绕使随机数序列分岔。

M. 点探测器

与其它减小方差技巧不同，点探测器不仅控制粒子随机行走，它本身就是一种记数类型，并且和在半径极小的DXTRAN 球面法线上记录面通量是等价的。在§3.2 D 中已对点探测器记数详细介绍，这里不在赘述。

第四章 数据卡

数据卡在INP文件的数据块中使用。常用的数据卡分以下几类：

§4.1 问题类型卡

格式: MODE X

X: N, 中子输运(卡可缺省)

N P, 中子-光子联合输运

P, 光子输运

此卡用于指定计算问题的类型，MODE N问题的源粒子是中子，MODE P的源粒子是光子，MODE N P的源粒子可以是中子也可以是光子，省略MODE卡表示MODE N。MCNP中，中子碰撞产生光子，光子碰撞不产生中子。

§4.2 栅元参数和曲面参数卡

以上卡中只有IMP卡是需要有的(使用权窗时除外)，多数卡可用于降低方差。栅元参数卡或曲面参数卡上的数据项必须与栅元描述卡或曲面描述卡在INP文件中的顺序一致、个数相等。

A. IMP 栅元重要性卡(Cell Importance Card)

格式: IMP:n X1 X2 ... X l

n: N, 中子; P, 光子

X i: 栅元i的重要性

l: 栅元总数

栅元重要性用于几何分裂和轮盘赌技巧，能使更多的粒子轨迹朝向几何位置上较为重要的区域运动。

粒子在0重要性栅元被杀死，对探测器和DXTRAN球的贡献也终止。

在MODE N P问题中，若省略IMP:P卡，则全部栅元光子重要性取1，除非中子重要性为0的栅元，此栅元光子重要性也取0。

B. VOL 栅元体积卡(Cell V olume Card)

C. AREA 曲面面积卡(Surface Area Card)格式: VOL X1 X2 ... X l

X i: 栅元i的体积

l: 栅元总数

格式: AREA X1 X2 ... X l

X i: 曲面i的面积

l: 栅元总数

通常MCNP可以自动计算栅元体积和曲面面积，并在OUTP文件中打印。为了某种记数的目的，用户可以用VOL 和AREA卡修改这些参数。有时MCNP不能计算全部的栅元或曲面参数，如果记数中用到它们，用户又没手工输入，则会出现致命错误。以上二卡可被分段除数卡SDn代替。

D. PWT 光子产生权重卡(Photon Weight Card)

格式: PWT W1 W2 ... W l

W i: 中子在栅元i内碰撞产生光子的权重的相对阈值

l: 栅元总数

仅在MODE N P问题中PWT卡有效，它的作用是控制中子碰撞时诱发光子的权重和数目。通过诱发光子的权重与PWT卡上指定的有关权重阈值之比来控制光子的产生及其权重，权重阈值与栅元中子重要性有关，若W i<0还与源中子权重有关。设中子从栅元s出发，在栅元i与原子核碰撞，两个栅元的中子重要性分别是I s和I i，该中子在栅元s内时的权重为W s，抽样出的诱发光子权重是W p，W i是PWT卡上栅元i的数据项：

1)若W i≥0：W p≥W i?I i/I s时，光子权重不变；W p

时，取W p=W i?I i/I s，用轮盘赌决定光子是否存活。

若W i=0，诱发光子权重不变，都存活。

2)若W i<0：用-W i?W s?I i/I s替代W i?I i/I s作为是否启动

轮盘赌和存活光子权重的判据。若W i=-1E6：关闭

栅元i的光子产生。

如果希望较少的诱发光子(携带较大权重)，应增大W i 的绝对值。如果缺省此卡，则所有栅元的W i均被置为-1。

E. EXT 指数变换卡(Exponential Transform Card)

F. VECT 矢量输入卡(Vector Input Card)

格式: EXT:n A1 A2 ... A l

n: N, 中子; P, 光子

A i: 对应栅元i的数据项, 格式为QVm, 其中Q是拉

伸长度, Vm是坐标轴或VECT卡定义的方向, 表

示拉伸方向

l: 栅元总数

格式: VECT Vm X m Y m Z m ... Vn X n Y n Z n

m,n: 标记矢量Vm, Vn的符号, 任意

X m Y m Z m: 定义Vm方向的末点坐标

指数变换方法是通过调整总截面来拉伸特定方向上碰撞之间的径迹长度，权重则相应调低。总截面调整为∑t* =∑t(1-pμ)，其中∑t*是调整后总截面，∑t是真实总截面，p是拉伸参数，μ是粒子方向和拉伸方向夹角余弦。

EXT卡数据项中Q可以是绝对值小于1的常数或符号S，前一种情况下该常数绝对值就是拉伸参数p，后一种情况下p=∑a/∑t，∑a是俘获截面。A i=0表示该栅元中不拉伸。

可以用VECT卡指定的Vm在EXT卡数据项中定义拉伸方向，即碰撞点到(X m Y m Z m)点的方向。或者用X、Y或Z表示拉伸方向与坐标轴平行。这两种情况下若Q为负值都表示拉伸方向与Vm相反。也可以在A i项中不设置V m，表示各项同性地拉伸，一般不推荐使用这个方法，除非想在飞行方向上做暗含俘获处理，这时Q填写符号S。

在强迫碰撞栅元中或在没有权窗控制的栅元中不应使用指数变换。仅当粒子通量按指数分布时，如强吸收问

题，使用指数变换才较为合适。

G. FCL 强迫碰撞卡(Forced Collision Card)

格式: FCL:n x1 x2 ... x l

n: N, 中子; P, 光子

x i: 第i栅元强迫碰撞控制参数, -1

i

<1

l: 栅元总数

FCL卡控制栅元内中子或光子的强迫碰撞，主要用于光疏介质中的粒子输运问题，对向探测器或DXTRAN球产生贡献也非常有用。粒子进入强迫碰撞栅元时在界面上不做权窗处理。

如果x i≠0，进入栅元i的全部粒子被分成碰撞部分和不碰撞部分，并配以相应的权重。如果|x i|≠1，对碰撞部分用存活几率|x i|做轮盘赌以防止径迹过多。如果一些强迫碰撞栅元彼此邻近，则建议选|x i|<1。

若x i<0，强迫碰撞仅用于进入栅元i的粒子。强迫碰撞之后，忽略权重截断并以通常的方式处理后续碰撞。不忽略权窗，并在对探测器或DXTRAN球作贡献之后使用权窗。

若x i>0，对进入栅元i的粒子和权截断或权窗处理后存活的粒子作强迫碰撞处理。粒子将继续被分成不碰撞和碰撞(概率|x i|)两部分，直到粒子被权截断或权窗杀死为止。H-J. 权窗卡

除重要性(IMP)卡和能量分裂(ESPLT)卡外，权窗是对指定的位置和能量提供分裂和轮盘赌依据的另一种方法。权窗的优点是：1) 权窗同时提供能量和位置的重要性函数，2) 权窗控制粒子的权重，3) 与其它减小方差的技巧不矛盾，如指数变换(EXT卡)，4) 粒子跨越界面或碰撞时都可以起作用，5) 能控制分裂或轮盘赌的程度，6) 可在指定的位置或能量区域关闭权窗，7) 可用权窗产生器自动生成权窗。权窗的缺点是：1) 权窗不如重要性或能量分裂直接明了，2) 当改变源的权重时，权窗必须重新归一。

H. WWE 权窗能量卡(Weight Window Energy Card)

I. WWN 权窗边界卡(Weight Window Bound Card)

J. WWP 权窗参数卡(Weight Window Parameter Card)

格式: WWE:n E1 E2 ... E k

n: N, 中子; P, 光子

E i: 能量间隔上限, E

i-1

k: 能量间隔数, 要求k<100

格式: WWNi:n W i1 W i2 ... W il

n: N, 中子; P, 光子

W ij: 在栅元j、能量区间为E

i-1

若没有WWE卡，i=1

l: 栅元总数

格式: WWP:n WUP WSURV MXSPL WHERE SWITCH n: N, 中子; P, 光子

WUP: 当权超过WUP与权窗下限之积时，粒子将被分裂，要求WUP≥2

WSURV: 对轮盘赌存活的粒子，它的权变为WSURV与权窗下限之积，要求1

MXSPL: 粒子分裂时, 分支数不超过MXSPL,

做轮盘赌时, 存活概率不小于1/MXSPL,

要求MXSPL>1

WHERE: 此参数用于指定在什么地方检查粒子权重,

=-1, 仅在碰撞点处检查,

=0, 在碰撞点和界面上检查,

=1, 仅在界面上检查

SWITCH: 若大于0, 权窗下限将被指定为SWITCH除以粒

子所在栅元重要性

WWE卡定义在WWN卡上指定的权窗边界的能量间隔。最低能量是0，在WWE卡上不填写。如果没有WWE卡而又使用了权窗，则对应于要计算的问题只建立一个能量区间。

WWN卡指定与位置及能量相关的权窗下限，它必须与WWP卡一起使用，如果权窗与能量相关，还应有WWE 卡。如果使用WWN:n卡，可以不用IMP:n卡，这里n是粒子类型。如果W ij<0，则进入栅元j的粒子全被杀死，这时必须保证该栅元所有能区的W ij都小于零。如果W ij>0，进入栅元j或在栅元j碰撞的粒子根据WWP卡的选择做分裂或轮盘赌。如果W ij=0，则关掉栅元j能区i的权窗处理，改用权截断处理。

WWP卡上的参数用于控制WWN卡上指定的权窗下限的应用。缺省参数值为WUP=5，WSURV=0.6?WUP, MXSPL=5，WHERE=0，SWITCH=0。如果不使用WWN卡，也可以使用IMP卡与WWP卡设置权窗下限。

K-L. 权窗产生器

权窗产生器估算由用户指定的位置-能量相空间中区域的重要性，然后按与该重要性成反比关系计算权窗下限。权窗产生器应该用于记数与角度没多大关系的问题。

权窗产生器估算相空间的平均重要性，如果区域太大，其内部重要性变化也较大，平均重要性不能很好地代表这个区域的特征。不适当的几何说明会导致栅元之间重要性相差很大，恰好权窗产生器能够提供重要性方面的信息。

K. WWG 权窗产生器卡(Weight Window Generation Card)

L. WWGE 权窗产生器能量卡(Weight Window Generation Energy Card)

格式: WWG I t I c W g X r Y r Z r J

I t: 问题的记数号(Fn卡上的n)

由TFn卡定义的对权窗产生器合适的特定记数

箱

I c: 参考栅元，一般为源栅元

W g: 指定I

c

栅元的权窗下限。当W g=0时，这个值将是源权重的一半。

X r Y r Z r: 粒子朝向点或离开点的坐标

J: 调试打印标记。=1打印，=0不打印

格式: WWGE:n E1 E2 ... E j

n: N, 中子; P, 光子

E i: 能量间隔上限, E

i-1

j: 能量间隔数, j<100

WWG卡用于产生I t记数的最佳重要性函数，在输出文件OUTP中做打印和评价，所产生的权窗用于以后运行中的WWP和WWNi卡。对许多问题，这种产生的重要性分布(权重下限)比有经验的用户在IMP卡上猜想的要好。

为了产生与能量有关的权窗，需使用WWGE卡。如果没有WWGE卡，又要使用权窗产生器，则根据问题的能量范围只产生一个能量区间。如果有WWGE卡，但卡上没有数据项，则产生10个能区，每个E i=10i-8MeV。

M. PDn 探测器贡献卡(Detector Contribution Card)

格式: PDn P1 P2 ... P l

n: 记数号

P i: 栅元i对探测器贡献的概率

l: 栅元总数

PDn卡减少对给定探测器相对不太重要的栅元对探测

器记数的贡献，节省了计算时间。在栅元i的每次碰撞，将以概率P i判断它是否参与探测器记数。为了保证结果无偏，当判定它参与记数时，权重除以P i，P i=0的栅元除外。用户可以对离探测器较远(与平均自由程比较)的栅元置以小于1的P i值，这样可以提高计算速度。也可以有选择地对某些栅元置P i=0，从而抑制这些栅元对探测器地贡献。

此卡缺省P i=1。可以建立缺省探测器贡献卡PD0。

N. DXC DXTRAN贡献卡(DXTRAN Contribution Card)

格式: DXC:n P1 P2 ... P l

n: N为中子, P为光子

P i: 栅元i对DXTRAN球贡献的概率

l: 栅元总数

此卡用法与PDn卡类似，但它控制对DXTRAN球的贡献。

O. 例子

例1 IMP:N .25 2M 2M 2M 1

问题中共有5个栅元,中子重要性分别是0.25、0.5、1、2和1。

例2 EXT:N 0 0 .7V2 S -SV2 -.6V9 0 .5V9

SZ -.4X

VECT V9 0 0 0 V2 1 1 1

EXT卡上10个数据项对应10个栅元，对栅元1、2、7

例3WWE:N E1 E2 E3

WWWN1:N W11 W12 W13 W14

WWWN2:N W21 W22 W23 W24

WWWN3:N W31 W32 W33 W34

这些卡对一个有四个栅元的中子问题定义了三个能量区间及相应的权窗边界。

例4WWN1:P W11 W12 W13

这个例子没有WWE卡，只对三个栅元定义与能量无关的光子权窗。

P. 思考题

1. 栅元重要性卡IMP有什么作用，如何填写。什么时候必须给出栅元重要性卡，什么时候可以省略IMP:P卡。

2. 栅元体积卡VOL和曲面面积卡

AREA有什么作用，如何填写。什么情况下使用它们，二者可被什么卡取代。

3. 什么问题中可以使用光子权重卡PWT，它有什么作用，如何填写。PWT卡上数据项为正、为负和为零各有什么含义，缺省为何值。判断是否产生光子时为什么要把W i乘以I i/I s。如何关掉某栅元的光子产生。

4. 什么情况下考虑使用指数变换卡EXT(及矢量输入卡VECT)。卡上数据项如何填写，数据项的正负号有什么含义。什么情况下考虑使用强迫碰撞卡FCL，卡上数据项的正负号有什么含义，数据项缺省为何值。不同时使用EXT卡和FCL卡的物理原因是什么。

5. 权窗卡是比较常用和比较重要的卡，它和重要性及能量劈裂卡相比有什么优点。如何用WWE卡设置权窗的能量区间。如何使用WWN卡设置某一栅元某一能区的权重下限，数据项为正为负各有什么含义。WWP卡上各参数有什么作用，当粒子权重低于权重下限或高于权重上限时MCNP如何处理。

6. 权窗产生器有什么作用。WWG卡上各项含义是什么。如何使用权窗产生器自动计算栅元权窗。

7. 探测器贡献卡PDn和DXTRAN贡献卡DXC各有什么作用。

§4.3 源的描述

MCNP3B提供4种描述源的方式：通用源(使用SDEF 卡)、曲面源(SSW和SSR卡)、临界源(KCODE和KSRC卡)和用户提供源(子程序SOURCE)。

曲面源描述方式是在初始运行MCNP时用SSW卡将以一定方向穿过某些曲面或进入某些栅元的粒子轨迹记录下来，接续计算时再以这些粒子作为源粒子，用SSR卡读出后继续计算。这里讲的曲面源描述方式和后面介绍的用SDEF卡在曲面上定义源的曲面源是两个概念。

临界源用于计算裂变问题。

如果INP文件中没有SDEF、SSR或KCODE卡，便采用用户提供的子程序SOURCE描述源。

由于我们常用的只是通用源描述方式，下面仅就这种源的描述方法进行介绍。描述通用源分布函数的相关卡还有SIn、SPn、SBn和DSn。

MCNP对所产生的每一个源粒子确定以下参量：

ERG 粒子能量

TME 初始时间

UUU,VVV,WWW 飞行方向

XXX,YYY,ZZZ 粒子位置

IPT 哪种粒子，中子=1，光子=2

WGT 粒子权重

ICL

发射粒子的栅元

JSU 发射曲面。若不在曲面上, =0 A. SDEF 通用源卡(General Source Define Card)

格式：SDEF 源变量1=说明源变量2=说明... ...

用途：使用通用源的问题需要此卡。

所有源变量都可以缺省，所有源变量都有缺省值。等号可以用空格替代。源变量的说明有三种形式：1) 显值；

2) D加分布号；3) 在另一个变量前加F, 后面跟着D和分布号。var=Dn意味着从分布n抽取源变量var的值，每个分布仅用于一个源变量。var Fvar’ Dn意味着从分布n抽取源变量var，分布n与变量var’有关，仅允许一级相关。上述描述形式可以理解为分三级：当源变量有显值或使用缺省值时仅存在第一级描述；当用概率分布给出源变量时出现第二级，这一级要求SI、SP卡或其一；当一个源变量与另一变量相关时出现第三级，这时要求有DS卡。

产生哪种源粒子(IPT值)由MODE卡确定，MODE N和MODE N P的源粒子是中子，MODE P的源粒子是光子。

SDEF卡的源变量和最终要确定的源粒子参量不完全是一一对应关系，下表给出源变量及其缺省值：

WGT和EFF必须是数值，不是分布。大部分源变量是标量，要求一个值，但VEC、POS和AXS是矢量，要求三个数为一组的说明，即矢量的三个分量。

源粒子初始位置的坐标由SUR、POS、RAD、EXT、AXS、X、Y、Z和CCC的各种组合来确定。用这些源变量可以指定三种不同类型的体分布和三种不同类型的曲面分布，这些分布的退化形式是线源和点源。

三个体分布是长方、球和圆柱体。在SDEF卡上如果有CEL变量，就可以用体分布在一个栅元内均匀抽样。

1) 用变量X、Y、Z指定一个长方体的体分布，栅元的面要垂直于坐标轴。如果这三个值都是常数，体源退化为点源。其它的退化形式是线源和平面矩形源。

2) 用变量POS和RAD指定一个球的体分布，这时不能使用X、Y、Z和AXS。POS定义球中心，RAD定义半径，在球面上均匀抽取粒子。如果RAD是一个a=2的幂函数分布(这是缺省情况)，则得到一个球体内的均匀抽样。

3) 用变量POS、AXS、RAD和EXT指定一个圆柱体分布。圆柱体的轴沿AXS方向通过POS。在以RAD为半径的圆周上均匀抽取粒子的位置，圆心在圆柱体轴上，该圆位于垂直于AXS离开POS距离EXT的面上。圆柱壳体内均匀分布的设置方法是，在SIn卡上填写圆柱两端离开POS的距离EXT、RAD的内半径和外半径，MCNP将按缺省值给出两个EXT之间的均匀抽样和两个RAD之间的幂抽样(a=1)。一个常用的退化情况是EXT=0，即平面圆形源。

对于曲面分布，SUR≠0。若指定X、Y、Z(必须确保其分布在指定的曲面上)，它们确定源粒子位置；若否，则在SUR指定的面上抽取这个位置，曲面的形状决定了抽取方法：

1) 如果是椭球面，则按面积均匀抽取粒子位置。

2) 如果是球面，若没有指定AXS，按面积均匀抽样；指定AXS后，EXT的值做为AXS方向和球的中心到粒子位置点矢量之间夹角的余弦。

3) 如果是平面，则按面积均匀抽样源粒子位置，除非还使用其它源变量描述。若指定POS值，必须保证该点在曲面上。如果再指定RAD，则在以POS为中心、RAD为半径的圆周上抽取粒子位置。如果RAD的分布是一个a=1的幂函数(缺省)，则在圆内按面积均匀抽样。

栅元取舍对栅元源和曲面源都适用。如CCC出现，按上述方法抽取的位置在CCC栅元内则接受，不在CCC栅元内的舍去并重新抽样。

源变量SUR、VEC、NRM和DIR确定源粒子初始飞行方向。若参照参考方向VEC抽取飞行方向，可从一个分布抽取极角余弦DIR值，方位角在0~360?之间均匀分布。若一个体分布(SUR=0)没有指定VEC和DIR，则为各向同性分布。如果对一个曲面源的分布没有指定VEC，按缺省使用

该曲面的由NRM指定符号的法线方向。如果对一个曲面源

的分布没有指定DIR，则按缺省使用余弦分布P(μ)=2μ,

0<μ<1。使用DIR的偏倚分布可使较多的源粒子朝着记数区

域发射，这时用指数分布(-31)一般较为合适。允许DIR取

离散值，DIR=1表示VEC方向上的单向源。

有效指标EFF适用于CCC或CEL的舍弃。如果在任何CEL或CCC指定的栅元中接受率太低，问题将被终止。终

止判据是: MAX(成功次数,10)/试验次数

B. SIn 源信息卡(Source Information Card)

C. SPn 源概率卡(Source Probability Card)

D. SBn 源偏倚卡(Source Bias Card)

格式: SIn option I1 . . . I k

n: 分布号, 1 ~ 999

option: I

i

的说明。允许值为:

空格或H –定义直方分布的箱边界，标量(缺省)

A –定义概率密度分布的一些点

L –定义离散的一些源变量值

S –定义一些分布号

I1 . . .I k : 源变量值或分布号

格式: SPn option P1 . . . P k

或者: SPn f a b

n: 分布号, 1 ~ 999

option: P

i

的说明。允许值为:

空格或D –定义SI卡上的H或L分布的概率。SI卡

上A分布的概率密度。(缺省)

C –定义SI卡上的H或L分布的累积概率

V –仅对栅元分布，定义概率与栅元体积成正比。

有P i时，要乘以P i

P1 . . .P k: 源变量概率

f: 内部函数的识别符

a b: 内部函数的输入参数

格式: SBn option B1 . . . B k

或者: SBn f a b

SBn卡上的n, option, f, a, b均和SPn卡一样，只有

一点不同，f只能取-21或-31。

B1 . . .B k : 源变量偏倚的概率

SP卡的第一种格式用于和相对应的SI卡联合定义一

个概率分布。SI卡的数据项是源变量的值或下一级分布号，SP卡的数据项是SI卡对应项的概率。

使用H选项时，SI卡上数据项的值应该单调上升，这

些数据值是箱边界。SP卡第一个数据项必须为0，后面的数

据项是相应箱的概率或累积概率，取决于使用D或C选项，

这些概率可不归一。抽样过程如下：先根据箱的概率抽取

一个箱，然后在这个箱内均匀抽样。

使用A选项时，SI卡上数据项的值应该单调上升，这

些数据值定义源变量的概率密度分布，最小和最大值定义

边界。SP卡数据项是对应SI卡数据的概率密度，这些概率

密度可不归一，第一项及最后一项数据一般为0。抽样过程中，在指定值之间的概率密度靠线性插值得到。

使用L选项时，SI卡上数据项是源变量离散值，如一

些栅元号或光子能谱中一些能量，不必单调上升。SP卡上

的数据项是离散点的概率或累积概率，取决于使用D或C选项。

S选项允许用另一个分布来描述本分布，MCNP能够

处理20层这种结构。此时每一个分布号前可加D，D也可以

省略。分布号为0表示使用这个源变量的缺省值。借助于S

选项，一个分布可以在多处出现，但一个完整的分布描述

只能用于一个源变量。

仅当源变量是CEL时才使用SP卡上的V选项。当栅元

网络工程设计教程课后答案beta版

《网络工程设计教程》 第一章网络工程设计概述 1. 网络工程的定义是什么？ 答：定义1：将系统化的、规范的、可度量的方法使用于网络系统的设计、建造和维护的过程，即将工程化思想使用于计算机网络系统中。 定义2：对定义1中所述方法的研究。 2.和网络工程有关的工作可以分为哪些阶段？每个阶段的主要任务是什么？ 答：1、选择系统集成商或设备供货商 网络系统的需求分析 逻辑网络设计 物理网络设计 系统安装和调试 系统测试和验收 用户培训和系统维护 4. 详细描述网络工程的系统集成模型。为何将该模型称为网络设计的系统集成模型？该模型具有哪些优点？为何要在实际工作中大量使用该模型？ 答：下图给出了网络工程的系统集成模型，该模型提出了设计和实现网络系统的系统化工程方法。虽然该模型支持带有反馈的循环，但若将该模型视为严格线性关系可能更易于处理。该模型从系统级开始，接着是用户需求分析、逻辑网络设计、物理网络设计和测试。由于在物理网络设计阶段，网络设计者通常是采用系统集成方法来设计实现网络的，因此将该模型称为网络工程的系统集成模型。 5. 简述系统集成的定义。试讨论系统集成主要有哪些好处？ 答：抽象地讲，系统是指为实现某一目标而使用的一组元素的有机结合，而系统本身又可作为一个元素单位（或称子系统或组件）参和多次组合，这种组合过程可概括为系统集成。 系统集成的好处： -质量水准较高：选择一流网络设备厂商的设备和系统，选择高水平的具有资质的系统集成商通常能够保证系统的质量水平，建造系统的风险较小。 -系统建设速度快：由多年从事系统集成工作的专家和配套的项目组进行集成，辅以畅通的国际厂商设备的进货渠道，及处理用户关系的丰富经验，能加快系统建设速度。 -交钥匙解决方案：系统集成商全权负责处理所有的工程事宜，而用户能够将注意力放在系统的使用要求上。 -标准化配置：由于系统集成商承担的系统存在共性，因此系统集成商会总结出它认为

《计算机网络》(第五版)期末复习资料

计算机网络复习要点 第一章 概述 1、计算机网络向用户提供的最重要的功能是连通性与共享性。（P1） 2、连接在因特网上的计算机称为主机（host ）（P3） 3、网络（network ）由若干结点（node ）和连接这些结点的链路（link ）组合。（P2） 4、简述Internet 和internet 的区别（P4） （1） internet （互联网或互连网）：是一个通用名词，它泛指由多个计算机网络互连而成的网络。在这些网络之间的通信协议（即通信规则）可以是任意的。 （2）Internet （因特网）：是一个专用名词，它采用 TCP/IP 协议族作为通信的规则。 区别：后者实际上是前者的双向应用。 5、因特网服务提供者是什么？（P4） ISP ，Internet Servvice Provider ，是一个进行商业活动的公司，又通常译为因特网服务提供商。 6、在网络边缘的端系统中运行的程序之间的通信方式通常可划分为两大类？（P8） 客户服务器方式（C/S 方式）和对等方式（P2P 方式） 7、路由器是实现分组交换的关键构件，其任务是转发收到的分组，这是网络核心部分最重要的功能。（P10） 8、分组交换网以分组作为数据传输单元。（P12） 9、定义：计算机网络是一些互相连接的、自制的计算机的集合[TANE03]？（P17） 10、速率即数据率（data rate ）或比特率（bit rate ），是计算机网络中最重要的一个性能指标，速率的单位为b/s(比特每秒)或bit/s ，有时也写为bps 等。(P18) 11、计算机网络中“带宽”是数字信道所能传送的“最高数据率”的同义语，单位是“ 比特/秒”（P19） 12、掌握传输（发送时延）发送数据时，数据块从结点进入到传输媒体所需要的时间公式：（P19~20） ） 发送速率（比特数据块长度（比特） 发送时延s / 13、协议是控制两个对等实体进行通信的规划的集合。(P30) 14、协议是“水平的”，控制对等实体之间的通信的规则。(P30) 15、服务是“垂直的”，是由下层向上层通过间接提供的。(P30) 16、同一系统的相邻两层的实体进行交互的地方，称为服务访问点（Service Access Point ）(P30) 第二章 物理层 1.单工，半双工，全双工通信的区别。P38 （1）单工：单向通信，即只能有一个方向的通信而没有反方向的交互。 （2）半双工：双向交替通信，即通信的双方都可以发送信息，但不能双方同时发送（当然 也就不能同时接收）。 （3）全双工通信：双向同时通信，即通信的双方可以同时发送和接收信息。 2. 最基本的二元制调制方法有哪几种。P38 （1）调幅(AM)：载波的振幅随基带数字信号而变化。

《实用写作》课程教学大纲

《实用写作》课程教学大纲 一、基本信息 课程编号: 12110460课程名称：实用写作 英文名称:Practical Writing 课程性质: 必修通识课程 总 学 时：36学 分：2 适用对象:非中文专业本科学生先修课程：中学语文，大学语文 二、编写说明 （一）课程的性质 《实用写作》是我校为非中文专业学生开设的必修通识课程，是提高学生应用文写作能力的基础课程，是一门工具性、实用性、实践性很强的课程。矚慫润厲钐瘗睞枥庑赖賃軔。 （二）课程教学目标和基本要求 《实用写作》课程旨在提高学生的写作素质和实际写作能力，以适应大学阶段的学习要求和毕业后就业、参加工作的需要。通过系统地学习应用文写作的基础理论，学习实践类、校园文化类、求职就业类、行政公文类等应用文文体知识，使学生系统掌握实用类文体的实际用途及其写作要领，使其实际写作水平得到切实的提高，以适应当前和今后在学习、工作、生活中的写作需要，为其总体素质和能力的提高提供必要的知识保证。聞創沟燴鐺險爱氇谴净祸測。 （三）课程的重点和难点 本课程的讲授为一个学期，分为四个部分，重点是学术论文、消息、求职信、求职简历、行政公文等实用文文体知识，难点是学习实用文写作的基础理论，掌握写作的基础知识与基本技能，掌握常用文体文本写作要领。残骛楼諍锩瀨濟溆塹籟婭骒。 （四）课程教学方法与手段 《实用写作》教学以案例教学法为主，注重常用文体的写作训练。在教学中整合运用讨论式、启发式、质疑式等教学方法，激发学生的学习兴趣。同时还要强化写作练习，以提高学生的能力写作。教学中充分利用多媒体辅助教学，提高教学效率，增强教学效果。酽锕极額閉镇桧 猪訣锥顧荭。 （五）实践环节 无。 （六）教学时数分配表 《实用写作》教学时数分配表 教学内容 各教学环节学时分配 采用何种多媒体教学手段 章节 主要内容 讲授 实验 讨论 习题 课 外 其它 小 计 1 绪论 1 1 1 3 多媒体 2 校园篇：请假条、证明信、调查报告、学术论 文、消息、通讯、计划、总结、 10 1 4 1 16 多媒体 3 求职篇：求职信、个人简历、演讲词、申论 5 1 2 1 9 多媒体 4 职场篇：广告文案、合同、国家行政机关公文 处理办法、通知、通告、通报、请示、报告、函、 8 1 2 1 12 多媒体

实用应用文写作习题答案

. . 练习题答案第一章应用文写作概述 1．什么是应用文？应用文的写作特点是什么？ 答：应用文是指国家机关、企事业单位、社会团体及人民群众在日常工作、生活中办理公务和个人事务时经常使用的具有惯用格式的文书。应用文的写作特点有实用性、时效性、真实性、规性和限制性。 2．应用文写作的基本要素包括那几个部分？ 答：应用文写作的基本要素包括主题、材料、结构和语言。 3．如何对应用文进行审查和修改？ 答：在审查应用文时可以自己审查，也可以请别人帮忙审查；在修改应用文时，应灵活掌握增、减、调、换四种修改方式的使用。 4．应用文的写作者应具备哪些能力？ 答：应用文的写作者应具备观察认知能力、思维能力、文字表达能力以及多媒体技术的运用能力。 第二章行政公文 一、判断题 1．×２．√３．√４．×５．×６．×７．√８．√９．× 二、问答题 1．行政公文作为行政机关行政的工具，具有哪些特点？ 答：行政公文作为行政机关行政的工具，具有政治性、法定刑、规性、强制性和规律性的特点。 2．通告的正文部分有哪些组成要素？ 答：通告的正文通常由引据、事项和结语三部分组成。 3．通知有哪些类型？各自的作用是什么？ 答：通知按其功能可分为指示性通知，事务性通知，批转、转发性通知和会议通知等类型，其中指示性通知用于上级机关向下级机关布置工作、做出相应指示；事务性通知用来传达、安排事务性工作；发布性通知用于颁布除国家

法律外的其他文件；批转、转发性通知用于将本机关收到的上级、下级或不相隶属机关的来文，转给所属下级机关或其他机关；会议通知用于召开比较重要的会议前，将会议的有关情况告知有关单位或人员。 4．在写作通报时应注意哪些事项？ 答：在写作通报时应注意事例要典型、评议要得当、叙述要全面简明、要掌握结构比例。 5．报告的主要作用是什么？具有哪些特点？ 答：报告主要适用于向上级机关汇报工作、反映情况，答复上级机关的询问，此外，报告还用于向上级机关报送材料。报告具有单向性和述性的特点。 6．请示最多有几个主送机关？其正文部分惯用的结语有哪些？ 答：请示只能有一个主送机关，不能多头请示。请示正文常用的结语有“当否，请审批”、“当否，请批示”、“以上请示，请批复”和“以上请示，盼复”等。 7．批复的标题有几种形式？各是由哪些元素组成的？ 答：批复的标题有5种形式，分别是发文机关+事由+文种、事由+文种、发文机关+事由+主送机关+文种、事由+主送机关+文种和复+原请示标题。 8．致函和复函的正文结构和惯用结语有何不同？ 答：致函的开头要开门见山，直接提出发函缘由，简述发函的目的、理由或依据。然后用简洁的语言将需要商洽、联系、询问、告知的事项明确地表达出来。复函是被动行文，开头要引叙来函的日期、标题、发文字号等，以体现其针对性。致函的结尾惯用“请大力协助为盼”、“望能同意”等做结语；复函的结尾通常使用“特此回复”、“特此复函”等惯用语做结语。 9．意见的主要作用是什么？具有哪些特点？ 答：意见适用于对重要问题提出见解和处理办法。意见具有广泛性、多向性和灵活性的特点。 第三章事务文书 一、判断题 1．√2．×3．√4．×5．√6．×7．√8．× 二、问答题 1．简报的作用是什么？具有哪些特点？ 答：简报是机关、团体、企事业单位部，或是某项中心工作、某次重要会议中，用以沟通信息、交流经验、反映情况、指导工作的一种期刊式常用文书。简报具有简、快、新的特点。 ·2·

计算机网络第五版答案完整版

计算机网络第五版答案完整版 计算机网络》课后习题答案 第一章概述 1-10 试在下列条件下比较电路交换和分组交换。要传送的报文共x( bit )，从源站到目的站共经过k 段 链路，每段链路的传播时延为d (s),数据率为C (bit/s )。在电路交换时电路的建立时间为s (s)。在分组交换时分组长度为p( bit )，且各结点的排队等待时间可忽略不计。问在怎样的条件下，分组交换的时延比电路交换的要小？答：对电路交换，当t=s时，链路建立；当t=s+x/C，发送完最后一bit ;当t=s+x/C+kd，所有的信息到达目的地。对分组交换,当t=x/C , 发送完最后一bit ；为到达目的地,最后一个分组需经过k-1 个分组交换机的转发，每次转发的时间为p/C,所以总的延迟=x/C+(k-1)p/C+kd所以当分组交换的时延小于电路交换x/C+(k-1)p/C+kd v s+x/C+kd 时，(k-1)p/C v s 1-11 在上题的分组交换网中，设报文长度和分组长度分别为x 和( p+h)( bit )，其中p 为分组的数据部分的长度，而h 为每个分组所带的控制信息固定长度，与p 的大小无关。通信的两端共经过k 段链路。链路的数据率为b( bit/s )，但传播时延和结点的排队时间均可忽略不计。若打算使总的时延为最小，问分组的数据部分长度p 应取为多大？答：分组个x/p ，传输的总比特数：(p+h)x/p 源发送时延：(p+h)x/pb 最后一个分组经过k-1 个分组交换机的转发，中间发送时延：(k-1)(p+h)/b 总发送时延D=源发送时延+中间发送时延D=(p+h)x/pb+(k-1)(p+h)/b 令其对p的导数等于0，求极值p=Vhx/(k-1) 1-17收发两端之间的传输距离为1000km信号在媒体上的传播速率为2.3 X 108。试计算以下两种情况的发送时延和传播时延：(1)数据长度为107bit，数据发送速率为100kbit/s，传播距离为1000km信号在媒体上的传播速率为 2X108m/s°(2)数据长度为103bit，数据发送速率为1Gbit/s，传输距离和信号在媒体上的传播速率同上。 答：(1)：发送延迟=107/(100X1000) =100s 传播延迟=1000X1000/(2X108) =5X10-3s=5ms (2)：发送延迟=103/ (109) =10-6s=1us传播延迟=1000X 1000/ (2X 108) =5X 10-3s=5ms 1-18、假设信号在媒体上的传播速率为2.3 X108m/s。媒体长度I分别为： (1) 10cm( 网卡)(2) 100m( 局域网)(3) 100km( 城域网)(4) 5000km( 广域网)

《应用文写作教程》教学大纲

《应用文写作教程》教学大纲(段轩如、高玲主编)人大 一、课程性质及教学目的要求 应用文写作,是指在日常工作、生活、生产、学习、科研中广泛应用的文书写作。《应用文写作》是面向在校大学生开设的一门公共课,也可以是相关专业的专业基础课。它既是一门研究应用写作规律与方法的学科,又是一门实践性、应用性较强的学科,具有综合性、工具性的特征。实用性强,适用面广,是本课程的特点。 通过应用文写作基础理论和各种应用知识的教学与写作训练,使学生掌握应用文体的格式规范并熟练地写作各类各种应用文体,以适应在学习、生活、工作以及科学研究中的写作需要。 课程具体教学要求： 本课程的理论性不是很强,但基本理论和基本知识的讲授应简明、准确； 应用类文体的规范性是至关紧要的,要引导学生多接触文章实际,加深对所学文体的全面的认识,补充的例文一定要选用规范的文本； 课堂讲授与讲授、训练相结合,并适当增加训练比重,把理论讲授的内容通过练习巩固提高,重在培养学生的实际写作能力和文字信息处理能力； 有条件的话,建议做应用写作模拟训练：老师把搜集大量的资料,并设计成合理有序的背景材料,创造一个逼真的写作环境,然后引导学生进入这个写作场景中,进入应有的角色定位,运用所提供的材料完成角色需要完成写作任务。作业设计方面,建议对基本知识的考察的同时,更要引导学生关注实际问题,培养学生的独立思考能力和分析研究能力。通过对作业的批改和课堂讲解,进一步强化学生的理论知识和写作能力。 学时分配：本课程讲授内容及详略以够用为度,教学总时数可因专业性质对应用写作的不同要求而定,一般为36—54学时,其中实训不低于总学时的三分之一。 二、教学内容与重点 第一章应用写作基本理论 【教学目的】 通过学习明确应用类文章的基本特点与作用等总体情况,掌握应用写作相关主题、材料、结构、语言等基础知识,把握文体的一般特征。了解应用写作的含义与种类及提高应用写作

常用写作文体

常用写作文体 一、应用文的概念 应用文是国家机关、企事业单位、社会团体和个人在日常生活、学习、工作中处理公私事务所使用的具有某种固定格式和直接实用价值的文章。 应用文也称“实用文”是人类在长期的社会实践活动中形成的带有实际应用性质的一种文体是人们用于传递信息、处理事务、交流感情的一种工具有的应用文还可作为凭证和依据。在各类文体中应用文的使用频率最高。党政机关、企事业单位用它上传下达办理公务普通人用它来交流信息和思想。 应用文与我们日常生活和工作有密切的关系。我们要了解天下大事就要阅读报刊、收听广播、收看电视以法治国要有各种法规文件召开会议要有会议文件党政机关指导工作要有许多公文机关、企事业单位要正常运转要有计划、总结、报告等事务文书人们礼尚往来常常借助于请柬、贺卡等等。这诸多的文字材料大部分都是应用文。可见大到整个国家小到某一个单位甚至个人要进行正常的活动都离不开应用文当然也就离不开应用写作。 随着社会的不断发展信息化特征日益显著人们在工作和生活中的交往越来越频繁事情也越来越复杂应用文的功能也就越来越多应用文写作的重要性愈加突出各种应用文书更成为各级机关、企业、团体开展各项公务、实施有效管理的一种重要工具。 写好应用文不管是对个人还是对单位及社团塑造自身形象处理

各种关系都起着重要作用。 因此应用文写作成为现代公民必备的素质之一。 二、应用文的种类 应用文种类繁多可以从不同的角度划分成不同的类别。 按照其处理事情的性质的不同可以分为行政机关、团体和企事业单位用来处理公务的公务类应用文亦称为公文和个人或集体用来处理个人事务的私务类应用文即个人日常应用文书两大类； 按照应用领域的不同又可分为： （一）公务应用文 公务应用文包括行政机关公文和行政事务文书。 1、行政机关公文（法定公文） 根据国务院办公厅年印发的枟国家行政机关公文处理办法枠的规定行政机关公文包括命令令、决定、公告、通告、通知、通报、议案、报告、请示、批复、意见、函、会议纪要十三类十三种公文。这类公文是行政机关在行政管理过程中形成的具有法定效力和规范体式的文书是依法行政和进行公务活动的重要工具主要用来传达和贯彻党和国家的政策法令、指导工作、提出要求、答复问题、通报情况、交流经验、传递信息。这类公文制作比较严格具有一定的法律效力。 2、行政事务文书 行政事务文书是国家机关、团体、企事业单位在处理单位内部日常事务或与具体部门进行工作联系中普遍使用的公务文件亦称为“常规文书”。行政事务文书包括简报、计划、总结、调查报告、规章制

计算机网络(第五版·谢希仁)知识点

第一章概要 1、网络是指“三网”，即电信网络、有线电视网络和计算机网络。 2、计算机网络向用户提供的最重要的功能： （1）连通性：计算机网络使上网用户之间都可以交换信息，好像这些用户的计算机都可以彼此直接连通一样。（2）共享：就是指资源共享，如：信息共享、软件共享、硬件共享。 3、因特网的组成： （1）边缘部分：有所有连接在因特网上的主机组成。用户直接使用 ·在网络边缘的端系统中运行的程序之间的通信方式：客户服务器方式（C/S方式）和对等方式（P2P方式）。·客户服务器方式主要特征：客户是服务请求方，服务器是服务提供方。 （2）核心部分：由大量网络和连接这些网络的路由器组成。为边缘部分提供服务 ·主机的用途是为用户进行信息处理的，并且可以和其他主机通过网络交换信息。路由器的用途则是用来转发分组的，即进行分组交换的。 ·路由器是实现分组的关键构件，其任务的转发收到的分组，这是网络核心部分最重要的功能。 4、电路交换：电路交换必定是面向连接的。 ·电路交换的三个阶段：建立连接（占用通信资源）、通信（一直占用通信资源）、释放连接（归还通信资源）。·电路交换的重要特点：在通话的全部时间内，通话的两个用户始终占用端到端的通信资源。 5、分组交换：（面向无连接的） ·分组交换采用存储转发技术。路由器收到一个分组，先暂时存储下来，在检查其首部，查找转发表，按照首部中的目的地址，找到合适的接口转发出去，把分组交给下一个路由器。当路由器满了以后采用丢弃处理。 ·主机和路由器都是计算机，但作用不一样。主机是为用户进行信息处理的。路由器则是用来转发分组的，即进行分组交换的。 6、三种交换方式在数据传送阶段的主要特点：（因特网中三种核心交换方式？概念、特点） （1）电路交换：整个报文的比特流连续地从源点直达终点，好像在一个管道中传送。 （2）报文交换：整个报文先传送到相邻结点，全部存储下来后查找转发表，转发到下一个结点。 （3）分组交换：单个分组（这只是整个报文的一部分）传送到相邻结点，存储下来后查找转发表，转发到下一个结点。 7、几种不同类别的网络：（这几种方式是什么？） ·不同作用范围的网络： （1）广域网W AN （2）城域网MAN （3）局域网LAN （4）个人区域网PAN ·不同使用者的网络： （1）公用网（public network）（2）专用网（private network） ·用来把用户接入到因特网的网络： 接入网AN (Access Network)，它又称为本地接入网或居民接入网。 注：由ISP 提供的接入网只是起到让用户能够与因特网连接的“桥梁”作用。 8、计算机网络的性能指标：（各自含义） （1）速率：指的是连接在计算机网络上的主机在数字信道上传送数据的速率，也称为数据率或比特率。 （2）带宽：用来表示网络的通信线路所能传送数据的能力。因此网络带宽表示在单位时间内从网络中某一点到另一点所能通过的“最高数据率”。 （3）吞吐量：表示在单位时间内通过某个网络（或信道、接口）的数据量。 （4）时延(delay 或latency)是指数据（一个报文、分组，甚至比特）从网络（或链路）的一端传送到另一端所需的时间。 ·发送时延（传输时延) ——是主机或路由器发送数据帧所需要的时间，也就是从发送数据帧的第一个比特算起，到该帧的最后一个比特发送完毕所需的时间。 ·传播时延——电磁波在信道中需要传播一定的距离需要花费的时间。 注：信号传输速率（即发送速率）和信号在信道上的传播速率是完全不同的概念。 ·处理时延——主机或路由器在收到分组时要花费一定的时间进行处理。（交换结点为存储转发而进行一些必要的处理所花费的时间。）

网文写作教程编辑说你小说节奏太快慢这里有药

网文写作教程16：编辑说你小说节奏太快慢这里有药 今天要讲的这个主题又是一个很玄妙很抽象的问题——小说的节奏。这种看起来比较高深的写作技巧，实难有什么高明见解，无非根据已有的见解看法，望能将此事说个清楚明白。一、节奏的定义 说到节奏，其实我们条件反射想到的应该是音乐的节奏，如果我们评价一首音乐很有节奏感，应该是觉得这首音乐的节拍快慢、强弱组合很有规律，配合得恰到好处，能给人一种美妙的享受。 再然后我们还能想到的最平常的关于节奏的场景，应该是生活节奏。我们经常觉得，像北上广深那些一线城市的生活节奏都太快了，这个语境下的节奏，我们大概可以理解为，一种很匆忙的状态。 而为什么会很匆忙呢我们就不去探讨什么人多交通拥堵等问题了，无非是时间少而事情又多，于是就要在一定时间内做好一定事情，需要速度。 于是，我们说小说中的节奏就包揽了以上两个场景的含义，一重含义指情节的高低起伏和轻重缓急，二重含义是指在单位字数内要叙述的情节数量多少，相同的情节下，用较多的篇幅去描写叙述则感觉节奏慢，反之则节奏快。 下面我们就具体讲以上两重含义及应用。

二、节奏——高低起伏，详略得当 在情节起伏的问题上，出现最多的就是如下两种情形： 没有起伏，流水账式记录。2.没有铺垫，处处是高潮，也就无所谓高潮。 第一种，流水账式小说，就像我们平常写日记一样，有什么事说什么事，很无聊的记录，或者一直在铺垫，无法挑起读者兴趣。 怎么解决最好的方法是有选择的记录，剪掉不必要的铺陈的平淡情节，而是选取能产生矛盾的的故事，要有冲突，文似看山不喜平。 第二种，高潮式小说，情节的安排持续给读者紧张的感觉，没有铺垫，没有缓冲，这种情况读者看着会非常累，也很容易写崩文。 解决方法除了要学会铺垫（讲起来太多了，有机会可以详讲），还应跟第一种情形学习，学会注水，多点轻松的情节描写。 这让我想起看那种悬疑类的影视剧，它并不是一直紧张到让观众无法喘息，会偶尔穿插个主角的生活场景，比如吃个饭，浇个花，跟朋友打个趣等，看上去无关紧要的情节，但却是非常必要的。三、男频与女频的节奏差异 网文中有一种说法，三章一小高潮，五章一大高潮，也就是说长篇小说中的高低起伏其实应该是循环反复的，我们认为

应用文写作复习答案

应用文写作重点考试资料 名词解释 1.应用文书：指日常应用、运用的一些文章体裁，其概念小于实用文体，但包含在实用文体之内。 2.规章制度：是用人单位制定的组织劳动过程和进行劳动管理的规则和制度的总和。也称为内部劳动规则，是企业内部的“法律”。 3.通知：是党政机关、企事业单位广泛使用的下行文文种；适用于批转下级公文、转发上级和不相隶属机关的公文；传达要求下级机关办理或执行的事项；任免人员等。 单项1.找出下面各项中不符合规章制度有强制性特点的一项规章_D公约。 2.找出不属于行政公文的一项 C制度_。 3.下边的公文文种属于下行文的一项是_ A通报_。 4.对函的特点论述不正确的一项是_D用语简洁性。 改病句1.机构调整问题的请示。加“关于” 2.公文质量的高低都有损发文单位的形象。将“损”改为“关” 3.对下级机关的公文请示性函需要答复的，应用“批复”。用“复函” 简答题 1.简述公文报告与请示的异同。 答：两种公文都属于上行文；不同点是：请示需上级批复；报告不许批复；请示是事前行文，报告则是事中和事后；请示需一文一事，报告则可一文多事。 2.简述计划与总结的写作要点。 答：总结的要点是：标题；前言；正文：成绩经验；问题或教训；今后努力的方向；计划的正文要点是：目的；方法措施；时间安排与步骤。 3.写出经济合同导语与正文的的内容要点。 答：合同的导语：双方经协商，本着互利互惠的原则，签订本合同，以便共同遵守；主体是：标的；数量和质量；价款与酬金；履行合同的时间地点和方式；违约责任等。 4.简述通知的种类和特点。 答：通知的种类：印发、批转、转发性通知；指示性通知；知照性通知；事务性通知；任免、聘用通知。通知的特点：知照性（通知的主要功能在于知照）广泛性（通知的广泛性表现在多方面）时效性（通知有一定的时效要求）

简单的网页制作教程-设计一个个人网站

题目：设计一个个人网站 一、要求： 1．使用Dreamweave网页工具制作一个个人网站； 2．包含至少四个网页: 包括首页、个人简介、个人相册等（可随意设计），网页之间用超链接相连。 3．网页中要有图片和文字内容，用表格进行页面布局； 4．添加至少两种行为，并为首页添加背景音乐。 5. 在网站中设计一个表单页面。 6. 首页必须包含页面标题，动态按钮导航栏。 首先新建一个文件夹，文件夹的名字不能为汉字，做网站所有的路径都必须用字母或者数字， 不能用汉字，我们就用名字吧，譬如说名字张三，那文件夹名字就是zs，如图 打开Dreamweaver软件，得到图 做网页要新建站点，关于站点配置服务器什么的，这里不讲了，只讲建立站点。 选择站点——新建站点。 我们建的文件夹就是站点根文件夹。

新建站点后得到这样一个界面 点选高级，得到界面 站点名称与我们建文件夹得名字相同，zs填进去就可以了本地根文件夹就是我们新建的那个文件夹zs， http地址为http://localhost/zs

接下来选择左侧栏里远程信息 点击无后面的那个三角，选择本地网络，远端文件夹同样选择我们新建的那个文件夹 接下来点选左面菜单里的测试服务器， 点选访问后面那个三角，选择本地网络，测试服务器文件夹也为我们建好的文件夹zs，在url前缀后面加上zs

然后点击确定就可以了得到这样一个界面。 下面看老师的第一条要求，是要至少四个网页，那我们就做四个 单击新建，然后单击 接下来，选择 然后单击创建，接下单击文件——保存，保存这个文件，保存在我们一开始建好的文件夹里面，保存名字不能是汉字，只能是字母或者数字，因为我们只坐四个网页，可以简单一点，把这四个网页命名为a、b、c、d，或者1、2、3、4，当然一个网站默认的索引首页名为index，这里也用index，

计算机网络教程第五版微课版谢钧谢希仁编著课后习题参考答案

】 计算机网络教程第五版（微课版）答案 第一章概述 1-01 计算机网络向用户可以提供那些服务 答：连通性和共享 1-02 简述分组交换的要点。 答：（1）报文分组，加首部（2）经路由器储存转发（3）在目的地合并 1-03 试从多个方面比较电路交换、报文交换和分组交换的主要优缺点。 ; 答：（1）电路交换：端对端通信质量因约定了通信资源获得可靠保障，对连续传送大量数据效率高。 （2）报文交换：无须预约传输带宽，动态逐段利用传输带宽对突发式数据通信效率高 ，通信迅速。 （3）分组交换：具有报文交换之高效、迅速的要点，且各分组小，路由灵活，网络生 存性能好。 1-04 为什么说因特网是自印刷术以来人类通信方面最大的变革 答：融合其他通信网络，在信息化过程中起核心作用，提供最好的连通性和信息共享

，第一次提供了各种媒体形式的实时交互能力。 [ 1-05 因特网的发展大致分为哪几个阶段请指出这几个阶段的主要特点。 答：从单个网络APPANET向互联网发展；TCP/IP协议的初步成型 建成三级结构的Internet；分为主干网、地区网和校园网； 形成多层次ISP结构的Internet；ISP首次出现。 1-06 简述因特网标准制定的几个阶段 答：（1）因特网草案(Internet Draft) ——在这个阶段还不是 RFC 文档。 （2）建议标准(Proposed Standard) ——从这个阶段开始就成为 RFC 文档。 （3）草案标准(Draft Standard) — （4）因特网标准(Internet Standard) 1-07小写和大写开头的英文名字 internet 和Internet在意思上有何重要区别 答：（1） internet（互联网或互连网）：通用名词，它泛指由多个计算机网络 互连而成的网络。；协议无特指 （2）Internet（因特网）：专用名词，特指采用 TCP/IP 协议的互联网络

写作课程笔记整理

第五章文学写作 第一节诗歌 一．.概述 1.诗歌是一种以简洁凝练的语言表达、分行分节的句式章法、流转和谐的节奏韵律、蕴含情思的意象组合，借助于联想想象书法内心情感的文学样式。 诗是心灵的歌，感情的火，思想的光。 2.诗歌—按内容（叙事诗、抒情诗、哲理诗） —按形式（格律诗、自由诗） —按发展（古体诗、近体诗、新体诗） —按杂交（散文诗、寓言诗、童话诗） 3.诗歌的特点：表现生活的抒情性、流转和谐的音乐性、蕴籍形象的精炼性 二．诗歌的意象、意境 1.意象组合 意象中的意是意念，指创作主体的主观情感；象是具象，指表现对象的客观物态 意象是主观情感与客观物态相互融合、相互渗透所形成的审美契合的产物 2.意境创造 意境指情景交融、物我一体、主客体统一的整体意象所显示的艺术境界 诗歌的意境有大境界与小境界、所造之境与所写之境、宏伟之境与优美之境、有我之境与无我之境等范畴（1）有我之境：诗人所描绘的具有强烈主观情感色彩的境界 （2）无我之境：诗人把主观情感含而不露地融化在所描绘的景物之中,营造出”一切景语皆情语”的境界三．诗歌的语言锤炼 1.独特搭配 2.欧化叠句 3.提行转句 4.省略跳跃 5.想象比喻 第二节散文 一．概述 1.散文是出诗歌以外的一切文学样式的源头和母体，是一种在字里行间流露真情、表抒真心，述说真意、传递真思，让读者在阅读时产生确有其人、确有其事、确有其情的心理认同的文学样式 2.散文—广义散文-与韵文骈文相对，之不押韵、不重排偶的文章，包括经书、史书 —次义散文-指诗歌、小说、戏剧以外的文学作品，包括报告文学、特写、杂文 —狭义散文-叙写“我”的经历、见闻、感受的文学作品，强调“自我亲感”的文体特征 散文分类—表达方式-记叙性、抒情性、议论性散文 —表现对象-感性、知性、审美、审智、都市、乡土、文化、旅游、闲适、纪实、网络、运动、生态、新生代散文

《经济应用文写作》课件——日常应用文的写作《经济应用文写作》课件第四章

第四章日常应用文的写作 本章内容主要讲授机关单位常用应用文，如计划、总结、调查报告等的知识、结构、写法及注意事项等。 这些章节，关键点在于内容和写法上的处理。学习时应该加以重视，最好能以形式多样的练习加以巩固。 因此，本章安排两次作业，一是调查报告，二是总结。 第一节工作计划 一、工作计划的含义、特性和订立工作计划的意义。 （一）概念 人们根据一定的原则，对未来一定时期要做什么和准备怎样去做所作的具体安排和打算。 （二）作用 凡事预则立，不预则废。好的计划能有效地集中人力、物力、财力，更好地达成某个目标。 二、工作计划的种类及其写作特点。 （一）按计划时间和内容要求的不同，分别有计划、规划、设想、要点、意见、安排、打算、方案等。 对以上不同叫法的计划，应该注意加以区分。如： 规划：带有全局性、方向性的宏观方面的计划，时间较长。 设想：一种初步的、粗线条的、时间长而又非正式的计划。 安排：时间较短、内容较单一、具体。 要点和意见：领导机关向所属单位布置一个阶段或一项重要工作时，交代有关政策，提出具体要求。 方案：非常具体、成熟的计划

练习：填上合适的叫法： 1、广东省十年消灭荒山。 2、田野公司第一季度工作。 3、广东金融学院第九周工作。 4、广州市工资与物价挂钩。 （二）、按内容涉及方面来分有综合计划、专题计划。 （三）按形式分有条文式计划、表格式计划、条文兼表格式计划。 （四）特征 1、预见性：从分析现状准确预测未来 2、指导性：对工作有实际指导意义 3、实践性：必须能实施，把计划变为现实 二、计划的结构、基本内容和写法。 （一）计划前的准备工作 1、结合实际、开拓创新 2、具体明确、重点突出 3、重视反馈、及时调整 （二）计划的内容 这是计划写作的一个关键。学生习作上最大的问题往往出在这里。要想写好计划，首先要搞清楚的是计划到底应该写什么内容。 1、依据——为什么做？ 2、任务——要做什么？ 3、措施或步骤——怎样去做？ 4、目标——做到什么程度？ 归纳起来，计划的内容包括三部分： 1、依据和目标 2、任务或指标 3、措施或步骤

网络课程设计题目

网络课程设计 题目一 基本要求：根据用户需求，设计网络，并完成相关文档和文件工作。要求通过查找资料，独立完成设计，全部图、表只能使用WORD或VISIO的相关工具来画，不得粘贴扫描的图片。路由器和交换机、PC机配置利用boson netsim或类似软件来辅助进行，防火墙、服务器配置用文档描述。 1、某高校要求设计一个校园网， 一、用户需求 （1）用户规模500台计算机。 （2）用户大致平均分散在4栋楼房内，4栋楼房排成前后两排，楼房之间各相距200米，楼房高4层。每栋楼的4楼用户构成两个VLAN。 （3）中心机房设在其中1栋楼房的1楼靠近另一栋楼房的一端。 （4）安装对外WWW、业务WWW、邮件、FTP、BBS、DNS、数据库七个服务器。提供匿名服务，但FTP仅对内部开放。 （5）提供LAN、WLAN接入。 （6）在业务WWW服务器上配备基于Web的业务应用系统，所有用户使用业务系统实现网上办公。 （7）要求出口带宽为1Gbps。 二、设计要求 （1）写出简要的可行性分析报告。 （2）设计网络结构，并给出解释。 （3）除用户计算机已购置外，其余全部设备和通信线路需要重新购买、安装。试具体给出全部主要设备的配置、型号或技术指标及其测算依据。 （4）给出工程预算（包括设备、线路等，不含施工费）及其计算依据。 题目二 设计一个中小企业网络规划与设计的方案： 一、用户需求 （1）公司有1000 台PC （2）公司共有7个部门，不同部门的相互访问要求有限制，公司有3个跨省的分公司。（3）公司有自己的内部网页与外部网站，公司能够提供匿名的FTP,邮件，WWW服务，但FTP只对内部员工开放。 （4）公司有自己的OA系统 （5）公司中的每台机能上互联网，每个部门的办公室联合构成一个VLAN。 （6）核心技术采用VPN。 二、设计要求

网络工程师第五版教材解读

2018年网络工程师第五版最新教材解读2018年网络工程师第5版教程已经上市，那么与第4版相比有哪些变动呢？网络工程师新旧教程对比变化大吗？ 与第4版教程内容相比，第5版结合新版考试大纲的要求，对考试内容进行了较大的调整，考生们在备考的时候可以根据最新版教程来进行复习。 至于网络工程师新旧教程的具体变化，想了解的考生可以在希赛网观看2018网络工程师教程对比公开课，希赛讲师为大家解读了新版教材与老版教材的变化，让考生们更加深入了解第5版教材的改版内容。 网络工程师新版教程具体目录如下： 第1章?计算机网络概论1 ?计算机网络的形成和发展1 ?计算机网络的分类和应用3 计算机网络的分类3 计算机网络的应用6 ?我国互联网的发展7 ?计算机网络体系结构9 计算机网络的功能特性9 开放系统互连参考模型的

基本概念11 ?几种商用网络的体系结构18 ?OSI协议集22 第2章?数据通信基础27 ?数据通信的基本概念27 ?信道特性28 信道带宽28 误码率30 信道延迟30 ?传输介质30 双绞线30 同轴电缆31 光缆33 无线信道34 ?数据编码35 ?数字调制技术38

?脉冲编码调制40 取样40 量化41 编码41 ?通信方式和交换方式41 数据通信方式41 交换方式43 ?多路复用技术46 频分多路复用46 时分多路复用47 波分多路复用48 数字传输系统48 同步数字系列50 ?差错控制50 检错码51 海明码51 循环冗余校验码53 第3章?广域通信网56 ?公共交换电话网56

电话系统的结构56 本地回路57 调制解调器61 ?公共数据网63 流量控制和差错控制64 协议69 协议75 ?帧中继网78 帧中继业务78 帧中继协议80 帧中继的应用82 ?ISDN和ATM83 综合业务数字网83 虚电路87 高层88 适配层89 通信管理90 第4章?局域网与城域网93 ?局域网技术概论93

科幻网文写作培训 点子与故事——纠结不休的暖昧

科幻联合写作班 第一期课程　 点子与故事——纠结不休的暖昧 讲师：说书人 本节课的内容主要是关于点子与构思之间的关系。 一般来说，科幻有时被戏称为点子文学，虽然这是一种非常片面的看法，但不能不承认点子对于科幻的巨大影响。 一般而言，非科幻的文学类型关注的主题是关注于人或处于社会中的人。而科幻小说将外部世界置于与人本身同等重要的地位。 这就引申出一个问题：究竟用什么方法来表现物质的世界与人类本身的交互？ 由之，点子应运而生了。 所谓点子，就是基于科学之可能（注重不是现实），来假设人类在特定的技术条件下，可能具有的戏剧性行为。更简单的说法，就是科技进步对于人类个体或者社会产生了巨大的冲击。科幻小说就是借助一个个具体的点子，来考察人在宇宙中的地位，人与宇宙的关系，在这样的宇宙中人与人之间的关系。 科学技术的的进步，在人与自然之间，产生了某种的紧张关系，打破了原有的平衡。而科学小说试图通过虚构去解释，这种平衡是如何被打破的，更重要的是又是如何被修复，回到某种新的平衡状态的。 第一节为什么会有点子 1、为什么科幻小说往往被称为点子文学，而其他类型的小说不享有如此“殊荣”。 2、点子文学的内在的成因：科学与人类的紧张关系。 这里简单提一下：点子就那种将物理背景与小说连结起来了的加固件。 要理清，点子与故事的关系。首先要搞懂什么是故事。故事本质上来是以戏剧化方式的描述人在世界中的行为。一般来说，人在世界上有各种行为，绝大部分是按部就班，是稀松平常的，只有其中少数突变性的偶然行为是具有新意的，故事就是去描述这样的偶然性。 这里要强调一点，正是因为故事只描述了少数的偶然性，所以故事比现实生活要简单的得多。千万不能认为故事比生活更复杂，恰恰相反，生活是一只薛定谔的猫，而故事是一个

新编应用文写作教程刘宏彬

新编应用文写作教?程刘宏彬 新编应?用文写作教程刘宏?彬 新编应用文写?作教程【完整版?794p】食品?安全法》的实施将?会极大提高群众对?于食品消费的积极?性因为放下了心所?以才敢放进了肚消?费信心建立起来消?费热情调动起来于?国于民都将有极大?益处2过去食品安?全有两个顽疾一个?是乱抽检乱收费另?外一个是乱发布信?息同一个企业产业?有些部门发布说安?全另外一些部门发?布信息却说不安全?由于监管部门的多?头执法重复监管也?无形中增加了企业?的负担例如农产品?的种植养殖由农业?部门管食品生产由?质监部门管超市商?场等食品流通环节?由工商部门管小餐?馆大饭店等餐饮服?务由食品药品监督?管理局管对食品安?全的监管机制人们?笑称五龙治水等于?无龙治水我国的监?管模式原来是食品?卫生部门一家负责?但是由于食品生产?的链条比较长从农?田到餐桌一个部门?管力不从心所以形?成了多部门管理的?监管体制但多部门?管理易出现管理空?白等问题以三鹿奶?粉事件为例奶粉生?产链条当中的奶站?是哪个部门也不监?管的真空地带教训?非常深刻监管的真?空地带怎么来弥补?在《食品安全法》?当中成立国家食品?安全委员会是最重?要的新措施新法?重新明确了各个部?门的监管职责确立?了分段监管体制主?要是卫生农业质检?工商和食药监各司?其职分别负责对食?品安全风险的评估?食品标准的制定对?初期农产品对食品?生产环节食品流通?环节和餐饮服务方?面的监管即从原料?到产品从生产到流?通餐饮的全程监管?在此基础上设立国?务院食品安全委员?会加强对各个监管?部门监管发挥协调?组织惩治监督等职?能作用最重要的是?要锁定责任实现无?缝对接纵要到底横?要有边不能留任何?真空设立食品安全?委员会是对分段监?管体制的很大推进?因为有了一个国家?层级的协调机构可?以对各个下属监管?部门的情况进行综?合分析协调处理把?安全漏洞补上3孩?子们一放学就跑到?学校周边的小吃摊?买零食可那里的卫?生根本没有保障吃?坏了肚子可咋办我?家孩子昨天中午在?学校边上买了几种?油炸食品回家后胃?就不舒服下午还发?生呕吐现象到医院?检查确诊是因为食?用不洁食物导致胃?肠道疾病打点滴后?症状才缓解这些路?边摊卖的黑食品太?坑人了有关部门是?不是该管一管了?3月17日中午?记者走访齐齐哈尔?部分中小学门前时?发现这

日常应用文

关于理工科高校实用文体写作课的思考 实用文体写作有很强的适用价值，它是各种文体中适用最广的一种文体。自古以来，它在人们的社会生活中占有极其重要的地位，是人类社会生活不可缺少的交际工具。刘勰在《文心雕龙》中称之为“虽艺术之末品，而政事之先务也。”在现代化经济建设的今天，以适用为目的的实用文体写作已成为每一个大学生必须具备的基本技能。 当今是信息、科技和知识经济时代，这个时代的特征之一，就是大量的人类未知领域被认识，被发现，被开发。这些科研成果如何才能传播出去，使之转化成生产力，这就要靠“文章”这个媒介和途径。这个时代的另一个特征就是重视和强调应用。任何发明创造只有转化为产品和商品时，才能最终体现科学技术的价值和威力，这也要靠“文章”来承担此重任。现在在校的大学生是祖国的未来和希望，他们将来要在各自的专业上有所造诣和发展。他们毕业后将从事经营、管理或科学研究工作，在这些领域中，只要有发现和创造，就必须要反映，要传播。试想：一个大学毕业生如果不能运用祖国的语言文字写出文理通顺的文章，即使在科技方面研究出成果，又怎能转化为生产力呢？从这个角度讲，实用文体写作显得至关重要。如果他们将来所从事的是企业管理工作，那么他们将经常使用通知、报告、请示、函、计划和总结等文种；一个大学毕业生不会写常用公文，不会写计划和总结，甚至连工作经验也不会总结，他又怎能担负起相应的工作职责，发挥一个大学生应有的作用呢？因此，新的时代给现在的理工科大学生的写作能力提出了新的要求———既懂专业，又会写文章。：“如果没有写作能力，那么人的思想感情就不能表达，经验就不能总结，不能记载下来也就不能交流。人类的科学文化的进步，很重要的一条，就是能够总结经验，并把经验记载下来，相互交流，传之后代。”“一个国家，如果人民缺乏写作能力，科学文化就不能提高。”周扬也曾说过：“从培养写作能力的角度看，不论是想当作家的人，还是不当作家的人，都应努力提高写作水平，这是关系到我们国家的科学文化水平的提高问题，一个有文化的国家，人人都应获得这种能力。”!"我国著名化学家、中科院院长卢嘉锡曾在有关会议上号召：“培养科学工作者的老师们，要教会年轻人学会表达，表达是很重要的。一个只会创造，不会表达的人，不能算是一个真正的合格的科学工作者。”从这些讲话中可以看出提高写作素养和文字表达能力历来都被重视，这是时代发展的需要，经济建设的需要，作好各项工作的需要。古往今来，凡在大事业上有所建树的科学家，他们大都具有渊博的知识、深厚的语文基础和很强的写作能力。如我国古代的科学家张衡、郦道元、祖冲之、李括、李时珍等，他们不仅在各自的科学领域里有着重要的发明和发现，而且他们还将这些发明和