aparler/ecp-benchmarks
1
0
Fork 0
forked from crp/ecp-benchmarks
Code Pull requests Activity

Updated SMR model per comments by @paulromano

Browse source
This commit is contained in:
Will Boyd 2017-02-09 12:49:10 -05:00
parent 060f4d912e
commit 34cb034c12
12 changed files with 357 additions and 371 deletions
Download patch file Download diff file Expand all files Collapse all files

0
smr/__init__.py
View file

21
smr/build-xml.py
View file

@ -19,7 +19,7 @@ materials.export_to_xml()
# Query the user on whether to use multipole cross sections
multipole = input('Use multipole cross sections? (y/n): ').lower()
multipole = True if multipole == 'y' else False
multipole = (multipole == 'y')
# Construct uniform initial source distribution over fissionable zones
lower_left = [-7.*lattice_pitch/2., -7.*lattice_pitch/2., bottom_fuel_stack]
@ -27,19 +27,18 @@ upper_right = [+7.*lattice_pitch/2., +7.*lattice_pitch/2., top_active_core]
source = openmc.source.Source(space=openmc.stats.Box(lower_left, upper_right))
source.space.only_fissionable = True
settings_file = openmc.Settings()
settings_file.batches = 200
settings_file.inactive = 100
settings_file.particles = 10000
settings_file.ptables = True
settings_file.output = {'tallies': False}
settings_file.source = source
settings_file.sourcepoint_write = False
settings = openmc.Settings()
settings.batches = 200
settings.inactive = 100
settings.particles = 10000
settings.output = {'tallies': False}
settings.source = source
settings.sourcepoint_write = False
if multipole:
settings_file.temperature = {'multipole': True, 'tolerance': 1000}
settings.temperature = {'multipole': True, 'tolerance': 1000}
settings_file.export_to_xml()
settings.export_to_xml()
#### Create OpenMC "plots.xml" file

334
smr/geometry.xml
View file

@ -6,11 +6,11 @@
<cell id="10003" material="10007" name="GT empty (last)" region="10005" universe="10001" />
<cell id="10004" material="10007" name="GT empty grid (top/bottom) (0)" region="-10004" universe="10002" />
<cell id="10005" material="10005" name="GT empty grid (top/bottom) (1)" region="10004 -10005" universe="10002" />
<cell id="10006" material="10007" name="GT empty grid (top/bottom) (last)" region="10005 (10019 -10020 10021 -10022)" universe="10002" />
<cell id="10006" material="10007" name="GT empty grid (top/bottom) (last)" region="10005 10019 -10020 10021 -10022" universe="10002" />
<cell id="10007" material="10005" name="GT empty grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10002" />
<cell id="10008" material="10007" name="GT empty grid (intermediate) (0)" region="-10004" universe="10003" />
<cell id="10009" material="10005" name="GT empty grid (intermediate) (1)" region="10004 -10005" universe="10003" />
<cell id="10010" material="10007" name="GT empty grid (intermediate) (last)" region="10005 (10023 -10024 10025 -10026)" universe="10003" />
<cell id="10010" material="10007" name="GT empty grid (intermediate) (last)" region="10005 10023 -10024 10025 -10026" universe="10003" />
<cell id="10011" material="10005" name="GT empty grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10003" />
<cell id="10012" material="10007" name="GT empty nozzle (0)" region="-10004" universe="10004" />
<cell id="10013" material="10005" name="GT empty nozzle (1)" region="10004 -10005" universe="10004" />
@ -20,11 +20,11 @@
<cell id="10017" material="10007" name="GT empty at dashpot (last)" region="10007" universe="10005" />
<cell id="10018" material="10007" name="GT empty at dashpot grid (top/bottom) (0)" region="-10006" universe="10006" />
<cell id="10019" material="10005" name="GT empty at dashpot grid (top/bottom) (1)" region="10006 -10007" universe="10006" />
<cell id="10020" material="10007" name="GT empty at dashpot grid (top/bottom) (last)" region="10007 (10019 -10020 10021 -10022)" universe="10006" />
<cell id="10020" material="10007" name="GT empty at dashpot grid (top/bottom) (last)" region="10007 10019 -10020 10021 -10022" universe="10006" />
<cell id="10021" material="10005" name="GT empty at dashpot grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10006" />
<cell id="10022" material="10007" name="GT empty at dashpot grid (intermediate) (0)" region="-10006" universe="10007" />
<cell id="10023" material="10005" name="GT empty at dashpot grid (intermediate) (1)" region="10006 -10007" universe="10007" />
<cell id="10024" material="10007" name="GT empty at dashpot grid (intermediate) (last)" region="10007 (10023 -10024 10025 -10026)" universe="10007" />
<cell id="10024" material="10007" name="GT empty at dashpot grid (intermediate) (last)" region="10007 10023 -10024 10025 -10026" universe="10007" />
<cell id="10025" material="10005" name="GT empty at dashpot grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10007" />
<cell id="10026" material="10007" name="GT empty nozzle (0)" region="-10006" universe="10008" />
<cell id="10027" material="10005" name="GT empty nozzle (1)" region="10006 -10007" universe="10008" />
@ -76,13 +76,13 @@
<cell id="10073" material="10005" name="IT grid (top/bottom) (1)" region="10015 -10016" universe="10012" />
<cell id="10074" material="10007" name="IT grid (top/bottom) (2)" region="10016 -10004" universe="10012" />
<cell id="10075" material="10005" name="IT grid (top/bottom) (3)" region="10004 -10005" universe="10012" />
<cell id="10076" material="10007" name="IT grid (top/bottom) (last)" region="10005 (10019 -10020 10021 -10022)" universe="10012" />
<cell id="10076" material="10007" name="IT grid (top/bottom) (last)" region="10005 10019 -10020 10021 -10022" universe="10012" />
<cell id="10077" material="10005" name="IT grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10012" />
<cell id="10078" material="10001" name="IT grid (intermediate) (0)" region="-10015" universe="10013" />
<cell id="10079" material="10005" name="IT grid (intermediate) (1)" region="10015 -10016" universe="10013" />
<cell id="10080" material="10007" name="IT grid (intermediate) (2)" region="10016 -10004" universe="10013" />
<cell id="10081" material="10005" name="IT grid (intermediate) (3)" region="10004 -10005" universe="10013" />
<cell id="10082" material="10007" name="IT grid (intermediate) (last)" region="10005 (10023 -10024 10025 -10026)" universe="10013" />
<cell id="10082" material="10007" name="IT grid (intermediate) (last)" region="10005 10023 -10024 10025 -10026" universe="10013" />
<cell id="10083" material="10005" name="IT grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10013" />
<cell id="10089" material="10001" name="IT dashpot (0)" region="-10015" universe="10015" />
<cell id="10090" material="10005" name="IT dashpot (1)" region="10015 -10016" universe="10015" />
@ -115,13 +115,13 @@
<cell id="10117" material="10001" name="CR grid (top/bottom) (1)" region="10008 -10009" universe="10018" />
<cell id="10118" material="10003" name="CR grid (top/bottom) (2)" region="10009 -10004" universe="10018" />
<cell id="10119" material="10007" name="CR grid (top/bottom) (3)" region="10004 -10005" universe="10018" />
<cell id="10120" material="10007" name="CR grid (top/bottom) (last)" region="10005 (10019 -10020 10021 -10022)" universe="10018" />
<cell id="10120" material="10007" name="CR grid (top/bottom) (last)" region="10005 10019 -10020 10021 -10022" universe="10018" />
<cell id="10121" material="10005" name="CR grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10018" />
<cell id="10122" material="10006" name="CR grid (intermediate) (0)" region="-10008" universe="10019" />
<cell id="10123" material="10001" name="CR grid (intermediate) (1)" region="10008 -10009" universe="10019" />
<cell id="10124" material="10003" name="CR grid (intermediate) (2)" region="10009 -10004" universe="10019" />
<cell id="10125" material="10007" name="CR grid (intermediate) (3)" region="10004 -10005" universe="10019" />
<cell id="10126" material="10007" name="CR grid (intermediate) (last)" region="10005 (10023 -10024 10025 -10026)" universe="10019" />
<cell id="10126" material="10007" name="CR grid (intermediate) (last)" region="10005 10023 -10024 10025 -10026" universe="10019" />
<cell id="10127" material="10005" name="CR grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10019" />
<cell id="10128" material="10006" name="CR nozzle (0)" region="-10008" universe="10020" />
<cell id="10129" material="10001" name="CR nozzle (1)" region="10008 -10009" universe="10020" />
@ -138,14 +138,14 @@
<cell id="10140" material="10003" name="CR blank grid (top/bottom) (2)" region="10009 -10010" universe="10022" />
<cell id="10141" material="10007" name="CR blank grid (top/bottom) (3)" region="10010 -10004" universe="10022" />
<cell id="10142" material="10005" name="CR blank grid (top/bottom) (4)" region="10004 -10005" universe="10022" />
<cell id="10143" material="10007" name="CR blank grid (top/bottom) (last)" region="10005 (10019 -10020 10021 -10022)" universe="10022" />
<cell id="10143" material="10007" name="CR blank grid (top/bottom) (last)" region="10005 10019 -10020 10021 -10022" universe="10022" />
<cell id="10144" material="10005" name="CR blank grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10022" />
<cell id="10145" material="10003" name="CR blank grid (intermediate) (0)" region="-10008" universe="10023" />
<cell id="10146" material="10001" name="CR blank grid (intermediate) (1)" region="10008 -10009" universe="10023" />
<cell id="10147" material="10003" name="CR blank grid (intermediate) (2)" region="10009 -10010" universe="10023" />
<cell id="10148" material="10007" name="CR blank grid (intermediate) (3)" region="10010 -10004" universe="10023" />
<cell id="10149" material="10005" name="CR blank grid (intermediate) (4)" region="10004 -10005" universe="10023" />
<cell id="10150" material="10007" name="CR blank grid (intermediate) (last)" region="10005 (10023 -10024 10025 -10026)" universe="10023" />
<cell id="10150" material="10007" name="CR blank grid (intermediate) (last)" region="10005 10023 -10024 10025 -10026" universe="10023" />
<cell id="10151" material="10005" name="CR blank grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10023" />
<cell id="10152" material="10003" name="CR blank nozzle (0)" region="-10008" universe="10024" />
<cell id="10153" material="10001" name="CR blank nozzle (1)" region="10008 -10009" universe="10024" />
@ -220,7 +220,7 @@
<cell id="10584" material="10003" name="BA grid (intermediate) (5)" region="10015 -10016" universe="10083" />
<cell id="10585" material="10007" name="BA grid (intermediate) (6)" region="10016 -10017" universe="10083" />
<cell id="10586" material="10005" name="BA grid (intermediate) (7)" region="10017 -10018" universe="10083" />
<cell id="10587" material="10007" name="BA grid (intermediate) (last)" region="10018 (10023 -10024 10025 -10026)" universe="10083" />
<cell id="10587" material="10007" name="BA grid (intermediate) (last)" region="10018 10023 -10024 10025 -10026" universe="10083" />
<cell id="10588" material="10005" name="BA grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10083" />
<cell id="10589" material="10001" name="BA dashpot (0)" region="-10011" universe="10084" />
<cell id="10590" material="10003" name="BA dashpot (1)" region="10011 -10012" universe="10084" />
@ -239,7 +239,7 @@
<cell id="10603" material="10003" name="BA dashpot grid (top/bottom) (5)" region="10015 -10016" universe="10085" />
<cell id="10604" material="10007" name="BA dashpot grid (top/bottom) (6)" region="10016 -10006" universe="10085" />
<cell id="10605" material="10005" name="BA dashpot grid (top/bottom) (7)" region="10006 -10007" universe="10085" />
<cell id="10606" material="10007" name="BA dashpot grid (top/bottom) (last)" region="10007 (10019 -10020 10021 -10022)" universe="10085" />
<cell id="10606" material="10007" name="BA dashpot grid (top/bottom) (last)" region="10007 10019 -10020 10021 -10022" universe="10085" />
<cell id="10607" material="10005" name="BA dashpot grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10085" />
<cell id="10618" material="10003" name="BA blank SS (0)" region="-10016" universe="10087" />
<cell id="10619" material="10007" name="BA blank SS (1)" region="10016 -10017" universe="10087" />
@ -278,151 +278,151 @@
<cell id="10652" material="10002" name="pin plenum grid (top/bottom) (0)" region="-10001" universe="10093" />
<cell id="10653" material="10000" name="pin plenum grid (top/bottom) (1)" region="10001 -10002" universe="10093" />
<cell id="10654" material="10005" name="pin plenum grid (top/bottom) (2)" region="10002 -10003" universe="10093" />
<cell id="10655" material="10007" name="pin plenum grid (top/bottom) (last)" region="10003 (10019 -10020 10021 -10022)" universe="10093" />
<cell id="10655" material="10007" name="pin plenum grid (top/bottom) (last)" region="10003 10019 -10020 10021 -10022" universe="10093" />
<cell id="10656" material="10005" name="pin plenum grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10093" />
<cell id="10657" material="10009" name="Fuel (1.6\%) (0)" region="-10000" universe="10094" />
<cell id="10658" material="10000" name="Fuel (1.6\%) (1)" region="10000 -10002" universe="10094" />
<cell id="10659" material="10005" name="Fuel (1.6\%) (2)" region="10002 -10003" universe="10094" />
<cell id="10660" material="10007" name="Fuel (1.6\%) (last)" region="10003" universe="10094" />
<cell id="10661" material="10009" name="Fuel (1.6\%) grid (top/bottom) (0)" region="-10000" universe="10095" />
<cell id="10662" material="10000" name="Fuel (1.6\%) grid (top/bottom) (1)" region="10000 -10002" universe="10095" />
<cell id="10663" material="10005" name="Fuel (1.6\%) grid (top/bottom) (2)" region="10002 -10003" universe="10095" />
<cell id="10664" material="10007" name="Fuel (1.6\%) grid (top/bottom) (last)" region="10003 (10019 -10020 10021 -10022)" universe="10095" />
<cell id="10665" material="10005" name="Fuel (1.6\%) grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10095" />
<cell id="10666" material="10009" name="Fuel (1.6\%) grid (intermediate) (0)" region="-10000" universe="10096" />
<cell id="10667" material="10000" name="Fuel (1.6\%) grid (intermediate) (1)" region="10000 -10002" universe="10096" />
<cell id="10668" material="10005" name="Fuel (1.6\%) grid (intermediate) (2)" region="10002 -10003" universe="10096" />
<cell id="10669" material="10007" name="Fuel (1.6\%) grid (intermediate) (last)" region="10003 (10023 -10024 10025 -10026)" universe="10096" />
<cell id="10670" material="10005" name="Fuel (1.6\%) grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10096" />
<cell fill="10000" id="10671" name="Fuel (1.6\%) stack (0)" region="-10035" universe="10097" />
<cell fill="10090" id="10672" name="Fuel (1.6\%) stack (1)" region="10035 -10036" universe="10097" />
<cell fill="10090" id="10673" name="Fuel (1.6\%) stack (2)" region="10036 -10037" universe="10097" />
<cell fill="10091" id="10674" name="Fuel (1.6\%) stack (3)" region="10037 -10038" universe="10097" />
<cell fill="10094" id="10675" name="Fuel (1.6\%) stack (4)" region="10038 -10041" universe="10097" />
<cell fill="10095" id="10676" name="Fuel (1.6\%) stack (5)" region="10041 -10042" universe="10097" />
<cell fill="10094" id="10677" name="Fuel (1.6\%) stack (6)" region="10042 -10043" universe="10097" />
<cell fill="10094" id="10678" name="Fuel (1.6\%) stack (7)" region="10043 -10044" universe="10097" />
<cell fill="10096" id="10679" name="Fuel (1.6\%) stack (8)" region="10044 -10045" universe="10097" />
<cell fill="10094" id="10680" name="Fuel (1.6\%) stack (9)" region="10045 -10046" universe="10097" />
<cell fill="10096" id="10681" name="Fuel (1.6\%) stack (10)" region="10046 -10047" universe="10097" />
<cell fill="10094" id="10682" name="Fuel (1.6\%) stack (11)" region="10047 -10048" universe="10097" />
<cell fill="10092" id="10683" name="Fuel (1.6\%) stack (12)" region="10048 -10049" universe="10097" />
<cell fill="10093" id="10684" name="Fuel (1.6\%) stack (13)" region="10049 -10039" universe="10097" />
<cell fill="10092" id="10685" name="Fuel (1.6\%) stack (14)" region="10039 -10050" universe="10097" />
<cell fill="10091" id="10686" name="Fuel (1.6\%) stack (15)" region="10050 -10051" universe="10097" />
<cell fill="10000" id="10687" name="Fuel (1.6\%) stack (16)" region="10051 -10052" universe="10097" />
<cell fill="10090" id="10688" name="Fuel (1.6\%) stack (17)" region="10052 -10053" universe="10097" />
<cell fill="10000" id="10689" name="Fuel (1.6\%) stack (last)" region="10053" universe="10097" />
<cell id="10690" material="10010" name="Fuel (2.4\%) (0)" region="-10000" universe="10098" />
<cell id="10691" material="10000" name="Fuel (2.4\%) (1)" region="10000 -10002" universe="10098" />
<cell id="10692" material="10005" name="Fuel (2.4\%) (2)" region="10002 -10003" universe="10098" />
<cell id="10693" material="10007" name="Fuel (2.4\%) (last)" region="10003" universe="10098" />
<cell id="10694" material="10010" name="Fuel (2.4\%) grid (top/bottom) (0)" region="-10000" universe="10099" />
<cell id="10695" material="10000" name="Fuel (2.4\%) grid (top/bottom) (1)" region="10000 -10002" universe="10099" />
<cell id="10696" material="10005" name="Fuel (2.4\%) grid (top/bottom) (2)" region="10002 -10003" universe="10099" />
<cell id="10697" material="10007" name="Fuel (2.4\%) grid (top/bottom) (last)" region="10003 (10019 -10020 10021 -10022)" universe="10099" />
<cell id="10698" material="10005" name="Fuel (2.4\%) grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10099" />
<cell id="10699" material="10010" name="Fuel (2.4\%) grid (intermediate) (0)" region="-10000" universe="10100" />
<cell id="10700" material="10000" name="Fuel (2.4\%) grid (intermediate) (1)" region="10000 -10002" universe="10100" />
<cell id="10701" material="10005" name="Fuel (2.4\%) grid (intermediate) (2)" region="10002 -10003" universe="10100" />
<cell id="10702" material="10007" name="Fuel (2.4\%) grid (intermediate) (last)" region="10003 (10023 -10024 10025 -10026)" universe="10100" />
<cell id="10703" material="10005" name="Fuel (2.4\%) grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10100" />
<cell fill="10000" id="10704" name="Fuel (2.4\%) stack (0)" region="-10035" universe="10101" />
<cell fill="10090" id="10705" name="Fuel (2.4\%) stack (1)" region="10035 -10036" universe="10101" />
<cell fill="10090" id="10706" name="Fuel (2.4\%) stack (2)" region="10036 -10037" universe="10101" />
<cell fill="10091" id="10707" name="Fuel (2.4\%) stack (3)" region="10037 -10038" universe="10101" />
<cell fill="10098" id="10708" name="Fuel (2.4\%) stack (4)" region="10038 -10041" universe="10101" />
<cell fill="10099" id="10709" name="Fuel (2.4\%) stack (5)" region="10041 -10042" universe="10101" />
<cell fill="10098" id="10710" name="Fuel (2.4\%) stack (6)" region="10042 -10043" universe="10101" />
<cell fill="10098" id="10711" name="Fuel (2.4\%) stack (7)" region="10043 -10044" universe="10101" />
<cell fill="10100" id="10712" name="Fuel (2.4\%) stack (8)" region="10044 -10045" universe="10101" />
<cell fill="10098" id="10713" name="Fuel (2.4\%) stack (9)" region="10045 -10046" universe="10101" />
<cell fill="10100" id="10714" name="Fuel (2.4\%) stack (10)" region="10046 -10047" universe="10101" />
<cell fill="10098" id="10715" name="Fuel (2.4\%) stack (11)" region="10047 -10048" universe="10101" />
<cell fill="10092" id="10716" name="Fuel (2.4\%) stack (12)" region="10048 -10049" universe="10101" />
<cell fill="10093" id="10717" name="Fuel (2.4\%) stack (13)" region="10049 -10039" universe="10101" />
<cell fill="10092" id="10718" name="Fuel (2.4\%) stack (14)" region="10039 -10050" universe="10101" />
<cell fill="10091" id="10719" name="Fuel (2.4\%) stack (15)" region="10050 -10051" universe="10101" />
<cell fill="10000" id="10720" name="Fuel (2.4\%) stack (16)" region="10051 -10052" universe="10101" />
<cell fill="10090" id="10721" name="Fuel (2.4\%) stack (17)" region="10052 -10053" universe="10101" />
<cell fill="10000" id="10722" name="Fuel (2.4\%) stack (last)" region="10053" universe="10101" />
<cell id="10723" material="10011" name="Fuel (3.1\%) (0)" region="-10000" universe="10102" />
<cell id="10724" material="10000" name="Fuel (3.1\%) (1)" region="10000 -10002" universe="10102" />
<cell id="10725" material="10005" name="Fuel (3.1\%) (2)" region="10002 -10003" universe="10102" />
<cell id="10726" material="10007" name="Fuel (3.1\%) (last)" region="10003" universe="10102" />
<cell id="10727" material="10011" name="Fuel (3.1\%) grid (top/bottom) (0)" region="-10000" universe="10103" />
<cell id="10728" material="10000" name="Fuel (3.1\%) grid (top/bottom) (1)" region="10000 -10002" universe="10103" />
<cell id="10729" material="10005" name="Fuel (3.1\%) grid (top/bottom) (2)" region="10002 -10003" universe="10103" />
<cell id="10730" material="10007" name="Fuel (3.1\%) grid (top/bottom) (last)" region="10003 (10019 -10020 10021 -10022)" universe="10103" />
<cell id="10731" material="10005" name="Fuel (3.1\%) grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10103" />
<cell id="10732" material="10011" name="Fuel (3.1\%) grid (intermediate) (0)" region="-10000" universe="10104" />
<cell id="10733" material="10000" name="Fuel (3.1\%) grid (intermediate) (1)" region="10000 -10002" universe="10104" />
<cell id="10734" material="10005" name="Fuel (3.1\%) grid (intermediate) (2)" region="10002 -10003" universe="10104" />
<cell id="10735" material="10007" name="Fuel (3.1\%) grid (intermediate) (last)" region="10003 (10023 -10024 10025 -10026)" universe="10104" />
<cell id="10736" material="10005" name="Fuel (3.1\%) grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10104" />
<cell fill="10000" id="10737" name="Fuel (3.1\%) stack (0)" region="-10035" universe="10105" />
<cell fill="10090" id="10738" name="Fuel (3.1\%) stack (1)" region="10035 -10036" universe="10105" />
<cell fill="10090" id="10739" name="Fuel (3.1\%) stack (2)" region="10036 -10037" universe="10105" />
<cell fill="10091" id="10740" name="Fuel (3.1\%) stack (3)" region="10037 -10038" universe="10105" />
<cell fill="10102" id="10741" name="Fuel (3.1\%) stack (4)" region="10038 -10041" universe="10105" />
<cell fill="10103" id="10742" name="Fuel (3.1\%) stack (5)" region="10041 -10042" universe="10105" />
<cell fill="10102" id="10743" name="Fuel (3.1\%) stack (6)" region="10042 -10043" universe="10105" />
<cell fill="10102" id="10744" name="Fuel (3.1\%) stack (7)" region="10043 -10044" universe="10105" />
<cell fill="10104" id="10745" name="Fuel (3.1\%) stack (8)" region="10044 -10045" universe="10105" />
<cell fill="10102" id="10746" name="Fuel (3.1\%) stack (9)" region="10045 -10046" universe="10105" />
<cell fill="10104" id="10747" name="Fuel (3.1\%) stack (10)" region="10046 -10047" universe="10105" />
<cell fill="10102" id="10748" name="Fuel (3.1\%) stack (11)" region="10047 -10048" universe="10105" />
<cell fill="10092" id="10749" name="Fuel (3.1\%) stack (12)" region="10048 -10049" universe="10105" />
<cell fill="10093" id="10750" name="Fuel (3.1\%) stack (13)" region="10049 -10039" universe="10105" />
<cell fill="10092" id="10751" name="Fuel (3.1\%) stack (14)" region="10039 -10050" universe="10105" />
<cell fill="10091" id="10752" name="Fuel (3.1\%) stack (15)" region="10050 -10051" universe="10105" />
<cell fill="10000" id="10753" name="Fuel (3.1\%) stack (16)" region="10051 -10052" universe="10105" />
<cell fill="10090" id="10754" name="Fuel (3.1\%) stack (17)" region="10052 -10053" universe="10105" />
<cell fill="10000" id="10755" name="Fuel (3.1\%) stack (last)" region="10053" universe="10105" />
<cell fill="10124" id="10855" name="Assembly (1.6\%) no BAs instr lattice" region="10027 -10028 10029 -10030" universe="10125" />
<cell fill="10000" id="10856" name="Assembly (1.6\%) no BAs instr lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10125" />
<cell id="10857" material="10007" name="Assembly (1.6\%) no BAs instr lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10041" universe="10125" />
<cell id="10858" material="10003" name="Assembly (1.6\%) no BAs instr lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10041 -10042" universe="10125" />
<cell id="10859" material="10007" name="Assembly (1.6\%) no BAs instr lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10044" universe="10125" />
<cell id="10860" material="10005" name="Assembly (1.6\%) no BAs instr lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10044 -10045" universe="10125" />
<cell id="10861" material="10007" name="Assembly (1.6\%) no BAs instr lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10125" />
<cell id="10862" material="10005" name="Assembly (1.6\%) no BAs instr lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10125" />
<cell id="10863" material="10007" name="Assembly (1.6\%) no BAs instr lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10125" />
<cell id="10864" material="10003" name="Assembly (1.6\%) no BAs instr lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10125" />
<cell id="10865" material="10007" name="Assembly (1.6\%) no BAs instr lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049" universe="10125" />
<cell fill="10144" id="10965" name="Assembly (2.4\%) CR D lattice" region="10027 -10028 10029 -10030" universe="10145" />
<cell fill="10000" id="10966" name="Assembly (2.4\%) CR D lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10145" />
<cell id="10967" material="10007" name="Assembly (2.4\%) CR D lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10041" universe="10145" />
<cell id="10968" material="10003" name="Assembly (2.4\%) CR D lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10041 -10042" universe="10145" />
<cell id="10969" material="10007" name="Assembly (2.4\%) CR D lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10044" universe="10145" />
<cell id="10970" material="10005" name="Assembly (2.4\%) CR D lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10044 -10045" universe="10145" />
<cell id="10971" material="10007" name="Assembly (2.4\%) CR D lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10145" />
<cell id="10972" material="10005" name="Assembly (2.4\%) CR D lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10145" />
<cell id="10973" material="10007" name="Assembly (2.4\%) CR D lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10145" />
<cell id="10974" material="10003" name="Assembly (2.4\%) CR D lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10145" />
<cell id="10975" material="10007" name="Assembly (2.4\%) CR D lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049" universe="10145" />
<cell fill="10164" id="11075" name="Assembly (3.1\%) 16BA lattice" region="10027 -10028 10029 -10030" universe="10165" />
<cell fill="10000" id="11076" name="Assembly (3.1\%) 16BA lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10165" />
<cell id="11077" material="10007" name="Assembly (3.1\%) 16BA lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10041" universe="10165" />
<cell id="11078" material="10003" name="Assembly (3.1\%) 16BA lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10041 -10042" universe="10165" />
<cell id="11079" material="10007" name="Assembly (3.1\%) 16BA lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10044" universe="10165" />
<cell id="11080" material="10005" name="Assembly (3.1\%) 16BA lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10044 -10045" universe="10165" />
<cell id="11081" material="10007" name="Assembly (3.1\%) 16BA lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10165" />
<cell id="11082" material="10005" name="Assembly (3.1\%) 16BA lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10165" />
<cell id="11083" material="10007" name="Assembly (3.1\%) 16BA lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10165" />
<cell id="11084" material="10003" name="Assembly (3.1\%) 16BA lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10165" />
<cell id="11085" material="10007" name="Assembly (3.1\%) 16BA lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049" universe="10165" />
<cell fill="10182" id="11174" name="Assembly (3.1\%) no BAs instr lattice" region="10027 -10028 10029 -10030" universe="10183" />
<cell fill="10000" id="11175" name="Assembly (3.1\%) no BAs instr lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10183" />
<cell id="11176" material="10007" name="Assembly (3.1\%) no BAs instr lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10041" universe="10183" />
<cell id="11177" material="10003" name="Assembly (3.1\%) no BAs instr lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10041 -10042" universe="10183" />
<cell id="11178" material="10007" name="Assembly (3.1\%) no BAs instr lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10044" universe="10183" />
<cell id="11179" material="10005" name="Assembly (3.1\%) no BAs instr lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10044 -10045" universe="10183" />
<cell id="11180" material="10007" name="Assembly (3.1\%) no BAs instr lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10183" />
<cell id="11181" material="10005" name="Assembly (3.1\%) no BAs instr lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10183" />
<cell id="11182" material="10007" name="Assembly (3.1\%) no BAs instr lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10183" />
<cell id="11183" material="10003" name="Assembly (3.1\%) no BAs instr lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10183" />
<cell id="11184" material="10007" name="Assembly (3.1\%) no BAs instr lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049" universe="10183" />
<cell id="10657" material="10009" name="Fuel (1.6%) (0)" region="-10000" universe="10094" />
<cell id="10658" material="10000" name="Fuel (1.6%) (1)" region="10000 -10002" universe="10094" />
<cell id="10659" material="10005" name="Fuel (1.6%) (2)" region="10002 -10003" universe="10094" />
<cell id="10660" material="10007" name="Fuel (1.6%) (last)" region="10003" universe="10094" />
<cell id="10661" material="10009" name="Fuel (1.6%) grid (top/bottom) (0)" region="-10000" universe="10095" />
<cell id="10662" material="10000" name="Fuel (1.6%) grid (top/bottom) (1)" region="10000 -10002" universe="10095" />
<cell id="10663" material="10005" name="Fuel (1.6%) grid (top/bottom) (2)" region="10002 -10003" universe="10095" />
<cell id="10664" material="10007" name="Fuel (1.6%) grid (top/bottom) (last)" region="10003 10019 -10020 10021 -10022" universe="10095" />
<cell id="10665" material="10005" name="Fuel (1.6%) grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10095" />
<cell id="10666" material="10009" name="Fuel (1.6%) grid (intermediate) (0)" region="-10000" universe="10096" />
<cell id="10667" material="10000" name="Fuel (1.6%) grid (intermediate) (1)" region="10000 -10002" universe="10096" />
<cell id="10668" material="10005" name="Fuel (1.6%) grid (intermediate) (2)" region="10002 -10003" universe="10096" />
<cell id="10669" material="10007" name="Fuel (1.6%) grid (intermediate) (last)" region="10003 10023 -10024 10025 -10026" universe="10096" />
<cell id="10670" material="10005" name="Fuel (1.6%) grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10096" />
<cell fill="10000" id="10671" name="Fuel (1.6%) stack (0)" region="-10035" universe="10097" />
<cell fill="10090" id="10672" name="Fuel (1.6%) stack (1)" region="10035 -10036" universe="10097" />
<cell fill="10090" id="10673" name="Fuel (1.6%) stack (2)" region="10036 -10037" universe="10097" />
<cell fill="10091" id="10674" name="Fuel (1.6%) stack (3)" region="10037 -10038" universe="10097" />
<cell fill="10094" id="10675" name="Fuel (1.6%) stack (4)" region="10038 -10041" universe="10097" />
<cell fill="10095" id="10676" name="Fuel (1.6%) stack (5)" region="10041 -10042" universe="10097" />
<cell fill="10094" id="10677" name="Fuel (1.6%) stack (6)" region="10042 -10043" universe="10097" />
<cell fill="10094" id="10678" name="Fuel (1.6%) stack (7)" region="10043 -10044" universe="10097" />
<cell fill="10096" id="10679" name="Fuel (1.6%) stack (8)" region="10044 -10045" universe="10097" />
<cell fill="10094" id="10680" name="Fuel (1.6%) stack (9)" region="10045 -10046" universe="10097" />
<cell fill="10096" id="10681" name="Fuel (1.6%) stack (10)" region="10046 -10047" universe="10097" />
<cell fill="10094" id="10682" name="Fuel (1.6%) stack (11)" region="10047 -10048" universe="10097" />
<cell fill="10092" id="10683" name="Fuel (1.6%) stack (12)" region="10048 -10049" universe="10097" />
<cell fill="10093" id="10684" name="Fuel (1.6%) stack (13)" region="10049 -10039" universe="10097" />
<cell fill="10092" id="10685" name="Fuel (1.6%) stack (14)" region="10039 -10050" universe="10097" />
<cell fill="10091" id="10686" name="Fuel (1.6%) stack (15)" region="10050 -10051" universe="10097" />
<cell fill="10000" id="10687" name="Fuel (1.6%) stack (16)" region="10051 -10052" universe="10097" />
<cell fill="10090" id="10688" name="Fuel (1.6%) stack (17)" region="10052 -10053" universe="10097" />
<cell fill="10000" id="10689" name="Fuel (1.6%) stack (last)" region="10053" universe="10097" />
<cell id="10690" material="10010" name="Fuel (2.4%) (0)" region="-10000" universe="10098" />
<cell id="10691" material="10000" name="Fuel (2.4%) (1)" region="10000 -10002" universe="10098" />
<cell id="10692" material="10005" name="Fuel (2.4%) (2)" region="10002 -10003" universe="10098" />
<cell id="10693" material="10007" name="Fuel (2.4%) (last)" region="10003" universe="10098" />
<cell id="10694" material="10010" name="Fuel (2.4%) grid (top/bottom) (0)" region="-10000" universe="10099" />
<cell id="10695" material="10000" name="Fuel (2.4%) grid (top/bottom) (1)" region="10000 -10002" universe="10099" />
<cell id="10696" material="10005" name="Fuel (2.4%) grid (top/bottom) (2)" region="10002 -10003" universe="10099" />
<cell id="10697" material="10007" name="Fuel (2.4%) grid (top/bottom) (last)" region="10003 10019 -10020 10021 -10022" universe="10099" />
<cell id="10698" material="10005" name="Fuel (2.4%) grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10099" />
<cell id="10699" material="10010" name="Fuel (2.4%) grid (intermediate) (0)" region="-10000" universe="10100" />
<cell id="10700" material="10000" name="Fuel (2.4%) grid (intermediate) (1)" region="10000 -10002" universe="10100" />
<cell id="10701" material="10005" name="Fuel (2.4%) grid (intermediate) (2)" region="10002 -10003" universe="10100" />
<cell id="10702" material="10007" name="Fuel (2.4%) grid (intermediate) (last)" region="10003 10023 -10024 10025 -10026" universe="10100" />
<cell id="10703" material="10005" name="Fuel (2.4%) grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10100" />
<cell fill="10000" id="10704" name="Fuel (2.4%) stack (0)" region="-10035" universe="10101" />
<cell fill="10090" id="10705" name="Fuel (2.4%) stack (1)" region="10035 -10036" universe="10101" />
<cell fill="10090" id="10706" name="Fuel (2.4%) stack (2)" region="10036 -10037" universe="10101" />
<cell fill="10091" id="10707" name="Fuel (2.4%) stack (3)" region="10037 -10038" universe="10101" />
<cell fill="10098" id="10708" name="Fuel (2.4%) stack (4)" region="10038 -10041" universe="10101" />
<cell fill="10099" id="10709" name="Fuel (2.4%) stack (5)" region="10041 -10042" universe="10101" />
<cell fill="10098" id="10710" name="Fuel (2.4%) stack (6)" region="10042 -10043" universe="10101" />
<cell fill="10098" id="10711" name="Fuel (2.4%) stack (7)" region="10043 -10044" universe="10101" />
<cell fill="10100" id="10712" name="Fuel (2.4%) stack (8)" region="10044 -10045" universe="10101" />
<cell fill="10098" id="10713" name="Fuel (2.4%) stack (9)" region="10045 -10046" universe="10101" />
<cell fill="10100" id="10714" name="Fuel (2.4%) stack (10)" region="10046 -10047" universe="10101" />
<cell fill="10098" id="10715" name="Fuel (2.4%) stack (11)" region="10047 -10048" universe="10101" />
<cell fill="10092" id="10716" name="Fuel (2.4%) stack (12)" region="10048 -10049" universe="10101" />
<cell fill="10093" id="10717" name="Fuel (2.4%) stack (13)" region="10049 -10039" universe="10101" />
<cell fill="10092" id="10718" name="Fuel (2.4%) stack (14)" region="10039 -10050" universe="10101" />
<cell fill="10091" id="10719" name="Fuel (2.4%) stack (15)" region="10050 -10051" universe="10101" />
<cell fill="10000" id="10720" name="Fuel (2.4%) stack (16)" region="10051 -10052" universe="10101" />
<cell fill="10090" id="10721" name="Fuel (2.4%) stack (17)" region="10052 -10053" universe="10101" />
<cell fill="10000" id="10722" name="Fuel (2.4%) stack (last)" region="10053" universe="10101" />
<cell id="10723" material="10011" name="Fuel (3.1%) (0)" region="-10000" universe="10102" />
<cell id="10724" material="10000" name="Fuel (3.1%) (1)" region="10000 -10002" universe="10102" />
<cell id="10725" material="10005" name="Fuel (3.1%) (2)" region="10002 -10003" universe="10102" />
<cell id="10726" material="10007" name="Fuel (3.1%) (last)" region="10003" universe="10102" />
<cell id="10727" material="10011" name="Fuel (3.1%) grid (top/bottom) (0)" region="-10000" universe="10103" />
<cell id="10728" material="10000" name="Fuel (3.1%) grid (top/bottom) (1)" region="10000 -10002" universe="10103" />
<cell id="10729" material="10005" name="Fuel (3.1%) grid (top/bottom) (2)" region="10002 -10003" universe="10103" />
<cell id="10730" material="10007" name="Fuel (3.1%) grid (top/bottom) (last)" region="10003 10019 -10020 10021 -10022" universe="10103" />
<cell id="10731" material="10005" name="Fuel (3.1%) grid (top/bottom) (grid)" region="~(10019 -10020 10021 -10022)" universe="10103" />
<cell id="10732" material="10011" name="Fuel (3.1%) grid (intermediate) (0)" region="-10000" universe="10104" />
<cell id="10733" material="10000" name="Fuel (3.1%) grid (intermediate) (1)" region="10000 -10002" universe="10104" />
<cell id="10734" material="10005" name="Fuel (3.1%) grid (intermediate) (2)" region="10002 -10003" universe="10104" />
<cell id="10735" material="10007" name="Fuel (3.1%) grid (intermediate) (last)" region="10003 10023 -10024 10025 -10026" universe="10104" />
<cell id="10736" material="10005" name="Fuel (3.1%) grid (intermediate) (grid)" region="~(10023 -10024 10025 -10026)" universe="10104" />
<cell fill="10000" id="10737" name="Fuel (3.1%) stack (0)" region="-10035" universe="10105" />
<cell fill="10090" id="10738" name="Fuel (3.1%) stack (1)" region="10035 -10036" universe="10105" />
<cell fill="10090" id="10739" name="Fuel (3.1%) stack (2)" region="10036 -10037" universe="10105" />
<cell fill="10091" id="10740" name="Fuel (3.1%) stack (3)" region="10037 -10038" universe="10105" />
<cell fill="10102" id="10741" name="Fuel (3.1%) stack (4)" region="10038 -10041" universe="10105" />
<cell fill="10103" id="10742" name="Fuel (3.1%) stack (5)" region="10041 -10042" universe="10105" />
<cell fill="10102" id="10743" name="Fuel (3.1%) stack (6)" region="10042 -10043" universe="10105" />
<cell fill="10102" id="10744" name="Fuel (3.1%) stack (7)" region="10043 -10044" universe="10105" />
<cell fill="10104" id="10745" name="Fuel (3.1%) stack (8)" region="10044 -10045" universe="10105" />
<cell fill="10102" id="10746" name="Fuel (3.1%) stack (9)" region="10045 -10046" universe="10105" />
<cell fill="10104" id="10747" name="Fuel (3.1%) stack (10)" region="10046 -10047" universe="10105" />
<cell fill="10102" id="10748" name="Fuel (3.1%) stack (11)" region="10047 -10048" universe="10105" />
<cell fill="10092" id="10749" name="Fuel (3.1%) stack (12)" region="10048 -10049" universe="10105" />
<cell fill="10093" id="10750" name="Fuel (3.1%) stack (13)" region="10049 -10039" universe="10105" />
<cell fill="10092" id="10751" name="Fuel (3.1%) stack (14)" region="10039 -10050" universe="10105" />
<cell fill="10091" id="10752" name="Fuel (3.1%) stack (15)" region="10050 -10051" universe="10105" />
<cell fill="10000" id="10753" name="Fuel (3.1%) stack (16)" region="10051 -10052" universe="10105" />
<cell fill="10090" id="10754" name="Fuel (3.1%) stack (17)" region="10052 -10053" universe="10105" />
<cell fill="10000" id="10755" name="Fuel (3.1%) stack (last)" region="10053" universe="10105" />
<cell fill="10124" id="10855" name="Assembly (1.6%) no BAs instr lattice" region="10027 -10028 10029 -10030" universe="10125" />
<cell fill="10000" id="10856" name="Assembly (1.6%) no BAs instr lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10125" />
<cell id="10857" material="10007" name="Assembly (1.6%) no BAs instr lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10041" universe="10125" />
<cell id="10858" material="10003" name="Assembly (1.6%) no BAs instr lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10041 -10042" universe="10125" />
<cell id="10859" material="10007" name="Assembly (1.6%) no BAs instr lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10044" universe="10125" />
<cell id="10860" material="10005" name="Assembly (1.6%) no BAs instr lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10044 -10045" universe="10125" />
<cell id="10861" material="10007" name="Assembly (1.6%) no BAs instr lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10125" />
<cell id="10862" material="10005" name="Assembly (1.6%) no BAs instr lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10125" />
<cell id="10863" material="10007" name="Assembly (1.6%) no BAs instr lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10125" />
<cell id="10864" material="10003" name="Assembly (1.6%) no BAs instr lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10125" />
<cell id="10865" material="10007" name="Assembly (1.6%) no BAs instr lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049" universe="10125" />
<cell fill="10144" id="10965" name="Assembly (2.4%) CR D lattice" region="10027 -10028 10029 -10030" universe="10145" />
<cell fill="10000" id="10966" name="Assembly (2.4%) CR D lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10145" />
<cell id="10967" material="10007" name="Assembly (2.4%) CR D lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10041" universe="10145" />
<cell id="10968" material="10003" name="Assembly (2.4%) CR D lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10041 -10042" universe="10145" />
<cell id="10969" material="10007" name="Assembly (2.4%) CR D lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10044" universe="10145" />
<cell id="10970" material="10005" name="Assembly (2.4%) CR D lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10044 -10045" universe="10145" />
<cell id="10971" material="10007" name="Assembly (2.4%) CR D lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10145" />
<cell id="10972" material="10005" name="Assembly (2.4%) CR D lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10145" />
<cell id="10973" material="10007" name="Assembly (2.4%) CR D lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10145" />
<cell id="10974" material="10003" name="Assembly (2.4%) CR D lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10145" />
<cell id="10975" material="10007" name="Assembly (2.4%) CR D lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049" universe="10145" />
<cell fill="10164" id="11075" name="Assembly (3.1%) 16BA lattice" region="10027 -10028 10029 -10030" universe="10165" />
<cell fill="10000" id="11076" name="Assembly (3.1%) 16BA lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10165" />
<cell id="11077" material="10007" name="Assembly (3.1%) 16BA lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10041" universe="10165" />
<cell id="11078" material="10003" name="Assembly (3.1%) 16BA lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10041 -10042" universe="10165" />
<cell id="11079" material="10007" name="Assembly (3.1%) 16BA lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10044" universe="10165" />
<cell id="11080" material="10005" name="Assembly (3.1%) 16BA lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10044 -10045" universe="10165" />
<cell id="11081" material="10007" name="Assembly (3.1%) 16BA lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10165" />
<cell id="11082" material="10005" name="Assembly (3.1%) 16BA lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10165" />
<cell id="11083" material="10007" name="Assembly (3.1%) 16BA lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10165" />
<cell id="11084" material="10003" name="Assembly (3.1%) 16BA lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10165" />
<cell id="11085" material="10007" name="Assembly (3.1%) 16BA lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049" universe="10165" />
<cell fill="10182" id="11174" name="Assembly (3.1%) no BAs instr lattice" region="10027 -10028 10029 -10030" universe="10183" />
<cell fill="10000" id="11175" name="Assembly (3.1%) no BAs instr lattice outer water" region="~(10031 -10032 10033 -10034)" universe="10183" />
<cell id="11176" material="10007" name="Assembly (3.1%) no BAs instr lattice axial (0)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) -10041" universe="10183" />
<cell id="11177" material="10003" name="Assembly (3.1%) no BAs instr lattice axial (1)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10041 -10042" universe="10183" />
<cell id="11178" material="10007" name="Assembly (3.1%) no BAs instr lattice axial (2)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10042 -10044" universe="10183" />
<cell id="11179" material="10005" name="Assembly (3.1%) no BAs instr lattice axial (3)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10044 -10045" universe="10183" />
<cell id="11180" material="10007" name="Assembly (3.1%) no BAs instr lattice axial (4)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10045 -10046" universe="10183" />
<cell id="11181" material="10005" name="Assembly (3.1%) no BAs instr lattice axial (5)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10046 -10047" universe="10183" />
<cell id="11182" material="10007" name="Assembly (3.1%) no BAs instr lattice axial (6)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10047 -10048" universe="10183" />
<cell id="11183" material="10003" name="Assembly (3.1%) no BAs instr lattice axial (7)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10048 -10049" universe="10183" />
<cell id="11184" material="10007" name="Assembly (3.1%) no BAs instr lattice axial (last)" region="10031 -10032 10033 -10034 ~(10027 -10028 10029 -10030) 10049" universe="10183" />
<cell id="11306" material="10003" name="baffle dummy SS" region="10083" universe="10206" />
<cell id="11307" material="10007" name="baffle south water" region="-10083" universe="10206" />
<cell fill="10207" id="11308" name="baffle south" universe="10208" />
@ -479,7 +479,7 @@
<cell id="11359" material="10007" name="neutron shield panel W" region="10073 -10074 -10075 -10077 10082 -10081" universe="0" />
<cell id="11360" material="10007" name="downcomer" region="10074 -10079 10082 -10081" universe="0" />
<cell id="11361" material="10004" name="reactor pressure vessel" region="10079 -10080 10082 -10081" universe="0" />
<lattice id="10124" name="Assembly (1.6\%) no BAs instr">
<lattice id="10124" name="Assembly (1.6%) no BAs instr">
<pitch>1.25984 1.25984</pitch>
<dimension>17 17</dimension>
<lower_left>-10.70864 -10.70864</lower_left>
@ -502,7 +502,7 @@
10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097
10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 </universes>
</lattice>
<lattice id="10144" name="Assembly (2.4\%) CR D">
<lattice id="10144" name="Assembly (2.4%) CR D">
<pitch>1.25984 1.25984</pitch>
<dimension>17 17</dimension>
<lower_left>-10.70864 -10.70864</lower_left>
@ -525,7 +525,7 @@
10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101
10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 10101 </universes>
</lattice>
<lattice id="10164" name="Assembly (3.1\%) 16BA">
<lattice id="10164" name="Assembly (3.1%) 16BA">
<pitch>1.25984 1.25984</pitch>
<dimension>17 17</dimension>
<lower_left>-10.70864 -10.70864</lower_left>
@ -548,7 +548,7 @@
10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105
10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 10105 </universes>
</lattice>
<lattice id="10182" name="Assembly (3.1\%) no BAs instr">
<lattice id="10182" name="Assembly (3.1%) no BAs instr">
<pitch>1.25984 1.25984</pitch>
<dimension>17 17</dimension>
<lower_left>-10.70864 -10.70864</lower_left>
@ -692,14 +692,14 @@
10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000
10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 </universes>
</lattice>
<surface coeffs="0.0 0.0 0.39218" id="10000" name="Pellet OR" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.405765" id="10000" name="Pellet OR" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.06459" id="10001" name="FR Plenum Spring OR" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.40005" id="10002" name="Clad IR" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.4572" id="10003" name="Clad OR" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.56134" id="10004" name="GT IR (above dashpot)" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.60198" id="10005" name="GT OR (above dashpot)" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.50419" id="10006" name="GT IR (at dashpot)" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.5461" id="10007" name="GT OR (at dashpot)" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.41402" id="10002" name="Clad IR" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.47498" id="10003" name="Clad OR" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.5715" id="10004" name="GT IR (above dashpot)" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.61214" id="10005" name="GT OR (above dashpot)" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.54019" id="10006" name="GT IR (at dashpot)" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.61214" id="10007" name="GT OR (at dashpot)" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.4331" id="10008" name="Control Poison OR" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.43688" id="10009" name="CR Clad IR" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.48387" id="10010" name="CR Clad OR" type="z-cylinder" />

8
smr/materials.xml
View file

@ -142,7 +142,7 @@
<nuclide name="Sn122" wo="0.0006894094798456136" />
<nuclide name="Sn124" wo="0.000876291604244001" />
</material>
<material id="10006" name="aic_rod">
<material id="10006" name="Ag-In-Cd">
<temperature>300</temperature>
<density units="g/cc" value="10.16" />
<nuclide name="Ag107" wo="0.4110094082785408" />
@ -181,7 +181,7 @@
<nuclide ao="0.007925516821765233" name="Si30" />
<nuclide ao="0.02423631977691551" name="Al27" />
</material>
<material id="10009" name="1.6\% Enr. UO2 Fuel">
<material id="10009" name="1.6% Enr. UO2 Fuel">
<temperature>300</temperature>
<density units="g/cc" value="10.31341" />
<nuclide ao="1.9998012531906832" name="O16" />
@ -190,7 +190,7 @@
<nuclide ao="0.016303176995310382" name="U235" />
<nuclide ao="0.9835658386532131" name="U238" />
</material>
<material id="10010" name="2.4\% Enr. UO2 Fuel">
<material id="10010" name="2.4% Enr. UO2 Fuel">
<temperature>300</temperature>
<density units="g/cc" value="10.29748" />
<nuclide ao="1.9998012531906837" name="O16" />
@ -199,7 +199,7 @@
<nuclide ao="0.024298776982208982" name="U235" />
<nuclide ao="0.9755059997465472" name="U238" />
</material>
<material id="10011" name="3.1\% Enr. UO2 Fuel">
<material id="10011" name="3.1% Enr. UO2 Fuel">
<temperature>300</temperature>
<density units="g/cc" value="10.30166" />
<nuclide ao="1.9998012531906824" name="O16" />

2
smr/plots.xml
View file

@ -22,7 +22,7 @@
<!--axial slice-->
<plot basis="xz" color="mat" filename="axial_xz_slice" id="10001" type="slice">
<origin>0.0 0.0 127.722</origin>
<width>127.722 127.722</width>
<width>270.0 255.444</width>
<pixels>1000 1000</pixels>
<background>255 255 255</background>
<col_spec id="10000" rgb="255 218 185" />

1
smr/settings.xml
View file

@ -13,7 +13,6 @@
<output>
<tallies>false</tallies>
</output>
<ptables>true</ptables>
<temperature_multipole>True</temperature_multipole>
<temperature_tolerance>1000</temperature_tolerance>
</settings>

115
smr/smr/assemblies.py
View file

@ -34,9 +34,8 @@ def make_assembly(name, universes):
lattice.pitch = [pin_pitch, pin_pitch]
lattice.universes = universes
# Create rectangular prism cylinder for lattice grid box
lat_grid_box = (surfs['lat grid box outer'] &
openmc.Complement(surfs['lat grid box inner']))
# Create rectangular prism for lattice grid box
lat_grid_box = (surfs['lat grid box outer'] & ~surfs['lat grid box inner'])
# Add lattice to bounding cell
univ_name = name + ' lattice'
@ -49,7 +48,7 @@ def make_assembly(name, universes):
# Add outer water cell
cell = openmc.Cell(name=univ_name + ' outer water')
cell.fill = univs['water pin']
cell.region = openmc.Complement(surfs['lat grid box outer'])
cell.region = ~surfs['lat grid box outer']
universe.add_cell(cell)
# Make bottom axial cell for outside of assembly (without sleeve)
@ -139,7 +138,7 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
# NO BURNABLE ABSORBERS
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (1.6\%) stack']
universes[:,:] = univs['Fuel (1.6%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, gtu,
gtu, gtu,
gtu, gtu, gtu, gtu, gtu,
@ -147,15 +146,15 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
gtu, gtu, gtu, gtu, gtu,
gtu, gtu,
gtu, gtu, gtu ]
univs['Assembly (1.6\%)' + comment] = \
make_assembly('Assembly (1.6\%) no BAs' + comment, universes)
univs['Assembly (1.6%)' + comment] = \
make_assembly('Assembly (1.6%) no BAs' + comment, universes)
# WITH EACH CONTROL ROD BANK
for bank, comment2 in [(crA, 'A'), (crB, 'B'), (crC, 'C'), (crD, 'D'),
(crSB, 'SB'), (crSC, 'SC'), (crSD, 'SD'), (crSE, 'SE')]:
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (1.6\%) stack']
universes[:,:] = univs['Fuel (1.6%) stack']
universes[nonfuel_y, nonfuel_x] = [ bank, bank, bank,
bank, bank,
bank, bank, bank, bank, bank,
@ -163,8 +162,8 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
bank, bank, bank, bank, bank,
bank, bank,
bank, bank, bank ]
univs['Assembly (1.6\%) CR {}'.format(comment2) + comment] = \
make_assembly('Assembly (1.6\%) CR {}'.format(comment2) + comment, universes)
univs['Assembly (1.6%) CR {}'.format(comment2) + comment] = \
make_assembly('Assembly (1.6%) CR {}'.format(comment2) + comment, universes)
#### 2.4% ENRICHED ASSEMBLIES
@ -173,7 +172,7 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
# NO BURNABLE ABSORBERS
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (2.4\%) stack']
universes[:,:] = univs['Fuel (2.4%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, gtu,
gtu, gtu,
gtu, gtu, gtu, gtu, gtu,
@ -181,12 +180,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
gtu, gtu, gtu, gtu, gtu,
gtu, gtu,
gtu, gtu, gtu ]
univs['Assembly (2.4\%) no BAs' + comment] = \
make_assembly('Assembly (2.4\%) no BAs' + comment, universes)
univs['Assembly (2.4%) no BAs' + comment] = \
make_assembly('Assembly (2.4%) no BAs' + comment, universes)
# WITH CONTROL ROD D BANK
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (2.4\%) stack']
universes[:,:] = univs['Fuel (2.4%) stack']
universes[nonfuel_y, nonfuel_x] = [ crD, crD, crD,
crD, crD,
crD, crD, crD, crD, crD,
@ -194,12 +193,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
crD, crD, crD, crD, crD,
crD, crD,
crD, crD, crD ]
univs['Assembly (2.4\%) CR D' + comment] = \
make_assembly('Assembly (2.4\%) CR D' + comment, universes)
univs['Assembly (2.4%) CR D' + comment] = \
make_assembly('Assembly (2.4%) CR D' + comment, universes)
# WITH 12 BURNABLE ABSORBERS
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (2.4\%) stack']
universes[:,:] = univs['Fuel (2.4%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, gtu, bas,
bas, bas,
bas, gtu, gtu, gtu, bas,
@ -207,12 +206,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
bas, gtu, gtu, gtu, bas,
bas, bas,
bas, gtu, bas ]
univs['Assembly (2.4\%) 12BA' + comment] = \
make_assembly('Assembly (2.4\%) 12BA' + comment, universes)
univs['Assembly (2.4%) 12BA' + comment] = \
make_assembly('Assembly (2.4%) 12BA' + comment, universes)
# WITH 16 BURNABLE ABSORBERS
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (2.4\%) stack']
universes[:,:] = univs['Fuel (2.4%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, bas, bas,
bas, bas,
bas, gtu, gtu, gtu, bas,
@ -220,8 +219,8 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
bas, gtu, gtu, gtu, bas,
bas, bas,
bas, bas, bas ]
univs['Assembly (2.4\%) 16BA' + comment] = \
make_assembly('Assembly (2.4\%) 16BA' + comment, universes)
univs['Assembly (2.4%) 16BA' + comment] = \
make_assembly('Assembly (2.4%) 16BA' + comment, universes)
#### 3.1% ENRICHED ASSEMBLIES
@ -230,7 +229,7 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
# NO BURNABLE ABSORBERS
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, gtu,
gtu, gtu,
gtu, gtu, gtu, gtu, gtu,
@ -238,12 +237,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
gtu, gtu, gtu, gtu, gtu,
gtu, gtu,
gtu, gtu, gtu ]
univs['Assembly (3.1\%)' + comment] = \
make_assembly('Assembly (3.1\%) no BAs' + comment, universes)
univs['Assembly (3.1%)' + comment] = \
make_assembly('Assembly (3.1%) no BAs' + comment, universes)
# WITH CONTROL ROD SA BANK
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ crSA, crSA, crSA,
crSA, crSA,
crSA, crSA, crSA, crSA, crSA,
@ -251,12 +250,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
crSA, crSA, crSA, crSA, crSA,
crSA, crSA,
crSA, crSA, crSA ]
univs['Assembly (3.1\%) CR SA' + comment] = \
make_assembly('Assembly (3.1\%) CR SA' + comment, universes)
univs['Assembly (3.1%) CR SA' + comment] = \
make_assembly('Assembly (3.1%) CR SA' + comment, universes)
# WITH 20 BURNABLE ABSORBERS
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, bas, bas,
bas, bas,
bas, bas, gtu, bas, bas,
@ -264,12 +263,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
bas, bas, gtu, bas, bas,
bas, bas,
bas, bas, bas ]
univs['Assembly (3.1\%) 20BA' + comment] = \
make_assembly('Assembly (3.1\%) 20BA' + comment, universes)
univs['Assembly (3.1%) 20BA' + comment] = \
make_assembly('Assembly (3.1%) 20BA' + comment, universes)
# WITH 16 BURNABLE ABSORBERS
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, bas, bas,
bas, bas,
bas, gtu, gtu, gtu, bas,
@ -277,12 +276,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
bas, gtu, gtu, gtu, bas,
bas, bas,
bas, bas, bas ]
univs['Assembly (3.1\%) 16BA' + comment] = \
make_assembly('Assembly (3.1\%) 16BA' + comment, universes)
univs['Assembly (3.1%) 16BA' + comment] = \
make_assembly('Assembly (3.1%) 16BA' + comment, universes)
# WITH 15 BURNABLE ABSORBERS NW
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, gtu,
gtu, gtu,
gtu, bas, bas, bas, bas,
@ -290,12 +289,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
gtu, bas, bas, bas, bas,
gtu, bas,
bas, bas, bas ]
univs['Assembly (3.1\%) 15BANW' + comment] = \
make_assembly('Assembly (3.1\%) 15BANW' + comment, universes)
univs['Assembly (3.1%) 15BANW' + comment] = \
make_assembly('Assembly (3.1%) 15BANW' + comment, universes)
# WITH 15 BURNABLE ABSORBERS NE
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, gtu,
gtu, gtu,
gtu, bas, bas, bas, gtu,
@ -303,12 +302,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
bas, bas, bas, bas, gtu,
bas, gtu,
bas, bas, bas ]
univs['Assembly (3.1\%) 15BANE' + comment] = \
make_assembly('Assembly (3.1\%) 15BANE' + comment, universes)
univs['Assembly (3.1%) 15BANE' + comment] = \
make_assembly('Assembly (3.1%) 15BANE' + comment, universes)
# WITH 15 BURNABLE ABSORBERS SW
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, bas, bas,
gtu, bas,
gtu, bas, bas, bas, bas,
@ -316,12 +315,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
gtu, bas, bas, bas, bas,
gtu, gtu,
gtu, gtu, gtu ]
univs['Assembly (3.1\%) 15BASW' + comment] = \
make_assembly('Assembly (3.1\%) 15BASW' + comment, universes)
univs['Assembly (3.1%) 15BASW' + comment] = \
make_assembly('Assembly (3.1%) 15BASW' + comment, universes)
# WITH 15 BURNABLE ABSORBERS SE
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, bas, bas,
bas, gtu,
bas, bas, bas, bas, gtu,
@ -329,12 +328,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
bas, bas, bas, bas, gtu,
gtu, gtu,
gtu, gtu, gtu ]
univs['Assembly (3.1\%) 15BASE' + comment] = \
make_assembly('Assembly (3.1\%) 15BASE' + comment, universes)
univs['Assembly (3.1%) 15BASE' + comment] = \
make_assembly('Assembly (3.1%) 15BASE' + comment, universes)
# WITH 6 BURNABLE ABSORBERS N
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, gtu,
gtu, gtu,
gtu, gtu, gtu, gtu, gtu,
@ -342,12 +341,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
bas, gtu, gtu, gtu, bas,
bas, bas,
bas, gtu, bas ]
univs['Assembly (3.1\%) 6BAN' + comment] = \
make_assembly('Assembly (3.1\%) 6BAN' + comment, universes)
univs['Assembly (3.1%) 6BAN' + comment] = \
make_assembly('Assembly (3.1%) 6BAN' + comment, universes)
# WITH 6 BURNABLE ABSORBERS S
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, gtu, bas,
bas, bas,
bas, gtu, gtu, gtu, bas,
@ -355,12 +354,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
gtu, gtu, gtu, gtu, gtu,
gtu, gtu,
gtu, gtu, gtu ]
univs['Assembly (3.1\%) 6BAS' + comment] = \
make_assembly('Assembly (3.1\%) 6BAS' + comment, universes)
univs['Assembly (3.1%) 6BAS' + comment] = \
make_assembly('Assembly (3.1%) 6BAS' + comment, universes)
# WITH 6 BURNABLE ABSORBERS W
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, bas,
gtu, bas,
gtu, gtu, gtu, gtu, bas,
@ -368,12 +367,12 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
gtu, gtu, gtu, gtu, bas,
gtu, bas,
gtu, gtu, bas ]
univs['Assembly (3.1\%) 6BAW' + comment] = \
make_assembly('Assembly (3.1\%) 6BAW' + comment, universes)
univs['Assembly (3.1%) 6BAW' + comment] = \
make_assembly('Assembly (3.1%) 6BAW' + comment, universes)
# WITH 6 BURNABLE ABSORBERS E
universes = np.empty((17,17), dtype=openmc.Universe)
universes[:,:] = univs['Fuel (3.1\%) stack']
universes[:,:] = univs['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, gtu, gtu,
bas, gtu,
bas, gtu, gtu, gtu, gtu,
@ -381,5 +380,5 @@ for cent, comment in [(gti, ''), (ins, ' instr')]:
bas, gtu, gtu, gtu, gtu,
bas, gtu,
bas, gtu, gtu ]
univs['Assembly (3.1\%) 6BAE' + comment] = \
make_assembly('Assembly (3.1\%) 6BAE' + comment, universes)
univs['Assembly (3.1%) 6BAE' + comment] = \
make_assembly('Assembly (3.1%) 6BAE' + comment, universes)

77
smr/smr/core.py
View file

@ -26,67 +26,67 @@ core.universes[5, 11] = univs['baffle northeast corner']
core.universes[6, 6] = univs['baffle northwest corner']
core.universes[6, 7] = univs['baffle northwest']
core.universes[6, 8] = univs['Assembly (3.1\%) instr']
core.universes[6, 9] = univs['Assembly (2.4\%) CR D']
core.universes[6, 10] = univs['Assembly (3.1\%) instr']
core.universes[6, 8] = univs['Assembly (3.1%) instr']
core.universes[6, 9] = univs['Assembly (2.4%) CR D']
core.universes[6, 10] = univs['Assembly (3.1%) instr']
core.universes[6, 11] = univs['baffle northeast']
core.universes[6, 12] = univs['baffle northeast corner']
core.universes[7, 5] = univs['baffle northwest corner']
core.universes[7, 6] = univs['baffle northwest']
core.universes[7, 7] = univs['Assembly (3.1\%) instr']
core.universes[7, 8] = univs['Assembly (2.4\%) CR D']
core.universes[7, 9] = univs['Assembly (3.1\%) 16BA']
core.universes[7, 10] = univs['Assembly (2.4\%) CR D']
core.universes[7, 11] = univs['Assembly (3.1\%) instr']
core.universes[7, 7] = univs['Assembly (3.1%) instr']
core.universes[7, 8] = univs['Assembly (2.4%) CR D']
core.universes[7, 9] = univs['Assembly (3.1%) 16BA']
core.universes[7, 10] = univs['Assembly (2.4%) CR D']
core.universes[7, 11] = univs['Assembly (3.1%) instr']
core.universes[7, 12] = univs['baffle northeast']
core.universes[7, 13] = univs['baffle northeast corner']
core.universes[8, 5] = univs['baffle west']
core.universes[8, 6] = univs['Assembly (3.1\%) instr']
core.universes[8, 7] = univs['Assembly (2.4\%) CR D']
core.universes[8, 8] = univs['Assembly (3.1\%) 16BA']
core.universes[8, 9] = univs['Assembly (2.4\%) CR D']
core.universes[8, 10] = univs['Assembly (3.1\%) 16BA']
core.universes[8, 11] = univs['Assembly (2.4\%) CR D']
core.universes[8, 12] = univs['Assembly (3.1\%) instr']
core.universes[8, 6] = univs['Assembly (3.1%) instr']
core.universes[8, 7] = univs['Assembly (2.4%) CR D']
core.universes[8, 8] = univs['Assembly (3.1%) 16BA']
core.universes[8, 9] = univs['Assembly (2.4%) CR D']
core.universes[8, 10] = univs['Assembly (3.1%) 16BA']
core.universes[8, 11] = univs['Assembly (2.4%) CR D']
core.universes[8, 12] = univs['Assembly (3.1%) instr']
core.universes[8, 13] = univs['baffle east']
core.universes[9, 5] = univs['baffle west']
core.universes[9, 6] = univs['Assembly (2.4\%) CR D']
core.universes[9, 7] = univs['Assembly (3.1\%) 16BA']
core.universes[9, 8] = univs['Assembly (2.4\%) CR D']
core.universes[9, 9] = univs['Assembly (1.6\%) instr']
core.universes[9, 10] = univs['Assembly (2.4\%) CR D']
core.universes[9, 11] = univs['Assembly (3.1\%) 16BA']
core.universes[9, 12] = univs['Assembly (2.4\%) CR D']
core.universes[9, 6] = univs['Assembly (2.4%) CR D']
core.universes[9, 7] = univs['Assembly (3.1%) 16BA']
core.universes[9, 8] = univs['Assembly (2.4%) CR D']
core.universes[9, 9] = univs['Assembly (1.6%) instr']
core.universes[9, 10] = univs['Assembly (2.4%) CR D']
core.universes[9, 11] = univs['Assembly (3.1%) 16BA']
core.universes[9, 12] = univs['Assembly (2.4%) CR D']
core.universes[9, 13] = univs['baffle east']
core.universes[10, 5] = univs['baffle west']
core.universes[10, 6] = univs['Assembly (3.1\%) instr']
core.universes[10, 7] = univs['Assembly (2.4\%) CR D']
core.universes[10, 8] = univs['Assembly (3.1\%) 16BA']
core.universes[10, 9] = univs['Assembly (2.4\%) CR D']
core.universes[10, 10] = univs['Assembly (3.1\%) 16BA']
core.universes[10, 11] = univs['Assembly (2.4\%) CR D']
core.universes[10, 12] = univs['Assembly (3.1\%) instr']
core.universes[10, 6] = univs['Assembly (3.1%) instr']
core.universes[10, 7] = univs['Assembly (2.4%) CR D']
core.universes[10, 8] = univs['Assembly (3.1%) 16BA']
core.universes[10, 9] = univs['Assembly (2.4%) CR D']
core.universes[10, 10] = univs['Assembly (3.1%) 16BA']
core.universes[10, 11] = univs['Assembly (2.4%) CR D']
core.universes[10, 12] = univs['Assembly (3.1%) instr']
core.universes[10, 13] = univs['baffle east']
core.universes[11, 5] = univs['baffle southwest corner']
core.universes[11, 6] = univs['baffle southwest']
core.universes[11, 7] = univs['Assembly (3.1\%) instr']
core.universes[11, 8] = univs['Assembly (2.4\%) CR D']
core.universes[11, 9] = univs['Assembly (3.1\%) 16BA']
core.universes[11, 10] = univs['Assembly (2.4\%) CR D']
core.universes[11, 11] = univs['Assembly (3.1\%) instr']
core.universes[11, 7] = univs['Assembly (3.1%) instr']
core.universes[11, 8] = univs['Assembly (2.4%) CR D']
core.universes[11, 9] = univs['Assembly (3.1%) 16BA']
core.universes[11, 10] = univs['Assembly (2.4%) CR D']
core.universes[11, 11] = univs['Assembly (3.1%) instr']
core.universes[11, 12] = univs['baffle southeast']
core.universes[11, 13] = univs['baffle southeast corner']
core.universes[12, 6] = univs['baffle southwest corner']
core.universes[12, 7] = univs['baffle southwest']
core.universes[12, 8] = univs['Assembly (3.1\%) instr']
core.universes[12, 9] = univs['Assembly (2.4\%) CR D']
core.universes[12, 10] = univs['Assembly (3.1\%) instr']
core.universes[12, 8] = univs['Assembly (3.1%) instr']
core.universes[12, 9] = univs['Assembly (2.4%) CR D']
core.universes[12, 10] = univs['Assembly (3.1%) instr']
core.universes[12, 11] = univs['baffle southeast']
core.universes[12, 12] = univs['baffle southeast corner']
@ -204,5 +204,4 @@ root_univ.add_cell(cell)
#### CONSTRUCT GEOMETRY
geometry = openmc.Geometry()
geometry.root_universe = root_univ
geometry = openmc.Geometry(root_univ)

13
smr/smr/materials.py
View file

@ -122,7 +122,7 @@ mats['Zr'].add_element('Zr', 0.98115, 'wo')
mats['Zr'].add_element('Sn', 0.0145, 'wo')
# Create Ag-In-Cd control rod material
mats['AIC'] = openmc.Material(name='aic_rod')
mats['AIC'] = openmc.Material(name='Ag-In-Cd')
mats['AIC'].temperature = 300
mats['AIC'].set_density('g/cc', 10.16)
mats['AIC'].add_element('Ag', 0.80, 'wo')
@ -167,7 +167,7 @@ mats['H2O'].set_density('g/cc', rho_Bh2o)
mats['H2O'].add_element('B', aB_Bh2o, 'ao')
mats['H2O'].add_element('H', ah_Bh2o, 'ao')
mats['H2O'].add_element('O', aho_Bh2o, 'ao')
mats['H2O'].add_s_alpha_beta(name='lwtr')
mats['H2O'].add_s_alpha_beta(name='c_H_in_H2O')
#### Borosilicate Glass
@ -211,7 +211,7 @@ mats['BSG'].add_nuclide('B11', aB11_B, 'ao')
# Create 1.6% enriched UO2 fuel material
a_U234, a_U235, a_U238, a_U, a_O = get_fuel_aos(0.0161006)
mats['UO2 1.6'] = openmc.Material(name='1.6\% Enr. UO2 Fuel')
mats['UO2 1.6'] = openmc.Material(name='1.6% Enr. UO2 Fuel')
mats['UO2 1.6'].temperature = 300
mats['UO2 1.6'].set_density('g/cc', 10.31341)
mats['UO2 1.6'].add_element('O', a_O, 'ao')
@ -219,7 +219,7 @@ mats['UO2 1.6'].add_element('U', a_U, 'ao', enrichment=1.61006)
# Create 2.4% enriched UO2 fuel material
a_U234, a_U235, a_U238, a_U, a_O = get_fuel_aos(0.0239993)
mats['UO2 2.4'] = openmc.Material(name='2.4\% Enr. UO2 Fuel')
mats['UO2 2.4'] = openmc.Material(name='2.4% Enr. UO2 Fuel')
mats['UO2 2.4'].temperature = 300
mats['UO2 2.4'].set_density('g/cc', 10.29748)
mats['UO2 2.4'].add_element('O', a_O, 'ao')
@ -227,7 +227,7 @@ mats['UO2 2.4'].add_element('U', a_U, 'ao', enrichment=2.39993)
# Create 3.1% enriched UO2 fuel material
a_U234, a_U235, a_U238, a_U, a_O = get_fuel_aos(0.0310221)
mats['UO2 3.1'] = openmc.Material(name='3.1\% Enr. UO2 Fuel')
mats['UO2 3.1'] = openmc.Material(name='3.1% Enr. UO2 Fuel')
mats['UO2 3.1'].temperature = 300
mats['UO2 3.1'].set_density('g/cc', 10.30166)
mats['UO2 3.1'].add_element('O', a_O, 'ao')
@ -236,5 +236,4 @@ mats['UO2 3.1'].add_element('U', a_U, 'ao', enrichment=3.10221)
# Construct a collection of Materials to export to XML
materials = openmc.Materials()
materials.add_materials(mats.values())
materials = openmc.Materials(mats.values())

125
smr/smr/pins.py
View file

@ -33,9 +33,7 @@ def make_pin(name, surfaces, materials, grid=None):
# Create cell for interior of innermost ZCylinder
cell_name = name + ' (0)'
cell = openmc.Cell(name=cell_name)
cell.fill = materials[0]
cell.region = -surfaces[0]
cell = openmc.Cell(name=cell_name, fill=materials[0], region=-surfaces[0])
universe.add_cell(cell)
# Create cells between two ZCylinders
@ -48,19 +46,16 @@ def make_pin(name, surfaces, materials, grid=None):
# Create cell for exterior of outermost ZCylinder
cell_name = name + ' (last)'
cell = openmc.Cell(name=cell_name)
cell.fill = materials[-1]
cell.region = +surfaces[-1]
cell = openmc.Cell(name=cell_name, fill=materials[-1], region=+surfaces[-1])
universe.add_cell(cell)
# Add spacer grid cells if specified
if grid:
grid_name = 'rod grid box {}'.format(grid)
cell.region = openmc.Intersection(cell.region, surfs[grid_name])
cell.region &= surfs[grid_name]
cell_name = name + ' (grid)'
cell = openmc.Cell(name=cell_name)
cell.region = openmc.Complement(surfs[grid_name])
cell = openmc.Cell(name=cell_name, region=~surfs[grid_name])
if grid == 'top/bottom':
cell.fill = mats['In']
@ -96,24 +91,20 @@ def make_stack(name, surfaces, universes):
# Create cell for interior of innermost ZCylinder
cell_name = name + ' (0)'
cell = openmc.Cell(name=cell_name)
cell.fill = universes[0]
cell.region = -surfaces[0]
cell = openmc.Cell(name=cell_name, fill=universes[0], region=-surfaces[0])
universe.add_cell(cell)
# Create cells between two ZCylinders
for i, (univ, surf) in enumerate(zip(universes[1:-1], surfaces[:-1])):
cell_name = name + ' ({})'.format(i+1)
cell = openmc.Cell(name=cell_name)
cell.fill = univ
cell.region = +surf & -surfaces[i+1]
cell = openmc.Cell(name=cell_name, fill=univ,
region=+surf & -surfaces[i+1])
universe.add_cell(cell)
# Create cell for exterior of outermost ZCylinder
cell_name = name + ' (last)'
cell = openmc.Cell(name=cell_name)
cell.fill = universes[-1]
cell.region = +surfaces[-1]
cell = openmc.Cell(name=cell_name, fill=universes[-1],
region=+surfaces[-1])
universe.add_cell(cell)
return universe
@ -321,7 +312,7 @@ univs['CR bare'] = make_pin(
# Stack all axial pieces of control rod tubes together for each bank
banks = ['A','B','C','D','SA','SB','SC','SD','SE']
banks = ['A', 'B', 'C', 'D', 'SA', 'SB', 'SC', 'SD', 'SE']
for b in banks:
@ -663,8 +654,8 @@ univs['pin plenum grid (top/bottom)'] = make_pin(
#### 1.6% ENRICHED FUEL PIN CELL
univs['Fuel (1.6\%)'] = make_pin(
'Fuel (1.6\%)',
univs['Fuel (1.6%)'] = make_pin(
'Fuel (1.6%)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
@ -673,8 +664,8 @@ univs['Fuel (1.6\%)'] = make_pin(
mats['Zr'],
mats['H2O']])
univs['Fuel (1.6\%) grid (top/bottom)'] = make_pin(
'Fuel (1.6\%) grid (top/bottom)',
univs['Fuel (1.6%) grid (top/bottom)'] = make_pin(
'Fuel (1.6%) grid (top/bottom)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
@ -684,8 +675,8 @@ univs['Fuel (1.6\%) grid (top/bottom)'] = make_pin(
mats['H2O']],
grid='(top/bottom)')
univs['Fuel (1.6\%) grid (intermediate)'] = make_pin(
'Fuel (1.6\%) grid (intermediate)',
univs['Fuel (1.6%) grid (intermediate)'] = make_pin(
'Fuel (1.6%) grid (intermediate)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
@ -697,21 +688,21 @@ univs['Fuel (1.6\%) grid (intermediate)'] = make_pin(
# Stack all axial pieces of 1.6% enriched fuel pin cell
univs['Fuel (1.6\%) stack'] = make_stack(
'Fuel (1.6\%) stack',
univs['Fuel (1.6%) stack'] = make_stack(
'Fuel (1.6%) stack',
surfaces=stack_surfs,
universes=[univs['water pin'],
univs['SS pin'],
univs['SS pin'],
univs['end plug'],
univs['Fuel (1.6\%)'],
univs['Fuel (1.6\%) grid (top/bottom)'],
univs['Fuel (1.6\%)'],
univs['Fuel (1.6\%)'],
univs['Fuel (1.6\%) grid (intermediate)'],
univs['Fuel (1.6\%)'],
univs['Fuel (1.6\%) grid (intermediate)'],
univs['Fuel (1.6\%)'],
univs['Fuel (1.6%)'],
univs['Fuel (1.6%) grid (top/bottom)'],
univs['Fuel (1.6%)'],
univs['Fuel (1.6%)'],
univs['Fuel (1.6%) grid (intermediate)'],
univs['Fuel (1.6%)'],
univs['Fuel (1.6%) grid (intermediate)'],
univs['Fuel (1.6%)'],
univs['pin plenum'],
univs['pin plenum grid (top/bottom)'],
univs['pin plenum'],
@ -723,8 +714,8 @@ univs['Fuel (1.6\%) stack'] = make_stack(
#### 2.4% ENRICHED FUEL PIN CELL
univs['Fuel (2.4\%)'] = make_pin(
'Fuel (2.4\%)',
univs['Fuel (2.4%)'] = make_pin(
'Fuel (2.4%)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
@ -733,8 +724,8 @@ univs['Fuel (2.4\%)'] = make_pin(
mats['Zr'],
mats['H2O']])
univs['Fuel (2.4\%) grid (top/bottom)'] = make_pin(
'Fuel (2.4\%) grid (top/bottom)',
univs['Fuel (2.4%) grid (top/bottom)'] = make_pin(
'Fuel (2.4%) grid (top/bottom)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
@ -744,8 +735,8 @@ univs['Fuel (2.4\%) grid (top/bottom)'] = make_pin(
mats['H2O']],
grid='(top/bottom)')
univs['Fuel (2.4\%) grid (intermediate)'] = make_pin(
'Fuel (2.4\%) grid (intermediate)',
univs['Fuel (2.4%) grid (intermediate)'] = make_pin(
'Fuel (2.4%) grid (intermediate)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
@ -757,21 +748,21 @@ univs['Fuel (2.4\%) grid (intermediate)'] = make_pin(
# Stack all axial pieces of 2.4% enriched fuel pin cell
univs['Fuel (2.4\%) stack'] = make_stack(
'Fuel (2.4\%) stack',
univs['Fuel (2.4%) stack'] = make_stack(
'Fuel (2.4%) stack',
surfaces=stack_surfs,
universes=[univs['water pin'],
univs['SS pin'],
univs['SS pin'],
univs['end plug'],
univs['Fuel (2.4\%)'],
univs['Fuel (2.4\%) grid (top/bottom)'],
univs['Fuel (2.4\%)'],
univs['Fuel (2.4\%)'],
univs['Fuel (2.4\%) grid (intermediate)'],
univs['Fuel (2.4\%)'],
univs['Fuel (2.4\%) grid (intermediate)'],
univs['Fuel (2.4\%)'],
univs['Fuel (2.4%)'],
univs['Fuel (2.4%) grid (top/bottom)'],
univs['Fuel (2.4%)'],
univs['Fuel (2.4%)'],
univs['Fuel (2.4%) grid (intermediate)'],
univs['Fuel (2.4%)'],
univs['Fuel (2.4%) grid (intermediate)'],
univs['Fuel (2.4%)'],
univs['pin plenum'],
univs['pin plenum grid (top/bottom)'],
univs['pin plenum'],
@ -783,8 +774,8 @@ univs['Fuel (2.4\%) stack'] = make_stack(
#### 3.1% ENRICHED FUEL PIN CELL
univs['Fuel (3.1\%)'] = make_pin(
'Fuel (3.1\%)',
univs['Fuel (3.1%)'] = make_pin(
'Fuel (3.1%)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
@ -793,8 +784,8 @@ univs['Fuel (3.1\%)'] = make_pin(
mats['Zr'],
mats['H2O']])
univs['Fuel (3.1\%) grid (top/bottom)'] = make_pin(
'Fuel (3.1\%) grid (top/bottom)',
univs['Fuel (3.1%) grid (top/bottom)'] = make_pin(
'Fuel (3.1%) grid (top/bottom)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
@ -804,8 +795,8 @@ univs['Fuel (3.1\%) grid (top/bottom)'] = make_pin(
mats['H2O']],
grid='(top/bottom)')
univs['Fuel (3.1\%) grid (intermediate)'] = make_pin(
'Fuel (3.1\%) grid (intermediate)',
univs['Fuel (3.1%) grid (intermediate)'] = make_pin(
'Fuel (3.1%) grid (intermediate)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
@ -817,21 +808,21 @@ univs['Fuel (3.1\%) grid (intermediate)'] = make_pin(
# Stack all axial pieces of 3.1% enriched fuel pin cell
univs['Fuel (3.1\%) stack'] = make_stack(
'Fuel (3.1\%) stack',
univs['Fuel (3.1%) stack'] = make_stack(
'Fuel (3.1%) stack',
surfaces=stack_surfs,
universes=[univs['water pin'],
univs['SS pin'],
univs['SS pin'],
univs['end plug'],
univs['Fuel (3.1\%)'],
univs['Fuel (3.1\%) grid (top/bottom)'],
univs['Fuel (3.1\%)'],
univs['Fuel (3.1\%)'],
univs['Fuel (3.1\%) grid (intermediate)'],
univs['Fuel (3.1\%)'],
univs['Fuel (3.1\%) grid (intermediate)'],
univs['Fuel (3.1\%)'],
univs['Fuel (3.1%)'],
univs['Fuel (3.1%) grid (top/bottom)'],
univs['Fuel (3.1%)'],
univs['Fuel (3.1%)'],
univs['Fuel (3.1%) grid (intermediate)'],
univs['Fuel (3.1%)'],
univs['Fuel (3.1%) grid (intermediate)'],
univs['Fuel (3.1%)'],
univs['pin plenum'],
univs['pin plenum grid (top/bottom)'],
univs['pin plenum'],

4
smr/smr/plots.py
View file

@ -2,7 +2,7 @@
import openmc
from .surfaces import lowest_extent, highest_extent, lattice_pitch
from .surfaces import lowest_extent, highest_extent, lattice_pitch, rpv_OR
from .materials import mats
@ -39,7 +39,7 @@ plot = openmc.Plot(name='axial slice')
plot.basis = 'xz'
plot.color = 'mat'
plot.origin = [0., 0., (highest_extent-lowest_extent)/2.]
plot.width = [(highest_extent-lowest_extent)/2.] * 2
plot.width = [rpv_OR*2., (highest_extent-lowest_extent)]
plot.filename = 'axial_xz_slice'
plot.col_spec = col_spec
plot.background = [255, 255, 255]

28
smr/smr/surfaces.py
View file

@ -4,7 +4,7 @@ The geometric parameters defining the core model are tabulated here.
"""
import copy
import math
from math import tan, pi
import openmc
@ -21,15 +21,15 @@ import openmc
# RPV: Reactor Pressure Vessel
# pin cell parameters
pellet_OR = 0.39218
clad_IR = 0.40005
clad_OR = 0.45720
pellet_OR = 0.405765
clad_IR = 0.41402
clad_OR = 0.47498
rod_grid_side_tb = 1.24416
rod_grid_side_i = 1.21962
guide_tube_IR = 0.56134
guide_tube_OR = 0.60198
guide_tube_dash_IR = 0.50419
guide_tube_dash_OR = 0.54610
guide_tube_IR = 0.5715
guide_tube_OR = 0.61214
guide_tube_dash_IR = 0.54019
guide_tube_dash_OR = 0.61214
control_poison_OR = 0.43310
control_rod_IR = 0.43688
control_rod_OR = 0.48387
@ -41,8 +41,8 @@ burn_abs_r5 = 0.43688
burn_abs_r6 = 0.48387
burn_abs_r7 = 0.56134
burn_abs_r8 = 0.60198
instr_tube_IR = 0.43688
instr_tube_OR = 0.48387
instr_tube_IR = 0.5715
instr_tube_OR = 0.61214
plenum_spring_OR = 0.06459
# grid spacer parameters
@ -99,10 +99,10 @@ bank_bot = 405.713
bank_step = 228.
bank_top = 766.348
neutron_shield_NWbot_SEtop = math.tan(math.pi/3)
neutron_shield_NWtop_SEbot = math.tan(math.pi/6)
neutron_shield_NEbot_SWtop = math.tan(-math.pi/3)
neutron_shield_NEtop_SWbot = math.tan(-math.pi/6)
neutron_shield_NWbot_SEtop = tan(pi/3)
neutron_shield_NWtop_SEbot = tan(pi/6)
neutron_shield_NEbot_SWtop = tan(-pi/3)
neutron_shield_NEtop_SWbot = tan(-pi/6)
surfs = {}