Page 1 of 1

Test case 1: BWR burnup calculation

Posted: Sun Nov 25, 2012 3:55 pm
by Jaakko Leppänen
The input files for the first test case are found at:

http://virtual.vtt.fi/virtual/montecarl ... WR_BURNUP/

This is a BWR burnup calculation based on a benchmark:

A. Yamamoto, et al. "Benchmark Problem Suite for Reactor Physics Study of LWR Next Generation Fuels." J. Nucl. Sci. Technol. 39 (2002) 900-912.

Geometry, material and depletion parameters are defined in input file "bwr.inp". Division into 5, 41, 113 and 436 depletion zones is carried out using the "div" card in files "bwr5", "bwr41", "bwr113" and "bwr436". The number of neutron histories is increased along with the number of depletion zones. Since the calculation of speed-up factors requries a single-CPU reference run, the fuel is burnt only to 5 MWd/kgU in order to complete the 436-zone calculation within a reasonable time. The build-up of plutonium and fission products increases absorption in fuel, which shortens the neutron histories. This will most likely affect the scalability, so it might be a good idea to run a few tests with extended burnup. The depletion history is defined at the end of "bwr.inp".

The download directory also contains Matlab / Octave script "test.m", which automatically calculates speed-up factors, etc. from the Serpent output. Before running the script, you need to read the "_res.m" output files in the order of increasing parallelization.

I will try to get the next Serpent update (2.1.10) distributed next week, but this case can probably be run with version 2.1.9 as well.

Re: Test case 1: BWR burnup calculation

Posted: Mon Jan 28, 2013 9:44 am
by hartanto
This is the result from KAIST for the test case 1 "bwr5" and "bwr41". No change was made in the input file. I will update the result for bwr113 and bwr436 later after the calculation has been done.

The CPU type is Intel(R) Xeon(R) CPU E5-2690 0 @ 2.90GHz.

Code: Select all

BWR5
Optimization mode: 4
OpenMP reproducibility: 0
Shared scoring buffer: 0

The three values in the columns are running time (min), speed-up factor and
fraction of total running time.

------------------------------------------------------------------------------------------
NODE MPI OMP PARA :  Total     :  Transport        :  Burnup          :  Processing      :
------------------------------------------------------------------------------------------
  1    1   1    1 : 34.0   1.0 :  32.0   1.0  0.94 :  1.5   1.0  0.04 :  0.5   1.0  0.02 : 
  1    1   2    2 : 16.6   2.1 :  15.4   2.1  0.93 :  0.9   1.7  0.05 :  0.3   2.0  0.02 : 
  1    2   1    2 : 16.5   2.1 :  15.3   2.1  0.93 :  0.8   1.8  0.05 :  0.3   1.7  0.02 : 
  1    1   4    4 :  8.5   4.0 :   7.7   4.1  0.91 :  0.6   2.6  0.07 :  0.2   3.4  0.02 : 
  1    2   2    4 :  8.7   3.9 :   7.9   4.0  0.91 :  0.6   2.6  0.07 :  0.2   2.9  0.02 : 
  1    4   1    4 :  9.0   3.8 :   8.2   3.9  0.91 :  0.6   2.6  0.06 :  0.2   2.4  0.03 : 
  1    1   8    8 :  4.4   7.7 :   3.9   8.1  0.90 :  0.3   4.5  0.08 :  0.1   4.4  0.03 : 
  1    2   4    8 :  4.5   7.6 :   4.0   7.9  0.90 :  0.3   4.6  0.07 :  0.1   4.7  0.03 : 
  1    4   2    8 :  4.6   7.4 :   4.1   7.8  0.90 :  0.3   4.7  0.07 :  0.1   3.8  0.03 : 
  1    8   1    8 :  4.5   7.5 :   4.0   7.9  0.89 :  0.3   4.6  0.07 :  0.2   3.5  0.03 : 
  1    1  16   16 :  2.5  13.6 :   2.1  15.6  0.82 :  0.3   4.6  0.13 :  0.1   4.1  0.05 : 
  1    2   8   16 :  2.5  13.5 :   2.1  15.2  0.83 :  0.3   4.6  0.13 :  0.1   5.7  0.04 : 
  1    4   4   16 :  2.5  13.4 :   2.1  15.1  0.84 :  0.3   4.8  0.12 :  0.1   5.4  0.04 : 
  1    8   2   16 :  2.6  13.2 :   2.2  14.8  0.83 :  0.3   4.7  0.12 :  0.1   5.1  0.04 : 
  1   16   1   16 :  Lack of memory
  2    2   8   16 :  2.4  13.9 :   2.0  15.7  0.83 :  0.3   4.9  0.13 :  0.1   5.1  0.04 : 
  4    4   4   16 :  2.4  14.0 :   2.0  15.9  0.83 :  0.3   4.9  0.12 :  0.1   4.7  0.05 : 
  2    2  16   32 :  1.5  22.0 :   1.1  28.7  0.72 :  0.3   4.9  0.20 :  0.1   4.5  0.08 :
  2   32   1   32 :  Lack of memory
  4    4   8   32 :  1.5  23.3 :   1.1  30.1  0.73 :  0.3   5.0  0.21 :  0.1   5.8  0.06 : 
  8    8   4   32 :  1.5  22.5 :   1.1  29.2  0.72 :  0.3   4.9  0.20 :  0.1   4.8  0.07 : 
------------------------------------------------------------------------------------------
BWR5.png
BWR5
BWR5.png (25.78 KiB) Viewed 2628 times

Code: Select all

BWR41
Optimization mode: 4
OpenMP reproducibility: 0
Shared scoring buffer: 0

The three values in the columns are running time (min), speed-up factor and
fraction of total running time.

-------------------------------------------------------------------------------------------
NODE MPI OMP PARA :  Total      :  Transport        :  Burnup          :  Processing      :
-------------------------------------------------------------------------------------------
  1   1   1    1  : 277.7   1.0 : 261.5   1.0  0.94 : 12.2   1.0  0.04 :  3.9   1.0  0.01 : 
  1   1   2    2  : 134.9   2.1 : 126.7   2.1  0.94 :  6.2   2.0  0.05 :  2.0   2.0  0.01 : 
  1   2   1    2  : 135.1   2.1 : 126.9   2.1  0.94 :  6.2   2.0  0.05 :  2.0   2.0  0.01 : 
  1   1   4    4  :  68.3   4.1 :  63.8   4.1  0.93 :  3.3   3.7  0.05 :  1.2   3.4  0.02 : 
  1   2   2    4  :  70.1   4.0 :  65.7   4.0  0.94 :  3.3   3.7  0.05 :  1.1   3.7  0.02 : 
  1   4   1    4  :  70.6   3.9 :  66.2   4.0  0.94 :  3.3   3.7  0.05 :  1.1   3.4  0.02 : 
  1   1   8    8  :  35.1   7.9 :  32.3   8.1  0.92 :  1.9   6.3  0.06 :  0.9   4.6  0.02 : 
  1   2   4    8  :  35.8   7.8 :  33.2   7.9  0.93 :  1.9   6.3  0.05 :  0.7   5.5  0.02 : 
  1   4   2    8  :  35.2   7.9 :  32.7   8.0  0.93 :  1.9   6.5  0.05 :  0.6   6.1  0.02 : 
  1   8   1    8  :  35.8   7.8 :  33.1   7.9  0.93 :  2.0   6.2  0.06 :  0.7   5.6  0.02 : 
  1   1  16   16  :  19.0  14.6 :  16.7  15.7  0.88 :  1.3   9.7  0.07 :  1.1   3.7  0.06 : 
  1   2   8   16  :  19.0  14.6 :  17.2  15.2  0.90 :  1.2  10.2  0.06 :  0.6   6.4  0.03 : 
  1   4   4   16  :  19.6  14.2 :  17.7  14.8  0.90 :  1.3   9.5  0.07 :  0.6   6.4  0.03 : 
  1   8   2   16  :  19.7  14.1 :  17.9  14.6  0.91 :  1.2   9.8  0.06 :  0.5   7.5  0.03 :
  1  16   1   16  :  Lack of memory
  2   2   8   16  :  18.0  15.5 :  16.3  16.1  0.91 :  1.0  12.3  0.06 :  0.7   5.8  0.04 : 
  4   4   4   16  :  17.9  15.5 :  16.3  16.0  0.91 :  0.9  13.2  0.05 :  0.7   5.7  0.04 : 
  2   2  16   32  :  10.0  27.8 :   8.5  30.9  0.85 :  0.7  16.3  0.07 :  0.8   5.0  0.08 :
  2  32   1   32  :  Lack of memory  
  4   4   8   32  :   9.5  29.2 :   8.3  31.5  0.87 :  0.6  19.0  0.07 :  0.6   7.1  0.06 : 
  8   8   4   32  :   9.5  29.2 :   8.3  31.4  0.88 :  0.6  19.9  0.06 :  0.6   6.8  0.06 : 
 ------------------------------------------------------------------------------------------
bwr41.png
BWR41
bwr41.png (14.75 KiB) Viewed 2615 times

Re: Test case 1: BWR burnup calculation

Posted: Mon Jan 28, 2013 2:36 pm
by Jaakko Leppänen
Thanks! The scalability of the transport simulation looks pretty good. The performance of the processing and burnup routines is limited by the fact that there are only 5 burnable materials, so when the number of CPU's is large, most of them are idling during this part of the calculation. There's clearly a lot of room for improvement there.