ChinaNet

标题: 请教曲面重复的问题 [打印本页]

作者: Evan    时间: 2015-4-8 10:30
标题: 请教曲面重复的问题
这几天用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





欢迎光临 ChinaNet (http://www.nuclear.net.cn/) Powered by Discuz! X3.1