diff --git a/smr/core.py b/smr/core.py index e12c97d..8befdf9 100644 --- a/smr/core.py +++ b/smr/core.py @@ -25,7 +25,7 @@ core.universes[5, 11] = univs['baffle southwest corner'] core.universes[6, 6] = univs['baffle southeast corner'] core.universes[6, 7] = univs['baffle southeast'] core.universes[6, 8] = univs['Assembly (3.1\%) instr'] -core.universes[6, 9] = univs['Assembly (1.6\%) CR B'] +core.universes[6, 9] = univs['Assembly (2.4\%) CR D'] core.universes[6, 10] = univs['Assembly (3.1\%) instr'] core.universes[6, 11] = univs['baffle southwest'] core.universes[6, 12] = univs['baffle southwest corner'] @@ -33,49 +33,49 @@ core.universes[6, 12] = univs['baffle southwest corner'] core.universes[7, 5] = univs['baffle southeast corner'] core.universes[7, 6] = univs['baffle southeast'] core.universes[7, 7] = univs['Assembly (3.1\%) instr'] -core.universes[7, 8] = univs['Assembly (1.6\%) CR B'] +core.universes[7, 8] = univs['Assembly (2.4\%) CR D'] core.universes[7, 9] = univs['Assembly (3.1\%) 16BA'] -core.universes[7, 10] = univs['Assembly (1.6\%) CR B'] +core.universes[7, 10] = univs['Assembly (2.4\%) CR D'] core.universes[7, 11] = univs['Assembly (3.1\%) instr'] core.universes[7, 12] = univs['baffle southwest'] core.universes[7, 13] = univs['baffle southwest corner'] core.universes[8, 5] = univs['baffle east'] core.universes[8, 6] = univs['Assembly (3.1\%) instr'] -core.universes[8, 7] = univs['Assembly (1.6\%) CR B'] +core.universes[8, 7] = univs['Assembly (2.4\%) CR D'] core.universes[8, 8] = univs['Assembly (3.1\%) 16BA'] core.universes[8, 9] = univs['Assembly (2.4\%) CR D'] core.universes[8, 10] = univs['Assembly (3.1\%) 16BA'] -core.universes[8, 11] = univs['Assembly (1.6\%) CR B'] +core.universes[8, 11] = univs['Assembly (2.4\%) CR D'] core.universes[8, 12] = univs['Assembly (3.1\%) instr'] core.universes[8, 13] = univs['baffle west'] core.universes[9, 5] = univs['baffle east'] -core.universes[9, 6] = univs['Assembly (1.6\%) CR B'] +core.universes[9, 6] = univs['Assembly (2.4\%) CR D'] core.universes[9, 7] = univs['Assembly (3.1\%) 16BA'] core.universes[9, 8] = univs['Assembly (2.4\%) CR D'] core.universes[9, 9] = univs['Assembly (1.6\%) instr'] core.universes[9, 10] = univs['Assembly (2.4\%) CR D'] core.universes[9, 11] = univs['Assembly (3.1\%) 16BA'] -core.universes[9, 12] = univs['Assembly (1.6\%) CR B'] +core.universes[9, 12] = univs['Assembly (2.4\%) CR D'] core.universes[9, 13] = univs['baffle west'] core.universes[10, 5] = univs['baffle east'] core.universes[10, 6] = univs['Assembly (3.1\%) instr'] -core.universes[10, 7] = univs['Assembly (1.6\%) CR B'] +core.universes[10, 7] = univs['Assembly (2.4\%) CR D'] core.universes[10, 8] = univs['Assembly (3.1\%) 16BA'] core.universes[10, 9] = univs['Assembly (2.4\%) CR D'] core.universes[10, 10] = univs['Assembly (3.1\%) 16BA'] -core.universes[10, 11] = univs['Assembly (1.6\%) CR B'] +core.universes[10, 11] = univs['Assembly (2.4\%) CR D'] core.universes[10, 12] = univs['Assembly (3.1\%) instr'] core.universes[10, 13] = univs['baffle west'] core.universes[11, 5] = univs['baffle northeast corner'] core.universes[11, 6] = univs['baffle northeast'] core.universes[11, 7] = univs['Assembly (3.1\%) instr'] -core.universes[11, 8] = univs['Assembly (1.6\%) CR B'] +core.universes[11, 8] = univs['Assembly (2.4\%) CR D'] core.universes[11, 9] = univs['Assembly (3.1\%) 16BA'] -core.universes[11, 10] = univs['Assembly (1.6\%) CR B'] +core.universes[11, 10] = univs['Assembly (2.4\%) CR D'] core.universes[11, 11] = univs['Assembly (3.1\%) instr'] core.universes[11, 12] = univs['baffle northwest'] core.universes[11, 13] = univs['baffle northwest corner'] @@ -83,7 +83,7 @@ core.universes[11, 13] = univs['baffle northwest corner'] core.universes[12, 6] = univs['baffle northeast corner'] core.universes[12, 7] = univs['baffle northeast'] core.universes[12, 8] = univs['Assembly (3.1\%) instr'] -core.universes[12, 9] = univs['Assembly (1.6\%) CR B'] +core.universes[12, 9] = univs['Assembly (2.4\%) CR D'] core.universes[12, 10] = univs['Assembly (3.1\%) instr'] core.universes[12, 11] = univs['baffle northwest'] core.universes[12, 12] = univs['baffle northwest corner'] diff --git a/smr/geometry.xml b/smr/geometry.xml index 7986082..d2eab5e 100644 --- a/smr/geometry.xml +++ b/smr/geometry.xml @@ -1,6 +1,6 @@ - + @@ -29,44 +29,44 @@ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + @@ -87,25 +87,25 @@ - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + @@ -159,94 +159,50 @@ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + @@ -291,26 +247,26 @@ - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + @@ -338,25 +294,25 @@ - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + @@ -371,25 +327,25 @@ - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + @@ -404,159 +360,125 @@ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + - - - - - - - - - + + + + + + + + + - - - - - - - - - + + + + + + + + + - - - - - - - - - + + + + + + + + + - - - - - - - - - - - + + + + + + + + + + + - - + + - - + + - - + + - - - + + + - - - + + + - - - + + + - - - + + + - - - + + + - - - + + + - - - + + + - - - + + + - - - - - - - - - - - - - - 1.25984 1.25984 - 17 17 - -10.70864 -10.70864 - -10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 -10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 -10097 10097 10097 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10097 10097 10097 -10097 10097 10097 10038 10097 10097 10097 10097 10097 10097 10097 10097 10097 10038 10097 10097 10097 -10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 -10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 -10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 -10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 -10097 10097 10038 10097 10097 10038 10097 10097 10010 10097 10097 10038 10097 10097 10038 10097 10097 -10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 -10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 -10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 -10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 -10097 10097 10097 10038 10097 10097 10097 10097 10097 10097 10097 10097 10097 10038 10097 10097 10097 -10097 10097 10097 10097 10097 10038 10097 10097 10038 10097 10097 10038 10097 10097 10097 10097 10097 -10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 -10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 10097 - + + + + + + + + + + + + 1.25984 1.25984 17 17 @@ -756,13 +678,13 @@ 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10238 10211 10211 10211 10241 10000 10000 10000 10000 10000 10000 10000 -10000 10000 10000 10000 10000 10000 10238 10229 10183 10111 10183 10226 10241 10000 10000 10000 10000 10000 10000 -10000 10000 10000 10000 10000 10238 10229 10183 10111 10165 10111 10183 10226 10241 10000 10000 10000 10000 10000 -10000 10000 10000 10000 10000 10214 10183 10111 10165 10145 10165 10111 10183 10217 10000 10000 10000 10000 10000 -10000 10000 10000 10000 10000 10214 10111 10165 10145 10125 10145 10165 10111 10217 10000 10000 10000 10000 10000 -10000 10000 10000 10000 10000 10214 10183 10111 10165 10145 10165 10111 10183 10217 10000 10000 10000 10000 10000 -10000 10000 10000 10000 10000 10235 10223 10183 10111 10165 10111 10183 10220 10232 10000 10000 10000 10000 10000 -10000 10000 10000 10000 10000 10000 10235 10223 10183 10111 10183 10220 10232 10000 10000 10000 10000 10000 10000 +10000 10000 10000 10000 10000 10000 10238 10229 10183 10145 10183 10226 10241 10000 10000 10000 10000 10000 10000 +10000 10000 10000 10000 10000 10238 10229 10183 10145 10165 10145 10183 10226 10241 10000 10000 10000 10000 10000 +10000 10000 10000 10000 10000 10214 10183 10145 10165 10145 10165 10145 10183 10217 10000 10000 10000 10000 10000 +10000 10000 10000 10000 10000 10214 10145 10165 10145 10125 10145 10165 10145 10217 10000 10000 10000 10000 10000 +10000 10000 10000 10000 10000 10214 10183 10145 10165 10145 10165 10145 10183 10217 10000 10000 10000 10000 10000 +10000 10000 10000 10000 10000 10235 10223 10183 10145 10165 10145 10183 10220 10232 10000 10000 10000 10000 10000 +10000 10000 10000 10000 10000 10000 10235 10223 10183 10145 10183 10220 10232 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10235 10208 10208 10208 10232 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 @@ -805,42 +727,40 @@ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/smr/pins.py b/smr/pins.py index 18b7dfa..9b9232a 100644 --- a/smr/pins.py +++ b/smr/pins.py @@ -80,10 +80,9 @@ def make_stack(name, surfaces, universes): return universe -# FIXME: Is this a good idea??? univs = {} -cell = openmc.Cell(name='water pin 1') +cell = openmc.Cell(name='water pin') cell.fill = mats['H2O'] univs['water pin'] = openmc.Universe(name='Empty water pin cell universe') univs['water pin'].add_cell(cell) diff --git a/smr/plots.py b/smr/plots.py index 8c48857..7236abb 100644 --- a/smr/plots.py +++ b/smr/plots.py @@ -25,7 +25,7 @@ plots = openmc.Plots() plot = openmc.Plot(name='radial core map ') plot.basis = 'xy' plot.color = 'mat' -plot.origin = [0., 0., (highest_extent-lowest_extent)/2.] +plot.origin = [0., 0., 200.] #(highest_extent-lowest_extent)/2.] plot.width = [25*lattice_pitch/2, 25*lattice_pitch/2.] plot.filename = 'radial_core_map' plot.col_spec = col_spec diff --git a/smr/plots.xml b/smr/plots.xml index 4b8cabf..191b2ea 100644 --- a/smr/plots.xml +++ b/smr/plots.xml @@ -2,7 +2,7 @@ - 0.0 0.0 127.722 + 0.0 0.0 200.0 268.7955 268.7955 1000 1000 255 255 255 diff --git a/smr/settings.xml b/smr/settings.xml index 48b2727..15b370e 100644 --- a/smr/settings.xml +++ b/smr/settings.xml @@ -14,6 +14,4 @@ false true - True - 1000 diff --git a/smr/surfaces.py b/smr/surfaces.py index 64db3f7..4c07023 100644 --- a/smr/surfaces.py +++ b/smr/surfaces.py @@ -141,8 +141,6 @@ surfs['lat grid box inner'] = \ surfs['lat grid box outer'] = \ openmc.get_rectangular_prism(grid_strap_side, grid_strap_side) -surfs['lowest extent'] = openmc.ZPlane( - z0=lowest_extent, name='lowest extent') surfs['bot support plate'] = openmc.ZPlane( z0=bottom_support_plate, name='bot support plate') surfs['top support plate'] = openmc.ZPlane( @@ -188,8 +186,6 @@ surfs['bot upper nozzle'] = openmc.ZPlane( z0=bottom_upper_nozzle, name='bottom upper nozzle') surfs['top upper nozzle'] = openmc.ZPlane( z0=top_upper_nozzle, name='top upper nozzle') -surfs['highest extent'] = openmc.ZPlane( - z0=highest_extent, name='highest extent') # Control rod bank surfaces for ARO configuration for bank in ['A','B','C','D','E',]: