常用蒙特卡罗程序介绍2.ppt

第八章 蒙特卡罗方法应用程序介绍 蒙特卡罗方法应用软件的特点 常用的通用蒙特卡罗程序简介 MCNP程序输入的描述 例子 IMP 栅元重要性卡 格式： IMP:n x1 x2 … xi … xI n ＝ 中子为 N，光子为 P，电子为 E。 N, P 、 P, E 或 N, P , E 也是允许的， 如果它们的值相同。 xi ＝ 栅元 i 的重要性，i ＝1，2，…，I I ＝ 问题中的栅元总数。 缺省：在一个MODE N P 问题中，若省略了IMP:P卡，则所有栅元的光子重要性都置 1，除非其中子重要性为零，这时其光子重要性也为零。 栅元重要性卡用于输入各个栅元的重要性。 零重要性可以用来终止粒子的历史。 IMP 卡是必须的，除非用了WWN卡。 例： IMP:N 1 2 2M 0 1 20R 表示栅元1的中子重要性为1，栅元2为2，栅元3为4，栅元4为0，栅元5至25为1。此时如果一个粒子从栅元2进入栅元3，则粒子分裂为两个粒子，每个粒子权重为原来的一半。若粒子从栅元3进入栅元2，则以50％的概率进行俄国轮盘赌终止粒子的历史，如果未能终止，则粒子的权重加倍。 如果粒子进入“0”重要性栅元，则粒子被杀死。 如果粒子进入真空栅元，则即不分裂也不终止，然而，当粒子离开真空进入一个非真空栅元时，则要根据其重要性比值的增减进行分裂或俄国轮盘赌。 源定义 对于点探测器或DXTRAN球还需定义其它一些变量。 SDEF 通用源卡 格式： SDEF 源变量=描述 等号“=”是可选的。 源变量和MCNP源必须设置的变量不完全一样，有许多是用来控制最终变量抽样值的中间量。所有源变量都有缺省值。源变量的描述有以下三种格式： 显值； 在一个分布号前加D； 在一个变量前面加F，后面跟一个D在前面的分布号。 Var=Dn 表示从分布 n 中抽取源变量Var。 Var Fvar Dn 则表示从分布 n 中抽取源变量Var，而分布 n 与变量 var 有关，只允许一级相关。 每个分布仅用于一个源变量。 源变量说明 源变量说明（续一） 源变量说明（续二） 源变量说明（续三） SIn 源信息卡 格式： SIn 选项 I1 … Ik n ＝ 分布号（n ＝1～999）。 选项＝ Ii 的说明，允许的值是： 省略或H：直方图分布的区间边界，仅用于标量。 L：离散的源变量值。 A：定义概率分布密度的点。 S：分布号。 I1 … Ik＝ 源变量值或分布号。 缺省： SIn H I1 … Ik SPn 源概率卡 格式： SPn 选项 P1 … Pk 或： SPn f a b n ＝ 分布号（n ＝1～999）。 选项＝ Pi 的说明，允许的值是： 省略：对H或L分布与D相同， 对A分布为概率密度。 D： H或L分布的各个区间的概率值(不用归一)。 C： H或L分布的各个区间的累积概率值。 V：仅用于栅元分布。概率与栅元体积成 正比（若有Pi ，要乘以Pi ）。 P1 … Pk＝ 源变量概率。 f ＝ 内部函数的标识符。 a，b ＝ 内部函数的输入参数。 缺省： SPn D P1 … Pk SBn 源偏移卡 格式： SBn 选项 B1 … Bk 或： SBn f a b n，选项， f ，a，b 均和SPn卡一样，只有一点不同，即内部函数 f 只允许－21和－31。 B1 … Bk＝ 源变量偏移概率。 缺省： SBn D B1 … Bk SPn卡的第一种格式的第一项是正数或非数值，表明该卡及它的SIn卡定义了一个概率分布函数。 SPn卡的第二种格式的第一项是负数，表示使用一个内部解析函数产生源变量的一个连续的概率密度函数。 内部函数说明 内部函数说明（续） DSn 相关源分布卡（复合分布） 格式： DSn 选项 J1 … Jk 或： DSn T I1 J1 … Ik Jk 或： DSn Q V1 S1 … Vk Sk n ＝ 分布号（n ＝1～999）。 选项＝ Ji 的说明，允许的值是： 省略或 H：连续分布的源变量值，仅用于标量。 L：离散的源变量值。 S：分布号。 T ＝ 独立变量值后面跟相关变量的值，