Commit 02a675d8 authored by Ronald Rahaman's avatar Ronald Rahaman
Browse files

Changed temperature settings and keff in OpenMC model

parent 6feae859
......@@ -69,7 +69,7 @@ rescaled_file = 'pebble_centers_rescaled.txt'
np.savetxt(rescaled_file, pebble_centers.reshape((-1,3)))
print("Saved rescaled centers to " + rescaled_file)
# Vessel
extra_thickness = 7.0 # change this value to set keff=1
extra_thickness = 8.125 # change this value to set keff=1
x_vessel = 0.0 # x position vessel
y_vessel = 0.0 # y position vessel
z1_vessel = -2.54 - extra_thickness # bottom plane vessel ; scaled by TAMU experiment factor
......@@ -262,6 +262,8 @@ settings.source = openmc.Source(space=openmc.stats.Box(*c_vessel.bounding_box
settings.particles = 10000
settings.inactive = 50
settings.batches = 150
settings.temperature = dict(default=573, method='interpolation',
multipole=True, range=(300.0, 1500.0), tolerance=1000.0)
# Fuel volume calculation
volume_fuel = openmc.VolumeCalculation([m_fuel], 10000000, *c_pebble_central.bounding_box)
settings.volume_calculations = [volume_fuel]
......
This diff is collapsed.
......@@ -5,40 +5,48 @@
<nuclide name="U234" wo="0.03565284" />
<nuclide name="U235" wo="17.790767159999998" />
<nuclide name="U238" wo="71.75358" />
<nuclide name="C0" wo="1.8" />
<nuclide name="O16" wo="8.616527990092724" />
<nuclide name="O17" wo="0.003472009907274255" />
<nuclide name="C12" wo="1.778412321154449" />
<nuclide name="C13" wo="0.021587678845551017" />
<nuclide name="O16" wo="8.597130043616898" />
<nuclide name="O17" wo="0.0034711398506764002" />
<nuclide name="O18" wo="0.019398816532423476" />
</material>
<material id="2" name="m_graphite_c_buffer - triso partciles">
<density units="g/cm3" value="1.0" />
<nuclide ao="1.0" name="C0" />
<nuclide ao="0.988922" name="C12" />
<nuclide ao="0.011078" name="C13" />
<sab name="c_Graphite" />
</material>
<material id="3" name="m_graphite_pyc - triso particles">
<density units="g/cm3" value="1.87" />
<nuclide ao="1.0" name="C0" />
<nuclide ao="0.988922" name="C12" />
<nuclide ao="0.011078" name="C13" />
<sab name="c_Graphite" />
</material>
<material id="4" name="sic - triso partciles">
<density units="g/cm3" value="3.2" />
<nuclide ao="1.0" name="C0" />
<nuclide ao="0.9222968" name="Si28" />
<nuclide ao="0.988922" name="C12" />
<nuclide ao="0.011078" name="C13" />
<nuclide ao="0.0468316" name="Si29" />
<nuclide ao="0.9222968" name="Si28" />
<nuclide ao="0.0308716" name="Si30" />
</material>
<material id="5" name="m_graphite_matrix - triso particles">
<density units="g/cm3" value="1.6" />
<nuclide ao="1.0" name="C0" />
<nuclide ao="0.988922" name="C12" />
<nuclide ao="0.011078" name="C13" />
<sab name="c_Graphite" />
</material>
<material id="6" name="m_graphite_inner - pebble inner zone">
<density units="g/cm3" value="1.59368" />
<nuclide ao="1.0" name="C0" />
<nuclide ao="0.988922" name="C12" />
<nuclide ao="0.011078" name="C13" />
<sab name="c_Graphite" />
</material>
<material id="7" name="m_graphite_outer - pebble outer zone">
<density units="g/cm3" value="1.74" />
<nuclide ao="1.0" name="C0" />
<nuclide ao="0.988922" name="C12" />
<nuclide ao="0.011078" name="C13" />
<sab name="c_Graphite" />
</material>
<material id="8" name="m_flibe - 2LiF-BeF2">
......
......@@ -2,7 +2,7 @@
<plots>
<plot basis="xy" color_by="material" filename="plot1" id="1" type="slice">
<origin>0 0 0</origin>
<width>28.0 28.0</width>
<width>30.25 30.25</width>
<pixels>1000 1000</pixels>
<color id="1" rgb="165 42 42" />
<color id="2" rgb="176 196 222" />
......@@ -15,7 +15,7 @@
</plot>
<plot basis="xy" color_by="material" filename="plot2" id="2" type="slice">
<origin>0 0 1</origin>
<width>28.0 28.0</width>
<width>30.25 30.25</width>
<pixels>1000 1000</pixels>
<color id="1" rgb="165 42 42" />
<color id="2" rgb="176 196 222" />
......@@ -28,7 +28,7 @@
</plot>
<plot basis="xy" color_by="material" filename="plot3" id="3" type="slice">
<origin>0 0 1</origin>
<width>28.0 28.0</width>
<width>30.25 30.25</width>
<pixels>1000 1000</pixels>
<color id="1" rgb="165 42 42" />
<color id="2" rgb="176 196 222" />
......@@ -41,7 +41,7 @@
</plot>
<plot basis="xz" color_by="material" filename="plot4" id="4" type="slice">
<origin>0 0 5.08</origin>
<width>29.24 29.24</width>
<width>31.49 31.49</width>
<pixels>1000 1000</pixels>
<color id="1" rgb="165 42 42" />
<color id="2" rgb="176 196 222" />
......
<?xml version='1.0' encoding='utf-8'?>
<settings>
<run_mode>eigenvalue</run_mode>
<particles>200000</particles>
<batches>100</batches>
<particles>10000</particles>
<batches>150</batches>
<inactive>50</inactive>
<source strength="1.0">
<space type="box">
<parameters>-14.0 -14.0 -9.54 14.0 14.0 19.7</parameters>
<space type="fission">
<parameters>-15.125 -15.125 -10.665 15.125 15.125 20.825</parameters>
</space>
</source>
<state_point>
<batches />
</state_point>
<temperature_default>573</temperature_default>
<temperature_method>interpolation</temperature_method>
<temperature_multipole>true</temperature_multipole>
<temperature_range>300.0 1500.0</temperature_range>
<temperature_tolerance>1000.0</temperature_tolerance>
<volume_calc>
<domain_type>material</domain_type>
<domain_ids>1</domain_ids>
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment