蒙特卡罗方法应用程序介绍课件.pptVIP

缺省值 MCNP的许多输入参数都有缺省值，因此用户不需要每次都给出各个输入参量的值。当缺省值符合用户要求时，便可不在输入文件中指定。当省略某张输入卡时，则该卡上的全部参数均使用缺省值。如果只想改变一张卡上的某一个特定参量时，则它前面的参量仍需指明，或者用nJ方式跳过前面那些使用缺省值的参量。 例如：光子截断卡 CUT:P 3J -.10 表示前3个参量使用缺省值，只改变第四项参量的值。 精品 输入错误信息 MCNP对输入文件出现的错误作广泛的检查，如果用户违反了输入说明的规定，将在终端上以及输出文件中打印致命错误信息，MCNP不再进行粒子输运计算，作业中断。 第一个出现的致命错误是真的，而后面的错误可能不一定是真的，这取决于前面出现的致命错误的情况。 若在MCNP运行行上指定FATAL项，则MCNP忽略致命错误，照常运行。 对于MCNP的警告信息，用户不应忽视，应搞清楚它们的含义。 精品 检查几何错误 MCNP在处理输入文件的数据时，不能检查一种非常重要的输入错误。即MCNP无法查出各栅元之间的重叠和空隙，只有当粒子丢失时，才会发现几何错误。即使如此，可能仍然无法准确判断错误性质。 精品 栅元描述卡 格式： j m d geom params 或： j LIKE n BUT list j 栅元号，1≤ j ≤99999，写在第 1～5 列上。 m 栅元材料号，与材料卡（Mm）中的序号对应。 m ＝0 为真空栅元。 d 栅元材料密度。正值为原子密度，负值为质量密度。对于真空栅元，该项缺省，不填写。 geom 栅元的几何说明。由一系列带符号的曲面号经过布尔运算组成。 params 任选的栅元参数说明。 n 另一个栅元的名字（编号）。 list 描述栅元j和栅元n之间差别的栅元参数。 精品 在栅元的几何说明中，关于曲面的指向是一个很重要的概念。假定曲面 S 的曲面方程为 f (x,y,z)＝0，则对于f (x,y,z)＞0的区域对于曲面 S 具有正的指向；而对于f (x,y,z)＜0的区域对于曲面 S 具有负的指向。正指向的区域用+S表示，“+”号可不写；负指向的区域用-S表示。栅元用各相关曲面的布尔运算表示，布尔算符包括交（用空格表示 ）、并（用冒号:表示）和非（用#表示）。缺省的运算顺序是先非，其次是交，最后是并 ，使用括号可控制布尔运算的次序。 精品 非运算有两种形式： #n，n是某个栅元号， #n表示一个由不在栅元n内的点组成的空间区域。 # ( ---)，括号内是对某一个栅元进行描述的曲面——栅元关系组，这一形式定义的几何区域由不属于括号内描述区域的点组成的空间。 例如：3 0 -1 2 -4 $ 定义栅元3 #3 $ 与下行相同 #(-1 2 -4) 精品 在栅元卡上可定义栅元参数以代替在输入文件中数据卡部分定义的栅元参数。格式为：关键词＝值。这儿允许的关键词是：带有粒子标识符的IMP、VOL、PWT、EXT、FCL、WWN、DXC、NONU、PD和TMP，以及关于重复结构的4个栅元参数卡：U卡、TRCL卡、LAT卡和FILL卡。 在LIKE n BUT格式中，还有两个关键词MAT和RHO，分别表示栅元的介质号和密度。 精品 例如：10 16 -4.2 1 -2 3 IMP:N=4 IMP:P=8 表示栅元10由曲面1的正面、曲面2的负面和 曲面3的正面的交集组成，填充质量密度为 4.2 克 / 厘米3 的16号材料。该栅元的中子重要性为4，光子重要性为8。 例如：2 3 -3.7 -1 IMP:N=2 IMP:P=4 3 LIKE 2 BUT TRCL=1 IMP:N=10 精品 曲面描述卡 由方程定义曲面 格式： j n a list j 曲面号，1≤ j ≤99999，写在第 1～5 列上。 如果曲面号前有*号，则该曲面为反射面。 n 对应坐标变换卡TRn，表示该曲面是在辅助坐标系下描述的，而该辅助坐标系与基本坐标系之间的关系由TRn卡给出。如果没有坐标变换，即曲面是在基本坐标系下描述的，则该项缺省。 a 曲面助记符。 list 曲面方程参数，1～10项，取决于曲面类型。 参见MCNP手册，表3.1。 精品 用点定义轴对称曲面 类型为X、Y或Z的曲面卡是用坐标点描述曲面而不是用方程系数描述。用这些卡描述的曲面必须是分别关于X、Y或Z轴对称的，并且如果该曲面是由多叶组成的，则指定的坐标点必须全