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

Remove assemblies with burnable absorber pins

Browse source
This commit is contained in:
Paul Romano 2019-12-19 15:00:03 -06:00
parent c375e6d9cc
commit 0b3546f23a
3 changed files with 17 additions and 367 deletions
Download patch file Download diff file Expand all files Collapse all files

157
smr/smr/assemblies.py
View file

@ -147,7 +147,6 @@ def assembly_universes(ring_radii, num_axial, depleted):
# commonly needed universes
gtu = pins['GT empty stack']
gti = pins['GT empty instr']
bas = pins['BA stack']
ins = pins['IT stack']
crA = pins['GT CR bank A']
crB = pins['GT CR bank B']
@ -232,33 +231,6 @@ def assembly_universes(ring_radii, num_axial, depleted):
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[:,:] = pins['Fuel (2.4%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, gtu, bas,
bas, bas,
bas, gtu, gtu, gtu, bas,
gtu, gtu, cent, gtu, gtu,
bas, gtu, gtu, gtu, bas,
bas, bas,
bas, gtu, bas ]
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[:,:] = pins['Fuel (2.4%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, bas, bas,
bas, bas,
bas, gtu, gtu, gtu, bas,
bas, gtu, cent, gtu, bas,
bas, gtu, gtu, gtu, bas,
bas, bas,
bas, bas, bas ]
univs['Assembly (2.4%) 16BA' + comment] = \
make_assembly('Assembly (2.4%) 16BA' + comment, universes)
#### 3.1% ENRICHED ASSEMBLIES
for cent, comment in [(gti, ''), (ins, ' instr')]:
@ -289,134 +261,5 @@ def assembly_universes(ring_radii, num_axial, depleted):
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[:,:] = pins['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, bas, bas,
bas, bas,
bas, bas, gtu, bas, bas,
bas, gtu, cent, gtu, bas,
bas, bas, gtu, bas, bas,
bas, bas,
bas, bas, bas ]
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[:,:] = pins['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, bas, bas,
bas, bas,
bas, gtu, gtu, gtu, bas,
bas, gtu, cent, gtu, bas,
bas, gtu, gtu, gtu, bas,
bas, bas,
bas, bas, bas ]
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[:,:] = pins['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, gtu,
gtu, gtu,
gtu, bas, bas, bas, bas,
gtu, bas, cent, bas, bas,
gtu, bas, bas, bas, bas,
gtu, bas,
bas, bas, bas ]
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[:,:] = pins['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, gtu,
gtu, gtu,
bas, bas, bas, bas, gtu,
bas, bas, cent, bas, gtu,
bas, bas, bas, bas, gtu,
bas, gtu,
bas, bas, bas ]
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[:,:] = pins['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, bas, bas,
gtu, bas,
gtu, bas, bas, bas, bas,
gtu, bas, cent, bas, bas,
gtu, bas, bas, bas, bas,
gtu, gtu,
gtu, gtu, gtu ]
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[:,:] = pins['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, bas, bas,
bas, gtu,
bas, bas, bas, bas, gtu,
bas, bas, cent, bas, gtu,
bas, bas, bas, bas, gtu,
gtu, gtu,
gtu, gtu, gtu ]
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[:,:] = pins['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, gtu,
gtu, gtu,
gtu, gtu, gtu, gtu, gtu,
gtu, gtu, cent, gtu, gtu,
bas, gtu, gtu, gtu, bas,
bas, bas,
bas, gtu, bas ]
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[:,:] = pins['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, gtu, bas,
bas, bas,
bas, gtu, gtu, gtu, bas,
gtu, gtu, cent, gtu, gtu,
gtu, gtu, gtu, gtu, gtu,
gtu, gtu,
gtu, gtu, gtu ]
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[:,:] = pins['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ gtu, gtu, bas,
gtu, bas,
gtu, gtu, gtu, gtu, bas,
gtu, gtu, cent, gtu, gtu,
gtu, gtu, gtu, gtu, bas,
gtu, bas,
gtu, gtu, bas ]
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[:,:] = pins['Fuel (3.1%) stack']
universes[nonfuel_y, nonfuel_x] = [ bas, gtu, gtu,
bas, gtu,
bas, gtu, gtu, gtu, gtu,
gtu, gtu, cent, gtu, gtu,
bas, gtu, gtu, gtu, gtu,
bas, gtu,
bas, gtu, gtu ]
univs['Assembly (3.1%) 6BAE' + comment] = \
make_assembly('Assembly (3.1%) 6BAE' + comment, universes)
return univs

34
smr/smr/core.py
View file

@ -46,9 +46,9 @@ def core_geometry(ring_radii, num_axial, depleted):
universes[1, 1] = reflector['1,1']
universes[1, 2] = reflector['NW']
universes[1, 3] = assembly['Assembly (3.1%) instr']
universes[1, 3] = assembly['Assembly (3.1%)']
universes[1, 4] = assembly['Assembly (2.4%) CR D']
universes[1, 5] = assembly['Assembly (3.1%) instr']
universes[1, 5] = assembly['Assembly (3.1%)']
universes[1, 6] = reflector['NE']
universes[1, 7] = reflector['1,7']
@ -56,47 +56,47 @@ def core_geometry(ring_radii, num_axial, depleted):
universes[2, 1] = reflector['NW']
universes[2, 2] = assembly['Assembly (3.1%) instr']
universes[2, 3] = assembly['Assembly (2.4%) CR D']
universes[2, 4] = assembly['Assembly (3.1%) 16BA']
universes[2, 4] = assembly['Assembly (3.1%) instr']
universes[2, 5] = assembly['Assembly (2.4%) CR D']
universes[2, 6] = assembly['Assembly (3.1%) instr']
universes[2, 7] = reflector['NE']
universes[2, 8] = reflector['2,8']
universes[3, 0] = reflector['3,0']
universes[3, 1] = assembly['Assembly (3.1%) instr']
universes[3, 1] = assembly['Assembly (3.1%)']
universes[3, 2] = assembly['Assembly (2.4%) CR D']
universes[3, 3] = assembly['Assembly (3.1%) 16BA']
universes[3, 3] = assembly['Assembly (3.1%) instr']
universes[3, 4] = assembly['Assembly (2.4%) CR D']
universes[3, 5] = assembly['Assembly (3.1%) 16BA']
universes[3, 5] = assembly['Assembly (3.1%) instr']
universes[3, 6] = assembly['Assembly (2.4%) CR D']
universes[3, 7] = assembly['Assembly (3.1%) instr']
universes[3, 7] = assembly['Assembly (3.1%)']
universes[3, 8] = reflector['3,8']
universes[4, 0] = reflector['4,0']
universes[4, 1] = assembly['Assembly (2.4%) CR D']
universes[4, 2] = assembly['Assembly (3.1%) 16BA']
universes[4, 2] = assembly['Assembly (3.1%) instr']
universes[4, 3] = assembly['Assembly (2.4%) CR D']
universes[4, 4] = assembly['Assembly (1.6%) instr']
universes[4, 4] = assembly['Assembly (1.6%)']
universes[4, 5] = assembly['Assembly (2.4%) CR D']
universes[4, 6] = assembly['Assembly (3.1%) 16BA']
universes[4, 6] = assembly['Assembly (3.1%) instr']
universes[4, 7] = assembly['Assembly (2.4%) CR D']
universes[4, 8] = reflector['4,8']
universes[5, 0] = reflector['5,0']
universes[5, 1] = assembly['Assembly (3.1%) instr']
universes[5, 1] = assembly['Assembly (3.1%)']
universes[5, 2] = assembly['Assembly (2.4%) CR D']
universes[5, 3] = assembly['Assembly (3.1%) 16BA']
universes[5, 3] = assembly['Assembly (3.1%) instr']
universes[5, 4] = assembly['Assembly (2.4%) CR D']
universes[5, 5] = assembly['Assembly (3.1%) 16BA']
universes[5, 5] = assembly['Assembly (3.1%) instr']
universes[5, 6] = assembly['Assembly (2.4%) CR D']
universes[5, 7] = assembly['Assembly (3.1%) instr']
universes[5, 7] = assembly['Assembly (3.1%)']
universes[5, 8] = reflector['5,8']
universes[6, 0] = reflector['6,0']
universes[6, 1] = reflector['SW']
universes[6, 2] = assembly['Assembly (3.1%) instr']
universes[6, 3] = assembly['Assembly (2.4%) CR D']
universes[6, 4] = assembly['Assembly (3.1%) 16BA']
universes[6, 4] = assembly['Assembly (3.1%) instr']
universes[6, 5] = assembly['Assembly (2.4%) CR D']
universes[6, 6] = assembly['Assembly (3.1%) instr']
universes[6, 7] = reflector['SE']
@ -104,9 +104,9 @@ def core_geometry(ring_radii, num_axial, depleted):
universes[7, 1] = reflector['7,1']
universes[7, 2] = reflector['SW']
universes[7, 3] = assembly['Assembly (3.1%) instr']
universes[7, 3] = assembly['Assembly (3.1%)']
universes[7, 4] = assembly['Assembly (2.4%) CR D']
universes[7, 5] = assembly['Assembly (3.1%) instr']
universes[7, 5] = assembly['Assembly (3.1%)']
universes[7, 6] = reflector['SE']
universes[7, 7] = reflector['7,7']

193
smr/smr/pins.py
View file

@ -473,199 +473,6 @@ def pin_universes(ring_radii=None, num_axial=196, depleted=False):
univs['GT CR bank {} dummy bare'.format(b)]])
#### BURNABLE ABSORBER PIN CELLS
univs['BA'] = make_pin(
'BA',
surfaces=[surfs['BA IR 1'],
surfs['BA IR 2'],
surfs['BA IR 3'],
surfs['BA IR 4'],
surfs['BA IR 5'],
surfs['BA IR 6'],
surfs['BA IR 7'],
surfs['BA IR 8']],
materials=[mats['Air'],
mats['SS'],
mats['Air'],
mats['BSG'],
mats['Air'],
mats['SS'],
mats['H2O'],
mats['Zr'],
mats['H2O']])
univs['BA grid (bottom)'] = make_pin(
'BA grid (bottom)',
surfaces=[surfs['BA IR 1'],
surfs['BA IR 2'],
surfs['BA IR 3'],
surfs['BA IR 4'],
surfs['BA IR 5'],
surfs['BA IR 6'],
surfs['BA IR 7'],
surfs['BA IR 8']],
materials=[mats['Air'],
mats['SS'],
mats['Air'],
mats['BSG'],
mats['Air'],
mats['SS'],
mats['H2O'],
mats['Zr'],
mats['H2O']],
grid='bottom')
univs['BA grid (intermediate)'] = make_pin(
'BA grid (intermediate)',
surfaces=[surfs['BA IR 1'],
surfs['BA IR 2'],
surfs['BA IR 3'],
surfs['BA IR 4'],
surfs['BA IR 5'],
surfs['BA IR 6'],
surfs['BA IR 7'],
surfs['BA IR 8']],
materials=[mats['Air'],
mats['SS'],
mats['Air'],
mats['BSG'],
mats['Air'],
mats['SS'],
mats['H2O'],
mats['Zr'],
mats['H2O']],
grid='intermediate')
univs['BA dashpot'] = make_pin(
'BA dashpot',
surfaces=[surfs['BA IR 1'],
surfs['BA IR 2'],
surfs['BA IR 3'],
surfs['BA IR 4'],
surfs['BA IR 5'],
surfs['BA IR 6'],
surfs['GT dashpot IR'],
surfs['GT dashpot OR']],
materials=[mats['Air'],
mats['SS'],
mats['Air'],
mats['BSG'],
mats['Air'],
mats['SS'],
mats['H2O'],
mats['Zr'],
mats['H2O']])
univs['BA dashpot grid (bottom)'] = make_pin(
'BA dashpot grid (bottom)',
surfaces=[surfs['BA IR 1'],
surfs['BA IR 2'],
surfs['BA IR 3'],
surfs['BA IR 4'],
surfs['BA IR 5'],
surfs['BA IR 6'],
surfs['GT dashpot IR'],
surfs['GT dashpot OR']],
materials=[mats['Air'],
mats['SS'],
mats['Air'],
mats['BSG'],
mats['Air'],
mats['SS'],
mats['H2O'],
mats['Zr'],
mats['H2O']],
grid='bottom')
univs['BA dashpot grid (intermediate)'] = make_pin(
'BA dashpot grid (intermediate)',
surfaces=[surfs['BA IR 1'],
surfs['BA IR 2'],
surfs['BA IR 3'],
surfs['BA IR 4'],
surfs['BA IR 5'],
surfs['BA IR 6'],
surfs['GT dashpot IR'],
surfs['GT dashpot OR']],
materials=[mats['Air'],
mats['SS'],
mats['Air'],
mats['BSG'],
mats['Air'],
mats['SS'],
mats['H2O'],
mats['Zr'],
mats['H2O']],
grid='intermediate')
univs['BA blank SS'] = make_pin(
'BA blank SS',
surfaces=[surfs['BA IR 6'],
surfs['BA IR 7'],
surfs['BA IR 8']],
materials=[mats['SS'],
mats['H2O'],
mats['Zr'],
mats['H2O']])
univs['BA blank SS bare'] = make_pin(
'BA blank SS bare',
surfaces=[surfs['BA IR 6']],
materials=[mats['SS'],
mats['H2O']])
stack_surfs_BA = [
surfs['bot support plate'],
surfs['top support plate'],
surfs['top lower nozzle'],
surfs['top lower thimble'],
surfs['grid1bot'],
surfs['BA bot'],
surfs['grid1top'],
surfs['dashpot top'],
surfs['grid2bot'],
surfs['grid2top'],
surfs['grid3bot'],
surfs['grid3top'],
surfs['grid4bot'],
surfs['grid4top'],
surfs['top active core'],
surfs['grid5bot'],
surfs['grid5top'],
surfs['top pin plenum'],
surfs['top FR'],
surfs['bot upper nozzle'],
surfs['top upper nozzle']]
# Stack all axial pieces of control rod tubes together for each bank
univs['BA stack'] = make_stack(
'BA stack', stack_surfs_BA,
universes=[univs['water pin'],
univs['water pin'],
univs['water pin'],
univs['GTd empty'],
univs['GTd empty'],
univs['GTd empty grid (bottom)'],
univs['BA dashpot grid (bottom)'],
univs['BA dashpot'],
univs['BA'],
univs['BA grid (intermediate)'],
univs['BA'],
univs['BA grid (intermediate)'],
univs['BA'],
univs['BA grid (intermediate)'],
univs['BA'],
univs['BA blank SS'],
univs['BA blank SS'],
univs['BA blank SS'],
univs['BA blank SS'],
univs['BA blank SS'],
univs['BA blank SS bare'],
univs['water pin']])
# Fuel pin cells
univs['SS pin'] = make_pin(
'SS pin',