forked from crp/ecp-benchmarks
Add reflector holes in corners, still more needed
This commit is contained in:
Paul Romano
parent
1ad4b3c60f
commit
573b51a62a
5 changed files with 178 additions and 445 deletions
|
|
@ -439,19 +439,36 @@
|
|||
<cell id="1215" material="8" name="Assembly (3.1%) no BAs instr lattice axial (6)" region="28 -29 30 -31 ~(24 -25 26 -27) 44 -45" universe="184" />
|
||||
<cell id="1216" material="4" name="Assembly (3.1%) no BAs instr lattice axial (7)" region="28 -29 30 -31 ~(24 -25 26 -27) 45 -46" universe="184" />
|
||||
<cell id="1217" material="8" name="Assembly (3.1%) no BAs instr lattice axial (last)" region="28 -29 30 -31 ~(24 -25 26 -27) 46" universe="184" />
|
||||
<cell id="1339" material="4" name="heavy reflector SS" universe="207" />
|
||||
<cell fill="244" id="1384" name="Main core" region="-71 81 -80" universe="0" />
|
||||
<cell id="1385" material="4" name="core barrel" region="71 -72 81 -80" universe="0" />
|
||||
<cell id="1386" material="4" name="neutron shield panel NW" region="72 -73 74 -75 81 -80" universe="0" />
|
||||
<cell id="1387" material="8" name="neutron shield panel N" region="72 -73 75 -77 81 -80" universe="0" />
|
||||
<cell id="1388" material="4" name="neutron shield panel SE" region="72 -73 -74 75 81 -80" universe="0" />
|
||||
<cell id="1389" material="8" name="neutron shield panel E" region="72 -73 74 76 81 -80" universe="0" />
|
||||
<cell id="1390" material="4" name="neutron shield panel NE" region="72 -73 76 -77 81 -80" universe="0" />
|
||||
<cell id="1391" material="8" name="neutron shield panel S" region="72 -73 -75 77 81 -80" universe="0" />
|
||||
<cell id="1392" material="4" name="neutron shield panel SW" region="72 -73 -76 77 81 -80" universe="0" />
|
||||
<cell id="1393" material="8" name="neutron shield panel W" region="72 -73 -74 -76 81 -80" universe="0" />
|
||||
<cell id="1394" material="8" name="downcomer" region="73 -78 81 -80" universe="0" />
|
||||
<cell id="1395" material="5" name="reactor pressure vessel" region="78 -79 81 -80" universe="0" />
|
||||
<cell id="1339" material="8" region="-82" universe="207" />
|
||||
<cell id="1340" material="8" region="-83" universe="207" />
|
||||
<cell id="1341" material="8" region="-84" universe="207" />
|
||||
<cell id="1342" material="8" region="-85" universe="207" />
|
||||
<cell id="1343" material="8" region="-86" universe="207" />
|
||||
<cell id="1344" material="8" region="-87" universe="207" />
|
||||
<cell id="1345" material="8" region="-88" universe="207" />
|
||||
<cell id="1346" material="8" region="-89" universe="207" />
|
||||
<cell id="1347" material="8" region="-90" universe="207" />
|
||||
<cell id="1348" material="8" region="-91" universe="207" />
|
||||
<cell id="1349" material="8" region="-92" universe="207" />
|
||||
<cell id="1350" material="8" region="-93" universe="207" />
|
||||
<cell id="1351" material="8" region="-94" universe="207" />
|
||||
<cell id="1352" material="4" name="heavy reflector SS" region="82 83 84 85 86 87 88 89 90 91 92 93 94" universe="207" />
|
||||
<cell fill="207" id="1353" name="heavy reflector NE" rotation="0 180 0" universe="208" />
|
||||
<cell fill="207" id="1354" name="heavy reflector SW" rotation="180 0 0" universe="209" />
|
||||
<cell fill="207" id="1355" name="heavy reflector SE" rotation="0 0 180" universe="210" />
|
||||
<cell id="1356" material="4" name="heavy reflector" universe="211" />
|
||||
<cell fill="212" id="1357" name="Main core" region="-71 81 -80" universe="0" />
|
||||
<cell id="1358" material="4" name="core barrel" region="71 -72 81 -80" universe="0" />
|
||||
<cell id="1359" material="4" name="neutron shield panel NW" region="72 -73 74 -75 81 -80" universe="0" />
|
||||
<cell id="1360" material="8" name="neutron shield panel N" region="72 -73 75 -77 81 -80" universe="0" />
|
||||
<cell id="1361" material="4" name="neutron shield panel SE" region="72 -73 -74 75 81 -80" universe="0" />
|
||||
<cell id="1362" material="8" name="neutron shield panel E" region="72 -73 74 76 81 -80" universe="0" />
|
||||
<cell id="1363" material="4" name="neutron shield panel NE" region="72 -73 76 -77 81 -80" universe="0" />
|
||||
<cell id="1364" material="8" name="neutron shield panel S" region="72 -73 -75 77 81 -80" universe="0" />
|
||||
<cell id="1365" material="4" name="neutron shield panel SW" region="72 -73 -76 77 81 -80" universe="0" />
|
||||
<cell id="1366" material="8" name="neutron shield panel W" region="72 -73 -74 -76 81 -80" universe="0" />
|
||||
<cell id="1367" material="8" name="downcomer" region="73 -78 81 -80" universe="0" />
|
||||
<cell id="1368" material="5" name="reactor pressure vessel" region="78 -79 81 -80" universe="0" />
|
||||
<lattice id="125" name="Assembly (1.6%) no BAs instr">
|
||||
<pitch>1.25984 1.25984</pitch>
|
||||
<dimension>17 17</dimension>
|
||||
|
|
@ -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 </universes>
|
||||
</lattice>
|
||||
<lattice id="244" name="Main core">
|
||||
<lattice id="212" name="Main core">
|
||||
<pitch>21.503639999999997 21.503639999999997</pitch>
|
||||
<dimension>9 9</dimension>
|
||||
<lower_left>-96.76637999999998 -96.76637999999998</lower_left>
|
||||
<universes>
|
||||
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 </universes>
|
||||
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 </universes>
|
||||
</lattice>
|
||||
<surface coeffs="0.0 0.0 0.405765" id="1" name="Pellet OR" type="z-cylinder" />
|
||||
<surface coeffs="0.0 0.0 0.06459" id="2" name="FR Plenum Spring OR" type="z-cylinder" />
|
||||
|
|
@ -624,4 +641,17 @@
|
|||
<surface boundary="vacuum" coeffs="0.0 0.0 135.0" id="79" name="RPV OR" type="z-cylinder" />
|
||||
<surface boundary="vacuum" coeffs="272.564" id="80" name="upper problem boundary" type="z-plane" />
|
||||
<surface boundary="vacuum" coeffs="0.0" id="81" name="lower problem boundary" type="z-plane" />
|
||||
<surface coeffs="-9.816879130434781 -6.155027391304348 0.506426304347826" id="82" type="z-cylinder" />
|
||||
<surface coeffs="-4.518880869565217 -6.155027391304348 0.506426304347826" id="83" type="z-cylinder" />
|
||||
<surface coeffs="0.9349408695652173 -6.155027391304348 0.506426304347826" id="84" type="z-cylinder" />
|
||||
<surface coeffs="6.155027391304348 -6.155027391304348 1.1686760869565216" id="85" type="z-cylinder" />
|
||||
<surface coeffs="6.155027391304348 -0.9349408695652173 0.506426304347826" id="86" type="z-cylinder" />
|
||||
<surface coeffs="6.155027391304348 4.518880869565216 0.506426304347826" id="87" type="z-cylinder" />
|
||||
<surface coeffs="6.155027391304348 9.816879130434781 0.506426304347826" id="88" type="z-cylinder" />
|
||||
<surface coeffs="2.5710873913043475 -2.5710873913043475 0.506426304347826" id="89" type="z-cylinder" />
|
||||
<surface coeffs="-2.0257052173913053 -1.246587826086957 0.506426304347826" id="90" type="z-cylinder" />
|
||||
<surface coeffs="1.2465878260869552 2.0257052173913035 0.506426304347826" id="91" type="z-cylinder" />
|
||||
<surface coeffs="-6.544586086956521 0.0 0.506426304347826" id="92" type="z-cylinder" />
|
||||
<surface coeffs="0.0 6.544586086956521 0.506426304347826" id="93" type="z-cylinder" />
|
||||
<surface coeffs="-9.816879130434781 8.49237956521739 0.506426304347826" id="94" type="z-cylinder" />
|
||||
</geometry>
|
||||
|
|
|
|||
|
|
@ -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)
|
||||
|
|
@ -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']
|
||||
|
|
|
|||
116
smr/smr/reflector.py
Normal file
116
smr/smr/reflector.py
Normal file
|
|
@ -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])
|
||||
|
|
@ -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')
|
||||
Loading…
Add table
Add a link
Reference in a new issue