crp/ecp-benchmarks
1
0
Fork 1
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Fix subdivided fuel (position of lower z-plane, add outer universe)

Browse source
This commit is contained in:
Paul Romano 2018-05-02 22:04:10 -05:00
parent 4af7da2ea7
commit 2b175b4411
2 changed files with 5 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

7
smr/smr/pins.py
View file

@ -7,7 +7,7 @@ import openmc
from openmc.model import subdivide
from .materials import mats
from .surfaces import surfs, pellet_OR, bottom_fuel_rod, top_active_core
from .surfaces import surfs, pellet_OR, bottom_fuel_stack, top_active_core
def make_pin(name, surfaces, materials, grid=None):
@ -723,10 +723,11 @@ def pin_universes(num_rings=10, num_axial=196, depleted=False):
# performance cost of neighbor lookups
if num_axial > 1:
zlattice = openmc.RectLattice()
zlattice.lower_left = (-100., -100., bottom_fuel_rod)
zlattice.pitch = (200., 200., (top_active_core - bottom_fuel_rod)/num_axial)
zlattice.lower_left = (-100., -100., bottom_fuel_stack)
zlattice.pitch = (200., 200., (top_active_core - bottom_fuel_stack)/num_axial)
zlattice.dimension = (1, 1, num_axial)
zlattice.universes = np.full((num_axial, 1, 1), univ)
zlattice.outer = univ
return zlattice
else:
return univ

2
smr/smr/surfaces.py
View file

@ -192,7 +192,7 @@ surfs['bot active core'] = openmc.ZPlane(
surfs['top active core'] = openmc.ZPlane(
z0=top_active_core, name='top active core')
surfs['top lower thimble'] = copy.deepcopy(surfs['bot active core'])
surfs['top lower thimble'] = surfs['bot active core']
surfs['BA bot'] = openmc.ZPlane(
z0=bot_burn_abs, name='bottom of BA')