GitLab maintenance scheduled form Friday, 2021-06-18 5:00pm to Satursday, 2021-06-19 10:00pm CT - Services will be unavailable during this time.

Commit fcbaafde authored by Elia Merzari's avatar Elia Merzari

Merge branch '10k_openmc_model' into 'master'

10k openmc model

See merge request !5
parents fb643297 478349ca
This diff is collapsed.
# Files
## dat_to_csv.py
Used to convert the initial data file for the pebble centers to a CSV format.
## openmc_model.py
Creates an OpenMC model using the pebble centers found in the pebble_centers.txt file.
import numpy as np
header_length = 8
with open('11145_pebble_yuhsiang_E3p5M_peb.11k.dat', 'r') as fh:
for i in range(header_length):
fh.readline()
centers = np.fromfile(fh, sep=' ', dtype=float)
centers = centers.reshape((centers.size // 3, 3))
centers *= 1.5
np.savetxt('centers.txt', centers)
This diff is collapsed.
<?xml version='1.0' encoding='utf-8'?>
<materials>
<material depletable="true" id="1" name="m_fuel - uranium oxycarbide - triso partciles">
<density units="g/cm3" value="10.5" />
<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" />
</material>
<material id="2" name="m_graphite_c_buffer - triso partciles">
<density units="g/cm3" value="1.0" />
<nuclide ao="1.0" name="C0" />
<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" />
<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.0468316" name="Si29" />
<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" />
<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" />
<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" />
<sab name="c_Graphite" />
</material>
<material id="8" name="m_flibe - 2LiF-BeF2">
<density units="kg/m3" value="1962.13354" />
<nuclide ao="1.9999" name="Li7" />
<nuclide ao="9.999999999998899e-05" name="Li6" />
<nuclide ao="1.0" name="Be9" />
<nuclide ao="4.0" name="F19" />
</material>
</materials>
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
<?xml version='1.0' encoding='utf-8'?>
<plots>
<plot basis="xy" color_by="material" filename="plot1" id="1" type="slice">
<origin>0 0 0</origin>
<width>153.81691779124347 153.81691779124347</width>
<pixels>1000 1000</pixels>
<color id="1" rgb="165 42 42" />
<color id="2" rgb="176 196 222" />
<color id="3" rgb="0 0 255" />
<color id="4" rgb="255 165 0" />
<color id="5" rgb="0 255 255" />
<color id="6" rgb="0 191 255" />
<color id="7" rgb="0 0 128" />
<color id="8" rgb="255 255 0" />
</plot>
<plot basis="xy" color_by="material" filename="plot2" id="2" type="slice">
<origin>0 0 1</origin>
<width>153.81691779124347 153.81691779124347</width>
<pixels>1000 1000</pixels>
<color id="1" rgb="165 42 42" />
<color id="2" rgb="176 196 222" />
<color id="3" rgb="0 0 255" />
<color id="4" rgb="255 165 0" />
<color id="5" rgb="0 255 255" />
<color id="6" rgb="0 191 255" />
<color id="7" rgb="0 0 128" />
<color id="8" rgb="255 255 0" />
</plot>
<plot basis="xy" color_by="material" filename="plot3" id="3" type="slice">
<origin>0 0 1</origin>
<width>153.81691779124347 153.81691779124347</width>
<pixels>1000 1000</pixels>
<color id="1" rgb="165 42 42" />
<color id="2" rgb="176 196 222" />
<color id="3" rgb="0 0 255" />
<color id="4" rgb="255 165 0" />
<color id="5" rgb="0 255 255" />
<color id="6" rgb="0 191 255" />
<color id="7" rgb="0 0 128" />
<color id="8" rgb="255 255 0" />
</plot>
<plot basis="xz" color_by="material" filename="plot4" id="4" type="slice">
<origin>0 0 0.02052855750000049</origin>
<width>153.81691779124347 153.81691779124347</width>
<pixels>1000 1000</pixels>
<color id="1" rgb="165 42 42" />
<color id="2" rgb="176 196 222" />
<color id="3" rgb="0 0 255" />
<color id="4" rgb="255 165 0" />
<color id="5" rgb="0 255 255" />
<color id="6" rgb="0 191 255" />
<color id="7" rgb="0 0 128" />
<color id="8" rgb="255 255 0" />
</plot>
</plots>
<?xml version='1.0' encoding='utf-8'?>
<settings>
<run_mode>eigenvalue</run_mode>
<particles>50000</particles>
<batches>250</batches>
<inactive>50</inactive>
<source strength="1.0">
<space type="box">
<parameters>-76.90845889562173 -76.90845889562173 -30.41567061 76.90845889562173 76.90845889562173 30.456727725</parameters>
</space>
</source>
<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>
<samples>10000000</samples>
<lower_left>-1.4 -1.4 -1.4</lower_left>
<upper_right>1.4 1.4 1.4</upper_right>
</volume_calc>
<material_cell_offsets>false</material_cell_offsets>
</settings>
This diff is collapsed.
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