From 3f18aabda2f09043d329ce364fcb0991d463ba3a Mon Sep 17 00:00:00 2001 From: drewejohnson <1drew.e.johnson@gmail.com> Date: Wed, 8 Mar 2017 15:58:09 -0500 Subject: [PATCH 1/3] Iterate until convergence; Testing files to testing directory Signed-off-by: drewejohnson <1drew.e.johnson@gmail.com> --- iter_tester_1.inp | 112 --- iter_tester_2.inp | 112 --- iter_tester_3.inp | 112 --- iter_tester_4.inp | 112 --- iterator.py | 8 +- outputs.py | 3 +- iter_tester.inp => testing/iter_tester.inp | 4 +- .../iter_tester_16.inp | 18 +- testing/iter_tester_16.out | 911 ++++++++++++++++++ param_tester.txt => testing/param_tester.txt | 2 +- 10 files changed, 930 insertions(+), 464 deletions(-) delete mode 100644 iter_tester_1.inp delete mode 100644 iter_tester_2.inp delete mode 100644 iter_tester_3.inp delete mode 100644 iter_tester_4.inp rename iter_tester.inp => testing/iter_tester.inp (98%) rename iter_tester_5.inp => testing/iter_tester_16.inp (88%) create mode 100644 testing/iter_tester_16.out rename param_tester.txt => testing/param_tester.txt (77%) diff --git a/iter_tester_1.inp b/iter_tester_1.inp deleted file mode 100644 index 4a57325..0000000 --- a/iter_tester_1.inp +++ /dev/null @@ -1,112 +0,0 @@ -=t-newt parm=(centrm) -MSR iterative reflector optimization procedure -v7-252 - -read comp - -' --- NaCl-UCl3-PuCl3' -' --- 50/25/25 dens=3.60' -' --- 19.75 wt% eniched uranium' - Na-23 1 0 5.3767E-03 900 end - Cl-35 1 0 1.5991E-02 900 end - Cl-37 1 0 5.5155E-03 900 end - U-235 1 0 5.2556E-04 900 end - U-238 1 0 2.1355E-03 900 end - Pu-239 1 0 2.4994E-03 900 end - Pu-240 1 0 1.8892E-04 900 end - -' --- NaCl-UCl3-PuCl3' -' --- 50/25/25 dens=3.60' -' --- 19.75 wt% eniched uranium' - Na-23 10 0 5.3767E-03 900 end - Cl-35 10 0 1.5991E-02 900 end - Cl-37 10 0 5.5155E-03 900 end - U-235 10 0 5.2556E-04 900 end - U-238 10 0 2.1355E-03 900 end - Pu-239 10 0 2.4994E-03 900 end - Pu-240 10 0 1.8892E-04 900 end -' --- Heat exchanger fluid -' -- Nitrogen gas a.t.m. - n-14 2 1.0 293.0 end -' --- Piping -' -- Stainless steel - ss304 3 1.0 end -' --- Piping -' -- Stainless steel - ss304 30 1.0 end -' -- Reflector - ss304 4 1.0 600.0 end - -end comp - - -read celldata - - multiregion slab right_bdy=vacuum left_bdy=reflected origin=0.0 end - 1 55 - 3 60 - 2 70 - 30 75 - 10 100 - 4 120 end zone - -end celldata - -read model -Preliminary model - -read parameters - -drawit=yes -sn=6 -echo=yes -prtflux=yes - -end parameters - -read materials - - mix=1 com="in fuel" end - mix=10 com="out fuel" end - mix=2 pn=2 com="HX" end - mix=3 com="pipe" end - mix=4 com="refl" pn=2 end - -end materials - -read geometry - -global unit 1 - cuboid 99 120 0 20 -10.00000 - cuboid 10 70 0.0 20 0.0 - cuboid 11 72 70 20 0 - cuboid 12 83 72 20 0 - cuboid 13 85 83 20 0 - cuboid 14 120 85 20 0 - media 1 1 10 - media 3 1 11 - media 2 1 12 - media 3 1 13 - media 10 1 14 - cuboid 20 5 0 0 -10.00000 - cuboid 21 70 5 0 -10.00000 - cuboid 22 72 70 0 -10.00000 - cuboid 23 83 72 0 -10.00000 - cuboid 24 85 83 0 -10.00000 - cuboid 25 120 85 0 -10.00000 - media 3 1 20 - media 4 1 21 - media 3 1 22 - media 2 1 23 - media 3 1 24 - media 4 1 25 - boundary 99 4 4 -end geometry - -read bounds - +x=vac +y=refl -x=refl -y=vac -end bounds -end model - - -end diff --git a/iter_tester_2.inp b/iter_tester_2.inp deleted file mode 100644 index 0ef4d6c..0000000 --- a/iter_tester_2.inp +++ /dev/null @@ -1,112 +0,0 @@ -=t-newt parm=(centrm) -MSR iterative reflector optimization procedure -v7-252 - -read comp - -' --- NaCl-UCl3-PuCl3' -' --- 50/25/25 dens=3.60' -' --- 19.75 wt% eniched uranium' - Na-23 1 0 5.3767E-03 900 end - Cl-35 1 0 1.5991E-02 900 end - Cl-37 1 0 5.5155E-03 900 end - U-235 1 0 5.2556E-04 900 end - U-238 1 0 2.1355E-03 900 end - Pu-239 1 0 2.4994E-03 900 end - Pu-240 1 0 1.8892E-04 900 end - -' --- NaCl-UCl3-PuCl3' -' --- 50/25/25 dens=3.60' -' --- 19.75 wt% eniched uranium' - Na-23 10 0 5.3767E-03 900 end - Cl-35 10 0 1.5991E-02 900 end - Cl-37 10 0 5.5155E-03 900 end - U-235 10 0 5.2556E-04 900 end - U-238 10 0 2.1355E-03 900 end - Pu-239 10 0 2.4994E-03 900 end - Pu-240 10 0 1.8892E-04 900 end -' --- Heat exchanger fluid -' -- Nitrogen gas a.t.m. - n-14 2 1.0 293.0 end -' --- Piping -' -- Stainless steel - ss304 3 1.0 end -' --- Piping -' -- Stainless steel - ss304 30 1.0 end -' -- Reflector - ss304 4 1.0 600.0 end - -end comp - - -read celldata - - multiregion slab right_bdy=vacuum left_bdy=reflected origin=0.0 end - 1 55 - 3 60 - 2 70 - 30 75 - 10 100 - 4 120 end zone - -end celldata - -read model -Preliminary model - -read parameters - -drawit=yes -sn=6 -echo=yes -prtflux=yes - -end parameters - -read materials - - mix=1 com="in fuel" end - mix=10 com="out fuel" end - mix=2 pn=2 com="HX" end - mix=3 com="pipe" end - mix=4 com="refl" pn=2 end - -end materials - -read geometry - -global unit 1 - cuboid 99 120 0 20 -9.13394 - cuboid 10 70 0.0 20 0.0 - cuboid 11 72 70 20 0 - cuboid 12 83 72 20 0 - cuboid 13 85 83 20 0 - cuboid 14 120 85 20 0 - media 1 1 10 - media 3 1 11 - media 2 1 12 - media 3 1 13 - media 10 1 14 - cuboid 20 5 0 0 -9.13394 - cuboid 21 70 5 0 -9.13394 - cuboid 22 72 70 0 -9.13394 - cuboid 23 83 72 0 -9.13394 - cuboid 24 85 83 0 -9.13394 - cuboid 25 120 85 0 -9.13394 - media 3 1 20 - media 4 1 21 - media 3 1 22 - media 2 1 23 - media 3 1 24 - media 4 1 25 - boundary 99 4 4 -end geometry - -read bounds - +x=vac +y=refl -x=refl -y=vac -end bounds -end model - - -end diff --git a/iter_tester_3.inp b/iter_tester_3.inp deleted file mode 100644 index 6933ce0..0000000 --- a/iter_tester_3.inp +++ /dev/null @@ -1,112 +0,0 @@ -=t-newt parm=(centrm) -MSR iterative reflector optimization procedure -v7-252 - -read comp - -' --- NaCl-UCl3-PuCl3' -' --- 50/25/25 dens=3.60' -' --- 19.75 wt% eniched uranium' - Na-23 1 0 5.3767E-03 900 end - Cl-35 1 0 1.5991E-02 900 end - Cl-37 1 0 5.5155E-03 900 end - U-235 1 0 5.2556E-04 900 end - U-238 1 0 2.1355E-03 900 end - Pu-239 1 0 2.4994E-03 900 end - Pu-240 1 0 1.8892E-04 900 end - -' --- NaCl-UCl3-PuCl3' -' --- 50/25/25 dens=3.60' -' --- 19.75 wt% eniched uranium' - Na-23 10 0 5.3767E-03 900 end - Cl-35 10 0 1.5991E-02 900 end - Cl-37 10 0 5.5155E-03 900 end - U-235 10 0 5.2556E-04 900 end - U-238 10 0 2.1355E-03 900 end - Pu-239 10 0 2.4994E-03 900 end - Pu-240 10 0 1.8892E-04 900 end -' --- Heat exchanger fluid -' -- Nitrogen gas a.t.m. - n-14 2 1.0 293.0 end -' --- Piping -' -- Stainless steel - ss304 3 1.0 end -' --- Piping -' -- Stainless steel - ss304 30 1.0 end -' -- Reflector - ss304 4 1.0 600.0 end - -end comp - - -read celldata - - multiregion slab right_bdy=vacuum left_bdy=reflected origin=0.0 end - 1 55 - 3 60 - 2 70 - 30 75 - 10 100 - 4 120 end zone - -end celldata - -read model -Preliminary model - -read parameters - -drawit=yes -sn=6 -echo=yes -prtflux=yes - -end parameters - -read materials - - mix=1 com="in fuel" end - mix=10 com="out fuel" end - mix=2 pn=2 com="HX" end - mix=3 com="pipe" end - mix=4 com="refl" pn=2 end - -end materials - -read geometry - -global unit 1 - cuboid 99 120 0 20 -8.48840 - cuboid 10 70 0.0 20 0.0 - cuboid 11 72 70 20 0 - cuboid 12 83 72 20 0 - cuboid 13 85 83 20 0 - cuboid 14 120 85 20 0 - media 1 1 10 - media 3 1 11 - media 2 1 12 - media 3 1 13 - media 10 1 14 - cuboid 20 5 0 0 -8.48840 - cuboid 21 70 5 0 -8.48840 - cuboid 22 72 70 0 -8.48840 - cuboid 23 83 72 0 -8.48840 - cuboid 24 85 83 0 -8.48840 - cuboid 25 120 85 0 -8.48840 - media 3 1 20 - media 4 1 21 - media 3 1 22 - media 2 1 23 - media 3 1 24 - media 4 1 25 - boundary 99 4 4 -end geometry - -read bounds - +x=vac +y=refl -x=refl -y=vac -end bounds -end model - - -end diff --git a/iter_tester_4.inp b/iter_tester_4.inp deleted file mode 100644 index a157147..0000000 --- a/iter_tester_4.inp +++ /dev/null @@ -1,112 +0,0 @@ -=t-newt parm=(centrm) -MSR iterative reflector optimization procedure -v7-252 - -read comp - -' --- NaCl-UCl3-PuCl3' -' --- 50/25/25 dens=3.60' -' --- 19.75 wt% eniched uranium' - Na-23 1 0 5.3767E-03 900 end - Cl-35 1 0 1.5991E-02 900 end - Cl-37 1 0 5.5155E-03 900 end - U-235 1 0 5.2556E-04 900 end - U-238 1 0 2.1355E-03 900 end - Pu-239 1 0 2.4994E-03 900 end - Pu-240 1 0 1.8892E-04 900 end - -' --- NaCl-UCl3-PuCl3' -' --- 50/25/25 dens=3.60' -' --- 19.75 wt% eniched uranium' - Na-23 10 0 5.3767E-03 900 end - Cl-35 10 0 1.5991E-02 900 end - Cl-37 10 0 5.5155E-03 900 end - U-235 10 0 5.2556E-04 900 end - U-238 10 0 2.1355E-03 900 end - Pu-239 10 0 2.4994E-03 900 end - Pu-240 10 0 1.8892E-04 900 end -' --- Heat exchanger fluid -' -- Nitrogen gas a.t.m. - n-14 2 1.0 293.0 end -' --- Piping -' -- Stainless steel - ss304 3 1.0 end -' --- Piping -' -- Stainless steel - ss304 30 1.0 end -' -- Reflector - ss304 4 1.0 600.0 end - -end comp - - -read celldata - - multiregion slab right_bdy=vacuum left_bdy=reflected origin=0.0 end - 1 55 - 3 60 - 2 70 - 30 75 - 10 100 - 4 120 end zone - -end celldata - -read model -Preliminary model - -read parameters - -drawit=yes -sn=6 -echo=yes -prtflux=yes - -end parameters - -read materials - - mix=1 com="in fuel" end - mix=10 com="out fuel" end - mix=2 pn=2 com="HX" end - mix=3 com="pipe" end - mix=4 com="refl" pn=2 end - -end materials - -read geometry - -global unit 1 - cuboid 99 120 0 20 -8.01292 - cuboid 10 70 0.0 20 0.0 - cuboid 11 72 70 20 0 - cuboid 12 83 72 20 0 - cuboid 13 85 83 20 0 - cuboid 14 120 85 20 0 - media 1 1 10 - media 3 1 11 - media 2 1 12 - media 3 1 13 - media 10 1 14 - cuboid 20 5 0 0 -8.01292 - cuboid 21 70 5 0 -8.01292 - cuboid 22 72 70 0 -8.01292 - cuboid 23 83 72 0 -8.01292 - cuboid 24 85 83 0 -8.01292 - cuboid 25 120 85 0 -8.01292 - media 3 1 20 - media 4 1 21 - media 3 1 22 - media 2 1 23 - media 3 1 24 - media 4 1 25 - boundary 99 4 4 -end geometry - -read bounds - +x=vac +y=refl -x=refl -y=vac -end bounds -end model - - -end diff --git a/iterator.py b/iterator.py index 4a8895f..6e0c59d 100644 --- a/iterator.py +++ b/iterator.py @@ -61,7 +61,7 @@ def update_itervar(): # Assumes only one iteration variable for now _var = list(gp.iter_vars.keys())[0] - _des = gp.iter_vecs[_var][-1] * gp.k_target / gp.k_vec[-1] + _des = gp.iter_vecs[_var][-1] * (gp.k_target / gp.k_vec[-1]) ** 2 # may square the ratio of k for small convergence if _des > gp.iter_vars[_var][2]: @@ -83,7 +83,7 @@ def itermain(): gp.k_vec.append(gp.k_guess) - utils.vprint("Starting the iteration procedure....\n") + print("Starting the iteration procedure....\n") conv_flag = False conv_type = None @@ -96,7 +96,7 @@ def itermain(): utils.vprint(' done') stat, _k = parse_scale_out_eig(_iter_file.replace('.inp', '.out')) if stat: # successful operation - utils.vprint(" {0:<3}: {1}".format(_n, _k)) + print(" {0:<3}: {1}".format(_n, _k)) gp.k_vec.append(_k) else: utils.error('Could not find value of k-eff for iteration file {0}.inp\n' @@ -120,4 +120,6 @@ def itermain(): if _n == gp.iter_lim and conv_type is None: conv_type = 2 + print(' done') + output_landing(conv_type) diff --git a/outputs.py b/outputs.py index 9492c72..ae330ad 100644 --- a/outputs.py +++ b/outputs.py @@ -47,7 +47,7 @@ def parse_scale_out_eig(_ofile): return _stat, _rK -def output_landing(_out_type, _itervar=""): +def output_landing(_out_type): """ Write the output file according to the type of output :param _out_type: Flag indicating the reason the program terminated @@ -65,3 +65,4 @@ def output_landing(_out_type, _itervar=""): 2: "Terminated due to exceeding the iteration limit {}".format(gp.iter_lim) } print('End of operation. Status: {}'.format(out_messages[_out_type])) + diff --git a/iter_tester.inp b/testing/iter_tester.inp similarity index 98% rename from iter_tester.inp rename to testing/iter_tester.inp index 9afcc7d..99cb53c 100644 --- a/iter_tester.inp +++ b/testing/iter_tester.inp @@ -57,10 +57,10 @@ Preliminary model read parameters -drawit=yes +drawit=no sn=6 echo=yes -prtflux=yes +prtflux=no end parameters diff --git a/iter_tester_5.inp b/testing/iter_tester_16.inp similarity index 88% rename from iter_tester_5.inp rename to testing/iter_tester_16.inp index f178e15..7945334 100644 --- a/iter_tester_5.inp +++ b/testing/iter_tester_16.inp @@ -57,10 +57,10 @@ Preliminary model read parameters -drawit=yes +drawit=no sn=6 echo=yes -prtflux=yes +prtflux=no end parameters @@ -77,7 +77,7 @@ end materials read geometry global unit 1 - cuboid 99 120 0 20 -7.66603 + cuboid 99 120 0 20 -6.78229 cuboid 10 70 0.0 20 0.0 cuboid 11 72 70 20 0 cuboid 12 83 72 20 0 @@ -88,12 +88,12 @@ global unit 1 media 2 1 12 media 3 1 13 media 10 1 14 - cuboid 20 5 0 0 -7.66603 - cuboid 21 70 5 0 -7.66603 - cuboid 22 72 70 0 -7.66603 - cuboid 23 83 72 0 -7.66603 - cuboid 24 85 83 0 -7.66603 - cuboid 25 120 85 0 -7.66603 + cuboid 20 5 0 0 -6.78229 + cuboid 21 70 5 0 -6.78229 + cuboid 22 72 70 0 -6.78229 + cuboid 23 83 72 0 -6.78229 + cuboid 24 85 83 0 -6.78229 + cuboid 25 120 85 0 -6.78229 media 3 1 20 media 4 1 21 media 3 1 22 diff --git a/testing/iter_tester_16.out b/testing/iter_tester_16.out new file mode 100644 index 0000000..6cf7e2e --- /dev/null +++ b/testing/iter_tester_16.out @@ -0,0 +1,911 @@ + ******************************************************************************* + * * + * SCALE 6.2.1 * + * ------------- * + * March 2017 * + * * + * SCALE: A Comprehensive Modeling and Simulation Suite * + * for Nuclear Safety Analysis and Design * + * * + * Reactor and Nuclear Systems Division * + * Oak Ridge National Laboratory * + * * + * http://scale.ornl.gov * + * scalehelp@ornl.gov * + * * + ******************************************************************************* + ******************************************************************************* + + Job Information + --------------- + Job started on ME04W0358GRD04 on Wed 08/03/2017 15:55:49 + Working directory: C:\Users\AJOHNS~1\AppData\Local\Temp\scale.ajohnson400.12048 + Input file name : C:\Users\ajohnson400\Documents\classes\msr-refl-iter\testing\iter_tester_16.inp + Output file name : C:\Users\ajohnson400\Documents\classes\msr-refl-iter\testing\iter_tester_16.out + SCALE executable : C:\SCALE-6.2.1/bin/scale + + ******************************************************************************* + + + module t-newt will be called on Wed Mar 08 15:55:49 2017. + sequence specification record:=t-newt parm=(centrm) +Input Data: + + MSR iterative reflector optimization procedure + v7-252 + + read comp + + ' --- NaCl-UCl3-PuCl3' + ' --- 50/25/25 dens=3.60' + ' --- 19.75 wt% eniched uranium' + Na-23 1 0 5.3767E-03 900 end + Cl-35 1 0 1.5991E-02 900 end + Cl-37 1 0 5.5155E-03 900 end + U-235 1 0 5.2556E-04 900 end + U-238 1 0 2.1355E-03 900 end + Pu-239 1 0 2.4994E-03 900 end + Pu-240 1 0 1.8892E-04 900 end + + ' --- NaCl-UCl3-PuCl3' + ' --- 50/25/25 dens=3.60' + ' --- 19.75 wt% eniched uranium' + Na-23 10 0 5.3767E-03 900 end + Cl-35 10 0 1.5991E-02 900 end + Cl-37 10 0 5.5155E-03 900 end + U-235 10 0 5.2556E-04 900 end + U-238 10 0 2.1355E-03 900 end + Pu-239 10 0 2.4994E-03 900 end + Pu-240 10 0 1.8892E-04 900 end + ' --- Heat exchanger fluid + ' -- Nitrogen gas a.t.m. + n-14 2 1.0 293.0 end + ' --- Piping + ' -- Stainless steel + ss304 3 1.0 end + ' --- Piping + ' -- Stainless steel + ss304 30 1.0 end + ' -- Reflector + ss304 4 1.0 600.0 end + + end comp + + + read celldata + + multiregion slab right_bdy=vacuum left_bdy=reflected origin=0.0 end + 1 55 + 3 60 + 2 70 + 30 75 + 10 100 + 4 120 end zone + + end celldata + + read model + Preliminary model + + read parameters + + drawit=no + sn=6 + echo=yes + prtflux=no + + end parameters + + read materials + + mix=1 com="in fuel" end + mix=10 com="out fuel" end + mix=2 pn=2 com="HX" end + mix=3 com="pipe" end + mix=4 com="refl" pn=2 end + + end materials + + read geometry + + global unit 1 + cuboid 99 120 0 20 -6.78229 + cuboid 10 70 0.0 20 0.0 + cuboid 11 72 70 20 0 + cuboid 12 83 72 20 0 + cuboid 13 85 83 20 0 + cuboid 14 120 85 20 0 + media 1 1 10 + media 3 1 11 + media 2 1 12 + media 3 1 13 + media 10 1 14 + cuboid 20 5 0 0 -6.78229 + cuboid 21 70 5 0 -6.78229 + cuboid 22 72 70 0 -6.78229 + cuboid 23 83 72 0 -6.78229 + cuboid 24 85 83 0 -6.78229 + cuboid 25 120 85 0 -6.78229 + media 3 1 20 + media 4 1 21 + media 3 1 22 + media 2 1 23 + media 3 1 24 + media 4 1 25 + boundary 99 4 4 + end geometry + + read bounds + +x=vac +y=refl -x=refl -y=vac + end bounds + end model + + +end + + + + +************************************************************************************************************************************ +************************************************************************************************************************************ +************************************************************************************************************************************ +***** ***** +***** program verification information ***** +***** ***** +***** code system: SCALE version: 6.2 ***** +***** ***** +************************************************************************************************************************************ +************************************************************************************************************************************ +***** ***** +***** ***** +***** program: t-newt ***** +***** ***** +***** version: 6.2.1 (602001) ***** +***** ***** +***** username: ajohnson400 ***** +***** ***** +***** hostname: ME04W0358GRD04 ***** +***** ***** +***** ***** +***** date of execution: 2017-03-08 ***** +***** ***** +***** time of execution: 15:55:49 ***** +***** ***** +************************************************************************************************************************************ +************************************************************************************************************************************ +************************************************************************************************************************************ + + + +------------------------------------------------------------------------------------------------------------------------------------ + + preliminary model + +------------------------------------------------------------------------------------------------------------------------------------ + + + +1 *********************** + *********************** + ** Control Options ** + *********************** + *********************** + + Run problem to completion/test input only run= yes + Weighted cross-section collapse/do not collapse collapse= no + Source acceleration on/off accel= yes + Restart from earlier case/begin new case restart= no + CMFD acceleration (0/1/2/3 = no/rect/yes/part) cmfd= 0 + Number of fine mesh base grid/coarse grid in X xcmfd= 1 + Number of fine mesh base grid/coarse grid in Y ycmfd= 1 + CMFD Multigroup solution only (no 2g) cmfd2g= no + Create restart file upon completion savrest= yes + Calculation type: eigenvalue + Forward/adjoint solution forward= yes + +1 ********************** + ********************** + ** Output Options ** + ********************** + ********************** + + Print micro. cross-sections for problem nuclides prtxsec= no + Print macro. cross-sections for mixtures prtmxsec= no + Print balance table prtbalnc= no + Print mixing table prtmxtab= no + Print postscript graphics file for fluxes prtflux= no + Print postscript grid/mixture geometry graphics drawit= no + Print broad-group fluxes prtbroad= no + Echo convergence tracking data to std. output echo= yes + +1 *************************** + *************************** + ** I/O Unit Assignments ** + *************************** + *************************** + + Unit no. for mixing table mixtab= 92 + Unit no. of working library xnlib= 4 + Unit no. for weighted broad-group library wtdlib= 30 + Unit no. for homogenized macroscopic x-sect. lib. hmoglib= 13 + +1 ************************* + ************************* + ** Convergence Control ** + ************************* + ************************* + + Quadrature order sn= 6 + Maximum no. of inner itrs./outer/energy group inners= 5 + Maximum no. of outer energy groups outers= 250 + Performing thermal group subouter iterations therm= yes + Maximum no. of thermal group subouter itrs./outer suboutrs= 2 + Convergence criterion on inner iterations epsinner= 1.00E-04 + Convergence criterion on subouter iterations epsthrm= 1.00E-04 + Convergence criterion on outer iterations epsouter= 1.00E-04 + Convergence criterion on eigenvalue epseigen= 1.00E-04 + First guess for eigenvalue (k-eff) kguess= 1.000 + +1 ***************************** + ***************************** + ** Geometry Specifications ** + ***************************** + ***************************** + + Left side x-coordinate xmin= 0.00000E+00 + Right side x-coordinate xmax= 1.20000E+02 + Bottom side y-coordinate ymin= -6.78229E+00 + Top side y-coordinate ymax= 2.00000E+01 + Number of subdivisions in x direction nx= 4 + Number of subdivisions in y direction ny= 4 + + Boundary Conditions (0/1/2/3=vacuum/reflective/white/periodic): + Left (West) wbc= 1 + Right (East) ebc= 0 + Bottom (South) sbc= 0 + Top (North) nbc= 1 + + Body Characteristics. Rotation angle in degrees counterclockwise. + Size parameters. cyl: A=radius, B=no. of sides. cub: A=width, B=height. hex: A(B)=major(minor) radius. + Body Input Zone Subdiv --Center-- Size Parameters Rotation Chord + Type Unit Number nx ny x y A B Angle Information + ---- ----- ----- --- --- ------ ------ ----- ------ -------- ----------- + cub 1 1 4 4 60.000 6.609 120.000 26.782 0.000 + cub 1 2 0 0 35.000 10.000 70.000 20.000 0.000 + cub 1 3 0 0 71.000 10.000 2.000 20.000 0.000 + cub 1 4 0 0 77.500 10.000 11.000 20.000 0.000 + cub 1 5 0 0 84.000 10.000 2.000 20.000 0.000 + cub 1 6 0 0 102.500 10.000 35.000 20.000 0.000 + cub 1 7 0 0 2.500 -3.391 5.000 6.782 0.000 + cub 1 8 0 0 37.500 -3.391 65.000 6.782 0.000 + cub 1 9 0 0 71.000 -3.391 2.000 6.782 0.000 + cub 1 10 0 0 77.500 -3.391 11.000 6.782 0.000 + cub 1 11 0 0 84.000 -3.391 2.000 6.782 0.000 + cub 1 12 0 0 102.500 -3.391 35.000 6.782 0.000 + +1 ***************************** + ***************************** + ** Material Specifications ** + ***************************** + ***************************** + + Mat'l No. SCALE ID Pn order External source ID Description (if any) + 1 1 1 0 in fuel + 2 10 1 0 out fuel + 3 2 2 0 hx + 4 3 1 0 pipe + 5 4 2 0 refl + +1 ************************ + ************************ + ** Derived Parameters ** + ************************ + ************************ + + Total number of energy groups numen= 252 + Number of neutron groups igm= 252 + Number of gamma groups ipm= 0 + First thermal energy group iftg= 156 + Mixing table length nmix= 53 + Number of mixtures created nummat= 6 + Number of computational cells numcells= 42 + Number of computational surfaces (sides) numlines= 98 + Combine small grid cells together combine=no + Average number of sides/cell avnosides= 4.02 + Minimum number of sides/cell mnnosides= 4 + Maximum number of sides/cell mxnosides= 5 + Number of angles in quadrature set ndir= 24 + Maximum order of Legendre scattering maxpn= 2 + Number of flux moments nmom 6 + +1 ***************************** + ***************************** + ** Energy Group Structures ** + ***************************** + ***************************** + + Energy Energy Lethargy + Group Boundaries Boundaries + No. (eV) + 1 2.00000E+07 -6.93147E-01 127 3.17500E+01 1.26602E+01 + 2 1.73300E+07 -5.49854E-01 128 3.12500E+01 1.26761E+01 + 3 1.56800E+07 -4.49801E-01 129 3.00000E+01 1.27169E+01 + 4 1.45500E+07 -3.75006E-01 130 2.75000E+01 1.28039E+01 + 5 1.38400E+07 -3.24978E-01 131 2.50000E+01 1.28992E+01 + 6 1.28400E+07 -2.49980E-01 132 2.25000E+01 1.30046E+01 + 7 1.00000E+07 0.00000E+00 133 2.17500E+01 1.30385E+01 + 8 8.18700E+06 2.00038E-01 134 2.12000E+01 1.30641E+01 + 9 6.43400E+06 4.40989E-01 135 2.05000E+01 1.30977E+01 + 10 4.80000E+06 7.33969E-01 136 2.00000E+01 1.31224E+01 + 11 4.30400E+06 8.43040E-01 137 1.94000E+01 1.31528E+01 + 12 3.00000E+06 1.20397E+00 138 1.85000E+01 1.32003E+01 + 13 2.47900E+06 1.39473E+00 139 1.70000E+01 1.32849E+01 + 14 2.35400E+06 1.44647E+00 140 1.60000E+01 1.33455E+01 + 15 1.85000E+06 1.68740E+00 141 1.44000E+01 1.34509E+01 + 16 1.50000E+06 1.89712E+00 142 1.29000E+01 1.35609E+01 + 17 1.40000E+06 1.96611E+00 143 1.19000E+01 1.36416E+01 + 18 1.35600E+06 1.99805E+00 144 1.15000E+01 1.36757E+01 + 19 1.31700E+06 2.02723E+00 145 1.00000E+01 1.38155E+01 + 20 1.25000E+06 2.07944E+00 146 9.10000E+00 1.39098E+01 + 21 1.20000E+06 2.12026E+00 147 8.10000E+00 1.40262E+01 + 22 1.10000E+06 2.20727E+00 148 7.15000E+00 1.41510E+01 + 23 1.01000E+06 2.29263E+00 149 7.00000E+00 1.41722E+01 + 24 9.20000E+05 2.38597E+00 150 6.87500E+00 1.41902E+01 + 25 9.00000E+05 2.40795E+00 151 6.75000E+00 1.42086E+01 + 26 8.75000E+05 2.43612E+00 152 6.50000E+00 1.42463E+01 + 27 8.61100E+05 2.45213E+00 153 6.25000E+00 1.42855E+01 + 28 8.20000E+05 2.50104E+00 154 6.00000E+00 1.43263E+01 + 29 7.50000E+05 2.59027E+00 155 5.40000E+00 1.44317E+01 + 30 6.79000E+05 2.68972E+00 156 5.00000E+00 1.45087E+01 + 31 6.70000E+05 2.70306E+00 157 4.70000E+00 1.45705E+01 + 32 6.00000E+05 2.81341E+00 158 4.10000E+00 1.47071E+01 + 33 5.73000E+05 2.85945E+00 159 3.73000E+00 1.48017E+01 + 34 5.50000E+05 2.90042E+00 160 3.50000E+00 1.48653E+01 + 35 4.92000E+05 3.01186E+00 161 3.20000E+00 1.49549E+01 + 36 4.70000E+05 3.05761E+00 162 3.10000E+00 1.49867E+01 + 37 4.40000E+05 3.12357E+00 163 3.00000E+00 1.50195E+01 + 38 4.20000E+05 3.17009E+00 164 2.97000E+00 1.50295E+01 + 39 4.00000E+05 3.21888E+00 165 2.87000E+00 1.50638E+01 + 40 3.30000E+05 3.41125E+00 166 2.77000E+00 1.50992E+01 + 41 2.70000E+05 3.61192E+00 167 2.67000E+00 1.51360E+01 + 42 2.00000E+05 3.91202E+00 168 2.57000E+00 1.51742E+01 + 43 1.49000E+05 4.20639E+00 169 2.47000E+00 1.52139E+01 + 44 1.28300E+05 4.35597E+00 170 2.38000E+00 1.52510E+01 + 45 1.00000E+05 4.60517E+00 171 2.30000E+00 1.52852E+01 + 46 8.50000E+04 4.76769E+00 172 2.21000E+00 1.53251E+01 + 47 8.20000E+04 4.80362E+00 173 2.12000E+00 1.53667E+01 + 48 7.50000E+04 4.89285E+00 174 2.00000E+00 1.54249E+01 + 49 7.30000E+04 4.91988E+00 175 1.94000E+00 1.54554E+01 + 50 6.00000E+04 5.11600E+00 176 1.86000E+00 1.54975E+01 + 51 5.20000E+04 5.25910E+00 177 1.77000E+00 1.55471E+01 + 52 5.00000E+04 5.29832E+00 178 1.68000E+00 1.55993E+01 + 53 4.50000E+04 5.40368E+00 179 1.59000E+00 1.56544E+01 + 54 3.00000E+04 5.80914E+00 180 1.50000E+00 1.57126E+01 + 55 2.00000E+04 6.21461E+00 181 1.45000E+00 1.57465E+01 + 56 1.70000E+04 6.37713E+00 182 1.40000E+00 1.57816E+01 + 57 1.30000E+04 6.64539E+00 183 1.35000E+00 1.58180E+01 + 58 9.50000E+03 6.95905E+00 184 1.30000E+00 1.58557E+01 + 59 8.03000E+03 7.12716E+00 185 1.25000E+00 1.58950E+01 + 60 5.70000E+03 7.46987E+00 186 1.22500E+00 1.59152E+01 + 61 3.90000E+03 7.84936E+00 187 1.20000E+00 1.59358E+01 + 62 3.74000E+03 7.89125E+00 188 1.17500E+00 1.59568E+01 + 63 3.00000E+03 8.11173E+00 189 1.15000E+00 1.59783E+01 + 64 2.50000E+03 8.29405E+00 190 1.14000E+00 1.59871E+01 + 65 2.25000E+03 8.39941E+00 191 1.13000E+00 1.59959E+01 + 66 2.20000E+03 8.42188E+00 192 1.12000E+00 1.60048E+01 + 67 1.80000E+03 8.62255E+00 193 1.11000E+00 1.60137E+01 + 68 1.55000E+03 8.77209E+00 194 1.10000E+00 1.60228E+01 + 69 1.50000E+03 8.80488E+00 195 1.09000E+00 1.60319E+01 + 70 1.15000E+03 9.07058E+00 196 1.08000E+00 1.60411E+01 + 71 9.50000E+02 9.26163E+00 197 1.07000E+00 1.60504E+01 + 72 6.83000E+02 9.59160E+00 198 1.06000E+00 1.60598E+01 + 73 6.70000E+02 9.61082E+00 199 1.05000E+00 1.60693E+01 + 74 5.50000E+02 9.80818E+00 200 1.04000E+00 1.60789E+01 + 75 3.05000E+02 1.03978E+01 201 1.03000E+00 1.60885E+01 + 76 2.85000E+02 1.04656E+01 202 1.02000E+00 1.60983E+01 + 77 2.40000E+02 1.06375E+01 203 1.01000E+00 1.61081E+01 + 78 2.20000E+02 1.07245E+01 204 1.00000E+00 1.61181E+01 + 79 2.09500E+02 1.07734E+01 205 9.75000E-01 1.61434E+01 + 80 2.07400E+02 1.07834E+01 206 9.50000E-01 1.61694E+01 + 81 2.02000E+02 1.08098E+01 207 9.25000E-01 1.61961E+01 + 82 1.93000E+02 1.08554E+01 208 9.00000E-01 1.62235E+01 + 83 1.91500E+02 1.08632E+01 209 8.50000E-01 1.62806E+01 + 84 1.88500E+02 1.08790E+01 210 8.00000E-01 1.63412E+01 + 85 1.87700E+02 1.08833E+01 211 7.50000E-01 1.64058E+01 + 86 1.80000E+02 1.09251E+01 212 7.00000E-01 1.64748E+01 + 87 1.70000E+02 1.09823E+01 213 6.50000E-01 1.65489E+01 + 88 1.43000E+02 1.11553E+01 214 6.25000E-01 1.65881E+01 + 89 1.22000E+02 1.13141E+01 215 6.00000E-01 1.66289E+01 + 90 1.19000E+02 1.13390E+01 216 5.50000E-01 1.67159E+01 + 91 1.17500E+02 1.13517E+01 217 5.00000E-01 1.68112E+01 + 92 1.16000E+02 1.13645E+01 218 4.50000E-01 1.69166E+01 + 93 1.13000E+02 1.13907E+01 219 4.00000E-01 1.70344E+01 + 94 1.08000E+02 1.14360E+01 220 3.75000E-01 1.70989E+01 + 95 1.05000E+02 1.14641E+01 221 3.50000E-01 1.71679E+01 + 96 1.01200E+02 1.15010E+01 222 3.25000E-01 1.72420E+01 + 97 9.70000E+01 1.15434E+01 223 3.00000E-01 1.73221E+01 + 98 9.00000E+01 1.16183E+01 224 2.75000E-01 1.74091E+01 + 99 8.17000E+01 1.17150E+01 225 2.50000E-01 1.75044E+01 + 100 8.00000E+01 1.17361E+01 226 2.25000E-01 1.76098E+01 + 101 7.60000E+01 1.17874E+01 227 2.00000E-01 1.77275E+01 + 102 7.20000E+01 1.18414E+01 228 1.75000E-01 1.78611E+01 + 103 6.75000E+01 1.19060E+01 229 1.50000E-01 1.80152E+01 + 104 6.50000E+01 1.19437E+01 230 1.25000E-01 1.81975E+01 + 105 6.30000E+01 1.19750E+01 231 1.00000E-01 1.84207E+01 + 106 6.10000E+01 1.20072E+01 232 9.00000E-02 1.85260E+01 + 107 5.80000E+01 1.20577E+01 233 8.00000E-02 1.86438E+01 + 108 5.34000E+01 1.21403E+01 234 7.00000E-02 1.87774E+01 + 109 5.06000E+01 1.21941E+01 235 6.00000E-02 1.89315E+01 + 110 4.83000E+01 1.22407E+01 236 5.00000E-02 1.91138E+01 + 111 4.52000E+01 1.23070E+01 237 4.00000E-02 1.93370E+01 + 112 4.40000E+01 1.23339E+01 238 3.00000E-02 1.96247E+01 + 113 4.24000E+01 1.23709E+01 239 2.53000E-02 1.97950E+01 + 114 4.10000E+01 1.24045E+01 240 1.00000E-02 2.07233E+01 + 115 3.96000E+01 1.24393E+01 241 7.50000E-03 2.10109E+01 + 116 3.91000E+01 1.24520E+01 242 5.00000E-03 2.14164E+01 + 117 3.80000E+01 1.24805E+01 243 4.00000E-03 2.16396E+01 + 118 3.76300E+01 1.24903E+01 244 3.00000E-03 2.19272E+01 + 119 3.72700E+01 1.24999E+01 245 2.50000E-03 2.21096E+01 + 120 3.71300E+01 1.25037E+01 246 2.00000E-03 2.23327E+01 + 121 3.70000E+01 1.25072E+01 247 1.50000E-03 2.26204E+01 + 122 3.60000E+01 1.25346E+01 248 1.20000E-03 2.28435E+01 + 123 3.55000E+01 1.25486E+01 249 1.00000E-03 2.30259E+01 + 124 3.50000E+01 1.25627E+01 250 7.50000E-04 2.33135E+01 + 125 3.37500E+01 1.25991E+01 251 5.00000E-04 2.37190E+01 + 126 3.32500E+01 1.26140E+01 252 1.00000E-04 2.53284E+01 + 253 1.00000E-05 2.76310E+01 + +1 *************************** + *************************** + ** Quadrature Parameters ** + *************************** + *************************** + + Quadrature angles and weights for first quadrant: + ( 6 angles/quadrant for S- 6 quadrature.) + + Angle No. Mu Eta Weight Angle in x-y Elevation + plane (degrees) (degrees) + 1 2.30091940D-01 2.30091940D-01 4.23616400D-02 45.00 71.01 + 2 6.88134320D-01 2.30091940D-01 4.09716925D-02 18.49 43.48 + 3 9.45576760D-01 2.30091940D-01 4.23616400D-02 13.68 13.30 + 4 2.30091940D-01 6.88134320D-01 4.09716925D-02 71.51 43.48 + 5 6.88134320D-01 6.88134320D-01 4.09716925D-02 45.00 13.30 + 6 2.30091940D-01 9.45576760D-01 4.23616400D-02 76.32 13.30 + + +1 ********************************** + ********************************** + ** Quadrature Reflection Angles ** + ********************************** + ********************************** + (Angles are in degrees) + + 12 primary (incoming) angles on each side. + 12 reflected (outgoing) angles on each side. + + Side No. (slope) Direction No. Primary angle (slope/elevation) Reflecting Angle (slope/elevation) + + 2 ( -90.00) 1 1 ( 45.00/ 71.01) 7 ( 135.00/ 71.01) + 2 2 ( 18.49/ 43.48) 8 ( 161.51/ 43.48) + 3 3 ( 13.68/ 13.30) 9 ( 166.32/ 13.30) + 4 4 ( 71.51/ 43.48) 10 ( 108.49/ 43.48) + 5 5 ( 45.00/ 13.30) 11 ( 135.00/ 13.30) + 6 6 ( 76.32/ 13.30) 12 ( 103.68/ 13.30) + 7 19 ( -45.00/ 71.01) 13 (-135.00/ 71.01) + 8 20 ( -18.49/ 43.48) 14 (-161.51/ 43.48) + 9 21 ( -13.68/ 13.30) 15 (-166.32/ 13.30) + 10 22 ( -71.51/ 43.48) 16 (-108.49/ 43.48) + 11 23 ( -45.00/ 13.30) 17 (-135.00/ 13.30) + 12 24 ( -76.32/ 13.30) 18 (-103.68/ 13.30) + + 4 ( 180.00) 1 13 (-135.00/ 71.01) 7 ( 135.00/ 71.01) + 2 14 (-161.51/ 43.48) 8 ( 161.51/ 43.48) + 3 15 (-166.32/ 13.30) 9 ( 166.32/ 13.30) + 4 16 (-108.49/ 43.48) 10 ( 108.49/ 43.48) + 5 17 (-135.00/ 13.30) 11 ( 135.00/ 13.30) + 6 18 (-103.68/ 13.30) 12 ( 103.68/ 13.30) + 7 19 ( -45.00/ 71.01) 1 ( 45.00/ 71.01) + 8 20 ( -18.49/ 43.48) 2 ( 18.49/ 43.48) + 9 21 ( -13.68/ 13.30) 3 ( 13.68/ 13.30) + 10 22 ( -71.51/ 43.48) 4 ( 71.51/ 43.48) + 11 23 ( -45.00/ 13.30) 5 ( 45.00/ 13.30) + 12 24 ( -76.32/ 13.30) 6 ( 76.32/ 13.30) + + + + + + Scattering constants for first quadrant + ( 6 moments for p-2 scattering) + Angle ------------------------- Moment No. ------------------------- + No. 0 1 2 3 4 5 + 1 1.00000E+00 3.98531E-01 -2.81804E-01 -9.40460E-01 -1.44988E-01 -1.15182E+00 + 2 1.00000E+00 1.19188E+00 -2.81804E-01 4.70230E-01 -4.33616E-01 -5.75912E-01 + 3 1.00000E+00 1.63779E+00 -2.81804E-01 1.88092E+00 -5.95839E-01 1.47310E-08 + 4 1.00000E+00 3.98531E-01 -8.42789E-01 -9.40460E-01 -4.33616E-01 -1.97499E-08 + 5 1.00000E+00 1.19188E+00 -8.42789E-01 4.70230E-01 -1.29681E+00 5.75912E-01 + 6 1.00000E+00 3.98531E-01 -1.15809E+00 -9.40460E-01 -5.95839E-01 1.15182E+00 + + +1 ************************** + ************************** + **** Mixture Volumes ***** + ************************** + ************************** + + + Mixture Volume Volume + ID (cc/unit length) Fraction + 0 0.00000E+00 0.00000E+00 + 1 1.40000E+03 4.35611E-01 + 10 7.00000E+02 2.17806E-01 + 2 2.94605E+02 9.16667E-02 + 3 1.41041E+02 4.38849E-02 + 4 6.78229E+02 2.11032E-01 + Total 3.21387E+03 1.00000E+00 + + + Memory size of attenuation factors(MB)= 5.86 + + + + ===================================================================================== + Outer iteration sweep begins. + + Outer Eigen- Eigenvalue Max Flux Max Flux Max Fuel Max Fuel Inners + It. # value Delta Delta Location(r,g) Delta Location(r,g) Cnvrged + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + 1 1.00000 0.000E+00 6.185E+07 ( 3,252) 1.000E+00 ( 3, 0) F + 2 0.69905 4.305E-01 1.041E+05 ( 4,248) 1.182E+00 ( 25, 0) F + 3 1.02592 3.186E-01 7.092E+03 ( 4,240) 2.490E-01 ( 4, 0) F + 4 0.91001 1.274E-01 3.168E+02 ( 4, 1) 3.571E-01 ( 3, 0) F + 5 1.01478 1.032E-01 2.193E+02 ( 4, 1) 1.614E-01 ( 17, 0) F + 6 0.97406 4.180E-02 1.180E+02 ( 34, 1) 2.062E-01 ( 3, 0) F + 7 1.00581 3.157E-02 1.092E+02 ( 3, 1) 1.171E-01 ( 17, 0) F + 8 0.99252 1.340E-02 1.063E+02 ( 3, 1) 1.057E-01 ( 3, 0) F + 9 1.00205 9.518E-03 1.013E+02 ( 3, 1) 6.473E-02 ( 17, 0) F + 10 0.99783 4.235E-03 1.004E+02 ( 2, 1) 5.093E-02 ( 17, 0) F + 11 1.00071 2.876E-03 1.740E+00 ( 1,235) 3.229E-02 ( 17, 0) F + 12 0.99938 1.329E-03 9.695E-01 ( 1,235) 2.389E-02 ( 17, 0) F + 13 1.00025 8.700E-04 3.899E-01 ( 1,234) 1.551E-02 ( 17, 0) F + 14 0.99984 4.051E-04 1.297E-01 ( 1,234) 1.113E-02 ( 17, 0) F + 15 1.00004 1.952E-04 4.070E-02 ( 1,234) 7.368E-03 ( 17, 0) F + 16 1.00000 3.398E-05 2.138E-02 ( 1, 76) 5.159E-03 ( 17, 0) F + 17 1.00001 3.594E-06 1.530E-02 ( 1, 76) 3.518E-03 ( 17, 0) F + 18 1.00001 4.401E-07 1.073E-02 ( 1, 76) 2.417E-03 ( 17, 0) F + 19 1.00001 7.860E-07 7.422E-03 ( 1, 88) 1.663E-03 ( 22, 0) F + 20 1.00001 1.253E-06 5.059E-03 ( 1, 88) 1.145E-03 ( 22, 0) F + 21 1.00000 1.320E-06 3.378E-03 ( 1, 88) 7.891E-04 ( 22, 0) F + 22 1.00000 1.197E-06 2.207E-03 ( 1, 88) 5.438E-04 ( 22, 0) F + 23 1.00000 1.002E-06 1.410E-03 ( 1, 88) 3.746E-04 ( 22, 0) F + 24 1.00000 8.019E-07 8.788E-04 ( 1, 88) 2.575E-04 ( 22, 0) F + 25 1.00000 6.181E-07 5.309E-04 ( 1, 88) 1.772E-04 ( 22, 0) F + 26 1.00000 4.578E-07 3.103E-04 ( 1, 88) 1.217E-04 ( 22, 0) F + 27 1.00000 3.064E-07 1.745E-04 ( 1, 88) 7.985E-05 ( 22, 0) F + 28 1.00000 1.363E-07 1.091E-04 ( 2, 74) 4.820E-05 ( 17, 0) F + k-eff = 1.00000028 + + + + Four-Factor Estimate of k-infinity. Fast/Thermal boundary: 0.6250 eV + Fiss. neutrons/thermal abs. (eta): 1.596735 + Thermal utilization (f): 0.327632 + Resonance Escape probability (p): 0.000148 + Fast-fission factor (epsilon): 28640.652159 + -------------- + Infinite neutron multiplication 2.222343 + + Six-Factor Estimate of k-effective. + Fast non-leakage prob. (Pf): 0.449976 + Thermal non-leakage prob. (Pt): 0.999993 + -------------- + Effective neutron multiplication: 0.999993 + + + Alternate 3-Group Formulation Terms. + Fast/Res. boundary: 0.8210 MeV. Res./Therm. boundary: 0.6250eV + Fiss. neutrons/thermal abs. (eta): 1.596735 + Thermal utilization (f): 0.327632 + Resonance Escape probability (p): 0.000278 + Fast Escape probability (pf): 0.533425 + Fast-fission factor (epsilon): 28640.652159 + Fast non-leakage prob. (Pf): 0.728642 + Resonance non-leakage prob. (Pr): 0.617554 + Thermal non-leakage prob. (Pt): 0.999993 + + + + + Critical Spectra Calculated with the B1 Method + Input buckling of 0.00000E+00 gives a corresponding k-eff of 2.08115E+00 + Equivalent core height is infinite. + + Group Upper Critical Inf.Med. Critical Inf.Med. + No Energy(eV) Flux Flux Adj.Flux Adj.Flux + + 1 2.0000E+07 0.000E+00 0.000E+00 3.552E-03 3.552E-03 + 2 1.7330E+07 3.568E-11 3.568E-11 5.029E-03 5.029E-03 + 3 1.5680E+07 1.429E-07 1.429E-07 6.317E-03 6.317E-03 + 4 1.4550E+07 5.240E-06 5.240E-06 8.934E-03 8.934E-03 + 5 1.3840E+07 1.861E-05 1.861E-05 5.622E-03 5.622E-03 + 6 1.2840E+07 7.993E-05 7.993E-05 1.528E-03 1.528E-03 + 7 1.0000E+07 2.948E-04 2.948E-04 1.811E-03 1.811E-03 + 8 8.1870E+06 8.265E-04 8.265E-04 1.427E-03 1.427E-03 + 9 6.4340E+06 2.021E-03 2.021E-03 1.141E-03 1.141E-03 + 10 4.8000E+06 3.293E-03 3.293E-03 3.122E-03 3.122E-03 + 11 4.3040E+06 5.320E-03 5.320E-03 9.614E-04 9.614E-04 + 12 3.0000E+06 8.209E-03 8.209E-03 1.863E-03 1.863E-03 + 13 2.4790E+06 9.457E-03 9.457E-03 6.921E-03 6.921E-03 + 14 2.3540E+06 1.137E-02 1.137E-02 1.488E-03 1.488E-03 + 15 1.8500E+06 1.393E-02 1.393E-02 1.711E-03 1.711E-03 + 16 1.5000E+06 1.495E-02 1.495E-02 5.224E-03 5.224E-03 + 17 1.4000E+06 1.615E-02 1.615E-02 1.133E-02 1.133E-02 + 18 1.3560E+06 1.611E-02 1.611E-02 1.246E-02 1.246E-02 + 19 1.3170E+06 1.709E-02 1.709E-02 6.997E-03 6.997E-03 + 20 1.2500E+06 1.788E-02 1.788E-02 8.983E-03 8.983E-03 + 21 1.2000E+06 1.918E-02 1.918E-02 4.250E-03 4.250E-03 + 22 1.1000E+06 2.034E-02 2.034E-02 4.325E-03 4.325E-03 + 23 1.0100E+06 2.345E-02 2.345E-02 3.945E-03 3.945E-03 + 24 9.2000E+05 2.169E-02 2.169E-02 1.678E-02 1.678E-02 + 25 9.0000E+05 2.322E-02 2.322E-02 1.307E-02 1.307E-02 + 26 8.7500E+05 2.078E-02 2.078E-02 2.291E-02 2.291E-02 + 27 8.6110E+05 2.472E-02 2.472E-02 7.544E-03 7.544E-03 + 28 8.2000E+05 2.314E-02 2.314E-02 4.105E-03 4.105E-03 + 29 7.5000E+05 2.587E-02 2.587E-02 3.652E-03 3.652E-03 + 30 6.7900E+05 2.596E-02 2.596E-02 2.707E-02 2.707E-02 + 31 6.7000E+05 3.058E-02 3.058E-02 3.260E-03 3.260E-03 + 32 6.0000E+05 3.244E-02 3.244E-02 7.795E-03 7.795E-03 + 33 5.7300E+05 3.276E-02 3.276E-02 8.724E-03 8.724E-03 + 34 5.5000E+05 3.503E-02 3.503E-02 3.182E-03 3.182E-03 + 35 4.9200E+05 3.007E-02 3.007E-02 7.646E-03 7.646E-03 + 36 4.7000E+05 3.192E-02 3.192E-02 5.277E-03 5.277E-03 + 37 4.4000E+05 3.038E-02 3.038E-02 7.423E-03 7.423E-03 + 38 4.2000E+05 2.537E-02 2.537E-02 7.026E-03 7.026E-03 + 39 4.0000E+05 3.164E-02 3.164E-02 1.767E-03 1.767E-03 + 40 3.3000E+05 3.656E-02 3.656E-02 1.654E-03 1.654E-03 + 41 2.7000E+05 2.823E-02 2.823E-02 1.071E-03 1.071E-03 + 42 2.0000E+05 2.897E-02 2.897E-02 1.052E-03 1.052E-03 + 43 1.4900E+05 2.663E-02 2.663E-02 1.984E-03 1.984E-03 + 44 1.2830E+05 2.187E-02 2.187E-02 1.162E-03 1.162E-03 + 45 1.0000E+05 1.936E-02 1.936E-02 1.716E-03 1.716E-03 + 46 8.5000E+04 2.058E-02 2.058E-02 7.597E-03 7.597E-03 + 47 8.2000E+04 2.009E-02 2.009E-02 3.031E-03 3.031E-03 + 48 7.5000E+04 1.631E-02 1.631E-02 9.873E-03 9.873E-03 + 49 7.3000E+04 1.713E-02 1.713E-02 1.347E-03 1.347E-03 + 50 6.0000E+04 1.054E-02 1.054E-02 1.783E-03 1.783E-03 + 51 5.2000E+04 1.040E-02 1.040E-02 6.439E-03 6.439E-03 + 52 5.0000E+04 1.369E-02 1.369E-02 2.379E-03 2.379E-03 + 53 4.5000E+04 8.834E-03 8.834E-03 6.064E-04 6.064E-04 + 54 3.0000E+04 9.879E-03 9.879E-03 5.782E-04 5.782E-04 + 55 2.0000E+04 7.048E-03 7.048E-03 1.372E-03 1.372E-03 + 56 1.7000E+04 5.109E-03 5.109E-03 8.112E-04 8.112E-04 + 57 1.3000E+04 4.310E-03 4.310E-03 6.780E-04 6.780E-04 + 58 9.5000E+03 2.643E-03 2.643E-03 1.229E-03 1.229E-03 + 59 8.0300E+03 2.589E-03 2.589E-03 5.961E-04 5.961E-04 + 60 5.7000E+03 2.028E-03 2.028E-03 5.301E-04 5.301E-04 + 61 3.9000E+03 1.648E-03 1.648E-03 4.770E-03 4.770E-03 + 62 3.7400E+03 1.377E-03 1.377E-03 9.014E-04 9.014E-04 + 63 3.0000E+03 9.725E-04 9.725E-04 1.085E-03 1.085E-03 + 64 2.5000E+03 1.239E-03 1.239E-03 1.877E-03 1.877E-03 + 65 2.2500E+03 1.541E-03 1.541E-03 8.770E-03 8.770E-03 + 66 2.2000E+03 1.884E-03 1.884E-03 9.857E-04 9.857E-04 + 67 1.8000E+03 1.898E-03 1.898E-03 1.309E-03 1.309E-03 + 68 1.5500E+03 1.787E-03 1.787E-03 5.975E-03 5.975E-03 + 69 1.5000E+03 1.450E-03 1.450E-03 7.319E-04 7.319E-04 + 70 1.1500E+03 1.031E-03 1.031E-03 9.830E-04 9.830E-04 + 71 9.5000E+02 9.211E-04 9.211E-04 5.652E-04 5.652E-04 + 72 6.8300E+02 8.488E-04 8.488E-04 9.507E-03 9.507E-03 + 73 6.7000E+02 7.875E-04 7.875E-04 9.190E-04 9.190E-04 + 74 5.5000E+02 5.632E-04 5.632E-04 3.052E-04 3.052E-04 + 75 3.0500E+02 4.941E-04 4.941E-04 2.603E-03 2.603E-03 + 76 2.8500E+02 4.917E-04 4.917E-04 1.033E-03 1.033E-03 + 77 2.4000E+02 4.763E-04 4.763E-04 1.968E-03 1.968E-03 + 78 2.2000E+02 4.778E-04 4.778E-04 3.529E-03 3.529E-03 + 79 2.0950E+02 4.131E-04 4.131E-04 1.642E-02 1.642E-02 + 80 2.0740E+02 4.638E-04 4.638E-04 6.744E-03 6.744E-03 + 81 2.0200E+02 4.401E-04 4.401E-04 3.842E-03 3.842E-03 + 82 1.9300E+02 4.430E-04 4.430E-04 2.151E-02 2.151E-02 + 83 1.9150E+02 3.951E-04 3.951E-04 1.044E-02 1.044E-02 + 84 1.8850E+02 5.118E-04 5.118E-04 3.920E-02 3.920E-02 + 85 1.8770E+02 4.593E-04 4.593E-04 4.061E-03 4.061E-03 + 86 1.8000E+02 4.402E-04 4.402E-04 2.919E-03 2.919E-03 + 87 1.7000E+02 4.146E-04 4.146E-04 9.652E-04 9.652E-04 + 88 1.4300E+02 3.872E-04 3.872E-04 1.038E-03 1.038E-03 + 89 1.2200E+02 3.577E-04 3.577E-04 6.464E-03 6.464E-03 + 90 1.1900E+02 3.474E-04 3.474E-04 1.273E-02 1.273E-02 + 91 1.1750E+02 3.260E-04 3.260E-04 1.212E-02 1.212E-02 + 92 1.1600E+02 3.720E-04 3.720E-04 6.112E-03 6.112E-03 + 93 1.1300E+02 3.679E-04 3.679E-04 3.455E-03 3.455E-03 + 94 1.0800E+02 3.297E-04 3.297E-04 5.520E-03 5.520E-03 + 95 1.0500E+02 3.301E-04 3.301E-04 4.200E-03 4.200E-03 + 96 1.0120E+02 3.465E-04 3.465E-04 3.813E-03 3.813E-03 + 97 9.7000E+01 3.200E-04 3.200E-04 2.085E-03 2.085E-03 + 98 9.0000E+01 3.033E-04 3.033E-04 1.636E-03 1.636E-03 + 99 8.1700E+01 2.934E-04 2.934E-04 6.939E-03 6.939E-03 + 100 8.0000E+01 3.049E-04 3.049E-04 2.846E-03 2.846E-03 + 101 7.6000E+01 2.784E-04 2.784E-04 2.743E-03 2.743E-03 + 102 7.2000E+01 2.911E-04 2.911E-04 2.197E-03 2.197E-03 + 103 6.7500E+01 2.623E-04 2.623E-04 3.693E-03 3.693E-03 + 104 6.5000E+01 2.764E-04 2.764E-04 4.526E-03 4.526E-03 + 105 6.3000E+01 2.696E-04 2.696E-04 4.436E-03 4.436E-03 + 106 6.1000E+01 2.541E-04 2.541E-04 2.851E-03 2.851E-03 + 107 5.8000E+01 2.501E-04 2.501E-04 1.591E-03 1.591E-03 + 108 5.3400E+01 2.423E-04 2.423E-04 2.300E-03 2.300E-03 + 109 5.0600E+01 2.374E-04 2.374E-04 2.702E-03 2.702E-03 + 110 4.8300E+01 2.387E-04 2.387E-04 1.873E-03 1.873E-03 + 111 4.5200E+01 2.237E-04 2.237E-04 4.336E-03 4.336E-03 + 112 4.4000E+01 2.372E-04 2.372E-04 3.251E-03 3.251E-03 + 113 4.2400E+01 2.192E-04 2.192E-04 3.560E-03 3.560E-03 + 114 4.1000E+01 2.302E-04 2.302E-04 3.482E-03 3.482E-03 + 115 3.9600E+01 2.248E-04 2.248E-04 9.646E-03 9.646E-03 + 116 3.9100E+01 2.192E-04 2.192E-04 4.119E-03 4.119E-03 + 117 3.8000E+01 2.251E-04 2.251E-04 1.227E-02 1.227E-02 + 118 3.7630E+01 2.199E-04 2.199E-04 1.247E-02 1.247E-02 + 119 3.7270E+01 2.110E-04 2.110E-04 3.140E-02 3.140E-02 + 120 3.7130E+01 2.022E-04 2.022E-04 3.316E-02 3.316E-02 + 121 3.7000E+01 2.031E-04 2.031E-04 4.219E-03 4.219E-03 + 122 3.6000E+01 2.212E-04 2.212E-04 8.825E-03 8.825E-03 + 123 3.5500E+01 2.072E-04 2.072E-04 9.033E-03 9.033E-03 + 124 3.5000E+01 2.108E-04 2.108E-04 3.456E-03 3.456E-03 + 125 3.3750E+01 2.089E-04 2.089E-04 8.328E-03 8.328E-03 + 126 3.3250E+01 2.093E-04 2.093E-04 2.723E-03 2.723E-03 + 127 3.1750E+01 2.126E-04 2.126E-04 7.702E-03 7.702E-03 + 128 3.1250E+01 2.122E-04 2.122E-04 2.999E-03 2.999E-03 + 129 3.0000E+01 2.131E-04 2.131E-04 1.414E-03 1.414E-03 + 130 2.7500E+01 1.987E-04 1.987E-04 1.310E-03 1.310E-03 + 131 2.5000E+01 1.913E-04 1.913E-04 1.161E-03 1.161E-03 + 132 2.2500E+01 1.702E-04 1.702E-04 3.739E-03 3.739E-03 + 133 2.1750E+01 1.775E-04 1.775E-04 4.400E-03 4.400E-03 + 134 2.1200E+01 1.590E-04 1.590E-04 3.192E-03 3.192E-03 + 135 2.0500E+01 1.681E-04 1.681E-04 4.660E-03 4.660E-03 + 136 2.0000E+01 1.679E-04 1.679E-04 3.914E-03 3.914E-03 + 137 1.9400E+01 1.638E-04 1.638E-04 2.521E-03 2.521E-03 + 138 1.8500E+01 1.577E-04 1.577E-04 1.416E-03 1.416E-03 + 139 1.7000E+01 1.564E-04 1.564E-04 1.918E-03 1.918E-03 + 140 1.6000E+01 1.441E-04 1.441E-04 1.126E-03 1.126E-03 + 141 1.4400E+01 1.387E-04 1.387E-04 9.724E-04 9.724E-04 + 142 1.2900E+01 1.336E-04 1.336E-04 1.249E-03 1.249E-03 + 143 1.1900E+01 1.250E-04 1.250E-04 3.015E-03 3.015E-03 + 144 1.1500E+01 1.211E-04 1.211E-04 7.249E-04 7.249E-04 + 145 1.0000E+01 1.209E-04 1.209E-04 9.269E-04 9.269E-04 + 146 9.1000E+00 1.157E-04 1.157E-04 7.469E-04 7.469E-04 + 147 8.1000E+00 1.052E-04 1.052E-04 6.868E-04 6.868E-04 + 148 7.1500E+00 1.034E-04 1.034E-04 3.566E-03 3.566E-03 + 149 7.0000E+00 1.019E-04 1.019E-04 4.091E-03 4.091E-03 + 150 6.8750E+00 9.662E-05 9.662E-05 3.843E-03 3.843E-03 + 151 6.7500E+00 9.382E-05 9.382E-05 1.888E-03 1.888E-03 + 152 6.5000E+00 9.558E-05 9.558E-05 1.960E-03 1.960E-03 + 153 6.2500E+00 9.558E-05 9.558E-05 1.915E-03 1.915E-03 + 154 6.0000E+00 9.546E-05 9.546E-05 7.433E-04 7.433E-04 + 155 5.4000E+00 9.434E-05 9.434E-05 1.019E-03 1.019E-03 + 156 5.0000E+00 9.521E-05 9.521E-05 1.264E-03 1.264E-03 + 157 4.7000E+00 9.114E-05 9.114E-05 5.779E-04 5.779E-04 + 158 4.1000E+00 8.608E-05 8.608E-05 8.370E-04 8.370E-04 + 159 3.7300E+00 8.089E-05 8.089E-05 1.267E-03 1.267E-03 + 160 3.5000E+00 7.937E-05 7.937E-05 8.798E-04 8.798E-04 + 161 3.2000E+00 7.698E-05 7.698E-05 2.490E-03 2.490E-03 + 162 3.1000E+00 7.779E-05 7.779E-05 2.390E-03 2.390E-03 + 163 3.0000E+00 7.938E-05 7.938E-05 7.734E-03 7.734E-03 + 164 2.9700E+00 7.552E-05 7.552E-05 2.253E-03 2.253E-03 + 165 2.8700E+00 7.401E-05 7.401E-05 2.161E-03 2.161E-03 + 166 2.7700E+00 7.271E-05 7.271E-05 2.072E-03 2.072E-03 + 167 2.6700E+00 7.130E-05 7.130E-05 1.990E-03 1.990E-03 + 168 2.5700E+00 6.981E-05 6.981E-05 1.911E-03 1.911E-03 + 169 2.4700E+00 6.831E-05 6.831E-05 2.038E-03 2.038E-03 + 170 2.3800E+00 6.690E-05 6.690E-05 2.207E-03 2.207E-03 + 171 2.3000E+00 6.541E-05 6.541E-05 1.884E-03 1.884E-03 + 172 2.2100E+00 6.372E-05 6.372E-05 1.802E-03 1.802E-03 + 173 2.1200E+00 6.143E-05 6.143E-05 1.278E-03 1.278E-03 + 174 2.0000E+00 5.984E-05 5.984E-05 2.434E-03 2.434E-03 + 175 1.9400E+00 5.857E-05 5.857E-05 1.754E-03 1.754E-03 + 176 1.8600E+00 5.687E-05 5.687E-05 1.482E-03 1.482E-03 + 177 1.7700E+00 5.498E-05 5.498E-05 1.399E-03 1.399E-03 + 178 1.6800E+00 5.301E-05 5.301E-05 1.317E-03 1.317E-03 + 179 1.5900E+00 5.095E-05 5.095E-05 1.235E-03 1.235E-03 + 180 1.5000E+00 4.926E-05 4.926E-05 2.110E-03 2.110E-03 + 181 1.4500E+00 4.794E-05 4.794E-05 2.030E-03 2.030E-03 + 182 1.4000E+00 4.656E-05 4.656E-05 1.951E-03 1.951E-03 + 183 1.3500E+00 4.509E-05 4.509E-05 1.872E-03 1.872E-03 + 184 1.3000E+00 4.351E-05 4.351E-05 1.793E-03 1.793E-03 + 185 1.2500E+00 4.224E-05 4.224E-05 3.472E-03 3.472E-03 + 186 1.2250E+00 4.133E-05 4.133E-05 3.394E-03 3.394E-03 + 187 1.2000E+00 4.041E-05 4.041E-05 3.314E-03 3.314E-03 + 188 1.1750E+00 3.947E-05 3.947E-05 3.228E-03 3.228E-03 + 189 1.1500E+00 3.881E-05 3.881E-05 7.912E-03 7.912E-03 + 190 1.1400E+00 3.844E-05 3.844E-05 7.825E-03 7.825E-03 + 191 1.1300E+00 3.809E-05 3.809E-05 7.748E-03 7.748E-03 + 192 1.1200E+00 3.777E-05 3.777E-05 7.684E-03 7.684E-03 + 193 1.1100E+00 3.748E-05 3.748E-05 7.634E-03 7.634E-03 + 194 1.1000E+00 3.721E-05 3.721E-05 7.595E-03 7.595E-03 + 195 1.0900E+00 3.695E-05 3.695E-05 7.566E-03 7.566E-03 + 196 1.0800E+00 3.671E-05 3.671E-05 7.543E-03 7.543E-03 + 197 1.0700E+00 3.647E-05 3.647E-05 7.524E-03 7.524E-03 + 198 1.0600E+00 3.624E-05 3.624E-05 7.507E-03 7.507E-03 + 199 1.0500E+00 3.601E-05 3.601E-05 7.491E-03 7.491E-03 + 200 1.0400E+00 3.579E-05 3.579E-05 7.476E-03 7.476E-03 + 201 1.0300E+00 3.557E-05 3.557E-05 7.463E-03 7.463E-03 + 202 1.0200E+00 3.536E-05 3.536E-05 7.452E-03 7.452E-03 + 203 1.0100E+00 3.516E-05 3.516E-05 7.444E-03 7.444E-03 + 204 1.0000E+00 3.484E-05 3.484E-05 2.977E-03 2.977E-03 + 205 9.7500E-01 3.442E-05 3.442E-05 2.978E-03 2.978E-03 + 206 9.5000E-01 3.400E-05 3.400E-05 2.970E-03 2.970E-03 + 207 9.2500E-01 3.353E-05 3.353E-05 2.948E-03 2.948E-03 + 208 9.0000E-01 3.275E-05 3.275E-05 1.446E-03 1.446E-03 + 209 8.5000E-01 3.162E-05 3.162E-05 1.396E-03 1.396E-03 + 210 8.0000E-01 3.039E-05 3.039E-05 1.335E-03 1.335E-03 + 211 7.5000E-01 2.906E-05 2.906E-05 1.264E-03 1.264E-03 + 212 7.0000E-01 2.763E-05 2.763E-05 1.185E-03 1.185E-03 + 213 6.5000E-01 2.652E-05 2.652E-05 2.245E-03 2.245E-03 + 214 6.2500E-01 2.570E-05 2.570E-05 2.157E-03 2.157E-03 + 215 6.0000E-01 2.443E-05 2.443E-05 1.009E-03 1.009E-03 + 216 5.5000E-01 2.271E-05 2.271E-05 9.125E-04 9.125E-04 + 217 5.0000E-01 2.090E-05 2.090E-05 8.118E-04 8.118E-04 + 218 4.5000E-01 1.906E-05 1.906E-05 7.075E-04 7.075E-04 + 219 4.0000E-01 1.770E-05 1.770E-05 1.256E-03 1.256E-03 + 220 3.7500E-01 1.673E-05 1.673E-05 1.147E-03 1.147E-03 + 221 3.5000E-01 1.573E-05 1.573E-05 1.040E-03 1.040E-03 + 222 3.2500E-01 1.471E-05 1.471E-05 9.335E-04 9.335E-04 + 223 3.0000E-01 1.365E-05 1.365E-05 8.303E-04 8.303E-04 + 224 2.7500E-01 1.257E-05 1.257E-05 7.307E-04 7.307E-04 + 225 2.5000E-01 1.144E-05 1.144E-05 6.352E-04 6.352E-04 + 226 2.2500E-01 1.026E-05 1.026E-05 5.443E-04 5.443E-04 + 227 2.0000E-01 9.011E-06 9.011E-06 4.584E-04 4.584E-04 + 228 1.7500E-01 7.688E-06 7.688E-06 3.776E-04 3.776E-04 + 229 1.5000E-01 6.282E-06 6.282E-06 3.018E-04 3.018E-04 + 230 1.2500E-01 4.941E-06 4.941E-06 2.310E-04 2.310E-04 + 231 1.0000E-01 4.075E-06 4.075E-06 4.669E-04 4.669E-04 + 232 9.0000E-02 3.485E-06 3.485E-06 4.020E-04 4.020E-04 + 233 8.0000E-02 2.893E-06 2.893E-06 3.395E-04 3.395E-04 + 234 7.0000E-02 2.310E-06 2.310E-06 2.798E-04 2.798E-04 + 235 6.0000E-02 1.750E-06 1.750E-06 2.232E-04 2.232E-04 + 236 5.0000E-02 1.229E-06 1.229E-06 1.701E-04 1.701E-04 + 237 4.0000E-02 7.694E-07 7.694E-07 1.209E-04 1.209E-04 + 238 3.0000E-02 4.885E-07 4.885E-07 1.879E-04 1.879E-04 + 239 2.5300E-02 1.867E-07 1.867E-07 3.013E-05 3.013E-05 + 240 1.0000E-02 4.799E-08 4.799E-08 8.010E-05 8.010E-05 + 241 7.5000E-03 2.319E-08 2.319E-08 5.114E-05 5.114E-05 + 242 5.0000E-03 1.155E-08 1.155E-08 8.438E-05 8.438E-05 + 243 4.0000E-03 6.726E-09 6.726E-09 6.125E-05 6.125E-05 + 244 3.0000E-03 4.033E-09 4.033E-09 9.117E-05 9.117E-05 + 245 2.5000E-03 2.630E-09 2.630E-09 7.134E-05 7.134E-05 + 246 2.0000E-03 1.541E-09 1.541E-09 5.276E-05 5.276E-05 + 247 1.5000E-03 8.953E-10 8.953E-10 6.518E-05 6.518E-05 + 248 1.2000E-03 5.841E-10 5.841E-10 7.753E-05 7.753E-05 + 249 1.0000E-03 3.599E-10 3.599E-10 4.783E-05 4.783E-05 + 250 7.5000E-04 1.770E-10 1.770E-10 3.286E-05 3.286E-05 + 251 5.0000E-04 3.222E-11 3.222E-11 7.860E-06 7.860E-06 + 252 1.0000E-04 5.613E-13 5.613E-13 5.157E-06 5.157E-06 + + NOTE: Flux values above are per unit-lethargy. + + + + Zone data + + Zone Mix Relative Relative Zone + No. No. Flux Fission Power Volume + ----------------------------------------------------------------- + 1 1 1.87787E-01 7.98654E-01 1.40000E+03 + 2 3 1.40417E-01 0.00000E+00 4.00000E+01 + 3 2 1.24511E-01 0.00000E+00 2.20000E+02 + 4 3 9.27423E-02 0.00000E+00 4.00000E+01 + 5 10 4.09596E-02 2.01346E-01 7.00000E+02 + 6 3 1.12383E-01 0.00000E+00 3.39115E+01 + 7 4 1.00545E-01 0.00000E+00 4.40849E+02 + 8 3 7.01993E-02 0.00000E+00 1.35646E+01 + 9 2 6.09558E-02 0.00000E+00 7.46052E+01 + 10 3 4.77119E-02 0.00000E+00 1.35646E+01 + 11 4 2.17874E-02 0.00000E+00 2.37380E+02 + + + +module xsproc is finished. completion code 0. +Sequence T-Newt is finished. completion code 0. + + SCALE is finished on Wed Mar 08 15:56:49 2017. +-------------------------- Summary -------------------------- +t-newt finished. used 59.7199 seconds. +SCALE driver required a maximum of 65.7 MiB of RAM. +------------------------ End Summary ------------------------ + diff --git a/param_tester.txt b/testing/param_tester.txt similarity index 77% rename from param_tester.txt rename to testing/param_tester.txt index 909f5b0..e7eb3bd 100644 --- a/param_tester.txt +++ b/testing/param_tester.txt @@ -1,3 +1,3 @@ var_char $ iter_var del_z 10 5.0 15.0 -iter_lim 5 \ No newline at end of file +iter_lim 50 \ No newline at end of file From e2c7aff501ba1e0d04a9e07f1f77354d6065b6b5 Mon Sep 17 00:00:00 2001 From: drewejohnson <1drew.e.johnson@gmail.com> Date: Mon, 13 Mar 2017 12:35:52 -0400 Subject: [PATCH 2/3] Display progression of iter_vars in output By use of pandas.DataFrame, outputs.py now displays the value of each iter_var over time, and the eigenvalue at each iteration at the end of the runs. Signed-off-by: drewejohnson <1drew.e.johnson@gmail.com> --- globalparams.py | 3 +- iterator.py | 8 +- outputs.py | 7 +- testing/iter_tester_16.inp | 112 ----- testing/iter_tester_16.out | 911 ------------------------------------- testing/param_tester.txt | 3 +- 6 files changed, 12 insertions(+), 1032 deletions(-) delete mode 100644 testing/iter_tester_16.inp delete mode 100644 testing/iter_tester_16.out diff --git a/globalparams.py b/globalparams.py index 70717b8..a1dab4a 100644 --- a/globalparams.py +++ b/globalparams.py @@ -55,13 +55,12 @@ # Iteration parameters iter_vars = {} # syntax: iter_vars['param'] = [initial_value, minimum_value, maximium_value] eps_k = 1E-6 -k_guess = 1.0 k_target = 1.0 inf = 1E30 tiny = 1E-16 iter_lim = 50 iter_ints = ('iter_lim',) -iter_floats = ('eps_k', 'inf', 'tiny', 'k_guess', "k_target") +iter_floats = ('eps_k', 'inf', 'tiny', "k_target") iter_vecs = {} k_vec = [] diff --git a/iterator.py b/iterator.py index be9e8b6..d69a7d3 100644 --- a/iterator.py +++ b/iterator.py @@ -81,9 +81,9 @@ def itermain(): for _var in gp.iter_vars.keys(): gp.iter_vecs[_var] = [gp.iter_vars[_var][0], ] - gp.k_vec.append(gp.k_guess) + # gp.k_vec.append(gp.k_guess) - print("Starting the iteration procedure....\n") + utils.oprint("Starting the iteration procedure....\n") conv_flag = False conv_type = None @@ -96,7 +96,7 @@ def itermain(): utils.vprint(' done') stat, _k = parse_scale_out_eig(_iter_file.replace('.inp', '.out')) if stat: # successful operation - print(" {0:<3}: {1}".format(_n, _k)) + utils.oprint(" {0:<3}: {1}".format(_n, _k)) gp.k_vec.append(_k) else: utils.error('Could not find value of k-eff for iteration file {0}.inp\n' @@ -120,6 +120,6 @@ def itermain(): if _n == gp.iter_lim and conv_type is None: conv_type = 2 - print(' done') + utils.oprint(' done') output_landing(conv_type) diff --git a/outputs.py b/outputs.py index ae330ad..a1b1c7f 100644 --- a/outputs.py +++ b/outputs.py @@ -14,6 +14,7 @@ Classes: """ +import pandas as pd import globalparams as gp import utils @@ -64,5 +65,7 @@ def output_landing(_out_type): list(gp.iter_vars.keys())[0], gp.iter_vars[list(gp.iter_vars.keys())[0]][2], ), 2: "Terminated due to exceeding the iteration limit {}".format(gp.iter_lim) } - print('End of operation. Status: {}'.format(out_messages[_out_type])) - + utils.oprint('End of operation. Status: {}'.format(out_messages[_out_type])) + var_df = pd.DataFrame({_var: gp.iter_vecs[_var] for _var in gp.iter_vars}) + var_df['k-eff'] = gp.k_vec + utils.oprint(var_df.to_string()) diff --git a/testing/iter_tester_16.inp b/testing/iter_tester_16.inp deleted file mode 100644 index 7945334..0000000 --- a/testing/iter_tester_16.inp +++ /dev/null @@ -1,112 +0,0 @@ -=t-newt parm=(centrm) -MSR iterative reflector optimization procedure -v7-252 - -read comp - -' --- NaCl-UCl3-PuCl3' -' --- 50/25/25 dens=3.60' -' --- 19.75 wt% eniched uranium' - Na-23 1 0 5.3767E-03 900 end - Cl-35 1 0 1.5991E-02 900 end - Cl-37 1 0 5.5155E-03 900 end - U-235 1 0 5.2556E-04 900 end - U-238 1 0 2.1355E-03 900 end - Pu-239 1 0 2.4994E-03 900 end - Pu-240 1 0 1.8892E-04 900 end - -' --- NaCl-UCl3-PuCl3' -' --- 50/25/25 dens=3.60' -' --- 19.75 wt% eniched uranium' - Na-23 10 0 5.3767E-03 900 end - Cl-35 10 0 1.5991E-02 900 end - Cl-37 10 0 5.5155E-03 900 end - U-235 10 0 5.2556E-04 900 end - U-238 10 0 2.1355E-03 900 end - Pu-239 10 0 2.4994E-03 900 end - Pu-240 10 0 1.8892E-04 900 end -' --- Heat exchanger fluid -' -- Nitrogen gas a.t.m. - n-14 2 1.0 293.0 end -' --- Piping -' -- Stainless steel - ss304 3 1.0 end -' --- Piping -' -- Stainless steel - ss304 30 1.0 end -' -- Reflector - ss304 4 1.0 600.0 end - -end comp - - -read celldata - - multiregion slab right_bdy=vacuum left_bdy=reflected origin=0.0 end - 1 55 - 3 60 - 2 70 - 30 75 - 10 100 - 4 120 end zone - -end celldata - -read model -Preliminary model - -read parameters - -drawit=no -sn=6 -echo=yes -prtflux=no - -end parameters - -read materials - - mix=1 com="in fuel" end - mix=10 com="out fuel" end - mix=2 pn=2 com="HX" end - mix=3 com="pipe" end - mix=4 com="refl" pn=2 end - -end materials - -read geometry - -global unit 1 - cuboid 99 120 0 20 -6.78229 - cuboid 10 70 0.0 20 0.0 - cuboid 11 72 70 20 0 - cuboid 12 83 72 20 0 - cuboid 13 85 83 20 0 - cuboid 14 120 85 20 0 - media 1 1 10 - media 3 1 11 - media 2 1 12 - media 3 1 13 - media 10 1 14 - cuboid 20 5 0 0 -6.78229 - cuboid 21 70 5 0 -6.78229 - cuboid 22 72 70 0 -6.78229 - cuboid 23 83 72 0 -6.78229 - cuboid 24 85 83 0 -6.78229 - cuboid 25 120 85 0 -6.78229 - media 3 1 20 - media 4 1 21 - media 3 1 22 - media 2 1 23 - media 3 1 24 - media 4 1 25 - boundary 99 4 4 -end geometry - -read bounds - +x=vac +y=refl -x=refl -y=vac -end bounds -end model - - -end diff --git a/testing/iter_tester_16.out b/testing/iter_tester_16.out deleted file mode 100644 index 6cf7e2e..0000000 --- a/testing/iter_tester_16.out +++ /dev/null @@ -1,911 +0,0 @@ - ******************************************************************************* - * * - * SCALE 6.2.1 * - * ------------- * - * March 2017 * - * * - * SCALE: A Comprehensive Modeling and Simulation Suite * - * for Nuclear Safety Analysis and Design * - * * - * Reactor and Nuclear Systems Division * - * Oak Ridge National Laboratory * - * * - * http://scale.ornl.gov * - * scalehelp@ornl.gov * - * * - ******************************************************************************* - ******************************************************************************* - - Job Information - --------------- - Job started on ME04W0358GRD04 on Wed 08/03/2017 15:55:49 - Working directory: C:\Users\AJOHNS~1\AppData\Local\Temp\scale.ajohnson400.12048 - Input file name : C:\Users\ajohnson400\Documents\classes\msr-refl-iter\testing\iter_tester_16.inp - Output file name : C:\Users\ajohnson400\Documents\classes\msr-refl-iter\testing\iter_tester_16.out - SCALE executable : C:\SCALE-6.2.1/bin/scale - - ******************************************************************************* - - - module t-newt will be called on Wed Mar 08 15:55:49 2017. - sequence specification record:=t-newt parm=(centrm) -Input Data: - - MSR iterative reflector optimization procedure - v7-252 - - read comp - - ' --- NaCl-UCl3-PuCl3' - ' --- 50/25/25 dens=3.60' - ' --- 19.75 wt% eniched uranium' - Na-23 1 0 5.3767E-03 900 end - Cl-35 1 0 1.5991E-02 900 end - Cl-37 1 0 5.5155E-03 900 end - U-235 1 0 5.2556E-04 900 end - U-238 1 0 2.1355E-03 900 end - Pu-239 1 0 2.4994E-03 900 end - Pu-240 1 0 1.8892E-04 900 end - - ' --- NaCl-UCl3-PuCl3' - ' --- 50/25/25 dens=3.60' - ' --- 19.75 wt% eniched uranium' - Na-23 10 0 5.3767E-03 900 end - Cl-35 10 0 1.5991E-02 900 end - Cl-37 10 0 5.5155E-03 900 end - U-235 10 0 5.2556E-04 900 end - U-238 10 0 2.1355E-03 900 end - Pu-239 10 0 2.4994E-03 900 end - Pu-240 10 0 1.8892E-04 900 end - ' --- Heat exchanger fluid - ' -- Nitrogen gas a.t.m. - n-14 2 1.0 293.0 end - ' --- Piping - ' -- Stainless steel - ss304 3 1.0 end - ' --- Piping - ' -- Stainless steel - ss304 30 1.0 end - ' -- Reflector - ss304 4 1.0 600.0 end - - end comp - - - read celldata - - multiregion slab right_bdy=vacuum left_bdy=reflected origin=0.0 end - 1 55 - 3 60 - 2 70 - 30 75 - 10 100 - 4 120 end zone - - end celldata - - read model - Preliminary model - - read parameters - - drawit=no - sn=6 - echo=yes - prtflux=no - - end parameters - - read materials - - mix=1 com="in fuel" end - mix=10 com="out fuel" end - mix=2 pn=2 com="HX" end - mix=3 com="pipe" end - mix=4 com="refl" pn=2 end - - end materials - - read geometry - - global unit 1 - cuboid 99 120 0 20 -6.78229 - cuboid 10 70 0.0 20 0.0 - cuboid 11 72 70 20 0 - cuboid 12 83 72 20 0 - cuboid 13 85 83 20 0 - cuboid 14 120 85 20 0 - media 1 1 10 - media 3 1 11 - media 2 1 12 - media 3 1 13 - media 10 1 14 - cuboid 20 5 0 0 -6.78229 - cuboid 21 70 5 0 -6.78229 - cuboid 22 72 70 0 -6.78229 - cuboid 23 83 72 0 -6.78229 - cuboid 24 85 83 0 -6.78229 - cuboid 25 120 85 0 -6.78229 - media 3 1 20 - media 4 1 21 - media 3 1 22 - media 2 1 23 - media 3 1 24 - media 4 1 25 - boundary 99 4 4 - end geometry - - read bounds - +x=vac +y=refl -x=refl -y=vac - end bounds - end model - - -end - - - - -************************************************************************************************************************************ -************************************************************************************************************************************ -************************************************************************************************************************************ -***** ***** -***** program verification information ***** -***** ***** -***** code system: SCALE version: 6.2 ***** -***** ***** -************************************************************************************************************************************ -************************************************************************************************************************************ -***** ***** -***** ***** -***** program: t-newt ***** -***** ***** -***** version: 6.2.1 (602001) ***** -***** ***** -***** username: ajohnson400 ***** -***** ***** -***** hostname: ME04W0358GRD04 ***** -***** ***** -***** ***** -***** date of execution: 2017-03-08 ***** -***** ***** -***** time of execution: 15:55:49 ***** -***** ***** -************************************************************************************************************************************ -************************************************************************************************************************************ -************************************************************************************************************************************ - - - ------------------------------------------------------------------------------------------------------------------------------------- - - preliminary model - ------------------------------------------------------------------------------------------------------------------------------------- - - - -1 *********************** - *********************** - ** Control Options ** - *********************** - *********************** - - Run problem to completion/test input only run= yes - Weighted cross-section collapse/do not collapse collapse= no - Source acceleration on/off accel= yes - Restart from earlier case/begin new case restart= no - CMFD acceleration (0/1/2/3 = no/rect/yes/part) cmfd= 0 - Number of fine mesh base grid/coarse grid in X xcmfd= 1 - Number of fine mesh base grid/coarse grid in Y ycmfd= 1 - CMFD Multigroup solution only (no 2g) cmfd2g= no - Create restart file upon completion savrest= yes - Calculation type: eigenvalue - Forward/adjoint solution forward= yes - -1 ********************** - ********************** - ** Output Options ** - ********************** - ********************** - - Print micro. cross-sections for problem nuclides prtxsec= no - Print macro. cross-sections for mixtures prtmxsec= no - Print balance table prtbalnc= no - Print mixing table prtmxtab= no - Print postscript graphics file for fluxes prtflux= no - Print postscript grid/mixture geometry graphics drawit= no - Print broad-group fluxes prtbroad= no - Echo convergence tracking data to std. output echo= yes - -1 *************************** - *************************** - ** I/O Unit Assignments ** - *************************** - *************************** - - Unit no. for mixing table mixtab= 92 - Unit no. of working library xnlib= 4 - Unit no. for weighted broad-group library wtdlib= 30 - Unit no. for homogenized macroscopic x-sect. lib. hmoglib= 13 - -1 ************************* - ************************* - ** Convergence Control ** - ************************* - ************************* - - Quadrature order sn= 6 - Maximum no. of inner itrs./outer/energy group inners= 5 - Maximum no. of outer energy groups outers= 250 - Performing thermal group subouter iterations therm= yes - Maximum no. of thermal group subouter itrs./outer suboutrs= 2 - Convergence criterion on inner iterations epsinner= 1.00E-04 - Convergence criterion on subouter iterations epsthrm= 1.00E-04 - Convergence criterion on outer iterations epsouter= 1.00E-04 - Convergence criterion on eigenvalue epseigen= 1.00E-04 - First guess for eigenvalue (k-eff) kguess= 1.000 - -1 ***************************** - ***************************** - ** Geometry Specifications ** - ***************************** - ***************************** - - Left side x-coordinate xmin= 0.00000E+00 - Right side x-coordinate xmax= 1.20000E+02 - Bottom side y-coordinate ymin= -6.78229E+00 - Top side y-coordinate ymax= 2.00000E+01 - Number of subdivisions in x direction nx= 4 - Number of subdivisions in y direction ny= 4 - - Boundary Conditions (0/1/2/3=vacuum/reflective/white/periodic): - Left (West) wbc= 1 - Right (East) ebc= 0 - Bottom (South) sbc= 0 - Top (North) nbc= 1 - - Body Characteristics. Rotation angle in degrees counterclockwise. - Size parameters. cyl: A=radius, B=no. of sides. cub: A=width, B=height. hex: A(B)=major(minor) radius. - Body Input Zone Subdiv --Center-- Size Parameters Rotation Chord - Type Unit Number nx ny x y A B Angle Information - ---- ----- ----- --- --- ------ ------ ----- ------ -------- ----------- - cub 1 1 4 4 60.000 6.609 120.000 26.782 0.000 - cub 1 2 0 0 35.000 10.000 70.000 20.000 0.000 - cub 1 3 0 0 71.000 10.000 2.000 20.000 0.000 - cub 1 4 0 0 77.500 10.000 11.000 20.000 0.000 - cub 1 5 0 0 84.000 10.000 2.000 20.000 0.000 - cub 1 6 0 0 102.500 10.000 35.000 20.000 0.000 - cub 1 7 0 0 2.500 -3.391 5.000 6.782 0.000 - cub 1 8 0 0 37.500 -3.391 65.000 6.782 0.000 - cub 1 9 0 0 71.000 -3.391 2.000 6.782 0.000 - cub 1 10 0 0 77.500 -3.391 11.000 6.782 0.000 - cub 1 11 0 0 84.000 -3.391 2.000 6.782 0.000 - cub 1 12 0 0 102.500 -3.391 35.000 6.782 0.000 - -1 ***************************** - ***************************** - ** Material Specifications ** - ***************************** - ***************************** - - Mat'l No. SCALE ID Pn order External source ID Description (if any) - 1 1 1 0 in fuel - 2 10 1 0 out fuel - 3 2 2 0 hx - 4 3 1 0 pipe - 5 4 2 0 refl - -1 ************************ - ************************ - ** Derived Parameters ** - ************************ - ************************ - - Total number of energy groups numen= 252 - Number of neutron groups igm= 252 - Number of gamma groups ipm= 0 - First thermal energy group iftg= 156 - Mixing table length nmix= 53 - Number of mixtures created nummat= 6 - Number of computational cells numcells= 42 - Number of computational surfaces (sides) numlines= 98 - Combine small grid cells together combine=no - Average number of sides/cell avnosides= 4.02 - Minimum number of sides/cell mnnosides= 4 - Maximum number of sides/cell mxnosides= 5 - Number of angles in quadrature set ndir= 24 - Maximum order of Legendre scattering maxpn= 2 - Number of flux moments nmom 6 - -1 ***************************** - ***************************** - ** Energy Group Structures ** - ***************************** - ***************************** - - Energy Energy Lethargy - Group Boundaries Boundaries - No. (eV) - 1 2.00000E+07 -6.93147E-01 127 3.17500E+01 1.26602E+01 - 2 1.73300E+07 -5.49854E-01 128 3.12500E+01 1.26761E+01 - 3 1.56800E+07 -4.49801E-01 129 3.00000E+01 1.27169E+01 - 4 1.45500E+07 -3.75006E-01 130 2.75000E+01 1.28039E+01 - 5 1.38400E+07 -3.24978E-01 131 2.50000E+01 1.28992E+01 - 6 1.28400E+07 -2.49980E-01 132 2.25000E+01 1.30046E+01 - 7 1.00000E+07 0.00000E+00 133 2.17500E+01 1.30385E+01 - 8 8.18700E+06 2.00038E-01 134 2.12000E+01 1.30641E+01 - 9 6.43400E+06 4.40989E-01 135 2.05000E+01 1.30977E+01 - 10 4.80000E+06 7.33969E-01 136 2.00000E+01 1.31224E+01 - 11 4.30400E+06 8.43040E-01 137 1.94000E+01 1.31528E+01 - 12 3.00000E+06 1.20397E+00 138 1.85000E+01 1.32003E+01 - 13 2.47900E+06 1.39473E+00 139 1.70000E+01 1.32849E+01 - 14 2.35400E+06 1.44647E+00 140 1.60000E+01 1.33455E+01 - 15 1.85000E+06 1.68740E+00 141 1.44000E+01 1.34509E+01 - 16 1.50000E+06 1.89712E+00 142 1.29000E+01 1.35609E+01 - 17 1.40000E+06 1.96611E+00 143 1.19000E+01 1.36416E+01 - 18 1.35600E+06 1.99805E+00 144 1.15000E+01 1.36757E+01 - 19 1.31700E+06 2.02723E+00 145 1.00000E+01 1.38155E+01 - 20 1.25000E+06 2.07944E+00 146 9.10000E+00 1.39098E+01 - 21 1.20000E+06 2.12026E+00 147 8.10000E+00 1.40262E+01 - 22 1.10000E+06 2.20727E+00 148 7.15000E+00 1.41510E+01 - 23 1.01000E+06 2.29263E+00 149 7.00000E+00 1.41722E+01 - 24 9.20000E+05 2.38597E+00 150 6.87500E+00 1.41902E+01 - 25 9.00000E+05 2.40795E+00 151 6.75000E+00 1.42086E+01 - 26 8.75000E+05 2.43612E+00 152 6.50000E+00 1.42463E+01 - 27 8.61100E+05 2.45213E+00 153 6.25000E+00 1.42855E+01 - 28 8.20000E+05 2.50104E+00 154 6.00000E+00 1.43263E+01 - 29 7.50000E+05 2.59027E+00 155 5.40000E+00 1.44317E+01 - 30 6.79000E+05 2.68972E+00 156 5.00000E+00 1.45087E+01 - 31 6.70000E+05 2.70306E+00 157 4.70000E+00 1.45705E+01 - 32 6.00000E+05 2.81341E+00 158 4.10000E+00 1.47071E+01 - 33 5.73000E+05 2.85945E+00 159 3.73000E+00 1.48017E+01 - 34 5.50000E+05 2.90042E+00 160 3.50000E+00 1.48653E+01 - 35 4.92000E+05 3.01186E+00 161 3.20000E+00 1.49549E+01 - 36 4.70000E+05 3.05761E+00 162 3.10000E+00 1.49867E+01 - 37 4.40000E+05 3.12357E+00 163 3.00000E+00 1.50195E+01 - 38 4.20000E+05 3.17009E+00 164 2.97000E+00 1.50295E+01 - 39 4.00000E+05 3.21888E+00 165 2.87000E+00 1.50638E+01 - 40 3.30000E+05 3.41125E+00 166 2.77000E+00 1.50992E+01 - 41 2.70000E+05 3.61192E+00 167 2.67000E+00 1.51360E+01 - 42 2.00000E+05 3.91202E+00 168 2.57000E+00 1.51742E+01 - 43 1.49000E+05 4.20639E+00 169 2.47000E+00 1.52139E+01 - 44 1.28300E+05 4.35597E+00 170 2.38000E+00 1.52510E+01 - 45 1.00000E+05 4.60517E+00 171 2.30000E+00 1.52852E+01 - 46 8.50000E+04 4.76769E+00 172 2.21000E+00 1.53251E+01 - 47 8.20000E+04 4.80362E+00 173 2.12000E+00 1.53667E+01 - 48 7.50000E+04 4.89285E+00 174 2.00000E+00 1.54249E+01 - 49 7.30000E+04 4.91988E+00 175 1.94000E+00 1.54554E+01 - 50 6.00000E+04 5.11600E+00 176 1.86000E+00 1.54975E+01 - 51 5.20000E+04 5.25910E+00 177 1.77000E+00 1.55471E+01 - 52 5.00000E+04 5.29832E+00 178 1.68000E+00 1.55993E+01 - 53 4.50000E+04 5.40368E+00 179 1.59000E+00 1.56544E+01 - 54 3.00000E+04 5.80914E+00 180 1.50000E+00 1.57126E+01 - 55 2.00000E+04 6.21461E+00 181 1.45000E+00 1.57465E+01 - 56 1.70000E+04 6.37713E+00 182 1.40000E+00 1.57816E+01 - 57 1.30000E+04 6.64539E+00 183 1.35000E+00 1.58180E+01 - 58 9.50000E+03 6.95905E+00 184 1.30000E+00 1.58557E+01 - 59 8.03000E+03 7.12716E+00 185 1.25000E+00 1.58950E+01 - 60 5.70000E+03 7.46987E+00 186 1.22500E+00 1.59152E+01 - 61 3.90000E+03 7.84936E+00 187 1.20000E+00 1.59358E+01 - 62 3.74000E+03 7.89125E+00 188 1.17500E+00 1.59568E+01 - 63 3.00000E+03 8.11173E+00 189 1.15000E+00 1.59783E+01 - 64 2.50000E+03 8.29405E+00 190 1.14000E+00 1.59871E+01 - 65 2.25000E+03 8.39941E+00 191 1.13000E+00 1.59959E+01 - 66 2.20000E+03 8.42188E+00 192 1.12000E+00 1.60048E+01 - 67 1.80000E+03 8.62255E+00 193 1.11000E+00 1.60137E+01 - 68 1.55000E+03 8.77209E+00 194 1.10000E+00 1.60228E+01 - 69 1.50000E+03 8.80488E+00 195 1.09000E+00 1.60319E+01 - 70 1.15000E+03 9.07058E+00 196 1.08000E+00 1.60411E+01 - 71 9.50000E+02 9.26163E+00 197 1.07000E+00 1.60504E+01 - 72 6.83000E+02 9.59160E+00 198 1.06000E+00 1.60598E+01 - 73 6.70000E+02 9.61082E+00 199 1.05000E+00 1.60693E+01 - 74 5.50000E+02 9.80818E+00 200 1.04000E+00 1.60789E+01 - 75 3.05000E+02 1.03978E+01 201 1.03000E+00 1.60885E+01 - 76 2.85000E+02 1.04656E+01 202 1.02000E+00 1.60983E+01 - 77 2.40000E+02 1.06375E+01 203 1.01000E+00 1.61081E+01 - 78 2.20000E+02 1.07245E+01 204 1.00000E+00 1.61181E+01 - 79 2.09500E+02 1.07734E+01 205 9.75000E-01 1.61434E+01 - 80 2.07400E+02 1.07834E+01 206 9.50000E-01 1.61694E+01 - 81 2.02000E+02 1.08098E+01 207 9.25000E-01 1.61961E+01 - 82 1.93000E+02 1.08554E+01 208 9.00000E-01 1.62235E+01 - 83 1.91500E+02 1.08632E+01 209 8.50000E-01 1.62806E+01 - 84 1.88500E+02 1.08790E+01 210 8.00000E-01 1.63412E+01 - 85 1.87700E+02 1.08833E+01 211 7.50000E-01 1.64058E+01 - 86 1.80000E+02 1.09251E+01 212 7.00000E-01 1.64748E+01 - 87 1.70000E+02 1.09823E+01 213 6.50000E-01 1.65489E+01 - 88 1.43000E+02 1.11553E+01 214 6.25000E-01 1.65881E+01 - 89 1.22000E+02 1.13141E+01 215 6.00000E-01 1.66289E+01 - 90 1.19000E+02 1.13390E+01 216 5.50000E-01 1.67159E+01 - 91 1.17500E+02 1.13517E+01 217 5.00000E-01 1.68112E+01 - 92 1.16000E+02 1.13645E+01 218 4.50000E-01 1.69166E+01 - 93 1.13000E+02 1.13907E+01 219 4.00000E-01 1.70344E+01 - 94 1.08000E+02 1.14360E+01 220 3.75000E-01 1.70989E+01 - 95 1.05000E+02 1.14641E+01 221 3.50000E-01 1.71679E+01 - 96 1.01200E+02 1.15010E+01 222 3.25000E-01 1.72420E+01 - 97 9.70000E+01 1.15434E+01 223 3.00000E-01 1.73221E+01 - 98 9.00000E+01 1.16183E+01 224 2.75000E-01 1.74091E+01 - 99 8.17000E+01 1.17150E+01 225 2.50000E-01 1.75044E+01 - 100 8.00000E+01 1.17361E+01 226 2.25000E-01 1.76098E+01 - 101 7.60000E+01 1.17874E+01 227 2.00000E-01 1.77275E+01 - 102 7.20000E+01 1.18414E+01 228 1.75000E-01 1.78611E+01 - 103 6.75000E+01 1.19060E+01 229 1.50000E-01 1.80152E+01 - 104 6.50000E+01 1.19437E+01 230 1.25000E-01 1.81975E+01 - 105 6.30000E+01 1.19750E+01 231 1.00000E-01 1.84207E+01 - 106 6.10000E+01 1.20072E+01 232 9.00000E-02 1.85260E+01 - 107 5.80000E+01 1.20577E+01 233 8.00000E-02 1.86438E+01 - 108 5.34000E+01 1.21403E+01 234 7.00000E-02 1.87774E+01 - 109 5.06000E+01 1.21941E+01 235 6.00000E-02 1.89315E+01 - 110 4.83000E+01 1.22407E+01 236 5.00000E-02 1.91138E+01 - 111 4.52000E+01 1.23070E+01 237 4.00000E-02 1.93370E+01 - 112 4.40000E+01 1.23339E+01 238 3.00000E-02 1.96247E+01 - 113 4.24000E+01 1.23709E+01 239 2.53000E-02 1.97950E+01 - 114 4.10000E+01 1.24045E+01 240 1.00000E-02 2.07233E+01 - 115 3.96000E+01 1.24393E+01 241 7.50000E-03 2.10109E+01 - 116 3.91000E+01 1.24520E+01 242 5.00000E-03 2.14164E+01 - 117 3.80000E+01 1.24805E+01 243 4.00000E-03 2.16396E+01 - 118 3.76300E+01 1.24903E+01 244 3.00000E-03 2.19272E+01 - 119 3.72700E+01 1.24999E+01 245 2.50000E-03 2.21096E+01 - 120 3.71300E+01 1.25037E+01 246 2.00000E-03 2.23327E+01 - 121 3.70000E+01 1.25072E+01 247 1.50000E-03 2.26204E+01 - 122 3.60000E+01 1.25346E+01 248 1.20000E-03 2.28435E+01 - 123 3.55000E+01 1.25486E+01 249 1.00000E-03 2.30259E+01 - 124 3.50000E+01 1.25627E+01 250 7.50000E-04 2.33135E+01 - 125 3.37500E+01 1.25991E+01 251 5.00000E-04 2.37190E+01 - 126 3.32500E+01 1.26140E+01 252 1.00000E-04 2.53284E+01 - 253 1.00000E-05 2.76310E+01 - -1 *************************** - *************************** - ** Quadrature Parameters ** - *************************** - *************************** - - Quadrature angles and weights for first quadrant: - ( 6 angles/quadrant for S- 6 quadrature.) - - Angle No. Mu Eta Weight Angle in x-y Elevation - plane (degrees) (degrees) - 1 2.30091940D-01 2.30091940D-01 4.23616400D-02 45.00 71.01 - 2 6.88134320D-01 2.30091940D-01 4.09716925D-02 18.49 43.48 - 3 9.45576760D-01 2.30091940D-01 4.23616400D-02 13.68 13.30 - 4 2.30091940D-01 6.88134320D-01 4.09716925D-02 71.51 43.48 - 5 6.88134320D-01 6.88134320D-01 4.09716925D-02 45.00 13.30 - 6 2.30091940D-01 9.45576760D-01 4.23616400D-02 76.32 13.30 - - -1 ********************************** - ********************************** - ** Quadrature Reflection Angles ** - ********************************** - ********************************** - (Angles are in degrees) - - 12 primary (incoming) angles on each side. - 12 reflected (outgoing) angles on each side. - - Side No. (slope) Direction No. Primary angle (slope/elevation) Reflecting Angle (slope/elevation) - - 2 ( -90.00) 1 1 ( 45.00/ 71.01) 7 ( 135.00/ 71.01) - 2 2 ( 18.49/ 43.48) 8 ( 161.51/ 43.48) - 3 3 ( 13.68/ 13.30) 9 ( 166.32/ 13.30) - 4 4 ( 71.51/ 43.48) 10 ( 108.49/ 43.48) - 5 5 ( 45.00/ 13.30) 11 ( 135.00/ 13.30) - 6 6 ( 76.32/ 13.30) 12 ( 103.68/ 13.30) - 7 19 ( -45.00/ 71.01) 13 (-135.00/ 71.01) - 8 20 ( -18.49/ 43.48) 14 (-161.51/ 43.48) - 9 21 ( -13.68/ 13.30) 15 (-166.32/ 13.30) - 10 22 ( -71.51/ 43.48) 16 (-108.49/ 43.48) - 11 23 ( -45.00/ 13.30) 17 (-135.00/ 13.30) - 12 24 ( -76.32/ 13.30) 18 (-103.68/ 13.30) - - 4 ( 180.00) 1 13 (-135.00/ 71.01) 7 ( 135.00/ 71.01) - 2 14 (-161.51/ 43.48) 8 ( 161.51/ 43.48) - 3 15 (-166.32/ 13.30) 9 ( 166.32/ 13.30) - 4 16 (-108.49/ 43.48) 10 ( 108.49/ 43.48) - 5 17 (-135.00/ 13.30) 11 ( 135.00/ 13.30) - 6 18 (-103.68/ 13.30) 12 ( 103.68/ 13.30) - 7 19 ( -45.00/ 71.01) 1 ( 45.00/ 71.01) - 8 20 ( -18.49/ 43.48) 2 ( 18.49/ 43.48) - 9 21 ( -13.68/ 13.30) 3 ( 13.68/ 13.30) - 10 22 ( -71.51/ 43.48) 4 ( 71.51/ 43.48) - 11 23 ( -45.00/ 13.30) 5 ( 45.00/ 13.30) - 12 24 ( -76.32/ 13.30) 6 ( 76.32/ 13.30) - - - - - - Scattering constants for first quadrant - ( 6 moments for p-2 scattering) - Angle ------------------------- Moment No. ------------------------- - No. 0 1 2 3 4 5 - 1 1.00000E+00 3.98531E-01 -2.81804E-01 -9.40460E-01 -1.44988E-01 -1.15182E+00 - 2 1.00000E+00 1.19188E+00 -2.81804E-01 4.70230E-01 -4.33616E-01 -5.75912E-01 - 3 1.00000E+00 1.63779E+00 -2.81804E-01 1.88092E+00 -5.95839E-01 1.47310E-08 - 4 1.00000E+00 3.98531E-01 -8.42789E-01 -9.40460E-01 -4.33616E-01 -1.97499E-08 - 5 1.00000E+00 1.19188E+00 -8.42789E-01 4.70230E-01 -1.29681E+00 5.75912E-01 - 6 1.00000E+00 3.98531E-01 -1.15809E+00 -9.40460E-01 -5.95839E-01 1.15182E+00 - - -1 ************************** - ************************** - **** Mixture Volumes ***** - ************************** - ************************** - - - Mixture Volume Volume - ID (cc/unit length) Fraction - 0 0.00000E+00 0.00000E+00 - 1 1.40000E+03 4.35611E-01 - 10 7.00000E+02 2.17806E-01 - 2 2.94605E+02 9.16667E-02 - 3 1.41041E+02 4.38849E-02 - 4 6.78229E+02 2.11032E-01 - Total 3.21387E+03 1.00000E+00 - - - Memory size of attenuation factors(MB)= 5.86 - - - - ===================================================================================== - Outer iteration sweep begins. - - Outer Eigen- Eigenvalue Max Flux Max Flux Max Fuel Max Fuel Inners - It. # value Delta Delta Location(r,g) Delta Location(r,g) Cnvrged - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 1 1.00000 0.000E+00 6.185E+07 ( 3,252) 1.000E+00 ( 3, 0) F - 2 0.69905 4.305E-01 1.041E+05 ( 4,248) 1.182E+00 ( 25, 0) F - 3 1.02592 3.186E-01 7.092E+03 ( 4,240) 2.490E-01 ( 4, 0) F - 4 0.91001 1.274E-01 3.168E+02 ( 4, 1) 3.571E-01 ( 3, 0) F - 5 1.01478 1.032E-01 2.193E+02 ( 4, 1) 1.614E-01 ( 17, 0) F - 6 0.97406 4.180E-02 1.180E+02 ( 34, 1) 2.062E-01 ( 3, 0) F - 7 1.00581 3.157E-02 1.092E+02 ( 3, 1) 1.171E-01 ( 17, 0) F - 8 0.99252 1.340E-02 1.063E+02 ( 3, 1) 1.057E-01 ( 3, 0) F - 9 1.00205 9.518E-03 1.013E+02 ( 3, 1) 6.473E-02 ( 17, 0) F - 10 0.99783 4.235E-03 1.004E+02 ( 2, 1) 5.093E-02 ( 17, 0) F - 11 1.00071 2.876E-03 1.740E+00 ( 1,235) 3.229E-02 ( 17, 0) F - 12 0.99938 1.329E-03 9.695E-01 ( 1,235) 2.389E-02 ( 17, 0) F - 13 1.00025 8.700E-04 3.899E-01 ( 1,234) 1.551E-02 ( 17, 0) F - 14 0.99984 4.051E-04 1.297E-01 ( 1,234) 1.113E-02 ( 17, 0) F - 15 1.00004 1.952E-04 4.070E-02 ( 1,234) 7.368E-03 ( 17, 0) F - 16 1.00000 3.398E-05 2.138E-02 ( 1, 76) 5.159E-03 ( 17, 0) F - 17 1.00001 3.594E-06 1.530E-02 ( 1, 76) 3.518E-03 ( 17, 0) F - 18 1.00001 4.401E-07 1.073E-02 ( 1, 76) 2.417E-03 ( 17, 0) F - 19 1.00001 7.860E-07 7.422E-03 ( 1, 88) 1.663E-03 ( 22, 0) F - 20 1.00001 1.253E-06 5.059E-03 ( 1, 88) 1.145E-03 ( 22, 0) F - 21 1.00000 1.320E-06 3.378E-03 ( 1, 88) 7.891E-04 ( 22, 0) F - 22 1.00000 1.197E-06 2.207E-03 ( 1, 88) 5.438E-04 ( 22, 0) F - 23 1.00000 1.002E-06 1.410E-03 ( 1, 88) 3.746E-04 ( 22, 0) F - 24 1.00000 8.019E-07 8.788E-04 ( 1, 88) 2.575E-04 ( 22, 0) F - 25 1.00000 6.181E-07 5.309E-04 ( 1, 88) 1.772E-04 ( 22, 0) F - 26 1.00000 4.578E-07 3.103E-04 ( 1, 88) 1.217E-04 ( 22, 0) F - 27 1.00000 3.064E-07 1.745E-04 ( 1, 88) 7.985E-05 ( 22, 0) F - 28 1.00000 1.363E-07 1.091E-04 ( 2, 74) 4.820E-05 ( 17, 0) F - k-eff = 1.00000028 - - - - Four-Factor Estimate of k-infinity. Fast/Thermal boundary: 0.6250 eV - Fiss. neutrons/thermal abs. (eta): 1.596735 - Thermal utilization (f): 0.327632 - Resonance Escape probability (p): 0.000148 - Fast-fission factor (epsilon): 28640.652159 - -------------- - Infinite neutron multiplication 2.222343 - - Six-Factor Estimate of k-effective. - Fast non-leakage prob. (Pf): 0.449976 - Thermal non-leakage prob. (Pt): 0.999993 - -------------- - Effective neutron multiplication: 0.999993 - - - Alternate 3-Group Formulation Terms. - Fast/Res. boundary: 0.8210 MeV. Res./Therm. boundary: 0.6250eV - Fiss. neutrons/thermal abs. (eta): 1.596735 - Thermal utilization (f): 0.327632 - Resonance Escape probability (p): 0.000278 - Fast Escape probability (pf): 0.533425 - Fast-fission factor (epsilon): 28640.652159 - Fast non-leakage prob. (Pf): 0.728642 - Resonance non-leakage prob. (Pr): 0.617554 - Thermal non-leakage prob. (Pt): 0.999993 - - - - - Critical Spectra Calculated with the B1 Method - Input buckling of 0.00000E+00 gives a corresponding k-eff of 2.08115E+00 - Equivalent core height is infinite. - - Group Upper Critical Inf.Med. Critical Inf.Med. - No Energy(eV) Flux Flux Adj.Flux Adj.Flux - - 1 2.0000E+07 0.000E+00 0.000E+00 3.552E-03 3.552E-03 - 2 1.7330E+07 3.568E-11 3.568E-11 5.029E-03 5.029E-03 - 3 1.5680E+07 1.429E-07 1.429E-07 6.317E-03 6.317E-03 - 4 1.4550E+07 5.240E-06 5.240E-06 8.934E-03 8.934E-03 - 5 1.3840E+07 1.861E-05 1.861E-05 5.622E-03 5.622E-03 - 6 1.2840E+07 7.993E-05 7.993E-05 1.528E-03 1.528E-03 - 7 1.0000E+07 2.948E-04 2.948E-04 1.811E-03 1.811E-03 - 8 8.1870E+06 8.265E-04 8.265E-04 1.427E-03 1.427E-03 - 9 6.4340E+06 2.021E-03 2.021E-03 1.141E-03 1.141E-03 - 10 4.8000E+06 3.293E-03 3.293E-03 3.122E-03 3.122E-03 - 11 4.3040E+06 5.320E-03 5.320E-03 9.614E-04 9.614E-04 - 12 3.0000E+06 8.209E-03 8.209E-03 1.863E-03 1.863E-03 - 13 2.4790E+06 9.457E-03 9.457E-03 6.921E-03 6.921E-03 - 14 2.3540E+06 1.137E-02 1.137E-02 1.488E-03 1.488E-03 - 15 1.8500E+06 1.393E-02 1.393E-02 1.711E-03 1.711E-03 - 16 1.5000E+06 1.495E-02 1.495E-02 5.224E-03 5.224E-03 - 17 1.4000E+06 1.615E-02 1.615E-02 1.133E-02 1.133E-02 - 18 1.3560E+06 1.611E-02 1.611E-02 1.246E-02 1.246E-02 - 19 1.3170E+06 1.709E-02 1.709E-02 6.997E-03 6.997E-03 - 20 1.2500E+06 1.788E-02 1.788E-02 8.983E-03 8.983E-03 - 21 1.2000E+06 1.918E-02 1.918E-02 4.250E-03 4.250E-03 - 22 1.1000E+06 2.034E-02 2.034E-02 4.325E-03 4.325E-03 - 23 1.0100E+06 2.345E-02 2.345E-02 3.945E-03 3.945E-03 - 24 9.2000E+05 2.169E-02 2.169E-02 1.678E-02 1.678E-02 - 25 9.0000E+05 2.322E-02 2.322E-02 1.307E-02 1.307E-02 - 26 8.7500E+05 2.078E-02 2.078E-02 2.291E-02 2.291E-02 - 27 8.6110E+05 2.472E-02 2.472E-02 7.544E-03 7.544E-03 - 28 8.2000E+05 2.314E-02 2.314E-02 4.105E-03 4.105E-03 - 29 7.5000E+05 2.587E-02 2.587E-02 3.652E-03 3.652E-03 - 30 6.7900E+05 2.596E-02 2.596E-02 2.707E-02 2.707E-02 - 31 6.7000E+05 3.058E-02 3.058E-02 3.260E-03 3.260E-03 - 32 6.0000E+05 3.244E-02 3.244E-02 7.795E-03 7.795E-03 - 33 5.7300E+05 3.276E-02 3.276E-02 8.724E-03 8.724E-03 - 34 5.5000E+05 3.503E-02 3.503E-02 3.182E-03 3.182E-03 - 35 4.9200E+05 3.007E-02 3.007E-02 7.646E-03 7.646E-03 - 36 4.7000E+05 3.192E-02 3.192E-02 5.277E-03 5.277E-03 - 37 4.4000E+05 3.038E-02 3.038E-02 7.423E-03 7.423E-03 - 38 4.2000E+05 2.537E-02 2.537E-02 7.026E-03 7.026E-03 - 39 4.0000E+05 3.164E-02 3.164E-02 1.767E-03 1.767E-03 - 40 3.3000E+05 3.656E-02 3.656E-02 1.654E-03 1.654E-03 - 41 2.7000E+05 2.823E-02 2.823E-02 1.071E-03 1.071E-03 - 42 2.0000E+05 2.897E-02 2.897E-02 1.052E-03 1.052E-03 - 43 1.4900E+05 2.663E-02 2.663E-02 1.984E-03 1.984E-03 - 44 1.2830E+05 2.187E-02 2.187E-02 1.162E-03 1.162E-03 - 45 1.0000E+05 1.936E-02 1.936E-02 1.716E-03 1.716E-03 - 46 8.5000E+04 2.058E-02 2.058E-02 7.597E-03 7.597E-03 - 47 8.2000E+04 2.009E-02 2.009E-02 3.031E-03 3.031E-03 - 48 7.5000E+04 1.631E-02 1.631E-02 9.873E-03 9.873E-03 - 49 7.3000E+04 1.713E-02 1.713E-02 1.347E-03 1.347E-03 - 50 6.0000E+04 1.054E-02 1.054E-02 1.783E-03 1.783E-03 - 51 5.2000E+04 1.040E-02 1.040E-02 6.439E-03 6.439E-03 - 52 5.0000E+04 1.369E-02 1.369E-02 2.379E-03 2.379E-03 - 53 4.5000E+04 8.834E-03 8.834E-03 6.064E-04 6.064E-04 - 54 3.0000E+04 9.879E-03 9.879E-03 5.782E-04 5.782E-04 - 55 2.0000E+04 7.048E-03 7.048E-03 1.372E-03 1.372E-03 - 56 1.7000E+04 5.109E-03 5.109E-03 8.112E-04 8.112E-04 - 57 1.3000E+04 4.310E-03 4.310E-03 6.780E-04 6.780E-04 - 58 9.5000E+03 2.643E-03 2.643E-03 1.229E-03 1.229E-03 - 59 8.0300E+03 2.589E-03 2.589E-03 5.961E-04 5.961E-04 - 60 5.7000E+03 2.028E-03 2.028E-03 5.301E-04 5.301E-04 - 61 3.9000E+03 1.648E-03 1.648E-03 4.770E-03 4.770E-03 - 62 3.7400E+03 1.377E-03 1.377E-03 9.014E-04 9.014E-04 - 63 3.0000E+03 9.725E-04 9.725E-04 1.085E-03 1.085E-03 - 64 2.5000E+03 1.239E-03 1.239E-03 1.877E-03 1.877E-03 - 65 2.2500E+03 1.541E-03 1.541E-03 8.770E-03 8.770E-03 - 66 2.2000E+03 1.884E-03 1.884E-03 9.857E-04 9.857E-04 - 67 1.8000E+03 1.898E-03 1.898E-03 1.309E-03 1.309E-03 - 68 1.5500E+03 1.787E-03 1.787E-03 5.975E-03 5.975E-03 - 69 1.5000E+03 1.450E-03 1.450E-03 7.319E-04 7.319E-04 - 70 1.1500E+03 1.031E-03 1.031E-03 9.830E-04 9.830E-04 - 71 9.5000E+02 9.211E-04 9.211E-04 5.652E-04 5.652E-04 - 72 6.8300E+02 8.488E-04 8.488E-04 9.507E-03 9.507E-03 - 73 6.7000E+02 7.875E-04 7.875E-04 9.190E-04 9.190E-04 - 74 5.5000E+02 5.632E-04 5.632E-04 3.052E-04 3.052E-04 - 75 3.0500E+02 4.941E-04 4.941E-04 2.603E-03 2.603E-03 - 76 2.8500E+02 4.917E-04 4.917E-04 1.033E-03 1.033E-03 - 77 2.4000E+02 4.763E-04 4.763E-04 1.968E-03 1.968E-03 - 78 2.2000E+02 4.778E-04 4.778E-04 3.529E-03 3.529E-03 - 79 2.0950E+02 4.131E-04 4.131E-04 1.642E-02 1.642E-02 - 80 2.0740E+02 4.638E-04 4.638E-04 6.744E-03 6.744E-03 - 81 2.0200E+02 4.401E-04 4.401E-04 3.842E-03 3.842E-03 - 82 1.9300E+02 4.430E-04 4.430E-04 2.151E-02 2.151E-02 - 83 1.9150E+02 3.951E-04 3.951E-04 1.044E-02 1.044E-02 - 84 1.8850E+02 5.118E-04 5.118E-04 3.920E-02 3.920E-02 - 85 1.8770E+02 4.593E-04 4.593E-04 4.061E-03 4.061E-03 - 86 1.8000E+02 4.402E-04 4.402E-04 2.919E-03 2.919E-03 - 87 1.7000E+02 4.146E-04 4.146E-04 9.652E-04 9.652E-04 - 88 1.4300E+02 3.872E-04 3.872E-04 1.038E-03 1.038E-03 - 89 1.2200E+02 3.577E-04 3.577E-04 6.464E-03 6.464E-03 - 90 1.1900E+02 3.474E-04 3.474E-04 1.273E-02 1.273E-02 - 91 1.1750E+02 3.260E-04 3.260E-04 1.212E-02 1.212E-02 - 92 1.1600E+02 3.720E-04 3.720E-04 6.112E-03 6.112E-03 - 93 1.1300E+02 3.679E-04 3.679E-04 3.455E-03 3.455E-03 - 94 1.0800E+02 3.297E-04 3.297E-04 5.520E-03 5.520E-03 - 95 1.0500E+02 3.301E-04 3.301E-04 4.200E-03 4.200E-03 - 96 1.0120E+02 3.465E-04 3.465E-04 3.813E-03 3.813E-03 - 97 9.7000E+01 3.200E-04 3.200E-04 2.085E-03 2.085E-03 - 98 9.0000E+01 3.033E-04 3.033E-04 1.636E-03 1.636E-03 - 99 8.1700E+01 2.934E-04 2.934E-04 6.939E-03 6.939E-03 - 100 8.0000E+01 3.049E-04 3.049E-04 2.846E-03 2.846E-03 - 101 7.6000E+01 2.784E-04 2.784E-04 2.743E-03 2.743E-03 - 102 7.2000E+01 2.911E-04 2.911E-04 2.197E-03 2.197E-03 - 103 6.7500E+01 2.623E-04 2.623E-04 3.693E-03 3.693E-03 - 104 6.5000E+01 2.764E-04 2.764E-04 4.526E-03 4.526E-03 - 105 6.3000E+01 2.696E-04 2.696E-04 4.436E-03 4.436E-03 - 106 6.1000E+01 2.541E-04 2.541E-04 2.851E-03 2.851E-03 - 107 5.8000E+01 2.501E-04 2.501E-04 1.591E-03 1.591E-03 - 108 5.3400E+01 2.423E-04 2.423E-04 2.300E-03 2.300E-03 - 109 5.0600E+01 2.374E-04 2.374E-04 2.702E-03 2.702E-03 - 110 4.8300E+01 2.387E-04 2.387E-04 1.873E-03 1.873E-03 - 111 4.5200E+01 2.237E-04 2.237E-04 4.336E-03 4.336E-03 - 112 4.4000E+01 2.372E-04 2.372E-04 3.251E-03 3.251E-03 - 113 4.2400E+01 2.192E-04 2.192E-04 3.560E-03 3.560E-03 - 114 4.1000E+01 2.302E-04 2.302E-04 3.482E-03 3.482E-03 - 115 3.9600E+01 2.248E-04 2.248E-04 9.646E-03 9.646E-03 - 116 3.9100E+01 2.192E-04 2.192E-04 4.119E-03 4.119E-03 - 117 3.8000E+01 2.251E-04 2.251E-04 1.227E-02 1.227E-02 - 118 3.7630E+01 2.199E-04 2.199E-04 1.247E-02 1.247E-02 - 119 3.7270E+01 2.110E-04 2.110E-04 3.140E-02 3.140E-02 - 120 3.7130E+01 2.022E-04 2.022E-04 3.316E-02 3.316E-02 - 121 3.7000E+01 2.031E-04 2.031E-04 4.219E-03 4.219E-03 - 122 3.6000E+01 2.212E-04 2.212E-04 8.825E-03 8.825E-03 - 123 3.5500E+01 2.072E-04 2.072E-04 9.033E-03 9.033E-03 - 124 3.5000E+01 2.108E-04 2.108E-04 3.456E-03 3.456E-03 - 125 3.3750E+01 2.089E-04 2.089E-04 8.328E-03 8.328E-03 - 126 3.3250E+01 2.093E-04 2.093E-04 2.723E-03 2.723E-03 - 127 3.1750E+01 2.126E-04 2.126E-04 7.702E-03 7.702E-03 - 128 3.1250E+01 2.122E-04 2.122E-04 2.999E-03 2.999E-03 - 129 3.0000E+01 2.131E-04 2.131E-04 1.414E-03 1.414E-03 - 130 2.7500E+01 1.987E-04 1.987E-04 1.310E-03 1.310E-03 - 131 2.5000E+01 1.913E-04 1.913E-04 1.161E-03 1.161E-03 - 132 2.2500E+01 1.702E-04 1.702E-04 3.739E-03 3.739E-03 - 133 2.1750E+01 1.775E-04 1.775E-04 4.400E-03 4.400E-03 - 134 2.1200E+01 1.590E-04 1.590E-04 3.192E-03 3.192E-03 - 135 2.0500E+01 1.681E-04 1.681E-04 4.660E-03 4.660E-03 - 136 2.0000E+01 1.679E-04 1.679E-04 3.914E-03 3.914E-03 - 137 1.9400E+01 1.638E-04 1.638E-04 2.521E-03 2.521E-03 - 138 1.8500E+01 1.577E-04 1.577E-04 1.416E-03 1.416E-03 - 139 1.7000E+01 1.564E-04 1.564E-04 1.918E-03 1.918E-03 - 140 1.6000E+01 1.441E-04 1.441E-04 1.126E-03 1.126E-03 - 141 1.4400E+01 1.387E-04 1.387E-04 9.724E-04 9.724E-04 - 142 1.2900E+01 1.336E-04 1.336E-04 1.249E-03 1.249E-03 - 143 1.1900E+01 1.250E-04 1.250E-04 3.015E-03 3.015E-03 - 144 1.1500E+01 1.211E-04 1.211E-04 7.249E-04 7.249E-04 - 145 1.0000E+01 1.209E-04 1.209E-04 9.269E-04 9.269E-04 - 146 9.1000E+00 1.157E-04 1.157E-04 7.469E-04 7.469E-04 - 147 8.1000E+00 1.052E-04 1.052E-04 6.868E-04 6.868E-04 - 148 7.1500E+00 1.034E-04 1.034E-04 3.566E-03 3.566E-03 - 149 7.0000E+00 1.019E-04 1.019E-04 4.091E-03 4.091E-03 - 150 6.8750E+00 9.662E-05 9.662E-05 3.843E-03 3.843E-03 - 151 6.7500E+00 9.382E-05 9.382E-05 1.888E-03 1.888E-03 - 152 6.5000E+00 9.558E-05 9.558E-05 1.960E-03 1.960E-03 - 153 6.2500E+00 9.558E-05 9.558E-05 1.915E-03 1.915E-03 - 154 6.0000E+00 9.546E-05 9.546E-05 7.433E-04 7.433E-04 - 155 5.4000E+00 9.434E-05 9.434E-05 1.019E-03 1.019E-03 - 156 5.0000E+00 9.521E-05 9.521E-05 1.264E-03 1.264E-03 - 157 4.7000E+00 9.114E-05 9.114E-05 5.779E-04 5.779E-04 - 158 4.1000E+00 8.608E-05 8.608E-05 8.370E-04 8.370E-04 - 159 3.7300E+00 8.089E-05 8.089E-05 1.267E-03 1.267E-03 - 160 3.5000E+00 7.937E-05 7.937E-05 8.798E-04 8.798E-04 - 161 3.2000E+00 7.698E-05 7.698E-05 2.490E-03 2.490E-03 - 162 3.1000E+00 7.779E-05 7.779E-05 2.390E-03 2.390E-03 - 163 3.0000E+00 7.938E-05 7.938E-05 7.734E-03 7.734E-03 - 164 2.9700E+00 7.552E-05 7.552E-05 2.253E-03 2.253E-03 - 165 2.8700E+00 7.401E-05 7.401E-05 2.161E-03 2.161E-03 - 166 2.7700E+00 7.271E-05 7.271E-05 2.072E-03 2.072E-03 - 167 2.6700E+00 7.130E-05 7.130E-05 1.990E-03 1.990E-03 - 168 2.5700E+00 6.981E-05 6.981E-05 1.911E-03 1.911E-03 - 169 2.4700E+00 6.831E-05 6.831E-05 2.038E-03 2.038E-03 - 170 2.3800E+00 6.690E-05 6.690E-05 2.207E-03 2.207E-03 - 171 2.3000E+00 6.541E-05 6.541E-05 1.884E-03 1.884E-03 - 172 2.2100E+00 6.372E-05 6.372E-05 1.802E-03 1.802E-03 - 173 2.1200E+00 6.143E-05 6.143E-05 1.278E-03 1.278E-03 - 174 2.0000E+00 5.984E-05 5.984E-05 2.434E-03 2.434E-03 - 175 1.9400E+00 5.857E-05 5.857E-05 1.754E-03 1.754E-03 - 176 1.8600E+00 5.687E-05 5.687E-05 1.482E-03 1.482E-03 - 177 1.7700E+00 5.498E-05 5.498E-05 1.399E-03 1.399E-03 - 178 1.6800E+00 5.301E-05 5.301E-05 1.317E-03 1.317E-03 - 179 1.5900E+00 5.095E-05 5.095E-05 1.235E-03 1.235E-03 - 180 1.5000E+00 4.926E-05 4.926E-05 2.110E-03 2.110E-03 - 181 1.4500E+00 4.794E-05 4.794E-05 2.030E-03 2.030E-03 - 182 1.4000E+00 4.656E-05 4.656E-05 1.951E-03 1.951E-03 - 183 1.3500E+00 4.509E-05 4.509E-05 1.872E-03 1.872E-03 - 184 1.3000E+00 4.351E-05 4.351E-05 1.793E-03 1.793E-03 - 185 1.2500E+00 4.224E-05 4.224E-05 3.472E-03 3.472E-03 - 186 1.2250E+00 4.133E-05 4.133E-05 3.394E-03 3.394E-03 - 187 1.2000E+00 4.041E-05 4.041E-05 3.314E-03 3.314E-03 - 188 1.1750E+00 3.947E-05 3.947E-05 3.228E-03 3.228E-03 - 189 1.1500E+00 3.881E-05 3.881E-05 7.912E-03 7.912E-03 - 190 1.1400E+00 3.844E-05 3.844E-05 7.825E-03 7.825E-03 - 191 1.1300E+00 3.809E-05 3.809E-05 7.748E-03 7.748E-03 - 192 1.1200E+00 3.777E-05 3.777E-05 7.684E-03 7.684E-03 - 193 1.1100E+00 3.748E-05 3.748E-05 7.634E-03 7.634E-03 - 194 1.1000E+00 3.721E-05 3.721E-05 7.595E-03 7.595E-03 - 195 1.0900E+00 3.695E-05 3.695E-05 7.566E-03 7.566E-03 - 196 1.0800E+00 3.671E-05 3.671E-05 7.543E-03 7.543E-03 - 197 1.0700E+00 3.647E-05 3.647E-05 7.524E-03 7.524E-03 - 198 1.0600E+00 3.624E-05 3.624E-05 7.507E-03 7.507E-03 - 199 1.0500E+00 3.601E-05 3.601E-05 7.491E-03 7.491E-03 - 200 1.0400E+00 3.579E-05 3.579E-05 7.476E-03 7.476E-03 - 201 1.0300E+00 3.557E-05 3.557E-05 7.463E-03 7.463E-03 - 202 1.0200E+00 3.536E-05 3.536E-05 7.452E-03 7.452E-03 - 203 1.0100E+00 3.516E-05 3.516E-05 7.444E-03 7.444E-03 - 204 1.0000E+00 3.484E-05 3.484E-05 2.977E-03 2.977E-03 - 205 9.7500E-01 3.442E-05 3.442E-05 2.978E-03 2.978E-03 - 206 9.5000E-01 3.400E-05 3.400E-05 2.970E-03 2.970E-03 - 207 9.2500E-01 3.353E-05 3.353E-05 2.948E-03 2.948E-03 - 208 9.0000E-01 3.275E-05 3.275E-05 1.446E-03 1.446E-03 - 209 8.5000E-01 3.162E-05 3.162E-05 1.396E-03 1.396E-03 - 210 8.0000E-01 3.039E-05 3.039E-05 1.335E-03 1.335E-03 - 211 7.5000E-01 2.906E-05 2.906E-05 1.264E-03 1.264E-03 - 212 7.0000E-01 2.763E-05 2.763E-05 1.185E-03 1.185E-03 - 213 6.5000E-01 2.652E-05 2.652E-05 2.245E-03 2.245E-03 - 214 6.2500E-01 2.570E-05 2.570E-05 2.157E-03 2.157E-03 - 215 6.0000E-01 2.443E-05 2.443E-05 1.009E-03 1.009E-03 - 216 5.5000E-01 2.271E-05 2.271E-05 9.125E-04 9.125E-04 - 217 5.0000E-01 2.090E-05 2.090E-05 8.118E-04 8.118E-04 - 218 4.5000E-01 1.906E-05 1.906E-05 7.075E-04 7.075E-04 - 219 4.0000E-01 1.770E-05 1.770E-05 1.256E-03 1.256E-03 - 220 3.7500E-01 1.673E-05 1.673E-05 1.147E-03 1.147E-03 - 221 3.5000E-01 1.573E-05 1.573E-05 1.040E-03 1.040E-03 - 222 3.2500E-01 1.471E-05 1.471E-05 9.335E-04 9.335E-04 - 223 3.0000E-01 1.365E-05 1.365E-05 8.303E-04 8.303E-04 - 224 2.7500E-01 1.257E-05 1.257E-05 7.307E-04 7.307E-04 - 225 2.5000E-01 1.144E-05 1.144E-05 6.352E-04 6.352E-04 - 226 2.2500E-01 1.026E-05 1.026E-05 5.443E-04 5.443E-04 - 227 2.0000E-01 9.011E-06 9.011E-06 4.584E-04 4.584E-04 - 228 1.7500E-01 7.688E-06 7.688E-06 3.776E-04 3.776E-04 - 229 1.5000E-01 6.282E-06 6.282E-06 3.018E-04 3.018E-04 - 230 1.2500E-01 4.941E-06 4.941E-06 2.310E-04 2.310E-04 - 231 1.0000E-01 4.075E-06 4.075E-06 4.669E-04 4.669E-04 - 232 9.0000E-02 3.485E-06 3.485E-06 4.020E-04 4.020E-04 - 233 8.0000E-02 2.893E-06 2.893E-06 3.395E-04 3.395E-04 - 234 7.0000E-02 2.310E-06 2.310E-06 2.798E-04 2.798E-04 - 235 6.0000E-02 1.750E-06 1.750E-06 2.232E-04 2.232E-04 - 236 5.0000E-02 1.229E-06 1.229E-06 1.701E-04 1.701E-04 - 237 4.0000E-02 7.694E-07 7.694E-07 1.209E-04 1.209E-04 - 238 3.0000E-02 4.885E-07 4.885E-07 1.879E-04 1.879E-04 - 239 2.5300E-02 1.867E-07 1.867E-07 3.013E-05 3.013E-05 - 240 1.0000E-02 4.799E-08 4.799E-08 8.010E-05 8.010E-05 - 241 7.5000E-03 2.319E-08 2.319E-08 5.114E-05 5.114E-05 - 242 5.0000E-03 1.155E-08 1.155E-08 8.438E-05 8.438E-05 - 243 4.0000E-03 6.726E-09 6.726E-09 6.125E-05 6.125E-05 - 244 3.0000E-03 4.033E-09 4.033E-09 9.117E-05 9.117E-05 - 245 2.5000E-03 2.630E-09 2.630E-09 7.134E-05 7.134E-05 - 246 2.0000E-03 1.541E-09 1.541E-09 5.276E-05 5.276E-05 - 247 1.5000E-03 8.953E-10 8.953E-10 6.518E-05 6.518E-05 - 248 1.2000E-03 5.841E-10 5.841E-10 7.753E-05 7.753E-05 - 249 1.0000E-03 3.599E-10 3.599E-10 4.783E-05 4.783E-05 - 250 7.5000E-04 1.770E-10 1.770E-10 3.286E-05 3.286E-05 - 251 5.0000E-04 3.222E-11 3.222E-11 7.860E-06 7.860E-06 - 252 1.0000E-04 5.613E-13 5.613E-13 5.157E-06 5.157E-06 - - NOTE: Flux values above are per unit-lethargy. - - - - Zone data - - Zone Mix Relative Relative Zone - No. No. Flux Fission Power Volume - ----------------------------------------------------------------- - 1 1 1.87787E-01 7.98654E-01 1.40000E+03 - 2 3 1.40417E-01 0.00000E+00 4.00000E+01 - 3 2 1.24511E-01 0.00000E+00 2.20000E+02 - 4 3 9.27423E-02 0.00000E+00 4.00000E+01 - 5 10 4.09596E-02 2.01346E-01 7.00000E+02 - 6 3 1.12383E-01 0.00000E+00 3.39115E+01 - 7 4 1.00545E-01 0.00000E+00 4.40849E+02 - 8 3 7.01993E-02 0.00000E+00 1.35646E+01 - 9 2 6.09558E-02 0.00000E+00 7.46052E+01 - 10 3 4.77119E-02 0.00000E+00 1.35646E+01 - 11 4 2.17874E-02 0.00000E+00 2.37380E+02 - - - -module xsproc is finished. completion code 0. -Sequence T-Newt is finished. completion code 0. - - SCALE is finished on Wed Mar 08 15:56:49 2017. --------------------------- Summary -------------------------- -t-newt finished. used 59.7199 seconds. -SCALE driver required a maximum of 65.7 MiB of RAM. ------------------------- End Summary ------------------------ - diff --git a/testing/param_tester.txt b/testing/param_tester.txt index e7eb3bd..9357d66 100644 --- a/testing/param_tester.txt +++ b/testing/param_tester.txt @@ -1,3 +1,4 @@ var_char $ iter_var del_z 10 5.0 15.0 -iter_lim 50 \ No newline at end of file +iter_lim 50 +eps_k 1E-4 \ No newline at end of file From ebbd512b06e087af7d158a9177282cd92172237f Mon Sep 17 00:00:00 2001 From: drewejohnson <1drew.e.johnson@gmail.com> Date: Mon, 13 Mar 2017 14:30:05 -0400 Subject: [PATCH 3/3] Cleaned up output options Signed-off-by: drewejohnson <1drew.e.johnson@gmail.com> --- globalparams.py | 4 ++-- iterator.py | 3 +-- 2 files changed, 3 insertions(+), 4 deletions(-) diff --git a/globalparams.py b/globalparams.py index a1dab4a..6070d8b 100644 --- a/globalparams.py +++ b/globalparams.py @@ -37,7 +37,7 @@ [1]: M. D. DeHart, and S. Bowman, "Reactor Physics Methods and Analysis \tCapabilities in SCALE," Nuclear Technology, Technical Paper -vol. 174, no.2, pp. 196-213, 2011. +\tvol. 174, no.2, pp. 196-213, 2011. """ + lineBreakLong @@ -66,7 +66,7 @@ # SCALE parameters var_char = '$' -exe_str = r'C:\SCALE-6.2.1\bin\scalerte' +exe_str = r'C:\SCALE-6.2.1\bin\scalerte.exe' template_file = None # Outputs diff --git a/iterator.py b/iterator.py index d69a7d3..6ed9f40 100644 --- a/iterator.py +++ b/iterator.py @@ -91,7 +91,7 @@ def itermain(): while _n < gp.iter_lim: _n += 1 _iter_file = makefile(gp.args.inp_file, _n) - utils.vprint('Running SCALE iteration number {}...'.format(_n)) + utils.vprint('Running SCALE iteration number {}...'.format(_n)) subprocess.call([gp.exe_str, _iter_file]) utils.vprint(' done') stat, _k = parse_scale_out_eig(_iter_file.replace('.inp', '.out')) @@ -102,7 +102,6 @@ def itermain(): utils.error('Could not find value of k-eff for iteration file {0}.inp\n' 'Check {0}.out for error message'.format(_iter_file.split('.')[0]), 'itermain() of iteration {}'.format(_n)) - # todo: Parse through output and find cause of error # check for convergance based on updated eigenvalue, and then a termination based on exceeding the specified # input range from the parameter file