diff --git a/smr/fresh/geometry.xml b/smr/fresh/geometry.xml
index 40f5504..68df767 100644
--- a/smr/fresh/geometry.xml
+++ b/smr/fresh/geometry.xml
@@ -439,62 +439,19 @@
|
|
|
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
+ |
1.25984 1.25984
17 17
@@ -587,126 +544,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
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-207 207
-1 1
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-1 1
-210 210
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-1 213
-1 213
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-216 1
-216 1
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-219 207
-216 1
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-207 222
-1 213
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-216 1
-225 210
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-1 213
-210 228
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-231 1
-1 1
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-1 234
-1 1
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-1 1
-1 237
-
-
- 10.751819999999999 10.751819999999999
- 2 2
- -10.751819999999999 -10.751819999999999
-
-1 1
-240 1
-
-
+
21.503639999999997 21.503639999999997
- 19 19
- -204.28457999999998 -204.28457999999998
+ 9 9
+ -96.76637999999998 -96.76637999999998
-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
-1 1 1 1 1 1 1 239 212 212 212 242 1 1 1 1 1 1 1
-1 1 1 1 1 1 239 230 184 146 184 227 242 1 1 1 1 1 1
-1 1 1 1 1 239 230 184 146 166 146 184 227 242 1 1 1 1 1
-1 1 1 1 1 215 184 146 166 146 166 146 184 218 1 1 1 1 1
-1 1 1 1 1 215 146 166 146 126 146 166 146 218 1 1 1 1 1
-1 1 1 1 1 215 184 146 166 146 166 146 184 218 1 1 1 1 1
-1 1 1 1 1 236 224 184 146 166 146 184 221 233 1 1 1 1 1
-1 1 1 1 1 1 236 224 184 146 184 221 233 1 1 1 1 1 1
-1 1 1 1 1 1 1 236 209 209 209 233 1 1 1 1 1 1 1
-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
+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
@@ -773,8 +624,4 @@
-
-
-
-
diff --git a/smr/smr/baffle.py b/smr/smr/baffle.py
index 793dc2f..9cc849e 100644
--- a/smr/smr/baffle.py
+++ b/smr/smr/baffle.py
@@ -7,6 +7,18 @@ 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')
diff --git a/smr/smr/core.py b/smr/smr/core.py
index bf5c336..5984985 100644
--- a/smr/smr/core.py
+++ b/smr/smr/core.py
@@ -12,90 +12,90 @@ from .baffle import univs
#### CONSTRUCT MAIN CORE LATTICE
core = openmc.RectLattice(name='Main core')
-core.lower_left = [-19.*lattice_pitch/2., -19.*lattice_pitch/2.]
+core.lower_left = [-9*lattice_pitch/2, -9*lattice_pitch/2]
core.pitch = [lattice_pitch, lattice_pitch]
-universes = np.empty((19,19), dtype=openmc.Universe)
-universes[:,:] = univs['water pin']
+universes = np.tile(univs['heavy reflector'], (9,9))
+
+universes[0, 2] = univs['heavy reflector']
+universes[0, 3] = univs['heavy reflector']
+universes[0, 4] = univs['heavy reflector']
+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, 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, 7] = univs['heavy reflector']
+
+universes[2, 0] = univs['heavy reflector']
+universes[2, 1] = univs['heavy reflector']
+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, 8] = univs['heavy reflector']
+
+universes[3, 0] = univs['heavy reflector']
+universes[3, 1] = univs['Assembly (3.1%) instr']
+universes[3, 2] = univs['Assembly (2.4%) CR D']
+universes[3, 3] = univs['Assembly (3.1%) 16BA']
+universes[3, 4] = univs['Assembly (2.4%) CR D']
+universes[3, 5] = univs['Assembly (3.1%) 16BA']
+universes[3, 6] = univs['Assembly (2.4%) CR D']
+universes[3, 7] = univs['Assembly (3.1%) instr']
+universes[3, 8] = univs['heavy reflector']
+
+universes[4, 0] = univs['heavy reflector']
+universes[4, 1] = univs['Assembly (2.4%) CR D']
+universes[4, 2] = univs['Assembly (3.1%) 16BA']
+universes[4, 3] = univs['Assembly (2.4%) CR D']
+universes[4, 4] = univs['Assembly (1.6%) instr']
+universes[4, 5] = univs['Assembly (2.4%) CR D']
+universes[4, 6] = univs['Assembly (3.1%) 16BA']
+universes[4, 7] = univs['Assembly (2.4%) CR D']
+universes[4, 8] = univs['heavy reflector']
+
+universes[5, 0] = univs['heavy reflector']
+universes[5, 1] = univs['Assembly (3.1%) instr']
+universes[5, 2] = univs['Assembly (2.4%) CR D']
+universes[5, 3] = univs['Assembly (3.1%) 16BA']
+universes[5, 4] = univs['Assembly (2.4%) CR D']
+universes[5, 5] = univs['Assembly (3.1%) 16BA']
+universes[5, 6] = univs['Assembly (2.4%) CR D']
+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, 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, 8] = univs['heavy reflector']
+
+universes[7, 1] = univs['heavy reflector']
+universes[7, 2] = univs['heavy reflector']
+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, 7] = univs['heavy reflector']
+
+universes[8, 2] = univs['heavy reflector']
+universes[8, 3] = univs['heavy reflector']
+universes[8, 4] = univs['heavy reflector']
+universes[8, 5] = univs['heavy reflector']
+universes[8, 6] = univs['heavy reflector']
+
core.universes = universes
-core.universes[5, 7] = univs['baffle northwest corner']
-core.universes[5, 8] = univs['baffle north']
-core.universes[5, 9] = univs['baffle north']
-core.universes[5, 10] = univs['baffle north']
-core.universes[5, 11] = univs['baffle northeast corner']
-
-core.universes[6, 6] = univs['baffle northwest corner']
-core.universes[6, 7] = univs['baffle northwest']
-core.universes[6, 8] = univs['Assembly (3.1%) instr']
-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 northeast']
-core.universes[6, 12] = univs['baffle northeast corner']
-
-core.universes[7, 5] = univs['baffle northwest corner']
-core.universes[7, 6] = univs['baffle northwest']
-core.universes[7, 7] = univs['Assembly (3.1%) instr']
-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 (2.4%) CR D']
-core.universes[7, 11] = univs['Assembly (3.1%) instr']
-core.universes[7, 12] = univs['baffle northeast']
-core.universes[7, 13] = univs['baffle northeast corner']
-
-core.universes[8, 5] = univs['baffle west']
-core.universes[8, 6] = univs['Assembly (3.1%) instr']
-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 (2.4%) CR D']
-core.universes[8, 12] = univs['Assembly (3.1%) instr']
-core.universes[8, 13] = univs['baffle east']
-
-core.universes[9, 5] = univs['baffle west']
-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 (2.4%) CR D']
-core.universes[9, 13] = univs['baffle east']
-
-core.universes[10, 5] = univs['baffle west']
-core.universes[10, 6] = univs['Assembly (3.1%) instr']
-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 (2.4%) CR D']
-core.universes[10, 12] = univs['Assembly (3.1%) instr']
-core.universes[10, 13] = univs['baffle east']
-
-core.universes[11, 5] = univs['baffle southwest corner']
-core.universes[11, 6] = univs['baffle southwest']
-core.universes[11, 7] = univs['Assembly (3.1%) instr']
-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 (2.4%) CR D']
-core.universes[11, 11] = univs['Assembly (3.1%) instr']
-core.universes[11, 12] = univs['baffle southeast']
-core.universes[11, 13] = univs['baffle southeast corner']
-
-core.universes[12, 6] = univs['baffle southwest corner']
-core.universes[12, 7] = univs['baffle southwest']
-core.universes[12, 8] = univs['Assembly (3.1%) instr']
-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 southeast']
-core.universes[12, 12] = univs['baffle southeast corner']
-
-core.universes[13, 7] = univs['baffle southwest corner']
-core.universes[13, 8] = univs['baffle south']
-core.universes[13, 9] = univs['baffle south']
-core.universes[13, 10] = univs['baffle south']
-core.universes[13, 11] = univs['baffle southeast corner']
-
#### CONSTRUCT ROOT UNIVERSE AND CELLS
@@ -204,4 +204,4 @@ root_univ.add_cell(cell)
#### CONSTRUCT GEOMETRY
-geometry = openmc.Geometry(root_univ)
\ No newline at end of file
+geometry = openmc.Geometry(root_univ)