Skip to content

Commit e860ea9

Browse files
pshriwisepshriwise
committed
Adding XDG unstructured mesh tallies mirroring those of the current unstructured mesh capabilities.
1 parent 612c891 commit e860ea9

15 files changed

Lines changed: 1068 additions & 0 deletions

‎tests/regression_tests/xdg/tallies/inputs_true0_libmesh.dat‎

Lines changed: 90 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,90 @@
1+
<?xml version='1.0' encoding='utf-8'?>
2+
<model>
3+
<materials>
4+
<material id="1" name="fuel" depletable="true">
5+
<density value="4.5" units="g/cc"/>
6+
<nuclide name="U235" ao="1.0"/>
7+
</material>
8+
<material id="2" name="zircaloy">
9+
<density value="5.77" units="g/cc"/>
10+
<nuclide name="Zr90" ao="0.5145"/>
11+
<nuclide name="Zr91" ao="0.1122"/>
12+
<nuclide name="Zr92" ao="0.1715"/>
13+
<nuclide name="Zr94" ao="0.1738"/>
14+
<nuclide name="Zr96" ao="0.028"/>
15+
</material>
16+
<material id="3" name="water">
17+
<density value="0.07416" units="atom/b-cm"/>
18+
<nuclide name="H1" ao="2.0"/>
19+
<nuclide name="O16" ao="1.0"/>
20+
</material>
21+
</materials>
22+
<geometry>
23+
<cell id="1" name="fuel" material="1" region="1 -2 3 -4 5 -6" universe="1"/>
24+
<cell id="2" name="clad" material="2" region="(-1 | 2 | -3 | 4 | -5 | 6) (7 -8 9 -10 11 -12)" universe="1"/>
25+
<cell id="3" name="water" material="3" region="(-7 | 8 | -9 | 10 | -11 | 12) (13 -14 15 -16 17 -18)" universe="1"/>
26+
<surface id="1" name="minimum x" type="x-plane" coeffs="-5.0"/>
27+
<surface id="2" name="maximum x" type="x-plane" coeffs="5.0"/>
28+
<surface id="3" name="minimum y" type="y-plane" coeffs="-5.0"/>
29+
<surface id="4" name="maximum y" type="y-plane" coeffs="5.0"/>
30+
<surface id="5" name="minimum z" type="z-plane" coeffs="-5.0"/>
31+
<surface id="6" name="maximum z" type="z-plane" coeffs="5.0"/>
32+
<surface id="7" name="minimum x" type="x-plane" coeffs="-6.0"/>
33+
<surface id="8" name="maximum x" type="x-plane" coeffs="6.0"/>
34+
<surface id="9" name="minimum y" type="y-plane" coeffs="-6.0"/>
35+
<surface id="10" name="maximum y" type="y-plane" coeffs="6.0"/>
36+
<surface id="11" name="minimum z" type="z-plane" coeffs="-6.0"/>
37+
<surface id="12" name="maximum z" type="z-plane" coeffs="6.0"/>
38+
<surface id="13" name="minimum x" type="x-plane" boundary="vacuum" coeffs="-10.0"/>
39+
<surface id="14" name="maximum x" type="x-plane" boundary="vacuum" coeffs="10.0"/>
40+
<surface id="15" name="minimum y" type="y-plane" boundary="vacuum" coeffs="-10.0"/>
41+
<surface id="16" name="maximum y" type="y-plane" boundary="vacuum" coeffs="10.0"/>
42+
<surface id="17" name="minimum z" type="z-plane" boundary="vacuum" coeffs="-10.0"/>
43+
<surface id="18" name="maximum z" type="z-plane" boundary="vacuum" coeffs="10.0"/>
44+
</geometry>
45+
<settings>
46+
<run_mode>fixed source</run_mode>
47+
<particles>1000</particles>
48+
<batches>10</batches>
49+
<source type="independent" strength="1.0" particle="neutron">
50+
<space type="spherical" origin="0.0 0.0 0.0">
51+
<r type="uniform" parameters="0.0 0.0"/>
52+
<cos_theta type="discrete">
53+
<parameters>1.0 1.0</parameters>
54+
</cos_theta>
55+
<phi type="discrete">
56+
<parameters>0.0 1.0</parameters>
57+
</phi>
58+
</space>
59+
<energy type="discrete">
60+
<parameters>15000000.0 1.0</parameters>
61+
</energy>
62+
</source>
63+
</settings>
64+
<tallies>
65+
<mesh id="1">
66+
<dimension>10 10 10</dimension>
67+
<lower_left>-10.0 -10.0 -10.0</lower_left>
68+
<upper_right>10.0 10.0 10.0</upper_right>
69+
</mesh>
70+
<mesh id="2" type="xdg" library="libmesh">
71+
<filename>test_mesh_tets.e</filename>
72+
</mesh>
73+
<filter id="1" type="mesh">
74+
<bins>1</bins>
75+
</filter>
76+
<filter id="2" type="mesh">
77+
<bins>2</bins>
78+
</filter>
79+
<tally id="1" name="regular mesh tally">
80+
<filters>1</filters>
81+
<scores>flux</scores>
82+
<estimator>collision</estimator>
83+
</tally>
84+
<tally id="2" name="xdg mesh tally">
85+
<filters>2</filters>
86+
<scores>flux</scores>
87+
<estimator>collision</estimator>
88+
</tally>
89+
</tallies>
90+
</model>

‎tests/regression_tests/xdg/tallies/inputs_true0_moab.dat‎

Lines changed: 90 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,90 @@
1+
<?xml version='1.0' encoding='utf-8'?>
2+
<model>
3+
<materials>
4+
<material id="1" name="fuel" depletable="true">
5+
<density value="4.5" units="g/cc"/>
6+
<nuclide name="U235" ao="1.0"/>
7+
</material>
8+
<material id="2" name="zircaloy">
9+
<density value="5.77" units="g/cc"/>
10+
<nuclide name="Zr90" ao="0.5145"/>
11+
<nuclide name="Zr91" ao="0.1122"/>
12+
<nuclide name="Zr92" ao="0.1715"/>
13+
<nuclide name="Zr94" ao="0.1738"/>
14+
<nuclide name="Zr96" ao="0.028"/>
15+
</material>
16+
<material id="3" name="water">
17+
<density value="0.07416" units="atom/b-cm"/>
18+
<nuclide name="H1" ao="2.0"/>
19+
<nuclide name="O16" ao="1.0"/>
20+
</material>
21+
</materials>
22+
<geometry>
23+
<cell id="1" name="fuel" material="1" region="1 -2 3 -4 5 -6" universe="1"/>
24+
<cell id="2" name="clad" material="2" region="(-1 | 2 | -3 | 4 | -5 | 6) (7 -8 9 -10 11 -12)" universe="1"/>
25+
<cell id="3" name="water" material="3" region="(-7 | 8 | -9 | 10 | -11 | 12) (13 -14 15 -16 17 -18)" universe="1"/>
26+
<surface id="1" name="minimum x" type="x-plane" coeffs="-5.0"/>
27+
<surface id="2" name="maximum x" type="x-plane" coeffs="5.0"/>
28+
<surface id="3" name="minimum y" type="y-plane" coeffs="-5.0"/>
29+
<surface id="4" name="maximum y" type="y-plane" coeffs="5.0"/>
30+
<surface id="5" name="minimum z" type="z-plane" coeffs="-5.0"/>
31+
<surface id="6" name="maximum z" type="z-plane" coeffs="5.0"/>
32+
<surface id="7" name="minimum x" type="x-plane" coeffs="-6.0"/>
33+
<surface id="8" name="maximum x" type="x-plane" coeffs="6.0"/>
34+
<surface id="9" name="minimum y" type="y-plane" coeffs="-6.0"/>
35+
<surface id="10" name="maximum y" type="y-plane" coeffs="6.0"/>
36+
<surface id="11" name="minimum z" type="z-plane" coeffs="-6.0"/>
37+
<surface id="12" name="maximum z" type="z-plane" coeffs="6.0"/>
38+
<surface id="13" name="minimum x" type="x-plane" boundary="vacuum" coeffs="-10.0"/>
39+
<surface id="14" name="maximum x" type="x-plane" boundary="vacuum" coeffs="10.0"/>
40+
<surface id="15" name="minimum y" type="y-plane" boundary="vacuum" coeffs="-10.0"/>
41+
<surface id="16" name="maximum y" type="y-plane" boundary="vacuum" coeffs="10.0"/>
42+
<surface id="17" name="minimum z" type="z-plane" boundary="vacuum" coeffs="-10.0"/>
43+
<surface id="18" name="maximum z" type="z-plane" boundary="vacuum" coeffs="10.0"/>
44+
</geometry>
45+
<settings>
46+
<run_mode>fixed source</run_mode>
47+
<particles>1000</particles>
48+
<batches>10</batches>
49+
<source type="independent" strength="1.0" particle="neutron">
50+
<space type="spherical" origin="0.0 0.0 0.0">
51+
<r type="uniform" parameters="0.0 0.0"/>
52+
<cos_theta type="discrete">
53+
<parameters>1.0 1.0</parameters>
54+
</cos_theta>
55+
<phi type="discrete">
56+
<parameters>0.0 1.0</parameters>
57+
</phi>
58+
</space>
59+
<energy type="discrete">
60+
<parameters>15000000.0 1.0</parameters>
61+
</energy>
62+
</source>
63+
</settings>
64+
<tallies>
65+
<mesh id="1">
66+
<dimension>10 10 10</dimension>
67+
<lower_left>-10.0 -10.0 -10.0</lower_left>
68+
<upper_right>10.0 10.0 10.0</upper_right>
69+
</mesh>
70+
<mesh id="2" type="xdg" library="moab">
71+
<filename>test_mesh_tets.e</filename>
72+
</mesh>
73+
<filter id="1" type="mesh">
74+
<bins>1</bins>
75+
</filter>
76+
<filter id="2" type="mesh">
77+
<bins>2</bins>
78+
</filter>
79+
<tally id="1" name="regular mesh tally">
80+
<filters>1</filters>
81+
<scores>flux</scores>
82+
<estimator>collision</estimator>
83+
</tally>
84+
<tally id="2" name="xdg mesh tally">
85+
<filters>2</filters>
86+
<scores>flux</scores>
87+
<estimator>collision</estimator>
88+
</tally>
89+
</tallies>
90+
</model>

‎tests/regression_tests/xdg/tallies/inputs_true1_libmesh.dat‎

Lines changed: 90 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,90 @@
1+
<?xml version='1.0' encoding='utf-8'?>
2+
<model>
3+
<materials>
4+
<material id="1" name="fuel" depletable="true">
5+
<density value="4.5" units="g/cc"/>
6+
<nuclide name="U235" ao="1.0"/>
7+
</material>
8+
<material id="2" name="zircaloy">
9+
<density value="5.77" units="g/cc"/>
10+
<nuclide name="Zr90" ao="0.5145"/>
11+
<nuclide name="Zr91" ao="0.1122"/>
12+
<nuclide name="Zr92" ao="0.1715"/>
13+
<nuclide name="Zr94" ao="0.1738"/>
14+
<nuclide name="Zr96" ao="0.028"/>
15+
</material>
16+
<material id="3" name="water">
17+
<density value="0.07416" units="atom/b-cm"/>
18+
<nuclide name="H1" ao="2.0"/>
19+
<nuclide name="O16" ao="1.0"/>
20+
</material>
21+
</materials>
22+
<geometry>
23+
<cell id="1" name="fuel" material="1" region="1 -2 3 -4 5 -6" universe="1"/>
24+
<cell id="2" name="clad" material="2" region="(-1 | 2 | -3 | 4 | -5 | 6) (7 -8 9 -10 11 -12)" universe="1"/>
25+
<cell id="3" name="water" material="3" region="(-7 | 8 | -9 | 10 | -11 | 12) (13 -14 15 -16 17 -18)" universe="1"/>
26+
<surface id="1" name="minimum x" type="x-plane" coeffs="-5.0"/>
27+
<surface id="2" name="maximum x" type="x-plane" coeffs="5.0"/>
28+
<surface id="3" name="minimum y" type="y-plane" coeffs="-5.0"/>
29+
<surface id="4" name="maximum y" type="y-plane" coeffs="5.0"/>
30+
<surface id="5" name="minimum z" type="z-plane" coeffs="-5.0"/>
31+
<surface id="6" name="maximum z" type="z-plane" coeffs="5.0"/>
32+
<surface id="7" name="minimum x" type="x-plane" coeffs="-6.0"/>
33+
<surface id="8" name="maximum x" type="x-plane" coeffs="6.0"/>
34+
<surface id="9" name="minimum y" type="y-plane" coeffs="-6.0"/>
35+
<surface id="10" name="maximum y" type="y-plane" coeffs="6.0"/>
36+
<surface id="11" name="minimum z" type="z-plane" coeffs="-6.0"/>
37+
<surface id="12" name="maximum z" type="z-plane" coeffs="6.0"/>
38+
<surface id="13" name="minimum x" type="x-plane" boundary="vacuum" coeffs="-10.0"/>
39+
<surface id="14" name="maximum x" type="x-plane" boundary="vacuum" coeffs="10.0"/>
40+
<surface id="15" name="minimum y" type="y-plane" boundary="vacuum" coeffs="-10.0"/>
41+
<surface id="16" name="maximum y" type="y-plane" boundary="vacuum" coeffs="10.0"/>
42+
<surface id="17" name="minimum z" type="z-plane" boundary="vacuum" coeffs="-10.0"/>
43+
<surface id="18" name="maximum z" type="z-plane" boundary="vacuum" coeffs="10.0"/>
44+
</geometry>
45+
<settings>
46+
<run_mode>fixed source</run_mode>
47+
<particles>1000</particles>
48+
<batches>10</batches>
49+
<source type="independent" strength="1.0" particle="neutron">
50+
<space type="spherical" origin="0.0 0.0 0.0">
51+
<r type="uniform" parameters="0.0 0.0"/>
52+
<cos_theta type="discrete">
53+
<parameters>1.0 1.0</parameters>
54+
</cos_theta>
55+
<phi type="discrete">
56+
<parameters>0.0 1.0</parameters>
57+
</phi>
58+
</space>
59+
<energy type="discrete">
60+
<parameters>15000000.0 1.0</parameters>
61+
</energy>
62+
</source>
63+
</settings>
64+
<tallies>
65+
<mesh id="1">
66+
<dimension>10 10 10</dimension>
67+
<lower_left>-10.0 -10.0 -10.0</lower_left>
68+
<upper_right>10.0 10.0 10.0</upper_right>
69+
</mesh>
70+
<mesh id="2" type="xdg" library="libmesh">
71+
<filename>test_mesh_tets.exo</filename>
72+
</mesh>
73+
<filter id="1" type="mesh">
74+
<bins>1</bins>
75+
</filter>
76+
<filter id="2" type="mesh">
77+
<bins>2</bins>
78+
</filter>
79+
<tally id="1" name="regular mesh tally">
80+
<filters>1</filters>
81+
<scores>flux</scores>
82+
<estimator>collision</estimator>
83+
</tally>
84+
<tally id="2" name="xdg mesh tally">
85+
<filters>2</filters>
86+
<scores>flux</scores>
87+
<estimator>collision</estimator>
88+
</tally>
89+
</tallies>
90+
</model>

‎tests/regression_tests/xdg/tallies/inputs_true2_libmesh.dat‎

Lines changed: 90 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,90 @@
1+
<?xml version='1.0' encoding='utf-8'?>
2+
<model>
3+
<materials>
4+
<material id="1" name="fuel" depletable="true">
5+
<density value="4.5" units="g/cc"/>
6+
<nuclide name="U235" ao="1.0"/>
7+
</material>
8+
<material id="2" name="zircaloy">
9+
<density value="5.77" units="g/cc"/>
10+
<nuclide name="Zr90" ao="0.5145"/>
11+
<nuclide name="Zr91" ao="0.1122"/>
12+
<nuclide name="Zr92" ao="0.1715"/>
13+
<nuclide name="Zr94" ao="0.1738"/>
14+
<nuclide name="Zr96" ao="0.028"/>
15+
</material>
16+
<material id="3" name="water">
17+
<density value="0.07416" units="atom/b-cm"/>
18+
<nuclide name="H1" ao="2.0"/>
19+
<nuclide name="O16" ao="1.0"/>
20+
</material>
21+
</materials>
22+
<geometry>
23+
<cell id="1" name="fuel" material="1" region="1 -2 3 -4 5 -6" universe="1"/>
24+
<cell id="2" name="clad" material="2" region="(-1 | 2 | -3 | 4 | -5 | 6) (7 -8 9 -10 11 -12)" universe="1"/>
25+
<cell id="3" name="water" material="3" region="(-7 | 8 | -9 | 10 | -11 | 12) (13 -14 15 -16 17 -18)" universe="1"/>
26+
<surface id="1" name="minimum x" type="x-plane" coeffs="-5.0"/>
27+
<surface id="2" name="maximum x" type="x-plane" coeffs="5.0"/>
28+
<surface id="3" name="minimum y" type="y-plane" coeffs="-5.0"/>
29+
<surface id="4" name="maximum y" type="y-plane" coeffs="5.0"/>
30+
<surface id="5" name="minimum z" type="z-plane" coeffs="-5.0"/>
31+
<surface id="6" name="maximum z" type="z-plane" coeffs="5.0"/>
32+
<surface id="7" name="minimum x" type="x-plane" coeffs="-6.0"/>
33+
<surface id="8" name="maximum x" type="x-plane" coeffs="6.0"/>
34+
<surface id="9" name="minimum y" type="y-plane" coeffs="-6.0"/>
35+
<surface id="10" name="maximum y" type="y-plane" coeffs="6.0"/>
36+
<surface id="11" name="minimum z" type="z-plane" coeffs="-6.0"/>
37+
<surface id="12" name="maximum z" type="z-plane" coeffs="6.0"/>
38+
<surface id="13" name="minimum x" type="x-plane" boundary="vacuum" coeffs="-10.0"/>
39+
<surface id="14" name="maximum x" type="x-plane" boundary="vacuum" coeffs="10.0"/>
40+
<surface id="15" name="minimum y" type="y-plane" boundary="vacuum" coeffs="-10.0"/>
41+
<surface id="16" name="maximum y" type="y-plane" boundary="vacuum" coeffs="10.0"/>
42+
<surface id="17" name="minimum z" type="z-plane" boundary="vacuum" coeffs="-10.0"/>
43+
<surface id="18" name="maximum z" type="z-plane" boundary="vacuum" coeffs="10.0"/>
44+
</geometry>
45+
<settings>
46+
<run_mode>fixed source</run_mode>
47+
<particles>1000</particles>
48+
<batches>10</batches>
49+
<source type="independent" strength="1.0" particle="neutron">
50+
<space type="spherical" origin="0.0 0.0 0.0">
51+
<r type="uniform" parameters="0.0 0.0"/>
52+
<cos_theta type="discrete">
53+
<parameters>1.0 1.0</parameters>
54+
</cos_theta>
55+
<phi type="discrete">
56+
<parameters>0.0 1.0</parameters>
57+
</phi>
58+
</space>
59+
<energy type="discrete">
60+
<parameters>15000000.0 1.0</parameters>
61+
</energy>
62+
</source>
63+
</settings>
64+
<tallies>
65+
<mesh id="1">
66+
<dimension>10 10 10</dimension>
67+
<lower_left>-10.0 -10.0 -10.0</lower_left>
68+
<upper_right>10.0 10.0 10.0</upper_right>
69+
</mesh>
70+
<mesh id="2" type="xdg" library="libmesh">
71+
<filename>test_mesh_tets_w_holes.e</filename>
72+
</mesh>
73+
<filter id="1" type="mesh">
74+
<bins>1</bins>
75+
</filter>
76+
<filter id="2" type="mesh">
77+
<bins>2</bins>
78+
</filter>
79+
<tally id="1" name="regular mesh tally">
80+
<filters>1</filters>
81+
<scores>flux</scores>
82+
<estimator>collision</estimator>
83+
</tally>
84+
<tally id="2" name="xdg mesh tally">
85+
<filters>2</filters>
86+
<scores>flux</scores>
87+
<estimator>collision</estimator>
88+
</tally>
89+
</tallies>
90+
</model>

0 commit comments

Comments
 (0)