Changed AP1000 XML input files to be created through the Python API
This commit is contained in:
parent
06300aefe3
commit
21854a7ab6
18 changed files with 45 additions and 859 deletions
|
|
@ -1,70 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<geometry>
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
<cell id="1" fill="5" region="1 -2 3 -4 10 -11" />
|
||||
<cell id="10" universe="1" material="1" region="-5" />
|
||||
<cell id="20" universe="1" material="2" region="5 -6 " />
|
||||
<cell id="30" universe="1" material="3" region="6 -7" />
|
||||
<cell id="40" universe="1" material="4" region="7" />
|
||||
|
||||
<cell id="50" universe="2" material="4" region="-8" />
|
||||
<cell id="60" universe="2" material="5" region="8 -9" />
|
||||
<cell id="70" universe="2" material="4" region="9" />
|
||||
|
||||
<lattice id="5" >
|
||||
<dimension>17 17</dimension>
|
||||
<lower_left>-10.659 -10.659</lower_left>
|
||||
<pitch>1.26 1.26</pitch>
|
||||
<universes>
|
||||
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 2 1 1 2 1 1 2 1 1 1 1 1
|
||||
1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
|
||||
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
|
||||
1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 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 2 1 1 2 1 1 2 1 1 2 1 1 2 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 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
|
||||
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
|
||||
1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
|
||||
1 1 1 1 1 2 1 1 2 1 1 2 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
|
||||
</universes>
|
||||
</lattice>
|
||||
|
||||
<surface id="1" type="x-plane" coeffs="-10.659" boundary="reflective" />
|
||||
<surface id="2" type="x-plane" coeffs="10.659" boundary="reflective" />
|
||||
<surface id="3" type="y-plane" coeffs="-10.659" boundary="reflective" />
|
||||
<surface id="4" type="y-plane" coeffs="10.659" boundary="reflective" />
|
||||
|
||||
<surface id="5" type="z-cylinder" coeffs="0. 0. 0.4095" />
|
||||
<surface id="6" type="z-cylinder" coeffs="0. 0. 0.4177" />
|
||||
<surface id="7" type="z-cylinder" coeffs="0. 0. 0.4749" />
|
||||
<surface id="8" type="z-cylinder" coeffs="0. 0. 0.5715" />
|
||||
<surface id="9" type="z-cylinder" coeffs="0. 0. 0.612" />
|
||||
|
||||
<surface id="10" type="z-plane" coeffs="-182.9" boundary="reflective" />
|
||||
<surface id="11" type="z-plane" coeffs="182.9" boundary="reflective" />
|
||||
|
||||
<surface id="12" type="x-plane" coeffs="-0.627" />
|
||||
<surface id="13" type="x-plane" coeffs="0.627" />
|
||||
<surface id="14" type="y-plane" coeffs="-0.627" />
|
||||
<surface id="15" type="y-plane" coeffs="0.627" />
|
||||
|
||||
</geometry>
|
||||
|
||||
|
||||
|
||||
|
|
@ -1,84 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<materials>
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
<cross_sections>/opt/xdata/endfb-vii.1-hdf5/cross_sections.xml</cross_sections>
|
||||
|
||||
<default_xs>71c</default_xs>
|
||||
|
||||
<!-- Fuel composed of U-238 O-16 and U-235 -->
|
||||
<material id="1">
|
||||
<density value="10.5" units="g/cc" />
|
||||
<nuclide name="U235" xs="71c" wo="0.0185105302348" />
|
||||
<nuclide name="U238" xs="71c" wo="0.862943290471" />
|
||||
<nuclide name="O16" xs="71c" wo="0.118546179294" />
|
||||
</material>
|
||||
|
||||
<!-- Helium for gap -->
|
||||
<material id="2">
|
||||
<density value="1.00000e-4" units="atom/b-cm" />
|
||||
<nuclide name="He4" xs="71c" ao="1.0" />
|
||||
</material>
|
||||
|
||||
<!-- Zircaloy 4 -->
|
||||
<material id="3">
|
||||
<density value="6.55" units="g/cc" />
|
||||
<nuclide name="O16" xs="71c" wo="0.00125" />
|
||||
<nuclide name="Nb93" xs="71c" wo="0.01" />
|
||||
<nuclide name="Sn112" xs="71c" wo="0.000097" />
|
||||
<nuclide name="Sn114" xs="71c" wo="0.000066" />
|
||||
<nuclide name="Sn115" xs="71c" wo="0.000034" />
|
||||
<nuclide name="Sn116" xs="71c" wo="0.001454" />
|
||||
<nuclide name="Sn117" xs="71c" wo="0.000768" />
|
||||
<nuclide name="Sn118" xs="71c" wo="0.002422 " />
|
||||
<nuclide name="Sn119" xs="71c" wo="0.000859" />
|
||||
<nuclide name="Sn120" xs="71c" wo="0.003258" />
|
||||
<nuclide name="Sn122" xs="71c" wo="0.000463" />
|
||||
<nuclide name="Sn124" xs="71c" wo="0.000579" />
|
||||
<nuclide name="Fe54" xs="71c" wo="0.00005845" />
|
||||
<nuclide name="Fe56" xs="71c" wo="0.00091754" />
|
||||
<nuclide name="Fe57" xs="71c" wo="0.00002119" />
|
||||
<nuclide name="Fe58" xs="71c" wo="0.00000282" />
|
||||
<nuclide name="Zr90" xs="71c" wo="0.503052375" />
|
||||
<nuclide name="Zr91" xs="71c" wo="0.10970355" />
|
||||
<nuclide name="Zr92" xs="71c" wo="0.167684125" />
|
||||
<nuclide name="Zr94" xs="71c" wo="0.16993295" />
|
||||
<nuclide name="Zr96" xs="71c" wo="0.027377" />
|
||||
</material>
|
||||
|
||||
<!-- Water -->
|
||||
<material id="4">
|
||||
<density value="0.73" units="g/cc" />
|
||||
<nuclide name="H1" xs="71c" ao="2.0" />
|
||||
<nuclide name="O16" xs="71c" ao="1.0" />
|
||||
<sab name="c_H_in_H2O" xs="04t" />
|
||||
</material>
|
||||
|
||||
<!-- Guide Tube -->
|
||||
<material id="5">
|
||||
<density value="7.92" units="g/cc" />
|
||||
<nuclide name="Ni58" xs="71c" wo="0.06467315" />
|
||||
<nuclide name="Ni60" xs="71c" wo="0.02492135" />
|
||||
<nuclide name="Ni61" xs="71c" wo="0.001082905" />
|
||||
<nuclide name="Ni62" xs="71c" wo="0.00345287" />
|
||||
<nuclide name="Ni64" xs="71c" wo="0.000879225" />
|
||||
<nuclide name="Fe54" xs="71c" wo="0.04062275" />
|
||||
<nuclide name="Fe56" xs="71c" wo="0.63762775" />
|
||||
<nuclide name="Fe57" xs="71c" wo="0.01472705" />
|
||||
<nuclide name="Fe58" xs="71c" wo="0.0019599" />
|
||||
<nuclide name="Mn55" xs="71c" wo="0.020" />
|
||||
<nuclide name="Cr50" xs="71c" wo="0.0082555" />
|
||||
<nuclide name="Cr52" xs="71c" wo="0.1591991" />
|
||||
<nuclide name="Cr53" xs="71c" wo="0.0180519" />
|
||||
<nuclide name="Cr54" xs="71c" wo="0.0044935" />
|
||||
</material>
|
||||
|
||||
</materials>
|
||||
|
|
@ -1,24 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<plots>
|
||||
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
|
||||
<plot id="1" color="mat" >
|
||||
<origin>0. 0. 0.</origin>
|
||||
<width>25.50 25.50 </width>
|
||||
<pixels>1000 1000 </pixels>
|
||||
<col_spec id="1" rgb="123 123 231" />
|
||||
|
||||
|
||||
</plot>
|
||||
|
||||
</plots>
|
||||
|
|
@ -1,29 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<settings>
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
<eigenvalue>
|
||||
<batches>600</batches>
|
||||
<inactive>50</inactive>
|
||||
<particles>10000</particles>
|
||||
</eigenvalue>
|
||||
|
||||
<source>
|
||||
<space>
|
||||
<type>box</type>
|
||||
<parameters>
|
||||
-10.659 -10.659 -182.9
|
||||
10.659 10.659 182.9
|
||||
</parameters>
|
||||
</space>
|
||||
</source>
|
||||
|
||||
</settings>
|
||||
|
|
@ -1,77 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<tallies>
|
||||
|
||||
<filter id="1" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="2" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="3" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="4" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="5" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="6" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="7" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="8" type="delayedgroup">
|
||||
<bins>1 2 3 4 5 6</bins>
|
||||
</filter>
|
||||
|
||||
<tally id="1">
|
||||
<filters>1</filters>
|
||||
<scores>flux</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="2">
|
||||
<filters>2</filters>
|
||||
<scores>nu-fission</scores>
|
||||
</tally>
|
||||
|
||||
|
||||
<tally id="3">
|
||||
<filters>3</filters>
|
||||
<scores>delayed-nu-fission</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="4">
|
||||
<filters>4</filters>
|
||||
<scores>fission</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="5">
|
||||
<filters>5</filters>
|
||||
<scores>absorption</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="6">
|
||||
<filters>6</filters>
|
||||
<scores>scatter</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="7">
|
||||
<filters>7</filters>
|
||||
<scores>total</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="8">
|
||||
<filters>8</filters>
|
||||
<scores>delayed-nu-fission</scores>
|
||||
</tally>
|
||||
|
||||
</tallies>
|
||||
|
|
@ -1,19 +1,22 @@
|
|||
|
||||
#!/usr/bin/env python3
|
||||
|
||||
"""
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
Conversion: Adam Parler
|
||||
Date: 12/03/2023
|
||||
====================================================
|
||||
"""
|
||||
|
||||
|
||||
|
||||
|
||||
import os
|
||||
import openmc
|
||||
|
||||
try:
|
||||
os.mkdir(os.path.dirname(os.path.realpath(__file__)) + '/problem')
|
||||
except OSError: pass
|
||||
|
||||
################## MATERIALS ##################
|
||||
|
||||
fuel_21 = openmc.Material()
|
||||
|
|
@ -179,7 +182,7 @@ plot_file.export_to_xml(path="problem/plots.xml")
|
|||
################## SETTINGS ##################
|
||||
|
||||
settings_file = openmc.Settings()
|
||||
settings_file.source = openmc.Source(space=openmc.stats.Box((-10.659, -10.659, -182.9),
|
||||
settings_file.source = openmc.IndependentSource(space=openmc.stats.Box((-10.659, -10.659, -182.9),
|
||||
(10.659, 10.659, 182.9)))
|
||||
settings_file.batches = 600
|
||||
settings_file.inactive = 50
|
||||
|
|
|
|||
|
|
@ -1,70 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<geometry>
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
<cell id="1" fill="5" region="1 -2 3 -4 10 -11" />
|
||||
<cell id="10" universe="1" material="1" region="-5" />
|
||||
<cell id="20" universe="1" material="2" region="5 -6 " />
|
||||
<cell id="30" universe="1" material="3" region="6 -7" />
|
||||
<cell id="40" universe="1" material="4" region="7" />
|
||||
|
||||
<cell id="50" universe="2" material="4" region="-8" />
|
||||
<cell id="60" universe="2" material="5" region="8 -9" />
|
||||
<cell id="70" universe="2" material="4" region="9" />
|
||||
|
||||
<lattice id="5" >
|
||||
<dimension>17 17</dimension>
|
||||
<lower_left>-10.659 -10.659</lower_left>
|
||||
<pitch>1.26 1.26</pitch>
|
||||
<universes>
|
||||
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 2 1 1 2 1 1 2 1 1 1 1 1
|
||||
1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
|
||||
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
|
||||
1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 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 2 1 1 2 1 1 2 1 1 2 1 1 2 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 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
|
||||
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
|
||||
1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
|
||||
1 1 1 1 1 2 1 1 2 1 1 2 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
|
||||
</universes>
|
||||
</lattice>
|
||||
|
||||
<surface id="1" type="x-plane" coeffs="-10.659" boundary="reflective" />
|
||||
<surface id="2" type="x-plane" coeffs="10.659" boundary="reflective" />
|
||||
<surface id="3" type="y-plane" coeffs="-10.659" boundary="reflective" />
|
||||
<surface id="4" type="y-plane" coeffs="10.659" boundary="reflective" />
|
||||
|
||||
<surface id="5" type="z-cylinder" coeffs="0. 0. 0.4095" />
|
||||
<surface id="6" type="z-cylinder" coeffs="0. 0. 0.4177" />
|
||||
<surface id="7" type="z-cylinder" coeffs="0. 0. 0.4749" />
|
||||
<surface id="8" type="z-cylinder" coeffs="0. 0. 0.5715" />
|
||||
<surface id="9" type="z-cylinder" coeffs="0. 0. 0.612" />
|
||||
|
||||
<surface id="10" type="z-plane" coeffs="-182.9" boundary="reflective" />
|
||||
<surface id="11" type="z-plane" coeffs="182.9" boundary="reflective" />
|
||||
|
||||
<surface id="12" type="x-plane" coeffs="-0.627" />
|
||||
<surface id="13" type="x-plane" coeffs="0.627" />
|
||||
<surface id="14" type="y-plane" coeffs="-0.627" />
|
||||
<surface id="15" type="y-plane" coeffs="0.627" />
|
||||
|
||||
</geometry>
|
||||
|
||||
|
||||
|
||||
|
|
@ -1,84 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<materials>
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
<cross_sections>/opt/xdata/endfb-vii.1-hdf5/cross_sections.xml</cross_sections>
|
||||
|
||||
<default_xs>71c</default_xs>
|
||||
|
||||
<!-- Fuel composed of U-238 O-16 and U-235 -->
|
||||
<material id="1">
|
||||
<density value="10.5" units="g/cc" />
|
||||
<nuclide name="U235" xs="71c" wo="0.0229176" />
|
||||
<nuclide name="U238" xs="71c" wo="0.858529" />
|
||||
<nuclide name="O16" xs="71c" wo="0.11855276" />
|
||||
</material>
|
||||
|
||||
<!-- Helium for gap -->
|
||||
<material id="2">
|
||||
<density value="1.00000e-4" units="atom/b-cm" />
|
||||
<nuclide name="He4" xs="71c" ao="1.0" />
|
||||
</material>
|
||||
|
||||
<!-- Zircaloy 4 -->
|
||||
<material id="3">
|
||||
<density value="6.55" units="g/cc" />
|
||||
<nuclide name="O16" xs="71c" wo="0.00125" />
|
||||
<nuclide name="Nb93" xs="71c" wo="0.01" />
|
||||
<nuclide name="Sn112" xs="71c" wo="0.000097" />
|
||||
<nuclide name="Sn114" xs="71c" wo="0.000066" />
|
||||
<nuclide name="Sn115" xs="71c" wo="0.000034" />
|
||||
<nuclide name="Sn116" xs="71c" wo="0.001454" />
|
||||
<nuclide name="Sn117" xs="71c" wo="0.000768" />
|
||||
<nuclide name="Sn118" xs="71c" wo="0.002422 " />
|
||||
<nuclide name="Sn119" xs="71c" wo="0.000859" />
|
||||
<nuclide name="Sn120" xs="71c" wo="0.003258" />
|
||||
<nuclide name="Sn122" xs="71c" wo="0.000463" />
|
||||
<nuclide name="Sn124" xs="71c" wo="0.000579" />
|
||||
<nuclide name="Fe54" xs="71c" wo="0.00005845" />
|
||||
<nuclide name="Fe56" xs="71c" wo="0.00091754" />
|
||||
<nuclide name="Fe57" xs="71c" wo="0.00002119" />
|
||||
<nuclide name="Fe58" xs="71c" wo="0.00000282" />
|
||||
<nuclide name="Zr90" xs="71c" wo="0.503052375" />
|
||||
<nuclide name="Zr91" xs="71c" wo="0.10970355" />
|
||||
<nuclide name="Zr92" xs="71c" wo="0.167684125" />
|
||||
<nuclide name="Zr94" xs="71c" wo="0.16993295" />
|
||||
<nuclide name="Zr96" xs="71c" wo="0.027377" />
|
||||
</material>
|
||||
|
||||
<!-- Water -->
|
||||
<material id="4">
|
||||
<density value="0.73" units="g/cc" />
|
||||
<nuclide name="H1" xs="71c" ao="2.0" />
|
||||
<nuclide name="O16" xs="71c" ao="1.0" />
|
||||
<sab name="c_H_in_H2O" xs="04t" />
|
||||
</material>
|
||||
|
||||
<!-- Guide Tube -->
|
||||
<material id="5">
|
||||
<density value="7.92" units="g/cc" />
|
||||
<nuclide name="Ni58" xs="71c" wo="0.06467315" />
|
||||
<nuclide name="Ni60" xs="71c" wo="0.02492135" />
|
||||
<nuclide name="Ni61" xs="71c" wo="0.001082905" />
|
||||
<nuclide name="Ni62" xs="71c" wo="0.00345287" />
|
||||
<nuclide name="Ni64" xs="71c" wo="0.000879225" />
|
||||
<nuclide name="Fe54" xs="71c" wo="0.04062275" />
|
||||
<nuclide name="Fe56" xs="71c" wo="0.63762775" />
|
||||
<nuclide name="Fe57" xs="71c" wo="0.01472705" />
|
||||
<nuclide name="Fe58" xs="71c" wo="0.0019599" />
|
||||
<nuclide name="Mn55" xs="71c" wo="0.020" />
|
||||
<nuclide name="Cr50" xs="71c" wo="0.0082555" />
|
||||
<nuclide name="Cr52" xs="71c" wo="0.1591991" />
|
||||
<nuclide name="Cr53" xs="71c" wo="0.0180519" />
|
||||
<nuclide name="Cr54" xs="71c" wo="0.0044935" />
|
||||
</material>
|
||||
|
||||
</materials>
|
||||
|
|
@ -1,24 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<plots>
|
||||
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
|
||||
<plot id="1" color="mat" >
|
||||
<origin>0. 0. 0.</origin>
|
||||
<width>25.50 25.50 </width>
|
||||
<pixels>1000 1000 </pixels>
|
||||
<col_spec id="1" rgb="123 123 231" />
|
||||
|
||||
|
||||
</plot>
|
||||
|
||||
</plots>
|
||||
|
|
@ -1,29 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<settings>
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
<eigenvalue>
|
||||
<batches>600</batches>
|
||||
<inactive>50</inactive>
|
||||
<particles>10000</particles>
|
||||
</eigenvalue>
|
||||
|
||||
<source>
|
||||
<space>
|
||||
<type>box</type>
|
||||
<parameters>
|
||||
-10.659 -10.659 -182.9
|
||||
10.659 10.659 182.9
|
||||
</parameters>
|
||||
</space>
|
||||
</source>
|
||||
|
||||
</settings>
|
||||
|
|
@ -1,77 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<tallies>
|
||||
|
||||
<filter id="1" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="2" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="3" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="4" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="5" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="6" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="7" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="8" type="delayedgroup">
|
||||
<bins>1 2 3 4 5 6</bins>
|
||||
</filter>
|
||||
|
||||
<tally id="1">
|
||||
<filters>1</filters>
|
||||
<scores>flux</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="2">
|
||||
<filters>2</filters>
|
||||
<scores>nu-fission</scores>
|
||||
</tally>
|
||||
|
||||
|
||||
<tally id="3">
|
||||
<filters>3</filters>
|
||||
<scores>delayed-nu-fission</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="4">
|
||||
<filters>4</filters>
|
||||
<scores>fission</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="5">
|
||||
<filters>5</filters>
|
||||
<scores>absorption</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="6">
|
||||
<filters>6</filters>
|
||||
<scores>scatter</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="7">
|
||||
<filters>7</filters>
|
||||
<scores>total</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="8">
|
||||
<filters>8</filters>
|
||||
<scores>delayed-nu-fission</scores>
|
||||
</tally>
|
||||
|
||||
</tallies>
|
||||
|
|
@ -1,5 +1,22 @@
|
|||
#!/usr/bin/env python3
|
||||
|
||||
"""
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Conversion: Adam Parler
|
||||
Date: 12/03/2023
|
||||
====================================================
|
||||
"""
|
||||
|
||||
import os
|
||||
import openmc
|
||||
|
||||
try:
|
||||
os.mkdir(os.path.dirname(os.path.realpath(__file__)) + '/problem')
|
||||
except OSError: pass
|
||||
|
||||
################## MATERIALS ##################
|
||||
|
||||
fuel_26 = openmc.Material()
|
||||
|
|
@ -165,7 +182,7 @@ plot_file.export_to_xml(path="problem/plots.xml")
|
|||
################## SETTINGS ##################
|
||||
|
||||
settings_file = openmc.Settings()
|
||||
settings_file.source = openmc.Source(space=openmc.stats.Box((-10.659, -10.659, -182.9),
|
||||
settings_file.source = openmc.IndependentSource(space=openmc.stats.Box((-10.659, -10.659, -182.9),
|
||||
(10.659, 10.659, 182.9)))
|
||||
settings_file.batches = 600
|
||||
settings_file.inactive = 50
|
||||
|
|
|
|||
|
|
@ -1,70 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<geometry>
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
<cell id="1" fill="5" region="1 -2 3 -4 10 -11" />
|
||||
<cell id="10" universe="1" material="1" region="-5" />
|
||||
<cell id="20" universe="1" material="2" region="5 -6 " />
|
||||
<cell id="30" universe="1" material="3" region="6 -7" />
|
||||
<cell id="40" universe="1" material="4" region="7" />
|
||||
|
||||
<cell id="50" universe="2" material="4" region="-8" />
|
||||
<cell id="60" universe="2" material="5" region="8 -9" />
|
||||
<cell id="70" universe="2" material="4" region="9" />
|
||||
|
||||
<lattice id="5" >
|
||||
<dimension>17 17</dimension>
|
||||
<lower_left>-10.659 -10.659</lower_left>
|
||||
<pitch>1.26 1.26</pitch>
|
||||
<universes>
|
||||
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 2 1 1 2 1 1 2 1 1 1 1 1
|
||||
1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
|
||||
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
|
||||
1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 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 2 1 1 2 1 1 2 1 1 2 1 1 2 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 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1
|
||||
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
|
||||
1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1
|
||||
1 1 1 1 1 2 1 1 2 1 1 2 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
|
||||
</universes>
|
||||
</lattice>
|
||||
|
||||
<surface id="1" type="x-plane" coeffs="-10.659" boundary="reflective" />
|
||||
<surface id="2" type="x-plane" coeffs="10.659" boundary="reflective" />
|
||||
<surface id="3" type="y-plane" coeffs="-10.659" boundary="reflective" />
|
||||
<surface id="4" type="y-plane" coeffs="10.659" boundary="reflective" />
|
||||
|
||||
<surface id="5" type="z-cylinder" coeffs="0. 0. 0.4095" />
|
||||
<surface id="6" type="z-cylinder" coeffs="0. 0. 0.4177" />
|
||||
<surface id="7" type="z-cylinder" coeffs="0. 0. 0.4749" />
|
||||
<surface id="8" type="z-cylinder" coeffs="0. 0. 0.5715" />
|
||||
<surface id="9" type="z-cylinder" coeffs="0. 0. 0.612" />
|
||||
|
||||
<surface id="10" type="z-plane" coeffs="-182.9" boundary="reflective" />
|
||||
<surface id="11" type="z-plane" coeffs="182.9" boundary="reflective" />
|
||||
|
||||
<surface id="12" type="x-plane" coeffs="-0.627" />
|
||||
<surface id="13" type="x-plane" coeffs="0.627" />
|
||||
<surface id="14" type="y-plane" coeffs="-0.627" />
|
||||
<surface id="15" type="y-plane" coeffs="0.627" />
|
||||
|
||||
</geometry>
|
||||
|
||||
|
||||
|
||||
|
|
@ -1,84 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<materials>
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
<cross_sections>/opt/xdata/endfb-vii.1-hdf5/cross_sections.xml</cross_sections>
|
||||
|
||||
<default_xs>71c</default_xs>
|
||||
|
||||
<!-- Fuel composed of U-238 O-16 and U-235 -->
|
||||
<material id="1">
|
||||
<density value="10.5" units="g/cc" />
|
||||
<nuclide name="U235" xs="71c" wo="0.0273246599755" />
|
||||
<nuclide name="U238" xs="71c" wo="0.854115984395" />
|
||||
<nuclide name="O16" xs="71c" wo="0.11855935563" />
|
||||
</material>
|
||||
|
||||
<!-- Helium for gap -->
|
||||
<material id="2">
|
||||
<density value="1.00000e-4" units="atom/b-cm" />
|
||||
<nuclide name="He4" xs="71c" ao="1.0" />
|
||||
</material>
|
||||
|
||||
<!-- Zircaloy 4 -->
|
||||
<material id="3">
|
||||
<density value="6.55" units="g/cc" />
|
||||
<nuclide name="O16" xs="71c" wo="0.00125" />
|
||||
<nuclide name="Nb93" xs="71c" wo="0.01" />
|
||||
<nuclide name="Sn112" xs="71c" wo="0.000097" />
|
||||
<nuclide name="Sn114" xs="71c" wo="0.000066" />
|
||||
<nuclide name="Sn115" xs="71c" wo="0.000034" />
|
||||
<nuclide name="Sn116" xs="71c" wo="0.001454" />
|
||||
<nuclide name="Sn117" xs="71c" wo="0.000768" />
|
||||
<nuclide name="Sn118" xs="71c" wo="0.002422 " />
|
||||
<nuclide name="Sn119" xs="71c" wo="0.000859" />
|
||||
<nuclide name="Sn120" xs="71c" wo="0.003258" />
|
||||
<nuclide name="Sn122" xs="71c" wo="0.000463" />
|
||||
<nuclide name="Sn124" xs="71c" wo="0.000579" />
|
||||
<nuclide name="Fe54" xs="71c" wo="0.00005845" />
|
||||
<nuclide name="Fe56" xs="71c" wo="0.00091754" />
|
||||
<nuclide name="Fe57" xs="71c" wo="0.00002119" />
|
||||
<nuclide name="Fe58" xs="71c" wo="0.00000282" />
|
||||
<nuclide name="Zr90" xs="71c" wo="0.503052375" />
|
||||
<nuclide name="Zr91" xs="71c" wo="0.10970355" />
|
||||
<nuclide name="Zr92" xs="71c" wo="0.167684125" />
|
||||
<nuclide name="Zr94" xs="71c" wo="0.16993295" />
|
||||
<nuclide name="Zr96" xs="71c" wo="0.027377" />
|
||||
</material>
|
||||
|
||||
<!-- Water -->
|
||||
<material id="4">
|
||||
<density value="0.73" units="g/cc" />
|
||||
<nuclide name="H1" xs="71c" ao="2.0" />
|
||||
<nuclide name="O16" xs="71c" ao="1.0" />
|
||||
<sab name="c_H_in_H2O" xs="04t" />
|
||||
</material>
|
||||
|
||||
<!-- Guide Tube -->
|
||||
<material id="5">
|
||||
<density value="7.92" units="g/cc" />
|
||||
<nuclide name="Ni58" xs="71c" wo="0.06467315" />
|
||||
<nuclide name="Ni60" xs="71c" wo="0.02492135" />
|
||||
<nuclide name="Ni61" xs="71c" wo="0.001082905" />
|
||||
<nuclide name="Ni62" xs="71c" wo="0.00345287" />
|
||||
<nuclide name="Ni64" xs="71c" wo="0.000879225" />
|
||||
<nuclide name="Fe54" xs="71c" wo="0.04062275" />
|
||||
<nuclide name="Fe56" xs="71c" wo="0.63762775" />
|
||||
<nuclide name="Fe57" xs="71c" wo="0.01472705" />
|
||||
<nuclide name="Fe58" xs="71c" wo="0.0019599" />
|
||||
<nuclide name="Mn55" xs="71c" wo="0.020" />
|
||||
<nuclide name="Cr50" xs="71c" wo="0.0082555" />
|
||||
<nuclide name="Cr52" xs="71c" wo="0.1591991" />
|
||||
<nuclide name="Cr53" xs="71c" wo="0.0180519" />
|
||||
<nuclide name="Cr54" xs="71c" wo="0.0044935" />
|
||||
</material>
|
||||
|
||||
</materials>
|
||||
|
|
@ -1,24 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<plots>
|
||||
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
|
||||
<plot id="1" color="mat" >
|
||||
<origin>0. 0. 0.</origin>
|
||||
<width>25.50 25.50 </width>
|
||||
<pixels>1000 1000 </pixels>
|
||||
<col_spec id="1" rgb="123 123 231" />
|
||||
|
||||
|
||||
</plot>
|
||||
|
||||
</plots>
|
||||
|
|
@ -1,29 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<settings>
|
||||
|
||||
<!--
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Date: 02/03/2017
|
||||
====================================================
|
||||
-->
|
||||
|
||||
<eigenvalue>
|
||||
<batches>600</batches>
|
||||
<inactive>50</inactive>
|
||||
<particles>10000</particles>
|
||||
</eigenvalue>
|
||||
|
||||
<source>
|
||||
<space>
|
||||
<type>box</type>
|
||||
<parameters>
|
||||
-10.659 -10.659 -182.9
|
||||
10.659 10.659 182.9
|
||||
</parameters>
|
||||
</space>
|
||||
</source>
|
||||
|
||||
</settings>
|
||||
|
|
@ -1,76 +0,0 @@
|
|||
<?xml version="1.0"?>
|
||||
<tallies>
|
||||
|
||||
<filter id="1" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="2" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="3" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="4" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="5" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="6" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="7" type="cell">
|
||||
<bins>1</bins>
|
||||
</filter>
|
||||
|
||||
<filter id="8" type="delayedgroup">
|
||||
<bins>1 2 3 4 5 6</bins>
|
||||
</filter>
|
||||
|
||||
<tally id="1">
|
||||
<filters>1</filters>
|
||||
<scores>flux</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="2">
|
||||
<filters>2</filters>
|
||||
<scores>nu-fission</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="3">
|
||||
<filters>3</filters>
|
||||
<scores>delayed-nu-fission</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="4">
|
||||
<filters>4</filters>
|
||||
<scores>fission</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="5">
|
||||
<filters>5</filters>
|
||||
<scores>absorption</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="6">
|
||||
<filters>6</filters>
|
||||
<scores>scatter</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="7">
|
||||
<filters>7</filters>
|
||||
<scores>total</scores>
|
||||
</tally>
|
||||
|
||||
<tally id="8">
|
||||
<filters>8</filters>
|
||||
<scores>delayed-nu-fission</scores>
|
||||
</tally>
|
||||
|
||||
</tallies>
|
||||
|
|
@ -1,5 +1,22 @@
|
|||
#!/usr/bin/env python3
|
||||
|
||||
"""
|
||||
====================================================
|
||||
Description: AP1000 PWR Assembly 17*17
|
||||
Case: Problem 2 (17*17 PWR Assembly)
|
||||
Written By: Ahmed K Madani
|
||||
Conversion: Adam Parler
|
||||
Date: 12/03/2023
|
||||
====================================================
|
||||
"""
|
||||
|
||||
import os
|
||||
import openmc
|
||||
|
||||
try:
|
||||
os.mkdir(os.path.dirname(os.path.realpath(__file__)) + '/problem')
|
||||
except OSError: pass
|
||||
|
||||
################## MATERIALS ##################
|
||||
|
||||
fuel_31 = openmc.Material()
|
||||
|
|
@ -165,7 +182,7 @@ plot_file.export_to_xml(path="problem/plots.xml")
|
|||
################## SETTINGS ##################
|
||||
|
||||
settings_file = openmc.Settings()
|
||||
settings_file.source = openmc.Source(space=openmc.stats.Box((-10.659, -10.659, -182.9),
|
||||
settings_file.source = openmc.IndependentSource(space=openmc.stats.Box((-10.659, -10.659, -182.9),
|
||||
(10.659, 10.659, 182.9)))
|
||||
settings_file.batches = 600
|
||||
settings_file.inactive = 50
|
||||
|
|
|
|||
Loading…
Add table
Add a link
Reference in a new issue