Solve MCNP Input File for AP1000: Get Help Now

In summary, the problem is that the MCNP code does not recognize the lattice shape for the material in 1/8 AP1000.
  • #1
zaidtaher
14
0
Hello Every body,

I hope all is well.

I have a problem with MCNP code for 1/8 AP1000, so please can anybody help me ?

Regrads
 
Engineering news on Phys.org
  • #2
Are we supposed to guess what the problem is ?
 
  • #3
haha,, no I just want to see if I can have a help.
the problem that I repeat a lattice but VISED didn't fill all of it in materials, and there is some red lines I don't know why.
 
  • #4
1 1 -10.97 -1 -8 9 u=2 $FUEL Type 1
2 2 -0.004 1 -2 -8 9 u=2 $GAP
3 3 -6.7 2 -3 -8 9 u=2 $CLAD
4 4 -1 3 -4 5 -6 7 -8 9 u=2 $MODIRATOR
5 0 4 :-5 :6 :-7 :8 :-9 u=2
6 0 -10 11 -12 13 -8 9 fill=1 u=10
7 0 10 :-11 :12 :-13 :8 :-9
8 0 -4 5 -6 7 u=1 lat=1 $ROW 1
fill=-8:8 -8:8 0:0
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 1
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 2
2 2 2 2 2 0 2 2 0 2 2 0 2 2 2 2 2 $ROW 3
2 2 2 0 2 2 2 2 2 2 2 2 2 0 2 2 2 $ROW 4
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 5
2 2 0 2 2 0 2 2 0 2 2 0 2 2 0 2 2 $ROW 6
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 7
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 8
2 2 0 2 2 0 2 2 0 2 2 0 2 2 0 2 2 $ROW 9
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 10
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 11
2 2 0 2 2 0 2 2 0 2 2 0 2 2 0 2 2 $ROW 12
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 13
2 2 2 0 2 2 2 2 2 2 2 2 2 0 2 2 2 $ROW 14
2 2 2 2 2 0 2 2 0 2 2 0 2 2 2 2 2 $ROW 15
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 16
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 17
9 7 -10.97 -1 -8 9 u=4 $FUEL Type 2
10 2 -0.004 1 -2 -8 9 u=4 $GAP
11 3 -6.7 2 -3 -8 9 u=4 $CLAD
12 4 -1 3 -4 5 -6 7 -8 9 u=4 $MODIRATOR
13 0 4 :-5 :6 :-7 :8 :-9 u=4
14 0 -10 11 -12 13 -8 9 fill=3 u=11
15 0 -4 5 -6 7 u=3 lat=1 $ROW 1
fill=-8:8 -8:8 0:0
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 1
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 2
4 4 4 4 4 0 4 4 0 4 4 0 4 4 4 4 4 $ROW 3
4 4 4 0 4 4 4 4 4 4 4 4 4 0 4 4 4 $ROW 4
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 5
4 4 0 4 4 0 4 4 0 4 4 0 4 4 0 4 4 $ROW 6
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 7
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 8
4 4 0 4 4 0 4 4 0 4 4 0 4 4 0 4 4 $ROW 9
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 10
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 11
4 4 0 4 4 0 4 4 0 4 4 0 4 4 0 4 4 $ROW 12
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 13
4 4 4 0 4 4 4 4 4 4 4 4 4 0 4 4 4 $ROW 14
4 4 4 4 4 0 4 4 0 4 4 0 4 4 4 4 4 $ROW 15
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 16
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 17
16 8 -10.97 -1 -8 9 u=6 $FUEL Type 3
17 2 -0.004 1 -2 -8 9 u=6 $GAP
18 3 -6.7 2 -3 -8 9 u=6 $CLAD
19 4 -1 3 -4 5 -6 7 -8 9 u=6 $MODIRATOR
20 0 4 :-5 :6 :-7 :8 :-9 u=6
21 0 -10 11 -12 13 -8 9 fill=5 u=12
23 0 -4 5 -6 7 u=5 lat=1 $ROW 1
fill=-8:8 -8:8 0:0
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 1
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 2
6 6 6 6 6 0 6 6 0 6 6 0 6 6 6 6 6 $ROW 3
6 6 6 0 6 6 6 6 6 6 6 6 6 0 6 6 6 $ROW 4
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 5
6 6 0 6 6 0 6 6 0 6 6 0 6 6 0 6 6 $ROW 6
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 7
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 8
6 6 0 6 6 0 6 6 0 6 6 0 6 6 0 6 6 $ROW 9
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 10
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 11
6 6 0 6 6 0 6 6 0 6 6 0 6 6 0 6 6 $ROW 12
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 13
6 6 6 0 6 6 6 6 6 6 6 6 6 0 6 6 6 $ROW 14
6 6 6 6 6 0 6 6 0 6 6 0 6 6 6 6 6 $ROW 15
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 16
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 17
24 6 -5 -1 -8 9 u=7 $Control Material(black)
25 2 -0.004 1 -2 -8 9 u=7 $Gap
26 3 -6.7 2 -18 -8 9 u=7 $Clad
27 4 -1 18 -4 5 -6 7 -8 9 u=7 $MODIRATOR
28 0 -10 11 -12 13 -8 9 u=9 lat=1 $ROW 1
fill=-4:4 -4:4 0:0
10 11 10 11 10 11 10 12 0 $ROW 1
0 10 11 10 11 10 12 12 0 $ROW 2
0 0 10 11 10 11 12 0 0 $ROW 3
0 0 0 10 11 12 12 0 0 $ROW 4
0 0 0 0 12 12 0 0 0 $ROW 5
0 0 0 0 0 0 0 0 0 $ROW 6
0 0 0 0 0 0 0 0 0 $ROW 7
0 0 0 0 0 0 0 0 0 $ROW 8
0 0 0 0 0 0 0 0 0 $ROW 9
29 0 -19 20 -21 22 -8 9 fill=9
30 0 19 :-20 :21 :-22 :8 :-9

1 cz 0.4095
2 cz 0.41775
3 cz 0.47475
4 px 0.63
5 px -0.63
6 py 0.63
7 py -0.63
8 pz 215
9 pz -215
10 px 10.71
11 px -10.71
12 py 10.71
13 py -10.71
18 cz 0.43625
19 px 96.39
20 px -96.39
21 py 96.39
22 py -96.39

mode n
kcode 500 1.000000 50 1500
ksrc 0.000000 0.000000 0.000000
m1 92238. 0.855 $type 1
92235. 0.045 8016. 0.1
m2 8016. 1 $He2
m3 40090.70c 0.038189 $Zirlo
40091.70c 0.008328 40092.70c 0.01273 40094.70c 0.0129
40096.70c 0.002078 50112.70c 7e-006 50114.70c 5e-006
50115.70c 3e-006 50116.70c 0.00011 50117.70c 5.8e-005
50118.70c 0.000183 50119.70c 6.5e-005 50120.70c 0.000247
50122.70c 3.5e-005 50124.70c 4.4e-005 41093.70c 0.000757
40090.71c 0.475264 40091.71c 0.101628 40092.71c 0.15534
40094.71c 0.157424 40096.71c 0.025362 50112.71c 9e-005
50114.71c 6.1e-005 50115.71c 3.1e-005 50116.71c 0.001344
50117.71c 0.00071 50118.71c 0.002239 50119.71c 0.000794
50120.71c 0.003011 50122.71c 0.000428 50124.71c 0.000535
41093.71c 0.009243
m4 8016. 0.66 $H2O
1001. 0.66
m5 6000.60c -0.0003 $Structure Material
14000.60c -0.005 15031.60c -0.000225 16000.60c -0.00015
24050.60c -0.00793 24052.60c -0.159031 24053.60c -0.018378
24054.60c -0.004661 25055.60c -0.01 26054.60c -0.039996
26056.60c -0.644764 26057.60c -0.015026 26058.60c -0.002039
28058.60c -0.06234 28060.60c -0.024654 28061.60c -0.001085
28062.60c -0.003504 28064.60c -0.000917
m6 5010. 0.8 $Control Material
6000. 0.2
m7 92235.66c 0.11 $type 2
92238.66c 0.22 8016.62c 0.67
m8 92235.66c 0.07755 $type 3
92238.66c 0.25245 8016.62c 0.67
imp:n 1 14r 0 1 11r $ 2, 30
f5:n 0 0 0 0.1
 
  • #5
zaidtaher,

You need to leave 5 position before each continuation cards,
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 1
change for
(5 spaces) 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 1
If you add the spaces in all the lattices and material continuation cards it runs.

Hernán
 
  • #6
I tried it, but it's already 6 spaces, but i loosed the format when I put it here, btw it's does not work
Thanks for help :)
 
  • #7
Zaidtaher,

You have a problem with importance definition, I just specify the importance for each cell, that give me no fatal error. You still have to check, volumes and materials.

AP-1000 Core
c -----------------------CELL CARD---------------------------------------------
c
1 1 -10.97 -1 -8 9 u=2 imp:n=1 $FUEL Type 1
2 2 -0.004 1 -2 -8 9 u=2 imp:n=1 $GAP
3 3 -6.7 2 -3 -8 9 u=2 imp:n=1 $CLAD
4 4 -1 3 -4 5 -6 7 -8 9 u=2 imp:n=1 $MODIRATOR
5 0 4 :-5 :6 :-7 :8 :-9 u=2 imp:n=1
6 0 -10 11 -12 13 -8 9 fill=1 u=10 imp:n=1
7 0 10 :-11 :12 :-13 :8 :-9 imp:n=1
8 0 -4 5 -6 7 u=1 lat=1 $ROW 1
fill=-8:8 -8:8 0:0
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 1
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 2
2 2 2 2 2 0 2 2 0 2 2 0 2 2 2 2 2 $ROW 3
2 2 2 0 2 2 2 2 2 2 2 2 2 0 2 2 2 $ROW 4
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 5
2 2 0 2 2 0 2 2 0 2 2 0 2 2 0 2 2 $ROW 6
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 7
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 8
2 2 0 2 2 0 2 2 0 2 2 0 2 2 0 2 2 $ROW 9
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 10
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 11
2 2 0 2 2 0 2 2 0 2 2 0 2 2 0 2 2 $ROW 12
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 13
2 2 2 0 2 2 2 2 2 2 2 2 2 0 2 2 2 $ROW 14
2 2 2 2 2 0 2 2 0 2 2 0 2 2 2 2 2 $ROW 15
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 16
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 $ROW 17
imp:n=1
9 7 -10.97 -1 -8 9 u=4 imp:n=1 $FUEL Type 2
10 2 -0.004 1 -2 -8 9 u=4 imp:n=1 $GAP
11 3 -6.7 2 -3 -8 9 u=4 imp:n=1 $CLAD
12 4 -1 3 -4 5 -6 7 -8 9 u=4 imp:n=1 $MODIRATOR
13 0 4 :-5 :6 :-7 :8 :-9 u=4 imp:n=1
14 0 -10 11 -12 13 -8 9 fill=3 u=11 imp:n=1
15 0 -4 5 -6 7 u=3 lat=1 $ROW 1
fill=-8:8 -8:8 0:0
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 1
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 2
4 4 4 4 4 0 4 4 0 4 4 0 4 4 4 4 4 $ROW 3
4 4 4 0 4 4 4 4 4 4 4 4 4 0 4 4 4 $ROW 4
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 5
4 4 0 4 4 0 4 4 0 4 4 0 4 4 0 4 4 $ROW 6
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 7
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 8
4 4 0 4 4 0 4 4 0 4 4 0 4 4 0 4 4 $ROW 9
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 10
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 11
4 4 0 4 4 0 4 4 0 4 4 0 4 4 0 4 4 $ROW 12
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 13
4 4 4 0 4 4 4 4 4 4 4 4 4 0 4 4 4 $ROW 14
4 4 4 4 4 0 4 4 0 4 4 0 4 4 4 4 4 $ROW 15
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 16
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 $ROW 17
imp:n=1
16 8 -10.97 -1 -8 9 u=6 imp:n=1 $FUEL Type 3
17 2 -0.004 1 -2 -8 9 u=6 imp:n=1 $GAP
18 3 -6.7 2 -3 -8 9 u=6 imp:n=1 $CLAD
19 4 -1 3 -4 5 -6 7 -8 9 u=6 imp:n=1 $MODIRATOR
20 0 4 :-5 :6 :-7 :8 :-9 u=6 imp:n=1
21 0 -10 11 -12 13 -8 9 fill=5 u=12 imp:n=1
23 0 -4 5 -6 7 u=5 lat=1 $ROW 1
fill=-8:8 -8:8 0:0
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 1
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 2
6 6 6 6 6 0 6 6 0 6 6 0 6 6 6 6 6 $ROW 3
6 6 6 0 6 6 6 6 6 6 6 6 6 0 6 6 6 $ROW 4
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 5
6 6 0 6 6 0 6 6 0 6 6 0 6 6 0 6 6 $ROW 6
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 7
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 8
6 6 0 6 6 0 6 6 0 6 6 0 6 6 0 6 6 $ROW 9
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 10
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 11
6 6 0 6 6 0 6 6 0 6 6 0 6 6 0 6 6 $ROW 12
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 13
6 6 6 0 6 6 6 6 6 6 6 6 6 0 6 6 6 $ROW 14
6 6 6 6 6 0 6 6 0 6 6 0 6 6 6 6 6 $ROW 15
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 16
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ROW 17
imp:n=1
24 6 -5 -1 -8 9 u=7 imp:n=1 $Control Material(black)
25 2 -0.004 1 -2 -8 9 u=7 imp:n=1 $Gap
26 3 -6.7 2 -18 -8 9 u=7 imp:n=1 $Clad
27 4 -1 18 -4 5 -6 7 -8 9 u=7 imp:n=1 $MODIRATOR
28 0 -10 11 -12 13 -8 9 u=9 lat=1 $ROW 1
fill=-4:4 -4:4 0:0
10 11 10 11 10 11 10 12 0 $ROW 1
0 10 11 10 11 10 12 12 0 $ROW 2
0 0 10 11 10 11 12 0 0 $ROW 3
0 0 0 10 11 12 12 0 0 $ROW 4
0 0 0 0 12 12 0 0 0 $ROW 5
0 0 0 0 0 0 0 0 0 $ROW 6
0 0 0 0 0 0 0 0 0 $ROW 7
0 0 0 0 0 0 0 0 0 $ROW 8
0 0 0 0 0 0 0 0 0 $ROW 9
imp:n=1
29 0 -19 20 -21 22 -8 9 fill=9 imp:n=1
30 0 19 :-20 :21 :-22 :8 :-9 imp:n=0

1 cz 0.4095
2 cz 0.41775
3 cz 0.47475
4 px 0.63
5 px -0.63
6 py 0.63
7 py -0.63
8 pz 215
9 pz -215
10 px 10.71
11 px -10.71
12 py 10.71
13 py -10.71
18 cz 0.43625
19 px 96.39
20 px -96.39
21 py 96.39
22 py -96.39

mode n
kcode 500 1.000000 50 1500
ksrc 0.000000 0.000000 0.000000
m1 92238. 0.855 $ type 1
92235. 0.045 8016. 0.1
m2 8016. 1 $ He2
m3 40090.70c 0.038189 $ Zirlo
40091.70c 0.008328 40092.70c 0.01273 40094.70c 0.0129
40096.70c 0.002078 50112.70c 7e-006 50114.70c 5e-006
50115.70c 3e-006 50116.70c 0.00011 50117.70c 5.8e-005
50118.70c 0.000183 50119.70c 6.5e-005 50120.70c 0.000247
50122.70c 3.5e-005 50124.70c 4.4e-005 41093.70c 0.000757
40090.71c 0.475264 40091.71c 0.101628 40092.71c 0.15534
40094.71c 0.157424 40096.71c 0.025362 50112.71c 9e-005
50114.71c 6.1e-005 50115.71c 3.1e-005 50116.71c 0.001344
50117.71c 0.00071 50118.71c 0.002239 50119.71c 0.000794
50120.71c 0.003011 50122.71c 0.000428 50124.71c 0.000535
41093.71c 0.009243
m4 8016. 0.66 $H2O
1001. 0.66
m5 6000.60c -0.0003 $ Structure Material
14000.60c -0.005 15031.60c -0.000225 16000.60c -0.00015
24050.60c -0.00793 24052.60c -0.159031 24053.60c -0.018378
24054.60c -0.004661 25055.60c -0.01 26054.60c -0.039996
26056.60c -0.644764 26057.60c -0.015026 26058.60c -0.002039
28058.60c -0.06234 28060.60c -0.024654 28061.60c -0.001085
28062.60c -0.003504 28064.60c -0.000917
m6 5010. 0.8 $ Control Material
6000. 0.2
m7 92235.66c 0.11 $ type 2
92238.66c 0.22 8016.62c 0.67
m8 92235.66c 0.07755 $ type 3
92238.66c 0.25245 8016.62c 0.67
f5:n 0 0 0 0.1
 
  • #8
Dear,

It didn't work either :/

Zaid
 
  • #9
if there a way to send the input file to you it's will be easier
 
  • #10
you just put it on VISED and will take the proper format.

Regards
Zaid
 

Attachments

  • Zaid MCNP 2.txt
    6.7 KB · Views: 652
  • #12
Thanks for the file, but the problem is that the VISED didn't color all of the assemblies, and so the code will not run
 
  • #13
you something I have been here for about month and I can't find what's wrong,, please help me
 
  • #14
Zaidther,

What you want us to check is not so easy, try to run the inner universes alone first and then then the filled ones and hopefully you will find where is the problem.

Hernán
 
  • Like
Likes 1 person
  • #15
But do u knw anybody can help me ? please
 
  • #16
please , I need help if u knw anyone ??
 

What is MCNP?

MCNP (Monte Carlo N-Particle) is a general-purpose Monte Carlo particle transport code used for simulating nuclear processes. It is widely used in the nuclear industry for design and analysis of nuclear systems.

What is an AP1000?

The AP1000 is a type of pressurized water reactor (PWR) used for electricity generation. It is known for its advanced safety features and simplified design, making it a popular choice for nuclear power plants.

What is an input file?

An input file is a text file that contains all the necessary information for running a simulation in MCNP. It includes details such as geometry, materials, particle sources, and simulation parameters.

How do I solve an MCNP input file for AP1000?

To solve an MCNP input file for AP1000, you will need to have MCNP installed on your computer. You can then open the input file in the MCNP software and run the simulation. You may need to adjust some parameters or make modifications to the input file based on your specific needs.

Where can I get help with solving an MCNP input file for AP1000?

There are many resources available for getting help with solving an MCNP input file for AP1000. You can consult the MCNP user manual, attend training courses or workshops, or seek assistance from experienced MCNP users or consultants.

Similar threads

Replies
8
Views
1K
Replies
19
Views
1K
Replies
3
Views
1K
Replies
5
Views
2K
Replies
6
Views
1K
  • Nuclear Engineering
Replies
2
Views
2K
Replies
2
Views
2K
  • Nuclear Engineering
Replies
5
Views
402
  • Nuclear Engineering
Replies
4
Views
2K
Back
Top