请教曲面重复的问题

Evan

这几天用MCNP时老是被这样的错误困扰：

fatal error.  surface      n   appears more than once in a chain.

surface      n   is in cells     A and     B in chain
     A  <    B

如果栅元A被栅元B包含，即A<B的关系。那么这两个栅元如果同时出现曲面n，就有可能警告甚至报错，如果曲面n是圆柱面，那么几乎一定会报错。

在很多情况下，曲面n对于A是可以省略的，直接把曲面A删掉就可以了，这还没事。

但是有些情况曲面n对A是不可以忽略的，同样对B也是不可忽略的，这样就比较麻烦了，不知道如何解决。比如说A被B包含的同时，也被C包含，而且C和B的形状不一样，那么A的曲面n就不能忽略。

请问版上的高手，有什么好的解决方法，我给了一个作为例子的代码，如下所示:
test
1 1 -1.0 -1 2 -3 u=1
2 0 #1 u=1
3 0 -6 4 -5 fill=1 trcl=(0 0 -4)
4 0 -1 4 -5 fill=1 trcl=(0 0 4)
5 0 #3 #4

1 cz 2
2 pz -1
3 pz 1
4 pz -2
5 pz 2
6 cz 3

mode n p e
imp:n 1 3r 0
imp:p 1 3r 0
imp:e 1 3r 0
M1 1001 2 8016 1
补充一下：其实我的输入文件是用程序（MCNP-MCL）自动生成的，所以很多场合下会有曲面重复的问题。

目前还不知道怎么从程序上解决，现阶段只能去改模型而避免这个现象，会造成一些不方便。

附件里是我出现这个错误的文件（不涉密，是一个MCNP-MCL的例程）。


在楼主的例子中：
问题出在了surface 1 上，栅元3与栅元1共用了曲面1，但是栅元1是填充栅元1的
正如zoid所说，填充的物质体积要比填充空间大一些
新添1个suface 7 cz 2.1
把cell 1改为：1 1 -1.0 -7 2 -3 u=1 就不会有问题。




转自52mc