Commit a7bc9040 authored by Poarch, Sandra's avatar Poarch, Sandra
Browse files

Upload New File

parent 90df9e85

06_MAVRIC/reactor.mavric.roof.inp

0 → 100644
+677 −0
Original line number Diff line number Diff line 
=shell

  cp  ${INPDIR}/reactor.csas.fission_source.3dmap fission_source.3dmap

  cp  ${INPDIR}/reactor.csas.kenoNuBar.txt  kenoNuBar.txt

end



=mt2msm

'fission_source.3dmap'

1

-1

1

'fission_source.msm'

end



=mavric

Megapower Model "A" core - full core

v7.1-28n19g

' 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 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 definitions



    distribution 1

      title="U-235 Watt fission spectrum"

      special="wattSpectrum"

      parameters 0.988 2.249 end

    end distribution



    energyBounds 1

        title='SCALE 200-group energy structure'

        200n

    end energyBounds



    gridGeometry 1

        xlinear 50  330.0 230.0

        zlinear 60  155.0  45.0

        ylinear 44  385.0 235.0



        xlinear 5 100.0 0.0

        xlinear 20 100.0 500.0

        xlinear 10 600.0 500.0

        xlinear 4 780.0 740.0

        xlinear 4 1200.0 1160.0



        ylinear 5 100.0 0.0

        ylinear 10 600.0 500.0



        zlinear 40  400.0    0.0

        zlinear 10    0.0 -100.0

        zlinear 5  -100.0 -200.0



        xlinear 22 1700.0 -500.0

        ylinear 16 1100.0 -500.0

        zlinear 10  850.0 -200.0

    end gridGeometry



    gridGeometry 2

        xlinear 220 1700.0 -500.0

        ylinear 160 1100.0 -500.0

        zlinear 100  850.0 -200.0

    end gridGeometry



    gridGeometry 3

        xlinear 50  330.0 230.0

        zlinear 60  155.0  45.0

        ylinear 44  385.0 235.0



        xplanes -330.5 -229.5  229.5 330.5 end

        yplanes 234.5 385.5 end

        zplanes 44.5 155.5 end



        xlinear 20   523 40

        ylinear 20   560 40

        zlinear 20   300 0

    end gridGeometry



    response 1

        title="Unit response"

        neutron

        bounds 20e7 1.0e-5 end

        values 1.0 1.0 end

    end response



end definitions



read sources



    src 40

       title="Fission neutrons"

       meshSourceFile="fission_source.msm"

'    strength equals to fissions * nubar (1.5605e17*2.501056)

       fissions=1.5605e17

       mixture=40

       allowResampling

    end src



'    src 40

'       title="Watt fission spectrum"

'       neutron

'    strength equals to fissions * nubar (1.5605e17*2.501)

'       strength=3.902834e+17

'       yCylinder 55.0  385.0 235.0 origin x=280 z=100

'       edistributionID=1 mixture=40

'    end src



'    meshSourceSaver

'        gridgeometryid=3

'        nenergyboundsid=1

'        penergyboundsid=1

'        maketotal

'        filename="source.msm"

'    end meshSourceSaver



end sources



read tallies



    regionTally 201

        title="Interior pump rotor"

        unit=200 region=1

        neutron

        energyBoundsID=1

    end regionTally

    regionTally 202

        title="Interior pump shaft"

        unit=200 region=2

        neutron

        energyBoundsID=1

    end regionTally



    regionTally 301

        title="Hallway pump rotor"

        unit=300 region=1

        neutron

        energyBoundsID=1

    end regionTally

    regionTally 302

        title="Hallway pump shaft"

        unit=300 region=2

        neutron

        energyBoundsID=1

    end regionTally



    regionTally 401

        title="Control pump rotor"

        unit=400 region=1

        neutron

        energyBoundsID=1

    end regionTally

    regionTally 402

        title="Control pump shaft"

        unit=400 region=2

        neutron

        energyBoundsID=1

    end regionTally



    regionTally 501

        title="External pump rotor"

        unit=500 region=1

        neutron

        energyBoundsID=1

    end regionTally

    regionTally 502

        title="External pump shaft"

        unit=500 region=2

        neutron

        energyBoundsID=1

    end regionTally



    regionTally 601

        title="Roof pump rotor"

        unit=600 region=1

        neutron

        energyBoundsID=1

    end regionTally

    regionTally 602

        title="Roof pump shaft"

        unit=600 region=2

        neutron

        energyBoundsID=1

    end regionTally



    regionTally 777

        title="Gold in the treasure chest"

        unit=777 region=3

        neutron

        energyBoundsID=1

    end regionTally





    meshTally 1

        title="Neutron Flux"

        neutron

        energyBoundsID=1

        gridGeometryID=2

        nogroupfluxes

    end meshTally



end tallies



'-------------------------------------------------------------------------------

' Parameters Block

'-------------------------------------------------------------------------------

read parameters

    celibrary="ce_v7.1_endf.xml"

    neutrons

    fissionmult=0 secondaryMult=2

    randomSeed=8655761690

    perBatch=1000000 batches=100

    maximumLost=100

    allowshortimpmap

end parameters



read importancemap

    gridGeometryID=1

'   roof pump

    adjointSource 600

        boundingbox 1120.0 1080.0 320.0 280.0 440.0 300.0

        responseID=1

    end adjointSource

    legendre=1

end importancemap



end data

end



=shell

    cp source.msm ${OUTBASE}.source.msm

end
 
 No newline at end of file