Adding test for hex lattice with coincident interior surface.

This commit is contained in:
Patrick Shriwise 2019-04-01 17:30:12 -05:00
parent da9a876185
commit 8673c924ac
5 changed files with 124 additions and 0 deletions

View file

@ -0,0 +1,34 @@
<?xml version='1.0' encoding='utf-8'?>
<geometry>
<cell id="1" material="1" region="-1 2 -3" universe="1" />
<cell id="2" material="3" region="1" universe="1" />
<cell id="3" material="6" name="coolant" region="-4 2 -3" universe="2" />
<cell id="4" material="7" name="zirconium_shell" region="4 -5 2 -3" universe="2" />
<cell id="5" material="4" name="lead_shell" region="5 -6 2 -3" universe="2" />
<cell id="6" material="3" name="matrix coolant surround" region="6 2 -3" universe="2" />
<cell id="7" material="3" universe="3" />
<cell fill="4" id="8" name="container cell" region="-7 8 -9 10 11 -12 2 -3" universe="5" />
<hex_lattice id="4" n_rings="2" name="regular fuel assembly">
<pitch>1.4</pitch>
<outer>3</outer>
<center>0.0 0.0</center>
<universes>
2
2 2
1
2 2
2</universes>
</hex_lattice>
<surface coeffs="0.0 0.0 0.7" id="1" type="z-cylinder" />
<surface boundary="reflective" coeffs="0.0" id="2" type="z-plane" />
<surface boundary="reflective" coeffs="10.0" id="3" type="z-plane" />
<surface coeffs="0.0 0.0 0.293" id="4" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.35" id="5" type="z-cylinder" />
<surface coeffs="0.0 0.0 0.352" id="6" type="z-cylinder" />
<surface boundary="reflective" coeffs="1.4" id="7" type="y-plane" />
<surface boundary="reflective" coeffs="-1.4" id="8" type="y-plane" />
<surface boundary="reflective" coeffs="1.7320508075688772 1.0 0.0 2.8" id="9" type="plane" />
<surface boundary="reflective" coeffs="-1.7320508075688772 1.0 0.0 -2.8" id="10" type="plane" />
<surface boundary="reflective" coeffs="1.7320508075688772 1.0 0.0 -2.8" id="11" type="plane" />
<surface boundary="reflective" coeffs="-1.7320508075688772 1.0 0.0 2.8" id="12" type="plane" />
</geometry>

View file

@ -0,0 +1,62 @@
<?xml version='1.0' encoding='utf-8'?>
<materials>
<material depletable="true" id="1" name="fuel1">
<density units="sum" />
<nuclide ao="0.0049817" name="U235" />
<nuclide ao="0.019989" name="U238" />
<nuclide ao="0.037372" name="O16" />
<nuclide ao="0.0099659" name="C0" />
</material>
<material id="2" name="matrix_pin (no boron)">
<density units="atom/b-cm" value="0.087742" />
<nuclide ao="1.0" name="C0" />
<sab name="c_Graphite" />
</material>
<material id="3" name="matrix (boron)">
<density units="atom/b-cm" value="0.087742" />
<nuclide ao="0.999985" name="C0" />
<nuclide ao="1.5e-05" name="B10" />
<sab name="c_Graphite" />
</material>
<material id="4" name="Lead">
<density units="g/cm3" value="10.32" />
<nuclide ao="0.014" name="Pb204" />
<nuclide ao="0.241" name="Pb206" />
<nuclide ao="0.221" name="Pb207" />
<nuclide ao="0.524" name="Pb208" />
</material>
<material id="5" name="intermodule">
<density units="atom/b-cm" value="5e-05" />
<nuclide ao="1.0" name="N14" />
</material>
<material id="6" name="coolant">
<density units="atom/b-cm" value="0.00054464" />
<nuclide ao="1.0" name="He4" />
</material>
<material id="7" name="Zirc4">
<density units="sum" />
<nuclide ao="8.549e-06" name="Fe54" />
<nuclide ao="0.0001294" name="Fe56" />
<nuclide ao="2.936e-06" name="Fe57" />
<nuclide ao="3.84e-07" name="Fe58" />
<nuclide ao="3.431e-06" name="Cr50" />
<nuclide ao="6.363e-05" name="Cr52" />
<nuclide ao="7.079e-06" name="Cr53" />
<nuclide ao="1.73e-06" name="Cr54" />
<nuclide ao="0.02217" name="Zr90" />
<nuclide ao="0.004781" name="Zr91" />
<nuclide ao="0.007228" name="Zr92" />
<nuclide ao="0.007169" name="Zr94" />
<nuclide ao="0.001131" name="Zr96" />
<nuclide ao="5.129e-06" name="Sn112" />
<nuclide ao="3.428e-06" name="Sn114" />
<nuclide ao="1.751e-06" name="Sn115" />
<nuclide ao="7.423e-05" name="Sn116" />
<nuclide ao="3.887e-05" name="Sn117" />
<nuclide ao="0.0001215" name="Sn118" />
<nuclide ao="4.274e-05" name="Sn119" />
<nuclide ao="0.0001608" name="Sn120" />
<nuclide ao="2.247e-05" name="Sn122" />
<nuclide ao="2.765e-05" name="Sn124" />
</material>
</materials>

View file

@ -0,0 +1,2 @@
k-combined:
1.328686E+00 1.538567E-02

View file

@ -0,0 +1,20 @@
<?xml version='1.0' encoding='utf-8'?>
<settings>
<run_mode>eigenvalue</run_mode>
<particles>1000</particles>
<batches>10</batches>
<inactive>5</inactive>
<source strength="1.0">
<space type="box">
<parameters>-0.9899494936611666 -0.9899494936611666 0.0 0.9899494936611666 0.9899494936611666 10.0</parameters>
</space>
</source>
<output>
<summary>false</summary>
</output>
<seed>46</seed>
<temperature_method>interpolation</temperature_method>
<resonance_scattering>
<enable>true</enable>
</resonance_scattering>
</settings>

View file

@ -0,0 +1,6 @@
from tests.testing_harness import TestHarness
def test_lattice_hex():
harness = TestHarness('statepoint.10.h5')
harness.main()