Commit 4e67b467 authored by Poarch, Sandra's avatar Poarch, Sandra
Browse files

Upload New File

parent b6a79e97
Loading
Loading
Loading
Loading
+492 −0
Original line number Diff line number Diff line
=csas6
Megapower Model "A" core - full core
ce_v7.1
' Megapower - design A
read comp
' HP vapor region (K)
  K-39  10 0 5.4191E-06 925.0 end
  K-40  10 0 6.7900E-10 925.0 end
  K-41  10 0 3.9109E-07 925.0 end
' HP liquid region
  K-39  11 0 9.8657E-03 925.0 end
  K-40  11 0 1.2377E-06 925.0 end
  K-41  11 0 7.1198E-04 925.0 end
' Clad
' 725 C; heat pipe clad
    Fe-54  15 0  3.2697E-03 925.0 end
    Fe-56  15 0  5.1327E-02 925.0 end
    Fe-57  15 0  1.1854E-03 925.0 end
    Fe-58  15 0  1.5775E-04 925.0 end
    Cr-50  15 0  6.7739E-04 925.0 end
    Cr-52  15 0  1.3063E-02 925.0 end
    Cr-53  15 0  1.4812E-03 925.0 end
    Cr-54  15 0  3.6870E-04 925.0 end
    Ni-58  15 0  6.6375E-03 925.0 end
    Ni-60  15 0  2.5568E-03 925.0 end
    Ni-61  15 0  1.1114E-04 925.0 end
    Ni-62  15 0  3.5436E-04 925.0 end
    Ni-64  15 0  9.0246E-05 925.0 end
    Mo-92  15 0  1.8402E-04 925.0 end
    Mo-94  15 0  1.1470E-04 925.0 end
    Mo-95  15 0  1.9741E-04 925.0 end
    Mo-96  15 0  2.0683E-04 925.0 end
    Mo-97  15 0  1.1842E-04 925.0 end
    Mo-98  15 0  2.9921E-04 925.0 end
    Mo-100 15 0  1.1941E-04 925.0 end
    Mn-55  15 0  1.7400E-03 925.0 end
    Si-28  15 0  1.5679E-03 925.0 end
    Si-29  15 0  7.9614E-05 925.0 end
    Si-30  15 0  5.2482E-05 925.0 end
' 725 C, inner clad
    Fe-54  16 0  3.2697E-03 925.0 end
    Fe-56  16 0  5.1327E-02 925.0 end
    Fe-57  16 0  1.1854E-03 925.0 end
    Fe-58  16 0  1.5775E-04 925.0 end
    Cr-50  16 0  6.7739E-04 925.0 end
    Cr-52  16 0  1.3063E-02 925.0 end
    Cr-53  16 0  1.4812E-03 925.0 end
    Cr-54  16 0  3.6870E-04 925.0 end
    Ni-58  16 0  6.6375E-03 925.0 end
    Ni-60  16 0  2.5568E-03 925.0 end
    Ni-61  16 0  1.1114E-04 925.0 end
    Ni-62  16 0  3.5436E-04 925.0 end
    Ni-64  16 0  9.0246E-05 925.0 end
    Mo-92  16 0  1.8402E-04 925.0 end
    Mo-94  16 0  1.1470E-04 925.0 end
    Mo-95  16 0  1.9741E-04 925.0 end
    Mo-96  16 0  2.0683E-04 925.0 end
    Mo-97  16 0  1.1842E-04 925.0 end
    Mo-98  16 0  2.9921E-04 925.0 end
    Mo-100 16 0  1.1941E-04 925.0 end
    Mn-55  16 0  1.7400E-03 925.0 end
    Si-28  16 0  1.5679E-03 925.0 end
    Si-29  16 0  7.9614E-05 925.0 end
    Si-30  16 0  5.2482E-05 925.0 end
' 750 C, outer clad
    Fe-54  17 0  3.2697E-03 1000.0 end
    Fe-56  17 0  5.1327E-02 1000.0 end
    Fe-57  17 0  1.1854E-03 1000.0 end
    Fe-58  17 0  1.5775E-04 1000.0 end
    Cr-50  17 0  6.7739E-04 1000.0 end
    Cr-52  17 0  1.3063E-02 1000.0 end
    Cr-53  17 0  1.4812E-03 1000.0 end
    Cr-54  17 0  3.6870E-04 1000.0 end
    Ni-58  17 0  6.6375E-03 1000.0 end
    Ni-60  17 0  2.5568E-03 1000.0 end
    Ni-61  17 0  1.1114E-04 1000.0 end
    Ni-62  17 0  3.5436E-04 1000.0 end
    Ni-64  17 0  9.0246E-05 1000.0 end
    Mo-92  17 0  1.8402E-04 1000.0 end
    Mo-94  17 0  1.1470E-04 1000.0 end
    Mo-95  17 0  1.9741E-04 1000.0 end
    Mo-96  17 0  2.0683E-04 1000.0 end
    Mo-97  17 0  1.1842E-04 1000.0 end
    Mo-98  17 0  2.9921E-04 1000.0 end
    Mo-100 17 0  1.1941E-04 1000.0 end
    Mn-55  17 0  1.7400E-03 1000.0 end
    Si-28  17 0  1.5679E-03 1000.0 end
    Si-29  17 0  7.9614E-05 1000.0 end
    Si-30  17 0  5.2482E-05 1000.0 end
' Gap regions
  he 20 0 1.0E-6  925.0 end
  he 21 0 1.0E-6  925.0 end
  he 22 0 1.0E-6  925.0 end
  he 23 0 1.0E-6  925.0 end
' Fuel
  uo2 40 den=10.4 1.0 1000.0 92234 0.15 92235 19.95 92238 79.9 end
'  Alumina (Al2O3) reflector
  Al2O3    30 DEN=3.9              end
' Control drums
  B4C 55 den=2.51 1.0 925.0  5010 90.0 5011 10.0 end
' BeO axial reflectors
  BeO 60 den=3.01 end
' Building materials
  orconcrete 2 1.0  293.0 end
  aluminum   3 1.0  293.0 end
  ss304      4 0.3  293.0 end
  dry-air    5 1.0  293.0 end
' Earth
 wtptEarth   6   1.52  10
                 8000 51.371535
                11000  0.613996
                12000  1.330392
                13000  6.856351
                14000 27.118001
                19000  1.432657
                20000  5.116630
                22000  0.460436
                25000  0.071548
                26000  5.628454
                    1 293.0 end
' Pump 80% carbon steel, 15% copper, and 5% cobalt by volume
  carbonsteel 100 0.80 293.0  end
  copper      100 0.15 293.0  end
  cobalt      100 0.05 293.0  end
' Pump shaft
  ss304       101 1.0  293.0 end
' Treasure chest
  redwood     200 1.0  293.0 end
  gold        201 1.0  293.0 end

end comp

read param
  sig=0.00005
  gen=530
  npg=100000
  nbk=125000
  nsk=30
  flx=yes
  htm=no
  cds=1
  gfx=2
  ngp=200
end param

read bounds  all=void  end bounds

read geometry
unit 10
com="Fuel unit cell"
' Heat pipe region
zcylinder 1 0.7100 2p100.0   com='K vapor'
zcylinder 2 0.7875 2p100.0   com='K liquid'
zcylinder 3 0.8875 2p100.0   com='HP clad'
zcylinder 4 0.8939 2p100.0   com='HP gap'
media 10 1  1                com='K vapor'
media 11 1  2 -1             com='K liquid'
media 15 1  3 -2             com='Clad'
media 20 1  4 -3             com='HP gap'

zcylinder 10 0.9339 2p100.0  com="Fuel clad (inner)"
zcylinder 11 0.9403 2p100.0  com="Fuel gap (inner)"
hexprism  12 1.2802 2p100.0  com="UO2 fuel"
hexprism  13 1.2866 2p100.0  com="Fuel gap (outer)"
hexprism  14 1.3866 2p100.0  com="Fuel clad (outer)"
hexprism  15 1.3930 2p100.0  com="Outer gap"

plane  20 zpl=1.0 con=-85.0 com="Top BeO"
plane  21 zpl=1.0 con=65.0  com="Bottom BeO"
plane  22 zpl=1.0 con=80.0  com="FP void bound"

media  16 1  10 -4  -20 21  com="Fuel clad (inner)"
media  21 1  11 -10 -20 21  com="Fuel gap (inner)"
media  40 1  12 -11 -20 21  com="UO2 fuel"
media  22 1  13 -12 -20 21  com="Fuel gap (outer)"
media  17 1  14 -13 -20 21  com="Fuel clad (outer)"
media  23 1  15 -14         com="Outer gap"

media  60 1  14  20 -4      com="BeO reflector (top)"
media  60 1  14 -21 22 -4   com="BeO reflector (bottom)"
media  23 1  14 -22 -4      com="FP plenum"
boundary 15

unit 11
com="Inner dummy cell"
hexprism  1 1.393  2p100.0
media 30 1 1
boundary 1

unit 12
com="Inner core void region"
zcylinder 1 9.05 2p100.0
hexprism 2 9.551 2p100.0
media 20 1 1
media 30 1 2 -1
boundary 2

unit 15
com="Radial control drums"
zcylinder 1  11.5 2p100.0 origin x=-1.0
zcylinder 2  12.5 2p100.0 origin x=0.0
media 30 1 1
media 55 1 2 -1
boundary 2


unit 100
hexprism  40 49.7 2p100.0  com="Core boundary"
hexprism  41 50.1 2p100.0  com="Outer core gap"
zcylinder 42 77.85 2p100.0 com="Alumina outer reflector"
' RPV
zcylinder 43 82.93 2p100.0 com="RPV core barrel"
' Neutron shield
zcylinder 44 98.17 2p100.0001 com="B4C Neutron shield"
' Fuel region
array 1 40 place 22 25 1 0.0 0.0 0.0
' Coolant channel
hole 12  origin x=0.0 y=0.0 z=0.0

' Control drums
hole 15 origin x=63.1 y=+14.35 com="CD #1, 0 degrees"
hole 15 origin x=63.1 y=-14.35 com="CD #2, 0 degrees"
hole 15 origin x=-63.1 y=+14.35 rotate a1=180 com="CD #1, 180 degrees"
hole 15 origin x=-63.1 y=-14.35 rotate a1=180 com="CD #2, 180 degrees"

hole 15 origin x=18.9926 y=+61.8962 rotate a1=60 com="CD #1, 60 degrees"
hole 15 origin x=44.1074 y=+47.3962 rotate a1=60 com="CD #2, 60 degrees"

hole 15 origin x=-44.1074 y=+47.3962 rotate a1=120 com="CD #1, 120 degrees"
hole 15 origin x=-18.9926 y=+61.8962 rotate a1=120 com="CD #2, 120 degrees"

hole 15 origin x=-18.9926 y=-61.8962 rotate a1=240 com="CD #1, 240 degrees"
hole 15 origin x=-44.1074 y=-47.3962 rotate a1=240 com="CD #2, 240 degrees"

hole 15 origin x=+44.1074 y=-47.3962 rotate a1=300 com="CD #1, 300 degrees"
hole 15 origin x=+18.9926 y=-61.8962 rotate a1=300 com="CD #2, 300 degrees"

media 20 1 41 -40
media 30 1 42 -41
media 17 1 43 -42
media 55 1 44 -43
boundary 44

unit 200
com="internal pump"
xcylinder 1 20.0 140.0  0.0
xcylinder 2 20.0  40.0  0.0
xcylinder 3 5.0  140.0 5.0
media  100 1  2 -3 vol=47516.6
media  101 1  3 vol=10602.9
media    2 1  1 -2 -3
boundary 1

unit 300
com="hallway pump"
xcylinder 1 20.0 140.0  0.0
xcylinder 2 20.0  40.0  0.0
xcylinder 3 5.0  140.0 5.0
media  100 1  2 -3 vol=47516.6
media  101 1  3 vol=10602.9
media    2 1  1 -2 -3
boundary 1

unit 400
com="control room pump"
xcylinder 1 20.0 140.0  0.0
xcylinder 2 20.0  40.0  0.0
xcylinder 3 5.0  140.0  5.0
xcylinder 4 20.0 140.0 100.0
media  100 1  2 -3 vol=47516.6
media  101 1  3 vol=10602.9
media    2 1  4 -3
media    5 1  1 -2 -3 -4
boundary 1

unit 500
com="external pump"
xcylinder 1 20.0 140.0  0.0
xcylinder 2 20.0  40.0  0.0
xcylinder 3 5.0  140.0  5.0
xcylinder 4 20.0 140.0 100.0
media  100 1  2 -3 vol=47516.6
media  101 1  3 vol=10602.9
media    2 1  4 -3
media    5 1  1 -2 -3 -4
boundary 1

unit 600
com="roof pump"
xcylinder 1 20.0 140.0  0.0
xcylinder 2 20.0  40.0  0.0
xcylinder 3 5.0  140.0 5.0
media  100 1  2 -3 vol=47516.6
media  101 1  3 vol=10602.9
media    2 1  1 -2 -3
boundary 1

unit 700
com="control panel"
cuboid   1  90 -90   100  0   180  0
cuboid   2  90 -90   100  5   180 80
cuboid   3  85 -85    95  5    75  5
cuboid   4  85 -85   100  5    75  0
com="screens and keybord"
cuboid  10  30 -30    35 30   130 95 origin x=38.125 rotate a1=15
cuboid  11  30 -30    35 30   130 95 origin x=-38.125 rotate a1=-15
cuboid  12  30 -30    10  0    2 0 origin y=50 z=81.75 rotate a2=-10
media  3 1 1 -2 -3 -4
media  5 1 2 -10 -11 -12
media  5 1 3 -4
media  5 1 4
media  2 1 10
media  2 1 11
media  2 1 12
boundary 1

unit 777
com="treasure chest"
cuboid 1 2p40.0 4p30.0
cuboid 2 2p35.0 4p25.0
cuboid 3 2p35.0 2p25.0 2p10.0 origin z=-15
media  200 1  1 -2
media    5 1  2 -3
media  201 1  3 vol=70000.0
boundary 1

global unit 1
com="site"
cuboid   1 1700 -500   1100 -500     850.0   -200.0
cuboid   2 1700 -500   1100 -500       0.0   -200.0

com="exterior of the building, roof, floor"
cuboid   10 1200 0    600 0   400.0  -100.0
cuboid   11 1200 0    600 0   400.0   300.0
cuboid   12 1200 0    600 0     0.0  -100.0

com="air space in building - two rooms and maze"
cuboid   20 1160 100   500 100    300 0
cuboid   21 1160 600   560  40    300 0

com="interior walls to form maze to prevent streaming"
cuboid   31  600 500  500 100   300 0
cuboid   32  780 740  560  40   300 0

com="first door"
cuboid   40  600 500  180  100   210 0
cuboid   41  600 595  180  100   210 0

com="second door"
cuboid   50  780 740  560 480   210 0
cuboid   51  780 775  560 480   210 0

com="exterior door"
cuboid   60 1200 1160  360 240   210 0
cuboid   61 1200 1195  360 240   210 0

com="walls, roof, floor"
media  2 1   10 -20 -11 -12 -21 -40 -50 -60
media  2 1   11
media  2 1   12

com="air space and maze walls"
media  5 1   20 -21 -31 -32 -11 -12
media  5 1   21 -31 -32     -11 -12
media  2 1   31 -40         -11 -12
media  2 1   32 -50         -11 -12

com="first door"
media  5 1   40 -41 -12
media  4 1   41 -12

com="second door"
media  5 1   50 -51 -12
media  4 1   51 -12

com="exterior door"
media  5 1   60 -61 -12
media  4 1   61 -12

com="reactor"
hole 100 origin x=280 y=320 z=100 rotate a2=90

com="interior pump"
hole 200 origin x=460 y=240 z=20

com="hallway pump"
hole 300 origin x=640 y=440 z=20 rotate a1=180

com="control room pump"
hole 400 origin x=880 y=300 z=20 rotate a1=180

com="external pump"
hole 500 origin x=1100 y=700 z=20 rotate a1=-90

com="roof pump"
hole 600 origin x=1100 y=300 z=440 rotate a2=90 a3=-90

com="control panel"
hole 700 origin x=1040 y=45 z=60

com="treasure chest"
hole 777 origin x=280 y=320 z=-150


com="site soil and air"
media 5 1 1 -2 -10
media 6 1 2 -10

boundary 1

end geometry
' Filler cells needs to be alumina
' Net of 1134 heat pipe (fuel) cells

read array
 ara=1 typ=shexagonal
  nux=43 nuy=49 nuz=1 com="Core lattice"
  fill
  43r11
  43r11
  20r11 3r10    20r11
  18r11 6r10    19r11
  17r11 9r10    17r11
  15r11 12r10   16r11
  14r11 15r10   14r11
  12r11 18r10   13r11
  11r11 21r10   11r11
  9r11  24r10   10r11
  8r11  27r10    8r11
  6r11  30r10    7r11
  5r11  33r10    5r11
  4r11  34r10    5r11
  4r11  35r10    4r11
  4r11  34r10    5r11
  4r11  35r10    4r11
  4r11  34r10    5r11
  4r11  35r10    4r11
  4r11  34r10    5r11
  4r11  16r10 3r11 16r10  4r11
  4r11  14r10 6r11 14r10  5r11
  4r11  14r10 7r11 14r10  4r11
  4r11  13r10 8r11 13r10  5r11
  4r11  14r10 7r11 14r10  4r11
  4r11  13r10 8r11 13r10  5r11
  4r11  14r10 7r11 14r10  4r11
  4r11  14r10 6r11 14r10  5r11
  4r11  16r10 3r11 16r10  4r11
  4r11  34r10    5r11
  4r11  35r10    4r11
  4r11  34r10    5r11
  4r11  35r10    4r11
  4r11  34r10    5r11
  4r11  35r10    4r11
  4r11  34r10    5r11
  5r11  33r10    5r11
  6r11  30r10    7r11
  8r11  27r10    8r11
  9r11  24r10   10r11
  11r11 21r10   11r11
  12r11 18r10   13r11
  14r11 15r10   14r11
  15r11 12r10   16r11
  17r11 9r10    17r11
  18r11 6r10    19r11
  20r11 3r10    20r11
  43r11
  43r11
 end fill
end array
read grid 1
  xlinear 20  330.0 230.0
  zlinear 30  155.0  45.0
  ylinear 22  385.0 235.0
end grid
read grid 2
  xlinear 40 1700.0 -500.0
  ylinear 40 1100.0 -500.0
  zlinear 40  850.0 -150.0
end grid
read start
  nst=0
  xsp=330.0 xsm=230.0
  ysp=385.0 ysm=235.0
  zsp=155.0 zsm=45.0
end start
'read volume
'  type=random
'  batches=100
'  points=10000
'end volume
end data
end