Fix axial segmentation of fuel pin

This commit is contained in:
Paul Romano 2019-06-04 22:25:45 -05:00
parent 99efb29643
commit b400c1ca21

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):
@ -703,7 +703,7 @@ def pin_universes(num_rings=10, num_axial=196, depleted=False):
if num_axial > 1:
# Determine z position between each fuel pellet, omitting the surfaces
# corresponding to the very bottom and top of the active fuel length
axial_splits = np.linspace(bottom_fuel_rod, top_active_core, num_axial + 1)[1:-1]
axial_splits = np.linspace(bottom_fuel_stack, top_active_core, num_axial + 1)[1:-1]
axial_surfs = [openmc.ZPlane(z0=z) for z in axial_splits]
if num_rings > 1: