aparler/ecp-benchmarks
1
0
Fork 0
forked from crp/ecp-benchmarks
Code Pull requests Activity

Add ten rings in fuel

Browse source
This commit is contained in:
Paul Romano 2017-10-20 07:54:57 -04:00
parent a0a7e58197
commit 6762ab0237
2 changed files with 39 additions and 65 deletions
Download patch file Download diff file Expand all files Collapse all files

92
smr/smr/pins.py
View file

@ -3,7 +3,7 @@
import openmc
from .materials import mats
from .surfaces import surfs
from .surfaces import surfs, n_rings
def make_pin(name, surfaces, materials, grid=None):
@ -669,36 +669,26 @@ univs['pin plenum grid (intermediate)'] = make_pin(
#### 1.6% ENRICHED FUEL PIN CELL
rings = [surfs['fuel ring {}'.format(i)] for i in range(1, n_rings)]
fuel_surfaces = rings + [surfs['pellet OR'], surfs['clad IR'],
surfs['clad OR']]
fuel_mats = [mats['UO2 1.6']]*10 + [mats['He'], mats['Zr'], mats['H2O']]
univs['Fuel (1.6%)'] = make_pin(
'Fuel (1.6%)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
materials=[mats['UO2 1.6'],
mats['He'],
mats['Zr'],
mats['H2O']])
surfaces=fuel_surfaces,
materials=fuel_mats)
univs['Fuel (1.6%) grid (bottom)'] = make_pin(
'Fuel (1.6%) grid (bottom)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
materials=[mats['UO2 1.6'],
mats['He'],
mats['Zr'],
mats['H2O']],
surfaces=fuel_surfaces,
materials=fuel_mats,
grid='bottom')
univs['Fuel (1.6%) grid (intermediate)'] = make_pin(
'Fuel (1.6%) grid (intermediate)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
materials=[mats['UO2 1.6'],
mats['He'],
mats['Zr'],
mats['H2O']],
surfaces=fuel_surfaces,
materials=fuel_mats,
grid='intermediate')
# Stack all axial pieces of 1.6% enriched fuel pin cell
@ -731,36 +721,23 @@ univs['Fuel (1.6%) stack'] = make_stack(
#### 2.4% ENRICHED FUEL PIN CELL
fuel_mats = [mats['UO2 2.4']]*10 + [mats['He'], mats['Zr'], mats['H2O']]
univs['Fuel (2.4%)'] = make_pin(
'Fuel (2.4%)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
materials=[mats['UO2 2.4'],
mats['He'],
mats['Zr'],
mats['H2O']])
surfaces=fuel_surfaces,
materials=fuel_mats)
univs['Fuel (2.4%) grid (bottom)'] = make_pin(
'Fuel (2.4%) grid (bottom)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
materials=[mats['UO2 2.4'],
mats['He'],
mats['Zr'],
mats['H2O']],
surfaces=fuel_surfaces,
materials=fuel_mats,
grid='bottom')
univs['Fuel (2.4%) grid (intermediate)'] = make_pin(
'Fuel (2.4%) grid (intermediate)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
materials=[mats['UO2 2.4'],
mats['He'],
mats['Zr'],
mats['H2O']],
surfaces=fuel_surfaces,
materials=fuel_mats,
grid='intermediate')
# Stack all axial pieces of 2.4% enriched fuel pin cell
@ -793,36 +770,23 @@ univs['Fuel (2.4%) stack'] = make_stack(
#### 3.1% ENRICHED FUEL PIN CELL
fuel_mats = [mats['UO2 3.1']]*10 + [mats['He'], mats['Zr'], mats['H2O']]
univs['Fuel (3.1%)'] = make_pin(
'Fuel (3.1%)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
materials=[mats['UO2 3.1'],
mats['He'],
mats['Zr'],
mats['H2O']])
surfaces=fuel_surfaces,
materials=fuel_mats)
univs['Fuel (3.1%) grid (bottom)'] = make_pin(
'Fuel (3.1%) grid (bottom)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
materials=[mats['UO2 3.1'],
mats['He'],
mats['Zr'],
mats['H2O']],
surfaces=fuel_surfaces,
materials=fuel_mats,
grid='bottom')
univs['Fuel (3.1%) grid (intermediate)'] = make_pin(
'Fuel (3.1%) grid (intermediate)',
surfaces=[surfs['pellet OR'],
surfs['clad IR'],
surfs['clad OR']],
materials=[mats['UO2 3.1'],
mats['He'],
mats['Zr'],
mats['H2O']],
surfaces=fuel_surfaces,
materials=fuel_mats,
grid='intermediate')
# Stack all axial pieces of 3.1% enriched fuel pin cell

12
smr/smr/surfaces.py
View file

@ -14,7 +14,7 @@ NuScale DC application, chapter 4: https://www.nrc.gov/docs/ML1701/ML17013A274.p
"""
import copy
from math import tan, pi
from math import tan, pi, sqrt
import numpy as np
import openmc
@ -126,6 +126,16 @@ neutron_shield_NEtop_SWbot = tan(-pi/6)
surfs = {}
# Create radii for rings in fuel
n_rings = 10
ring_radii = []
for i in range(1, n_rings):
R = sqrt(i*pellet_OR**2/n_rings)
ring_radii.append(R)
name = 'fuel ring {}'.format(i)
surfs[name] = openmc.ZCylinder(R=R, name=name)
surfs['pellet OR'] = openmc.ZCylinder(
R=pellet_OR, name='Pellet OR')
surfs['plenum spring OR'] = openmc.ZCylinder(