crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-26 21:25:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Add regression tests for cmfd restart, ref_d, rolling_window, expanding_window

Browse source
This commit is contained in:
Shikhar Kumar 2019-04-17 18:41:45 -04:00
parent 61616eedde
commit 5f1f37bd76
28 changed files with 2514 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

0
tests/regression_tests/cmfd_feed_expanding_window/__init__.py Normal file
View file

43
tests/regression_tests/cmfd_feed_expanding_window/geometry.xml Normal file
View file

@ -0,0 +1,43 @@
<?xml version="1.0" encoding="UTF-8"?>
<geometry>
<!-- Definition of Cells -->
<cell id="1">
<universe>0</universe>
<region>-1 2 -3 4 -5 6</region>
<material>1</material>
</cell>
<!-- Defition of Surfaces -->
<surface id="1">
<type>x-plane</type>
<coeffs>10</coeffs>
<boundary> vacuum </boundary>
</surface>
<surface id="2">
<type>x-plane</type>
<coeffs>-10</coeffs>
<boundary> vacuum </boundary>
</surface>
<surface id="3">
<type>y-plane</type>
<coeffs>1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="4">
<type>y-plane</type>
<coeffs>-1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="5">
<type>z-plane</type>
<coeffs>1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="6">
<type>z-plane</type>
<coeffs>-1</coeffs>
<boundary>reflective</boundary>
</surface>
</geometry>

12
tests/regression_tests/cmfd_feed_expanding_window/materials.xml Normal file
View file

@ -0,0 +1,12 @@
<?xml version="1.0"?>
<materials>
<!-- Definition of materials -->
<material id="1">
<density value="19" units="g/cc" />
<nuclide name="U235" wo="0.21" />
<nuclide name="U238" wo="0.68" />
<nuclide name="O16" wo="0.11" />
</material>
</materials>

488
tests/regression_tests/cmfd_feed_expanding_window/results_true.dat Normal file
View file

@ -0,0 +1,488 @@
k-combined:
1.165403E+00 1.129918E-02
tally 1:
1.141961E+01
1.307497E+01
2.075941E+01
4.316095E+01
2.834539E+01
8.063510E+01
3.513869E+01
1.238515E+02
3.699858E+01
1.372516E+02
3.612668E+01
1.310763E+02
3.373368E+01
1.140458E+02
2.836844E+01
8.075759E+01
2.158837E+01
4.674848E+01
1.141868E+01
1.310648E+01
tally 2:
1.126626E+00
1.269287E+00
7.749042E-01
6.004765E-01
1.991608E+00
3.966503E+00
1.411775E+00
1.993108E+00
2.978264E+00
8.870057E+00
2.130066E+00
4.537180E+00
3.708857E+00
1.375562E+01
2.698010E+00
7.279259E+00
4.233192E+00
1.791991E+01
3.023990E+00
9.144518E+00
4.062113E+00
1.650076E+01
2.910472E+00
8.470848E+00
3.506681E+00
1.229681E+01
2.497751E+00
6.238760E+00
2.974675E+00
8.848690E+00
2.121163E+00
4.499332E+00
2.329248E+00
5.425398E+00
1.658233E+00
2.749736E+00
1.158552E+00
1.342242E+00
8.282698E-01
6.860309E-01
tally 3:
7.459301E-01
5.564117E-01
4.877614E-02
2.379112E-03
1.356702E+00
1.840641E+00
8.361624E-02
6.991675E-03
2.047332E+00
4.191570E+00
1.521351E-01
2.314509E-02
2.610287E+00
6.813596E+00
1.718778E-01
2.954199E-02
2.911379E+00
8.476127E+00
1.823299E-01
3.324418E-02
2.806920E+00
7.878801E+00
2.090406E-01
4.369797E-02
2.415727E+00
5.835736E+00
1.486511E-01
2.209715E-02
2.051064E+00
4.206862E+00
1.196177E-01
1.430839E-02
1.593746E+00
2.540026E+00
1.126497E-01
1.268994E-02
7.994935E-01
6.391898E-01
5.806683E-02
3.371757E-03
tally 4:
1.341719E-01
1.800210E-02
0.000000E+00
0.000000E+00
1.420707E-01
2.018410E-02
2.633150E-01
6.933479E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.633150E-01
6.933479E-02
1.420707E-01
2.018410E-02
2.628613E-01
6.909608E-02
3.590382E-01
1.289084E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.590382E-01
1.289084E-01
2.628613E-01
6.909608E-02
3.851599E-01
1.483482E-01
4.595809E-01
2.112146E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.595809E-01
2.112146E-01
3.851599E-01
1.483482E-01
4.679967E-01
2.190209E-01
4.970187E-01
2.470276E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.970187E-01
2.470276E-01
4.679967E-01
2.190209E-01
4.909661E-01
2.410477E-01
4.838785E-01
2.341384E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.838785E-01
2.341384E-01
4.909661E-01
2.410477E-01
5.152690E-01
2.655021E-01
4.788649E-01
2.293116E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.788649E-01
2.293116E-01
5.152690E-01
2.655021E-01
4.266876E-01
1.820623E-01
3.401342E-01
1.156913E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.401342E-01
1.156913E-01
4.266876E-01
1.820623E-01
3.724186E-01
1.386956E-01
2.560013E-01
6.553669E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.560013E-01
6.553669E-02
3.724186E-01
1.386956E-01
2.892134E-01
8.364439E-02
1.556176E-01
2.421685E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.556176E-01
2.421685E-02
2.892134E-01
8.364439E-02
1.497744E-01
2.243237E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
tally 5:
7.459301E-01
5.564117E-01
1.149820E-01
1.322085E-02
1.356702E+00
1.840641E+00
1.801102E-01
3.243968E-02
2.047332E+00
4.191570E+00
2.627497E-01
6.903740E-02
2.609217E+00
6.808012E+00
3.360314E-01
1.129171E-01
2.909457E+00
8.464941E+00
4.095451E-01
1.677272E-01
2.805131E+00
7.868760E+00
3.957627E-01
1.566281E-01
2.415727E+00
5.835736E+00
3.196414E-01
1.021706E-01
2.049859E+00
4.201923E+00
3.162399E-01
1.000077E-01
1.593746E+00
2.540026E+00
2.169985E-01
4.708836E-02
7.994935E-01
6.391898E-01
1.002481E-01
1.004969E-02
cmfd indices
1.000000E+01
1.000000E+00
1.000000E+00
1.000000E+00
k cmfd
1.143597E+00
1.163387E+00
1.173384E+00
1.171035E+00
1.147196E+00
1.122260E+00
1.106380E+00
1.124693E+00
1.133192E+00
1.134435E+00
1.142380E+00
1.132168E+00
1.132560E+00
1.153781E+00
1.170308E+00
1.184540E+00
cmfd entropy
3.212002E+00
3.206393E+00
3.223984E+00
3.222764E+00
3.232123E+00
3.242083E+00
3.246067E+00
3.238869E+00
3.235585E+00
3.234611E+00
3.233575E+00
3.236922E+00
3.229545E+00
3.225232E+00
3.221485E+00
3.219108E+00
cmfd balance
8.26180E-03
4.27338E-03
2.22686E-03
1.93026E-03
1.96979E-03
2.13756E-03
2.01479E-03
1.74519E-03
2.09248E-03
1.25545E-03
1.48370E-03
1.75963E-03
1.98194E-03
1.87306E-03
1.24780E-03
1.15560E-03
cmfd dominance ratio
5.404E-01
5.406E-01
5.449E-01
5.473E-01
5.534E-01
5.623E-01
5.738E-01
5.611E-01
5.569E-01
5.556E-01
5.555E-01
5.583E-01
5.544E-01
5.486E-01
5.412E-01
5.383E-01
cmfd openmc source comparison
1.575499E-02
1.293688E-02
3.531746E-03
8.281178E-03
5.771681E-03
7.459013E-03
5.012869E-03
1.770224E-03
5.242540E-03
3.888027E-03
6.653433E-03
8.839928E-03
5.456904E-03
5.668412E-03
4.016377E-03
4.179381E-03
cmfd source
4.116210E-02
7.797480E-02
1.062822E-01
1.333719E-01
1.481091E-01
1.370422E-01
1.299765E-01
9.887040E-02
8.228753E-02
4.492330E-02

26
tests/regression_tests/cmfd_feed_expanding_window/settings.xml Normal file
View file

@ -0,0 +1,26 @@
<?xml version="1.0" encoding="UTF-8"?>
<settings>
<!-- Parameters for criticality calculation -->
<run_mode>eigenvalue</run_mode>
<batches>20</batches>
<inactive>10</inactive>
<particles>1000</particles>
<!-- Starting source -->
<source>
<space>
<type>box</type>
<parameters>-10 -1 -1 10 1 1</parameters>
</space>
</source>
<!-- Shannon Entropy -->
<mesh id="10">
<dimension> 10 1 1 </dimension>
<lower_left> -10.0 -1.0 -1.0 </lower_left>
<upper_right> 10.0 1.0 1.0 </upper_right>
</mesh>
<entropy_mesh>10</entropy_mesh>
</settings>

21
tests/regression_tests/cmfd_feed_expanding_window/tallies.xml Normal file
View file

@ -0,0 +1,21 @@
<?xml version="1.0"?>
<tallies>
<mesh id="1">
<type>regular</type>
<lower_left>-10 -1 -1 </lower_left>
<upper_right>10 1 1</upper_right>
<dimension>10 1 1</dimension>
</mesh>
<filter id="1">
<type>mesh</type>
<bins>1</bins>
</filter>
<tally id="1">
<filters>1</filters>
<scores>flux</scores>
</tally>
</tallies>

26
tests/regression_tests/cmfd_feed_expanding_window/test.py Normal file
View file

@ -0,0 +1,26 @@
from tests.testing_harness import CMFDTestHarness
from openmc import cmfd
def test_cmfd_feed_rolling_window():
"""Test 1 group CMFD solver with CMFD feedback"""
# Initialize and set CMFD mesh
cmfd_mesh = cmfd.CMFDMesh()
cmfd_mesh.lower_left = -10.0, -1.0, -1.0
cmfd_mesh.upper_right = 10.0, 1.0, 1.0
cmfd_mesh.dimension = 10, 1, 1
cmfd_mesh.albedo = 0.0, 0.0, 1.0, 1.0, 1.0, 1.0
# Initialize and run CMFDRun object
cmfd_run = cmfd.CMFDRun()
cmfd_run.mesh = cmfd_mesh
cmfd_run.tally_begin = 5
cmfd_run.feedback_begin = 10
cmfd_run.feedback = True
cmfd_run.gauss_seidel_tolerance = [1.e-15, 1.e-20]
cmfd_run.window_type = 'expanding'
cmfd_run.run()
# Initialize and run CMFD test harness
harness = CMFDTestHarness('statepoint.20.h5', cmfd_run)
harness.main()

0
tests/regression_tests/cmfd_feed_ref_d/__init__.py Normal file
View file

43
tests/regression_tests/cmfd_feed_ref_d/geometry.xml Normal file
View file

@ -0,0 +1,43 @@
<?xml version="1.0" encoding="UTF-8"?>
<geometry>
<!-- Definition of Cells -->
<cell id="1">
<universe>0</universe>
<region>-1 2 -3 4 -5 6</region>
<material>1</material>
</cell>
<!-- Defition of Surfaces -->
<surface id="1">
<type>x-plane</type>
<coeffs>10</coeffs>
<boundary> vacuum </boundary>
</surface>
<surface id="2">
<type>x-plane</type>
<coeffs>-10</coeffs>
<boundary> vacuum </boundary>
</surface>
<surface id="3">
<type>y-plane</type>
<coeffs>1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="4">
<type>y-plane</type>
<coeffs>-1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="5">
<type>z-plane</type>
<coeffs>1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="6">
<type>z-plane</type>
<coeffs>-1</coeffs>
<boundary>reflective</boundary>
</surface>
</geometry>

12
tests/regression_tests/cmfd_feed_ref_d/materials.xml Normal file
View file

@ -0,0 +1,12 @@
<?xml version="1.0"?>
<materials>
<!-- Definition of materials -->
<material id="1">
<density value="19" units="g/cc" />
<nuclide name="U235" wo="0.21" />
<nuclide name="U238" wo="0.68" />
<nuclide name="O16" wo="0.11" />
</material>
</materials>

488
tests/regression_tests/cmfd_feed_ref_d/results_true.dat Normal file
View file

@ -0,0 +1,488 @@
k-combined:
1.165408E+00 1.127320E-02
tally 1:
1.141009E+01
1.305436E+01
2.074878E+01
4.311543E+01
2.834824E+01
8.065180E+01
3.513901E+01
1.238542E+02
3.699846E+01
1.372507E+02
3.612638E+01
1.310741E+02
3.373329E+01
1.140431E+02
2.836806E+01
8.075546E+01
2.158796E+01
4.674672E+01
1.141841E+01
1.310586E+01
tally 2:
1.126662E+00
1.269368E+00
7.749275E-01
6.005127E-01
1.991651E+00
3.966673E+00
1.411802E+00
1.993185E+00
2.978294E+00
8.870233E+00
2.130084E+00
4.537258E+00
3.708845E+00
1.375553E+01
2.698001E+00
7.279210E+00
4.233143E+00
1.791950E+01
3.023958E+00
9.144324E+00
4.062060E+00
1.650033E+01
2.910435E+00
8.470632E+00
3.506628E+00
1.229644E+01
2.497713E+00
6.238573E+00
2.974637E+00
8.848463E+00
2.121135E+00
4.499215E+00
2.329228E+00
5.425303E+00
1.658218E+00
2.749687E+00
1.158542E+00
1.342220E+00
8.282626E-01
6.860190E-01
tally 3:
7.459525E-01
5.564451E-01
4.877161E-02
2.378670E-03
1.356729E+00
1.840713E+00
8.360848E-02
6.990377E-03
2.047350E+00
4.191641E+00
1.521210E-01
2.314079E-02
2.610278E+00
6.813552E+00
1.718619E-01
2.953650E-02
2.911348E+00
8.475950E+00
1.823129E-01
3.323800E-02
2.806884E+00
7.878600E+00
2.090212E-01
4.368986E-02
2.415690E+00
5.835559E+00
1.486373E-01
2.209304E-02
2.051037E+00
4.206752E+00
1.196066E-01
1.430573E-02
1.593732E+00
2.539980E+00
1.126392E-01
1.268759E-02
7.994865E-01
6.391786E-01
5.806144E-02
3.371131E-03
tally 4:
1.341763E-01
1.800329E-02
0.000000E+00
0.000000E+00
1.420749E-01
2.018527E-02
2.633206E-01
6.933776E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.633206E-01
6.933776E-02
1.420749E-01
2.018527E-02
2.628684E-01
6.909981E-02
3.590428E-01
1.289117E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.590428E-01
1.289117E-01
2.628684E-01
6.909981E-02
3.851638E-01
1.483512E-01
4.595776E-01
2.112116E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.595776E-01
2.112116E-01
3.851638E-01
1.483512E-01
4.679953E-01
2.190196E-01
4.970127E-01
2.470216E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.970127E-01
2.470216E-01
4.679953E-01
2.190196E-01
4.909606E-01
2.410423E-01
4.838733E-01
2.341333E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.838733E-01
2.341333E-01
4.909606E-01
2.410423E-01
5.152620E-01
2.654949E-01
4.788580E-01
2.293049E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.788580E-01
2.293049E-01
5.152620E-01
2.654949E-01
4.266812E-01
1.820568E-01
3.401298E-01
1.156883E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.401298E-01
1.156883E-01
4.266812E-01
1.820568E-01
3.724142E-01
1.386923E-01
2.559996E-01
6.553579E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.559996E-01
6.553579E-02
3.724142E-01
1.386923E-01
2.892105E-01
8.364271E-02
1.556166E-01
2.421652E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.556166E-01
2.421652E-02
2.892105E-01
8.364271E-02
1.497731E-01
2.243198E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
tally 5:
7.459525E-01
5.564451E-01
1.149859E-01
1.322177E-02
1.356729E+00
1.840713E+00
1.801152E-01
3.244149E-02
2.047350E+00
4.191641E+00
2.627538E-01
6.903954E-02
2.609208E+00
6.807967E+00
3.360312E-01
1.129170E-01
2.909427E+00
8.464765E+00
4.095398E-01
1.677229E-01
2.805095E+00
7.868559E+00
3.957580E-01
1.566244E-01
2.415690E+00
5.835559E+00
3.196366E-01
1.021676E-01
2.049833E+00
4.201813E+00
3.162366E-01
1.000056E-01
1.593732E+00
2.539980E+00
2.169968E-01
4.708761E-02
7.994865E-01
6.391786E-01
1.002479E-01
1.004964E-02
cmfd indices
1.000000E+01
1.000000E+00
1.000000E+00
1.000000E+00
k cmfd
1.143785E+00
1.163460E+00
1.173453E+00
1.171056E+00
1.147214E+00
1.122230E+00
1.106385E+00
1.124706E+00
1.133207E+00
1.134453E+00
1.142355E+00
1.132121E+00
1.132479E+00
1.153657E+00
1.170183E+00
1.184408E+00
cmfd entropy
3.211758E+00
3.206356E+00
3.223933E+00
3.222754E+00
3.232110E+00
3.242098E+00
3.246062E+00
3.238858E+00
3.235577E+00
3.234604E+00
3.233582E+00
3.236922E+00
3.229563E+00
3.225265E+00
3.221498E+00
3.219126E+00
cmfd balance
8.26180E-03
4.27338E-03
2.22686E-03
1.93026E-03
1.96979E-03
2.13756E-03
2.01521E-03
1.74538E-03
2.09240E-03
1.25495E-03
1.49542E-03
1.76968E-03
1.99174E-03
1.87753E-03
1.24170E-03
1.15645E-03
cmfd dominance ratio
5.408E-01
5.374E-01
5.420E-01
5.444E-01
5.513E-01
5.613E-01
5.722E-01
5.592E-01
5.542E-01
5.537E-01
5.541E-01
5.579E-01
5.543E-01
5.487E-01
5.413E-01
5.381E-01
cmfd openmc source comparison
1.597982E-02
1.276493E-02
3.495229E-03
8.157807E-03
5.715330E-03
7.433111E-03
5.006211E-03
1.766072E-03
5.184425E-03
3.864321E-03
6.617152E-03
8.841210E-03
5.454747E-03
5.652690E-03
4.006767E-03
4.167617E-03
cmfd source
4.116746E-02
7.798354E-02
1.062881E-01
1.333681E-01
1.481008E-01
1.370408E-01
1.299726E-01
9.886814E-02
8.228698E-02
4.492338E-02

26
tests/regression_tests/cmfd_feed_ref_d/settings.xml Normal file
View file

@ -0,0 +1,26 @@
<?xml version="1.0" encoding="UTF-8"?>
<settings>
<!-- Parameters for criticality calculation -->
<run_mode>eigenvalue</run_mode>
<batches>20</batches>
<inactive>10</inactive>
<particles>1000</particles>
<!-- Starting source -->
<source>
<space>
<type>box</type>
<parameters>-10 -1 -1 10 1 1</parameters>
</space>
</source>
<!-- Shannon Entropy -->
<mesh id="10">
<dimension> 10 1 1 </dimension>
<lower_left> -10.0 -1.0 -1.0 </lower_left>
<upper_right> 10.0 1.0 1.0 </upper_right>
</mesh>
<entropy_mesh>10</entropy_mesh>
</settings>

21
tests/regression_tests/cmfd_feed_ref_d/tallies.xml Normal file
View file

@ -0,0 +1,21 @@
<?xml version="1.0"?>
<tallies>
<mesh id="1">
<type>regular</type>
<lower_left>-10 -1 -1 </lower_left>
<upper_right>10 1 1</upper_right>
<dimension>10 1 1</dimension>
</mesh>
<filter id="1">
<type>mesh</type>
<bins>1</bins>
</filter>
<tally id="1">
<filters>1</filters>
<scores>flux</scores>
</tally>
</tallies>

27
tests/regression_tests/cmfd_feed_ref_d/test.py Normal file
View file

@ -0,0 +1,27 @@
from tests.testing_harness import CMFDTestHarness
from openmc import cmfd
def test_cmfd_feed_rolling_window():
"""Test 1 group CMFD solver with CMFD feedback"""
# Initialize and set CMFD mesh
cmfd_mesh = cmfd.CMFDMesh()
cmfd_mesh.lower_left = -10.0, -1.0, -1.0
cmfd_mesh.upper_right = 10.0, 1.0, 1.0
cmfd_mesh.dimension = 10, 1, 1
cmfd_mesh.albedo = 0.0, 0.0, 1.0, 1.0, 1.0, 1.0
# Initialize and run CMFDRun object
cmfd_run = cmfd.CMFDRun()
cmfd_run.mesh = cmfd_mesh
cmfd_run.tally_begin = 5
cmfd_run.feedback_begin = 10
cmfd_run.feedback = True
cmfd_run.gauss_seidel_tolerance = [1.e-15, 1.e-20]
cmfd_run.window_type = 'expanding'
cmfd_run.ref_d = [0.542]
cmfd_run.run()
# Initialize and run CMFD test harness
harness = CMFDTestHarness('statepoint.20.h5', cmfd_run)
harness.main()

0
tests/regression_tests/cmfd_feed_rolling_window/__init__.py Normal file
View file

43
tests/regression_tests/cmfd_feed_rolling_window/geometry.xml Normal file
View file

@ -0,0 +1,43 @@
<?xml version="1.0" encoding="UTF-8"?>
<geometry>
<!-- Definition of Cells -->
<cell id="1">
<universe>0</universe>
<region>-1 2 -3 4 -5 6</region>
<material>1</material>
</cell>
<!-- Defition of Surfaces -->
<surface id="1">
<type>x-plane</type>
<coeffs>10</coeffs>
<boundary> vacuum </boundary>
</surface>
<surface id="2">
<type>x-plane</type>
<coeffs>-10</coeffs>
<boundary> vacuum </boundary>
</surface>
<surface id="3">
<type>y-plane</type>
<coeffs>1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="4">
<type>y-plane</type>
<coeffs>-1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="5">
<type>z-plane</type>
<coeffs>1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="6">
<type>z-plane</type>
<coeffs>-1</coeffs>
<boundary>reflective</boundary>
</surface>
</geometry>

12
tests/regression_tests/cmfd_feed_rolling_window/materials.xml Normal file
View file

@ -0,0 +1,12 @@
<?xml version="1.0"?>
<materials>
<!-- Definition of materials -->
<material id="1">
<density value="19" units="g/cc" />
<nuclide name="U235" wo="0.21" />
<nuclide name="U238" wo="0.68" />
<nuclide name="O16" wo="0.11" />
</material>
</materials>

488
tests/regression_tests/cmfd_feed_rolling_window/results_true.dat Normal file
View file

@ -0,0 +1,488 @@
k-combined:
1.172120E+00 9.761693E-03
tally 1:
1.131066E+01
1.287315E+01
2.035268E+01
4.162146E+01
2.839961E+01
8.098764E+01
3.405187E+01
1.164951E+02
3.712390E+01
1.380899E+02
3.687780E+01
1.363090E+02
3.380903E+01
1.147695E+02
2.903657E+01
8.489183E+01
2.121421E+01
4.519765E+01
1.112628E+01
1.245629E+01
tally 2:
1.114845E+00
1.242879E+00
8.000235E-01
6.400375E-01
1.852141E+00
3.430427E+00
1.324803E+00
1.755104E+00
2.585263E+00
6.683587E+00
1.840969E+00
3.389167E+00
3.737263E+00
1.396714E+01
2.674345E+00
7.152122E+00
3.989056E+00
1.591257E+01
2.843382E+00
8.084820E+00
3.951710E+00
1.561601E+01
2.838216E+00
8.055468E+00
3.734018E+00
1.394289E+01
2.676071E+00
7.161356E+00
3.241209E+00
1.050543E+01
2.315684E+00
5.362392E+00
2.459493E+00
6.049108E+00
1.734344E+00
3.007949E+00
1.079306E+00
1.164902E+00
7.536940E-01
5.680547E-01
tally 3:
7.728275E-01
5.972624E-01
5.550061E-02
3.080318E-03
1.272034E+00
1.618070E+00
7.284456E-02
5.306329E-03
1.779821E+00
3.167762E+00
1.110012E-01
1.232127E-02
2.586020E+00
6.687502E+00
1.618768E-01
2.620410E-02
2.726552E+00
7.434084E+00
1.965647E-01
3.863767E-02
2.741747E+00
7.517178E+00
1.907834E-01
3.639829E-02
2.564084E+00
6.574527E+00
1.503142E-01
2.259435E-02
2.233122E+00
4.986832E+00
1.456891E-01
2.122532E-02
1.682929E+00
2.832250E+00
9.828234E-02
9.659418E-03
7.232788E-01
5.231322E-01
6.128193E-02
3.755475E-03
tally 4:
1.326662E-01
1.760031E-02
0.000000E+00
0.000000E+00
1.369313E-01
1.875018E-02
2.569758E-01
6.603657E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.569758E-01
6.603657E-02
1.369313E-01
1.875018E-02
2.494513E-01
6.222596E-02
3.551002E-01
1.260962E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.551002E-01
1.260962E-01
2.494513E-01
6.222596E-02
3.610027E-01
1.303230E-01
4.252566E-01
1.808432E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.252566E-01
1.808432E-01
3.610027E-01
1.303230E-01
4.673502E-01
2.184162E-01
5.073957E-01
2.574504E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.073957E-01
2.574504E-01
4.673502E-01
2.184162E-01
4.818855E-01
2.322136E-01
4.949753E-01
2.450006E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.949753E-01
2.450006E-01
4.818855E-01
2.322136E-01
5.176611E-01
2.679730E-01
4.756654E-01
2.262576E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.756654E-01
2.262576E-01
5.176611E-01
2.679730E-01
4.591191E-01
2.107904E-01
3.930251E-01
1.544688E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.930251E-01
1.544688E-01
4.591191E-01
2.107904E-01
4.201532E-01
1.765287E-01
3.188250E-01
1.016494E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.188250E-01
1.016494E-01
4.201532E-01
1.765287E-01
2.898817E-01
8.403142E-02
1.509185E-01
2.277639E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.509185E-01
2.277639E-02
2.898817E-01
8.403142E-02
1.482175E-01
2.196844E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
tally 5:
7.719619E-01
5.959252E-01
8.168322E-02
6.672148E-03
1.272034E+00
1.618070E+00
1.458271E-01
2.126554E-02
1.779821E+00
3.167762E+00
2.705917E-01
7.321987E-02
2.584069E+00
6.677414E+00
3.176497E-01
1.009013E-01
2.726552E+00
7.434084E+00
3.562714E-01
1.269293E-01
2.740788E+00
7.511919E+00
3.690107E-01
1.361689E-01
2.563145E+00
6.569715E+00
3.296244E-01
1.086523E-01
2.233122E+00
4.986832E+00
3.106254E-01
9.648813E-02
1.682929E+00
2.832250E+00
2.709632E-01
7.342104E-02
7.232788E-01
5.231322E-01
9.343486E-02
8.730073E-03
cmfd indices
1.000000E+01
1.000000E+00
1.000000E+00
1.000000E+00
k cmfd
1.143597E+00
1.163387E+00
1.162391E+00
1.163351E+00
1.145721E+00
1.134785E+00
1.119048E+00
1.116124E+00
1.109085E+00
1.126581E+00
1.158559E+00
1.163719E+00
1.162162E+00
1.160856E+00
1.166709E+00
1.167223E+00
cmfd entropy
3.212002E+00
3.206393E+00
3.222109E+00
3.221996E+00
3.229978E+00
3.234615E+00
3.246512E+00
3.244634E+00
3.244312E+00
3.236922E+00
3.232693E+00
3.221849E+00
3.215716E+00
3.215968E+00
3.203758E+00
3.201798E+00
cmfd balance
8.26180E-03
4.27338E-03
2.62159E-03
2.40301E-03
2.08484E-03
1.58351E-03
1.59196E-03
1.87591E-03
1.94451E-03
1.91803E-03
1.90044E-03
1.87968E-03
2.55363E-03
2.39932E-03
2.25515E-03
1.53613E-03
cmfd dominance ratio
5.404E-01
5.406E-01
5.432E-01
5.454E-01
5.507E-01
5.578E-01
5.679E-01
5.671E-01
5.690E-01
5.637E-01
5.575E-01
5.474E-01
5.480E-01
5.467E-01
5.374E-01
5.321E-01
cmfd openmc source comparison
1.575499E-02
1.293688E-02
5.920734E-03
9.746850E-03
7.183896E-03
7.693485E-03
4.158805E-03
2.962505E-03
4.415044E-03
2.304495E-03
7.921580E-03
8.609203E-03
8.945198E-03
8.054204E-03
1.164189E-02
5.945645E-03
cmfd source
4.077779E-02
6.659143E-02
1.017198E-01
1.172269E-01
1.458410E-01
1.559801E-01
1.337001E-01
1.137262E-01
8.397376E-02
4.046290E-02

26
tests/regression_tests/cmfd_feed_rolling_window/settings.xml Normal file
View file

@ -0,0 +1,26 @@
<?xml version="1.0" encoding="UTF-8"?>
<settings>
<!-- Parameters for criticality calculation -->
<run_mode>eigenvalue</run_mode>
<batches>20</batches>
<inactive>10</inactive>
<particles>1000</particles>
<!-- Starting source -->
<source>
<space>
<type>box</type>
<parameters>-10 -1 -1 10 1 1</parameters>
</space>
</source>
<!-- Shannon Entropy -->
<mesh id="10">
<dimension> 10 1 1 </dimension>
<lower_left> -10.0 -1.0 -1.0 </lower_left>
<upper_right> 10.0 1.0 1.0 </upper_right>
</mesh>
<entropy_mesh>10</entropy_mesh>
</settings>

21
tests/regression_tests/cmfd_feed_rolling_window/tallies.xml Normal file
View file

@ -0,0 +1,21 @@
<?xml version="1.0"?>
<tallies>
<mesh id="1">
<type>regular</type>
<lower_left>-10 -1 -1 </lower_left>
<upper_right>10 1 1</upper_right>
<dimension>10 1 1</dimension>
</mesh>
<filter id="1">
<type>mesh</type>
<bins>1</bins>
</filter>
<tally id="1">
<filters>1</filters>
<scores>flux</scores>
</tally>
</tallies>

27
tests/regression_tests/cmfd_feed_rolling_window/test.py Normal file
View file

@ -0,0 +1,27 @@
from tests.testing_harness import CMFDTestHarness
from openmc import cmfd
def test_cmfd_feed_rolling_window():
"""Test 1 group CMFD solver with CMFD feedback"""
# Initialize and set CMFD mesh
cmfd_mesh = cmfd.CMFDMesh()
cmfd_mesh.lower_left = -10.0, -1.0, -1.0
cmfd_mesh.upper_right = 10.0, 1.0, 1.0
cmfd_mesh.dimension = 10, 1, 1
cmfd_mesh.albedo = 0.0, 0.0, 1.0, 1.0, 1.0, 1.0
# Initialize and run CMFDRun object
cmfd_run = cmfd.CMFDRun()
cmfd_run.mesh = cmfd_mesh
cmfd_run.tally_begin = 5
cmfd_run.feedback_begin = 10
cmfd_run.feedback = True
cmfd_run.gauss_seidel_tolerance = [1.e-15, 1.e-20]
cmfd_run.window_type = 'rolling'
cmfd_run.window_size = 5
cmfd_run.run()
# Initialize and run CMFD test harness
harness = CMFDTestHarness('statepoint.20.h5', cmfd_run)
harness.main()

0
tests/regression_tests/cmfd_restart/__init__.py Normal file
View file

43
tests/regression_tests/cmfd_restart/geometry.xml Normal file
View file

@ -0,0 +1,43 @@
<?xml version="1.0" encoding="UTF-8"?>
<geometry>
<!-- Definition of Cells -->
<cell id="1">
<universe>0</universe>
<region>-1 2 -3 4 -5 6</region>
<material>1</material>
</cell>
<!-- Defition of Surfaces -->
<surface id="1">
<type>x-plane</type>
<coeffs>10</coeffs>
<boundary> vacuum </boundary>
</surface>
<surface id="2">
<type>x-plane</type>
<coeffs>-10</coeffs>
<boundary> vacuum </boundary>
</surface>
<surface id="3">
<type>y-plane</type>
<coeffs>1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="4">
<type>y-plane</type>
<coeffs>-1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="5">
<type>z-plane</type>
<coeffs>1</coeffs>
<boundary>reflective</boundary>
</surface>
<surface id="6">
<type>z-plane</type>
<coeffs>-1</coeffs>
<boundary>reflective</boundary>
</surface>
</geometry>

12
tests/regression_tests/cmfd_restart/materials.xml Normal file
View file

@ -0,0 +1,12 @@
<?xml version="1.0"?>
<materials>
<!-- Definition of materials -->
<material id="1">
<density value="19" units="g/cc" />
<nuclide name="U235" wo="0.21" />
<nuclide name="U238" wo="0.68" />
<nuclide name="O16" wo="0.11" />
</material>
</materials>

488
tests/regression_tests/cmfd_restart/results_true.dat Normal file
View file

@ -0,0 +1,488 @@
k-combined:
1.169891E+00 6.289489E-03
tally 1:
1.173922E+01
1.385461E+01
2.164076E+01
4.699368E+01
2.906462E+01
8.464937E+01
3.382312E+01
1.147095E+02
3.632006E+01
1.323878E+02
3.655413E+01
1.341064E+02
3.347757E+01
1.124264E+02
2.931336E+01
8.607239E+01
2.182947E+01
4.789565E+01
1.147668E+01
1.325716E+01
tally 2:
2.298190E+01
2.667071E+01
1.600292E+01
1.293670E+01
4.268506E+01
9.161216E+01
3.022909E+01
4.598915E+01
5.680399E+01
1.623879E+02
4.033805E+01
8.196263E+01
6.814742E+01
2.331778E+02
4.851618E+01
1.182330E+02
7.392923E+01
2.740255E+02
5.253586E+01
1.384152E+02
7.332860E+01
2.698608E+02
5.227405E+01
1.371810E+02
6.830172E+01
2.340687E+02
4.867159E+01
1.188724E+02
5.885634E+01
1.736180E+02
4.170434E+01
8.719622E+01
4.371848E+01
9.592893E+01
3.106403E+01
4.844308E+01
2.338413E+01
2.752467E+01
1.636713E+01
1.347770E+01
tally 3:
1.538752E+01
1.196478E+01
1.079685E+00
6.010786E-02
2.911906E+01
4.269070E+01
1.822657E+00
1.671850E-01
3.885421E+01
7.608218E+01
2.541516E+00
3.262451E-01
4.673300E+01
1.097036E+02
2.885307E+00
4.214444E-01
5.059247E+01
1.283984E+02
3.222796E+00
5.237329E-01
5.034856E+01
1.272538E+02
3.230225E+00
5.273424E-01
4.688476E+01
1.103152E+02
2.941287E+00
4.363749E-01
4.013746E+01
8.077506E+01
2.634234E+00
3.520270E-01
2.996887E+01
4.509953E+01
1.946504E+00
1.919104E-01
1.575260E+01
1.248707E+01
1.020705E+00
5.413569E-02
tally 4:
3.049469E+00
4.677325E-01
0.000000E+00
0.000000E+00
2.770358E+00
3.879191E-01
5.514939E+00
1.528899E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.514939E+00
1.528899E+00
2.770358E+00
3.879191E-01
5.032131E+00
1.275040E+00
7.294002E+00
2.675589E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
7.294002E+00
2.675589E+00
5.032131E+00
1.275040E+00
7.036008E+00
2.490718E+00
8.668860E+00
3.776102E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
8.668860E+00
3.776102E+00
7.036008E+00
2.490718E+00
8.352414E+00
3.501945E+00
9.345868E+00
4.380719E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
9.345868E+00
4.380719E+00
8.352414E+00
3.501945E+00
9.093766E+00
4.158282E+00
9.223771E+00
4.270120E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
9.223771E+00
4.270120E+00
9.093766E+00
4.158282E+00
9.219150E+00
4.264346E+00
8.530966E+00
3.651778E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
8.530966E+00
3.651778E+00
9.219150E+00
4.264346E+00
8.690373E+00
3.785262E+00
7.204424E+00
2.604203E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
7.204424E+00
2.604203E+00
8.690373E+00
3.785262E+00
7.513640E+00
2.833028E+00
5.326721E+00
1.426975E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.326721E+00
1.426975E+00
7.513640E+00
2.833028E+00
5.662215E+00
1.607757E+00
2.848381E+00
4.093396E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.848381E+00
4.093396E-01
5.662215E+00
1.607757E+00
3.025812E+00
4.597241E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
tally 5:
1.538652E+01
1.196332E+01
2.252427E+00
2.605738E-01
2.911344E+01
4.267319E+01
3.873926E+00
7.615035E-01
3.884516E+01
7.604619E+01
5.280610E+00
1.414008E+00
4.672391E+01
1.096625E+02
6.261805E+00
1.983205E+00
5.058447E+01
1.283588E+02
6.733810E+00
2.278242E+00
5.033589E+01
1.271898E+02
6.714658E+00
2.273652E+00
4.687563E+01
1.102719E+02
6.215002E+00
1.956978E+00
4.013134E+01
8.075062E+01
5.253064E+00
1.396224E+00
2.996497E+01
4.508839E+01
3.818076E+00
7.509442E-01
1.574994E+01
1.248291E+01
2.219928E+00
2.515492E-01
cmfd indices
1.000000E+01
1.000000E+00
1.000000E+00
1.000000E+00
k cmfd
1.170416E+00
1.172966E+00
1.165537E+00
1.170979E+00
1.161922E+00
1.157523E+00
1.158873E+00
1.162877E+00
1.167101E+00
1.168130E+00
1.170570E+00
1.168115E+00
1.174081E+00
1.169458E+00
1.167848E+00
1.165116E+00
cmfd entropy
3.203643E+00
3.207943E+00
3.213367E+00
3.214360E+00
3.219634E+00
3.222232E+00
3.221744E+00
3.224544E+00
3.225990E+00
3.227769E+00
3.227417E+00
3.230728E+00
3.231662E+00
3.233316E+00
3.233193E+00
3.232564E+00
cmfd balance
4.00906E-03
4.43177E-03
3.15267E-03
3.51038E-03
2.05209E-03
2.06865E-03
1.50243E-03
1.58983E-03
1.56602E-03
1.17001E-03
9.50759E-04
9.07259E-04
1.00900E-03
1.06470E-03
1.16361E-03
9.75631E-04
cmfd dominance ratio
5.397E-01
5.425E-01
5.481E-01
5.473E-01
5.503E-01
5.502E-01
5.483E-01
5.520E-01
5.505E-01
3.216E-01
5.373E-01
5.517E-01
5.508E-01
5.524E-01
5.524E-01
5.523E-01
cmfd openmc source comparison
6.959834E-03
5.655657E-03
3.886186E-03
4.035116E-03
3.043277E-03
5.455474E-03
4.515310E-03
2.439840E-03
2.114032E-03
2.673132E-03
2.431749E-03
4.330928E-03
3.404647E-03
3.680298E-03
3.309620E-03
3.705541E-03
cmfd source
4.697085E-02
7.920706E-02
1.107968E-01
1.250932E-01
1.383930E-01
1.380648E-01
1.246874E-01
1.113705E-01
8.203754E-02
4.337882E-02

28
tests/regression_tests/cmfd_restart/settings.xml Normal file
View file

@ -0,0 +1,28 @@
<?xml version="1.0" encoding="UTF-8"?>
<settings>
<!-- Parameters for criticality calculation -->
<run_mode>eigenvalue</run_mode>
<batches>20</batches>
<inactive>10</inactive>
<particles>1000</particles>
<!-- Starting source -->
<source>
<space>
<type>box</type>
<parameters>-10 -1 -1 10 1 1</parameters>
</space>
</source>
<!-- Shannon Entropy -->
<mesh id="10">
<dimension> 10 1 1 </dimension>
<lower_left> -10.0 -1.0 -1.0 </lower_left>
<upper_right> 10.0 1.0 1.0 </upper_right>
</mesh>
<entropy_mesh>10</entropy_mesh>
<state_point batches="15 20" />
</settings>

21
tests/regression_tests/cmfd_restart/tallies.xml Normal file
View file

@ -0,0 +1,21 @@
<?xml version="1.0"?>
<tallies>
<mesh id="1">
<type>regular</type>
<lower_left>-10 -1 -1 </lower_left>
<upper_right>10 1 1</upper_right>
<dimension>10 1 1</dimension>
</mesh>
<filter id="1">
<type>mesh</type>
<bins>1</bins>
</filter>
<tally id="1">
<filters>1</filters>
<scores>flux</scores>
</tally>
</tallies>

72
tests/regression_tests/cmfd_restart/test.py Normal file
View file

@ -0,0 +1,72 @@
import glob
import os
import copy
from tests.testing_harness import CMFDTestHarness
from openmc import cmfd
import numpy as np
class CMFDRestartTestHarness(CMFDTestHarness):
def __init__(self, final_sp, restart_sp, cmfd_run1, cmfd_run2):
super().__init__(final_sp, cmfd_run1)
self._cmfd_restart_run = cmfd_run2
self._restart_sp = restart_sp
def execute_test(self):
try:
# Compare results from first CMFD run
self._test_output_created()
results = self._get_results()
results += self._cmfdrun_results
self._write_results(results)
self._compare_results()
# Run CMFD from restart file
statepoint = glob.glob(os.path.join(os.getcwd(), self._restart_sp))
assert len(statepoint) == 1
statepoint = statepoint[0]
self._cmfd_restart_run.run(args=['-r', statepoint])
# Compare results from second CMFD run
self._test_output_created()
self._create_cmfd_result_str(self._cmfd_restart_run)
results = self._get_results()
results += self._cmfdrun_results
self._write_results(results)
self._compare_results()
finally:
self._cleanup()
def test_cmfd_restart():
"""Test 1 group CMFD solver with restart run"""
# Initialize and set CMFD mesh, create a copy for second run
cmfd_mesh = cmfd.CMFDMesh()
cmfd_mesh.lower_left = -10.0, -1.0, -1.0
cmfd_mesh.upper_right = 10.0, 1.0, 1.0
cmfd_mesh.dimension = 10, 1, 1
cmfd_mesh.albedo = 0.0, 0.0, 1.0, 1.0, 1.0, 1.0
cmfd_mesh2 = copy.deepcopy(cmfd_mesh)
# Initialize and run first CMFDRun object
cmfd_run = cmfd.CMFDRun()
cmfd_run.mesh = cmfd_mesh
cmfd_run.tally_begin = 5
cmfd_run.feedback_begin = 5
cmfd_run.feedback = True
cmfd_run.gauss_seidel_tolerance = [1.e-15, 1.e-20]
cmfd_run.run()
# Initialize second CMFDRun object which will be run from restart file
cmfd_run2 = cmfd.CMFDRun()
cmfd_run2.mesh = cmfd_mesh2
cmfd_run2.tally_begin = 5
cmfd_run2.feedback_begin = 5
cmfd_run2.feedback = True
cmfd_run2.gauss_seidel_tolerance = [1.e-15, 1.e-20]
# Initialize and run CMFD restart test harness
harness = CMFDRestartTestHarness('statepoint.20.h5', 'statepoint.15.h5',
cmfd_run, cmfd_run2)
harness.main()