diff --git a/smr/fresh/geometry.xml b/smr/fresh/geometry.xml index 68df767..44e70df 100644 --- a/smr/fresh/geometry.xml +++ b/smr/fresh/geometry.xml @@ -439,19 +439,36 @@ - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + 1.25984 1.25984 17 17 @@ -544,20 +561,20 @@ 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 106 - + 21.503639999999997 21.503639999999997 9 9 -96.76637999999998 -96.76637999999998 -207 207 207 207 207 207 207 207 207 -207 207 207 184 146 184 207 207 207 -207 207 184 146 166 146 184 207 207 -207 184 146 166 146 166 146 184 207 -207 146 166 146 126 146 166 146 207 -207 184 146 166 146 166 146 184 207 -207 207 184 146 166 146 184 207 207 -207 207 207 184 146 184 207 207 207 -207 207 207 207 207 207 207 207 207 +211 211 211 211 211 211 211 211 211 +211 211 207 184 146 184 208 211 211 +211 207 184 146 166 146 184 208 211 +211 184 146 166 146 166 146 184 211 +211 146 166 146 126 146 166 146 211 +211 184 146 166 146 166 146 184 211 +211 209 184 146 166 146 184 210 211 +211 211 209 184 146 184 210 211 211 +211 211 211 211 211 211 211 211 211 @@ -624,4 +641,17 @@ + + + + + + + + + + + + + diff --git a/smr/smr/baffle.py b/smr/smr/baffle.py deleted file mode 100644 index 9cc849e..0000000 --- a/smr/smr/baffle.py +++ /dev/null @@ -1,403 +0,0 @@ -"""Instantiate components to use for the stainless steel baffle.""" - -import openmc - -from .materials import mats -from .surfaces import surfs, lattice_pitch -from .assemblies import univs - - -#### HEAVY REFLECTOR - -cell = openmc.Cell( - name='heavy reflector SS', - fill=mats['SS'] -) - -univs['heavy reflector'] = openmc.Universe( - name='heavy reflector', - cells=[cell] -) - -#### SOUTH BAFFLE - -univs['baffle south dummy'] = openmc.Universe(name='baffle south dummy') - -cell = openmc.Cell(name='baffle dummy SS') -cell.fill = mats['SS'] -cell.region = +surfs['baffle south'] -univs['baffle south dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle south water') -cell.fill = mats['H2O'] -cell.region = -surfs['baffle south'] -univs['baffle south dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle south') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['baffle south dummy'], univs['baffle south dummy']], - [univs['water pin'], univs['water pin']]] - -univs['baffle south'] = openmc.Universe(name='baffle south') - -cell = openmc.Cell(name='baffle south') -cell.fill = lattice -univs['baffle south'].add_cell(cell) - - -#### NORTH BAFFLE - -univs['baffle north dummy'] = openmc.Universe(name='baffle north dummy') - -cell = openmc.Cell(name='baffle dummy SS') -cell.fill = mats['H2O'] -cell.region = +surfs['baffle north'] -univs['baffle north dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle north water') -cell.fill = mats['SS'] -cell.region = -surfs['baffle north'] -univs['baffle north dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle north') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['water pin'], univs['water pin']], - [univs['baffle north dummy'], univs['baffle north dummy']]] - -univs['baffle north'] = openmc.Universe(name='baffle north') - -cell = openmc.Cell(name='baffle north') -cell.fill = lattice -univs['baffle north'].add_cell(cell) - - -#### WEST BAFFLE - -univs['baffle west dummy'] = openmc.Universe(name='baffle west dummy') - -cell = openmc.Cell(name='baffle dummy SS') -cell.fill = mats['SS'] -cell.region = +surfs['baffle west'] -univs['baffle west dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle west water') -cell.fill = mats['H2O'] -cell.region = -surfs['baffle west'] -univs['baffle west dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle west') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['water pin'], univs['baffle west dummy']], - [univs['water pin'], univs['baffle west dummy']]] - -univs['baffle west'] = openmc.Universe(name='baffle west') - -cell = openmc.Cell(name='baffle west') -cell.fill = lattice -univs['baffle west'].add_cell(cell) - - -#### EAST BAFFLE - -univs['baffle east dummy'] = openmc.Universe(name='baffle east dummy') - -cell = openmc.Cell(name='baffle dummy SS') -cell.fill = mats['H2O'] -cell.region = +surfs['baffle east'] -univs['baffle east dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle east water') -cell.fill = mats['SS'] -cell.region = -surfs['baffle east'] -univs['baffle east dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle east') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['baffle east dummy'], univs['water pin']], - [univs['baffle east dummy'], univs['water pin']]] - -univs['baffle east'] = openmc.Universe(name='baffle east') - -cell = openmc.Cell(name='baffle east') -cell.fill = lattice -univs['baffle east'].add_cell(cell) - - -#### SOUTHEAST BAFFLE - -univs['baffle southeast dummy'] = openmc.Universe(name='baffle southeast dummy') - -cell = openmc.Cell(name='baffle southeast dummy 1') -cell.fill = mats['H2O'] -cell.region = +surfs['baffle east'] & -surfs['baffle south'] -univs['baffle southeast dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle southeast dummy 2') -cell.fill = mats['SS'] -cell.region = +surfs['baffle east'] & +surfs['baffle south'] -univs['baffle southeast dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle southeast dummy 3') -cell.fill = mats['SS'] -cell.region = -surfs['baffle east'] -univs['baffle southeast dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle southeast') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['baffle southeast dummy'], univs['baffle south dummy']], - [univs['baffle east dummy'], univs['water pin']]] - -univs['baffle southeast'] = openmc.Universe(name='baffle southeast') - -cell = openmc.Cell(name='baffle southeast') -cell.fill = lattice -univs['baffle southeast'].add_cell(cell) - - -#### SOUTHWEST BAFFLE - -univs['baffle southwest dummy'] = openmc.Universe(name='baffle southwest dummy') - -cell = openmc.Cell(name='baffle southwest dummy 1') -cell.fill = mats['H2O'] -cell.region = -surfs['baffle south'] & -surfs['baffle west'] -univs['baffle southwest dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle southwest dummy 2') -cell.fill = mats['SS'] -cell.region = +surfs['baffle south'] & -surfs['baffle west'] -univs['baffle southwest dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle southwest dummy 3') -cell.fill = mats['SS'] -cell.region = +surfs['baffle west'] -univs['baffle southwest dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle southwest') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['baffle south dummy'], univs['baffle southwest dummy']], - [univs['water pin'], univs['baffle west dummy']]] - -univs['baffle southwest'] = openmc.Universe(name='baffle southwest') - -cell = openmc.Cell(name='baffle southwest') -cell.fill = lattice -univs['baffle southwest'].add_cell(cell) - - -#### NORTHEAST BAFFLE - -univs['baffle northeast dummy'] = openmc.Universe(name='baffle northeast dummy') - -cell = openmc.Cell(name='baffle northeast dummy 1') -cell.fill = mats['H2O'] -cell.region = +surfs['baffle north'] & +surfs['baffle east'] -univs['baffle northeast dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle northeast dummy 2') -cell.fill = mats['SS'] -cell.region = -surfs['baffle north'] & +surfs['baffle east'] -univs['baffle northeast dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle northeast dummy 3') -cell.fill = mats['SS'] -cell.region = -surfs['baffle east'] -univs['baffle northeast dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle northeast') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['baffle east dummy'], univs['water pin']], - [univs['baffle northeast dummy'], univs['baffle north dummy']]] - -univs['baffle northeast'] = openmc.Universe(name='baffle northeast') - -cell = openmc.Cell(name='baffle northeast') -cell.fill = lattice -univs['baffle northeast'].add_cell(cell) - - -#### NORTHWEST BAFFLE - -univs['baffle northwest dummy'] = openmc.Universe(name='baffle northwest dummy') - -cell = openmc.Cell(name='baffle northwest dummy 1') -cell.fill = mats['H2O'] -cell.region = +surfs['baffle north'] & -surfs['baffle west'] -univs['baffle northwest dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle northwest dummy 2') -cell.fill = mats['SS'] -cell.region = -surfs['baffle north'] & -surfs['baffle west'] -univs['baffle northwest dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle northwest dummy 3') -cell.fill = mats['SS'] -cell.region = +surfs['baffle west'] -univs['baffle northwest dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle northwest') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['water pin'], univs['baffle west dummy']], - [univs['baffle north dummy'], univs['baffle northwest dummy']]] - -univs['baffle northwest'] = openmc.Universe(name='baffle northwest') - -cell = openmc.Cell(name='baffle northwest') -cell.fill = lattice -univs['baffle northwest'].add_cell(cell) - - -#### SOUTHEAST CORNER BAFFLE - -univs['baffle southeast corner dummy'] = \ - openmc.Universe(name='baffle southeast corner dummy') - -cell = openmc.Cell(name='baffle southeast corner dummy 1') -cell.fill = mats['H2O'] -cell.region = -surfs['baffle east'] & -surfs['baffle south'] -univs['baffle southeast corner dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle southeast corner dummy 2') -cell.fill = mats['H2O'] -cell.region = +surfs['baffle east'] -univs['baffle southeast corner dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle southeast corner dummy 3') -cell.fill = mats['SS'] -cell.region = +surfs['baffle south'] & -surfs['baffle east'] -univs['baffle southeast corner dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle southeast corner') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['baffle southeast corner dummy'], univs['water pin']], - [univs['water pin'], univs['water pin']]] - -univs['baffle southeast corner'] = \ - openmc.Universe(name='baffle southeast corner') - -cell = openmc.Cell(name='baffle southeast corner') -cell.fill = lattice -univs['baffle southeast corner'].add_cell(cell) - - -#### SOUTHWEST CORNER BAFFLE - -univs['baffle southwest corner dummy'] = \ - openmc.Universe(name='baffle southwest corner dummy') - -cell = openmc.Cell(name='baffle southwest corner dummy 1') -cell.fill = mats['H2O'] -cell.region = +surfs['baffle west'] & -surfs['baffle south'] -univs['baffle southwest corner dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle southwest corner dummy 2') -cell.fill = mats['H2O'] -cell.region = -surfs['baffle west'] -univs['baffle southwest corner dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle southwest corner dummy 3') -cell.fill = mats['SS'] -cell.region = +surfs['baffle south'] & +surfs['baffle west'] -univs['baffle southwest corner dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle southwest corner') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['water pin'], univs['baffle southwest corner dummy']], - [univs['water pin'], univs['water pin']]] - -univs['baffle southwest corner'] = \ - openmc.Universe(name='baffle southwest corner') - -cell = openmc.Cell(name='baffle southwest corner') -cell.fill = lattice -univs['baffle southwest corner'].add_cell(cell) - - -#### NORTHWEST CORNER BAFFLE - -univs['baffle northwest corner dummy'] = \ - openmc.Universe(name='baffle northwest corner dummy') - -cell = openmc.Cell(name='baffle northwest corner dummy 1') -cell.fill = mats['H2O'] -cell.region = +surfs['baffle west'] & +surfs['baffle north'] -univs['baffle northwest corner dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle northwest corner dummy 2') -cell.fill = mats['H2O'] -cell.region = -surfs['baffle west'] -univs['baffle northwest corner dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle northwest corner dummy 3') -cell.fill = mats['SS'] -cell.region = -surfs['baffle north'] & +surfs['baffle west'] -univs['baffle northwest corner dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle northwest corner') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['water pin'], univs['water pin']], - [univs['water pin'], univs['baffle northwest corner dummy']]] - -univs['baffle northwest corner'] = \ - openmc.Universe(name='baffle northwest corner') - -cell = openmc.Cell(name='baffle northwest corner') -cell.fill = lattice -univs['baffle northwest corner'].add_cell(cell) - - -#### NORTHEAST CORNER BAFFLE - -univs['baffle northeast corner dummy'] = \ - openmc.Universe(name='baffle northeast corner dummy') - -cell = openmc.Cell(name='baffle northeast corner dummy 1') -cell.fill = mats['H2O'] -cell.region = -surfs['baffle east'] & +surfs['baffle north'] -univs['baffle northeast corner dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle northeast corner dummy 2') -cell.fill = mats['H2O'] -cell.region = +surfs['baffle east'] -univs['baffle northeast corner dummy'].add_cell(cell) - -cell = openmc.Cell(name='baffle northeast corner dummy 3') -cell.fill = mats['SS'] -cell.region = -surfs['baffle north'] & -surfs['baffle east'] -univs['baffle northeast corner dummy'].add_cell(cell) - -lattice = openmc.RectLattice(name='baffle northeast corner') -lattice.lower_left = [-lattice_pitch/2., -lattice_pitch/2.] -lattice.pitch = [lattice_pitch/2., lattice_pitch/2.] -lattice.universes = [ - [univs['water pin'], univs['water pin']], - [univs['baffle northeast corner dummy'], univs['water pin']]] - -univs['baffle northeast corner'] = \ - openmc.Universe(name='baffle northeast corner') - -cell = openmc.Cell(name='baffle northeast corner') -cell.fill = lattice -univs['baffle northeast corner'].add_cell(cell) diff --git a/smr/smr/core.py b/smr/smr/core.py index 5984985..1515553 100644 --- a/smr/smr/core.py +++ b/smr/smr/core.py @@ -6,7 +6,7 @@ import openmc from .materials import mats from .surfaces import surfs, lattice_pitch -from .baffle import univs +from .reflector import univs #### CONSTRUCT MAIN CORE LATTICE @@ -23,21 +23,21 @@ universes[0, 5] = univs['heavy reflector'] universes[0, 6] = univs['heavy reflector'] universes[1, 1] = univs['heavy reflector'] -universes[1, 2] = univs['heavy reflector'] +universes[1, 2] = univs['heavy reflector NW'] universes[1, 3] = univs['Assembly (3.1%) instr'] universes[1, 4] = univs['Assembly (2.4%) CR D'] universes[1, 5] = univs['Assembly (3.1%) instr'] -universes[1, 6] = univs['heavy reflector'] +universes[1, 6] = univs['heavy reflector NE'] universes[1, 7] = univs['heavy reflector'] universes[2, 0] = univs['heavy reflector'] -universes[2, 1] = univs['heavy reflector'] +universes[2, 1] = univs['heavy reflector NW'] universes[2, 2] = univs['Assembly (3.1%) instr'] universes[2, 3] = univs['Assembly (2.4%) CR D'] universes[2, 4] = univs['Assembly (3.1%) 16BA'] universes[2, 5] = univs['Assembly (2.4%) CR D'] universes[2, 6] = univs['Assembly (3.1%) instr'] -universes[2, 7] = univs['heavy reflector'] +universes[2, 7] = univs['heavy reflector NE'] universes[2, 8] = univs['heavy reflector'] universes[3, 0] = univs['heavy reflector'] @@ -71,21 +71,21 @@ universes[5, 7] = univs['Assembly (3.1%) instr'] universes[5, 8] = univs['heavy reflector'] universes[6, 0] = univs['heavy reflector'] -universes[6, 1] = univs['heavy reflector'] +universes[6, 1] = univs['heavy reflector SW'] universes[6, 2] = univs['Assembly (3.1%) instr'] universes[6, 3] = univs['Assembly (2.4%) CR D'] universes[6, 4] = univs['Assembly (3.1%) 16BA'] universes[6, 5] = univs['Assembly (2.4%) CR D'] universes[6, 6] = univs['Assembly (3.1%) instr'] -universes[6, 7] = univs['heavy reflector'] +universes[6, 7] = univs['heavy reflector SE'] universes[6, 8] = univs['heavy reflector'] universes[7, 1] = univs['heavy reflector'] -universes[7, 2] = univs['heavy reflector'] +universes[7, 2] = univs['heavy reflector SW'] universes[7, 3] = univs['Assembly (3.1%) instr'] universes[7, 4] = univs['Assembly (2.4%) CR D'] universes[7, 5] = univs['Assembly (3.1%) instr'] -universes[7, 6] = univs['heavy reflector'] +universes[7, 6] = univs['heavy reflector SE'] universes[7, 7] = univs['heavy reflector'] universes[8, 2] = univs['heavy reflector'] diff --git a/smr/smr/reflector.py b/smr/smr/reflector.py new file mode 100644 index 0000000..4fa9ae0 --- /dev/null +++ b/smr/smr/reflector.py @@ -0,0 +1,116 @@ +"""Instantiate components to use for the stainless steel heavy reflector.""" + +import openmc + +from .materials import mats +from .surfaces import surfs, lattice_pitch +from .assemblies import univs + + +#### HEAVY REFLECTOR + +# All the dimensions of the water holes in the heavy reflectors were eyeballed +# from DC, Figure 4.3-25 + +# Pixel widths from the Figure -- we'll convert these to actual coordinates +width = 276 +p1 = 59 +p2 = 126 +p3 = 196 +p4 = 264 + +p5 = 105 + +p6 = 122 +p7 = 164 + +p8 = 138 +p9 = 222 + +p10 = 247 + +d_small = 13 +d_large = 30 + +scale = lattice_pitch/width + +# Physical positions +x1 = -lattice_pitch/2 + scale*(width - p4) +x2 = -lattice_pitch/2 + scale*(width - p3) +x3 = -lattice_pitch/2 + scale*(width - p2) +x4 = -lattice_pitch/2 + scale*(width - p1) +y1 = -lattice_pitch/2 + scale*p1 +y2 = -lattice_pitch/2 + scale*p2 +y3 = -lattice_pitch/2 + scale*p3 +y4 = -lattice_pitch/2 + scale*p4 + +x5 = -lattice_pitch/2 + scale*(width - p5) +y5 = -lattice_pitch/2 + scale*p5 + +x6 = -lattice_pitch/2 + scale*(width - p7) +y6 = -lattice_pitch/2 + scale*p6 + +x7 = -lattice_pitch/2 + scale*(width - p6) +y7 = -lattice_pitch/2 + scale*p7 + +x8 = -lattice_pitch/2 + scale*(width - p9) +y8 = -lattice_pitch/2 + scale*p8 + +x9 = -lattice_pitch/2 + scale*(width - p8) +y9 = -lattice_pitch/2 + scale*p9 + +y10 = -lattice_pitch/2 + scale*p10 + +r1 = scale*d_small/2 +r2 = scale*d_large/2 + +dims = [ + (x1, y1, r1), (x2, y1, r1), (x3, y1, r1), (x4, y1, r2), + (x4, y2, r1), (x4, y3, r1), (x4, y4, r1), (x5, y5, r1), + (x6, y6, r1), (x7, y7, r1), (x8, y8, r1), (x9, y9, r1), + (x1, y10, r1) +] + +water_holes = [] +for x, y, r in dims: + zcyl = openmc.ZCylinder(x0=x, y0=y, R=r) + hole = openmc.Cell(fill=mats['H2O'], region=-zcyl) + water_holes.append(hole) + +ss_region = openmc.Intersection(~c.region for c in water_holes) + +ss_cell = openmc.Cell(name='heavy reflector SS', fill=mats['SS'], + region=ss_region) + +univs['heavy reflector NW'] = openmc.Universe(name='heavy reflector NW') +univs['heavy reflector NW'].add_cells(water_holes) +univs['heavy reflector NW'].add_cell(ss_cell) + + +# NE corner + +cell = openmc.Cell(name='heavy reflector NE', fill=univs['heavy reflector NW']) +cell.rotation = (0, 180, 0) +univs['heavy reflector NE'] = openmc.Universe(name='heavy reflector NE') +univs['heavy reflector NE'].add_cell(cell) + +# SW corner + +cell = openmc.Cell(name='heavy reflector SW', fill=univs['heavy reflector NW']) +cell.rotation = (180, 0, 0) +univs['heavy reflector SW'] = openmc.Universe(name='heavy reflector SW') +univs['heavy reflector SW'].add_cell(cell) + +# SE corner + +cell = openmc.Cell(name='heavy reflector SE', fill=univs['heavy reflector NW']) +cell.rotation = (0, 0, 180) +univs['heavy reflector SE'] = openmc.Universe(name='heavy reflector SE') +univs['heavy reflector SE'].add_cell(cell) + + +# Solid stainless steel universe + +all_ss = openmc.Cell(name='heavy reflector', fill=mats['SS']) +univs['heavy reflector'] = openmc.Universe( + name='heavy reflector', cells=[all_ss]) diff --git a/smr/smr/surfaces.py b/smr/smr/surfaces.py index 2b36c18..908b302 100644 --- a/smr/smr/surfaces.py +++ b/smr/smr/surfaces.py @@ -284,13 +284,3 @@ surfs['upper bound'] = openmc.ZPlane( z0=highest_extent, name='upper problem boundary', boundary_type='vacuum') surfs['lower bound'] = openmc.ZPlane( z0=lowest_extent, name='lower problem boundary', boundary_type='vacuum') - -# baffle surfaces -surfs['baffle south'] = openmc.YPlane( - y0=(lattice_pitch/4. - baffle_width), name='baffle north') -surfs['baffle north'] = openmc.YPlane( - y0=(baffle_width - lattice_pitch/4.), name='baffle south') -surfs['baffle west'] = openmc.XPlane( - x0=(lattice_pitch/4. - baffle_width), name='baffle east') -surfs['baffle east'] = openmc.XPlane( - x0=(baffle_width - lattice_pitch/4), name='baffle west') \ No newline at end of file