crp/ecp-benchmarks
1
0
Fork 1
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Fix number of pellets in SMR

Browse source
This commit is contained in:
Paul Romano 2017-11-10 11:49:36 -06:00
parent 186c525c8a
commit 7e735bec85
4 changed files with 4435 additions and 4403 deletions
Download patch file Download diff file Expand all files Collapse all files

8811
smr/fresh/geometry.xml
View file

File diff suppressed because it is too large Load diff

6
smr/fresh/settings.xml
View file

@ -1,9 +1,9 @@
<?xml version='1.0' encoding='utf-8'?>
<settings>
<run_mode>eigenvalue</run_mode>
<particles>1000</particles>
<batches>20</batches>
<inactive>10</inactive>
<particles>10000</particles>
<batches>200</batches>
<inactive>100</inactive>
<source strength="1.0">
<space type="fission">
<parameters>-75.26274 -75.26274 36.007 75.26274 75.26274 236.0066</parameters>

18
smr/fresh/tallies.xml
View file

@ -1,31 +1,31 @@
<?xml version='1.0' encoding='utf-8'?>
<tallies>
<filter id="1" type="distribcell">
<bins>2634</bins>
<bins>2644</bins>
</filter>
<filter id="2" type="distribcell">
<bins>2638</bins>
<bins>2648</bins>
</filter>
<filter id="3" type="distribcell">
<bins>2643</bins>
<bins>2653</bins>
</filter>
<filter id="4" type="distribcell">
<bins>4619</bins>
<bins>4639</bins>
</filter>
<filter id="5" type="distribcell">
<bins>4623</bins>
<bins>4643</bins>
</filter>
<filter id="6" type="distribcell">
<bins>4628</bins>
<bins>4648</bins>
</filter>
<filter id="7" type="distribcell">
<bins>6604</bins>
<bins>6634</bins>
</filter>
<filter id="8" type="distribcell">
<bins>6608</bins>
<bins>6638</bins>
</filter>
<filter id="9" type="distribcell">
<bins>6613</bins>
<bins>6643</bins>
</filter>
<tally id="1" name="depletion tally">
<filters>1</filters>

3
smr/smr/pins.py
View file

@ -669,7 +669,8 @@ univs['pin plenum grid (intermediate)'] = make_pin(
# Determine z position between each fuel pellet, omitting the surfaces
# corresponding to the very bottom and top of the active fuel length
axial_splits = np.linspace(bottom_fuel_rod, top_active_core, 196)[1:-1]
n_pellets = 196
axial_splits = np.linspace(bottom_fuel_rod, top_active_core, n_pellets + 1)[1:-1]
axial_surfs = [openmc.ZPlane(z0=z) for z in axial_splits]
# Get z-cylinder surfaces for each ring