From e3161efa571581d458a41eebda3e6aa5b5a42763 Mon Sep 17 00:00:00 2001 From: yardasol Date: Wed, 16 Nov 2022 14:51:30 -0600 Subject: [PATCH] IA elements are actually IB elements --- examples/msbr/core_elements.py | 66 +++++++++++++++++----------------- 1 file changed, 33 insertions(+), 33 deletions(-) diff --git a/examples/msbr/core_elements.py b/examples/msbr/core_elements.py index 985adc25d..464e026a6 100644 --- a/examples/msbr/core_elements.py +++ b/examples/msbr/core_elements.py @@ -2,7 +2,7 @@ import numpy as np def _bound_zone_cells(cells_tuples, levels): - """Helper function that moves Zone IA and Zone IIA cells to the + """Helper function that moves Zone IB and Zone IIA cells to the appropriate height.""" cell_list = [] n_levels = len(cells_tuples) @@ -26,14 +26,14 @@ def _bound_zone_cells(cells_tuples, levels): cell_list.append(cell) return cell_list -def zoneIA(gr_sq_neg, +def zoneIB(gr_sq_neg, gr_extra_regions, inter_elem_channel, gr_round_4, moder, fuel, hast): - """Create universe for Zone IA element. Based on specification in + """Create universe for Zone IB element. Based on specification in Robertson, 1971 Fig 3.4 (p. 17) Parameters @@ -58,74 +58,74 @@ def zoneIA(gr_sq_neg, Returns ------- - ia : openmc.Universe - Universe for Zone IA element. + ib : openmc.Universe + Universe for Zone IB element. """ elem_levels = [22.86, 419.10, 438.15, 445.135] level_bounds = [] for level in elem_levels: level_bounds.append(openmc.ZPlane(z0=level)) - s1 = openmc.ZCylinder(r=4.953, name='ia_gr_round_1') - s2 = openmc.ZCylinder(r=1.71069, name='ia_fuel_hole') + s1 = openmc.ZCylinder(r=4.953, name='ib_gr_round_1') + s2 = openmc.ZCylinder(r=1.71069, name='ib_fuel_hole') h = 12.66 theta = np.arctan(4.953 / h) r2 = (1 / np.cos(theta))**2 - 1 s3 = openmc.ZCone(z0=h + elem_levels[2], r2=r2, name='cone_i') - c1 = openmc.Cell(fill=fuel, region=(-s2), name='ia_fuel_inner_1') - c2 = openmc.Cell(fill=moder, region=(+s2 & -s1), name='ia_moderator_1') - c3 = openmc.Cell(fill=fuel, region=(+s1), name='ia_fuel_outer_1') - ia1 = (c1, c2, c3) + c1 = openmc.Cell(fill=fuel, region=(-s2), name='ib_fuel_inner_1') + c2 = openmc.Cell(fill=moder, region=(+s2 & -s1), name='ib_moderator_1') + c3 = openmc.Cell(fill=fuel, region=(+s1), name='ib_fuel_outer_1') + ib1 = (c1, c2, c3) # I-A main (lower 2) c4 = c1.clone(clone_materials=False) - c4.name = 'ia_fuel_inner_main' + c4.name = 'ib_fuel_inner_main' c5 = openmc.Cell(fill=moder, region=(+s2 & gr_sq_neg & inter_elem_channel), - name='ia_moderator_main') + name='ib_moderator_main') c6 = openmc.Cell(fill=fuel, region=(~gr_sq_neg & inter_elem_channel), - name='ia_fuel_outer_main') - iam = [c4, c5, c6] + name='ib_fuel_outer_main') + ibm = [c4, c5, c6] for (reg, name) in gr_extra_regions: - iam.append(openmc.Cell(fill=moder, region=reg, - name=f'ia_moderator_main_{name}')) + ibm.append(openmc.Cell(fill=moder, region=reg, + name=f'ib_moderator_main_{name}')) # I-A 2 (upper 1) c7 = c1.clone(clone_materials=False) c8 = c2.clone(clone_materials=False) c9 = c3.clone(clone_materials=False) - c7.name = 'ia_fuel_inner_2' - c8.name = 'ia_moderator_2' - c9.name = 'ia_fuel_outer_2' + c7.name = 'ib_fuel_inner_2' + c8.name = 'ib_moderator_2' + c9.name = 'ib_fuel_outer_2' - ia2 = (c7, c8, c9) + ib2 = (c7, c8, c9) # I-A 3 (upper 2) c10 = c1.clone(clone_materials=False) - c10.name = 'ia_fuel_inner_3' - c11 = openmc.Cell(fill=moder, region=(+s2 & -s3), name='ia_moderator_3') - c12 = openmc.Cell(fill=fuel, region=(+s3), name='ia_fuel_outer_3') - ia3 = (c10, c11, c12) + c10.name = 'ib_fuel_inner_3' + c11 = openmc.Cell(fill=moder, region=(+s2 & -s3), name='ib_moderator_3') + c12 = openmc.Cell(fill=fuel, region=(+s3), name='ib_fuel_outer_3') + ib3 = (c10, c11, c12) # I-A 4 (upper 3) c13 = openmc.Cell(fill=hast, region=(-s2), - name='ia_hast') + name='ib_hast') c14 = openmc.Cell(fill=moder, region=(+s2 & -gr_round_4), - name='ia_moderator_4') + name='ib_moderator_4') c15 = openmc.Cell(fill=fuel, region=(+gr_round_4), - name='ia_fuel_outer_4') - ia4 = (c13, c14, c15) + name='ib_fuel_outer_4') + ib4 = (c13, c14, c15) - elem_cells = [ia1, iam, ia2, ia3, ia4] + elem_cells = [ib1, ibm, ib2, ib3, ib4] # universe_id=1 - ia = openmc.Universe(name='zone_ia') - ia.add_cells(_bound_zone_cells(elem_cells, level_bounds)) - return ia + ib = openmc.Universe(name='zone_ib') + ib.add_cells(_bound_zone_cells(elem_cells, level_bounds)) + return ib def zoneIIA(gr_sq_neg, gr_extra_regions,