## set power normalization

Discussion on physics, models and calculation methods
aleja311
Posts: 10
Joined: Tue Jun 04, 2013 10:29 am
Security question 1: No
Security question 2: 92

### set power normalization

I am attempting to calculate the leakage probability in my core. The model is composed of a TRU cylinder surrounded by a HT9 reflector.
My input also utilizes burnup and sets a power level. This then automatically normalizes my detector results for surface current.
For the normalization, is it taking the surface current in neutrons/cm^2 per source neutron and multiplying by a neutron source in neutrons/second?
If I divide my result by a neutron source I get a similar answer to the result if I did not set a power level so I am assuming this is correct.
However, if I want to divide by a neutron source I need the Q value and eta. Do I have to calculate these additionally or is this stated somewhere?

Also in regards to current surface output without a set power. The answer should be in neutrons/cm^2 per source neutron. Can I pretty much assume this to be the leakage probability? It seems correct because it is calculating surface current for one neutron which is pretty much the fraction of leakage but I want to double check.

Thank you!
Alejandra

Jaakko Leppänen
Posts: 2442
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

### Re: set power normalization

The units for the surface current tally are in neutron/s. The value is not divided by the surface area, because there is no safe way to calculate that within the code.

In the _res.m output there is a separate section for normalized total reaction rates:

Code: Select all

``````% Normalized total reaction rates:

TOT_POWER                 (idx, [1:   2]) = [  1.21101E-11 0.00055 ];
TOT_GENRATE               (idx, [1:   2]) = [  9.18263E-01 0.00055 ];
TOT_FISSRATE              (idx, [1:   2]) = [  3.73240E-01 0.00056 ];
TOT_ABSRATE               (idx, [1:   2]) = [  1.00000E+00 0.00000 ];
TOT_LEAKRATE              (idx, [1:   2]) = [  0.00000E+00 0.00000 ];
TOT_LOSSRATE              (idx, [1:   2]) = [  1.00000E+00 0.00000 ];
TOT_FLUX                  (idx, [1:   2]) = [  1.76238E+02 0.00062 ];
TOT_RR                    (idx, [1:   2]) = [  5.89023E+01 0.00062 ];
``````
By default the normalization is set to total loss rate = 1 (for compatibility with MCNP, If I remember correctly). In a criticality source simulation, the total neutron generation rate (TOT_GENRATE) divided by the analog estimator of k-eff (ANA_KEFF) should give you the number of fission neutrons emitted per second.

Leakage rate calculated by the detector gives you the total number of neutrons escaped from the system per second. Whether or not this is equal to the leakage probability depends on what you mean by the term. Some small fraction of the escaped neutrons were born in (n,2n) reactions, and in theory, if you have completely non-absorbing system, your leakage rate can exceed the rate at which neutrons are generated in fissions. In any case, the value given by the detector should be equal to the total leakage rate (TOT_LEAKRATE) in the _res.m output (at least within statistics, there may be some small numerical difference because the estimator is different).
- Jaakko

subterfug3x
Posts: 33
Joined: Thu Oct 04, 2018 7:53 am
Security question 1: No
Security question 2: 96
Location: Tomsk, Russia

### Re: set power normalization

Jaakko Leppänen wrote:
Mon Jul 01, 2013 12:56 pm
The units for the surface current tally are in neutron/s. The value is not divided by the surface area, because there is no safe way to calculate that within the code.

In the _res.m output there is a separate section for normalized total reaction rates:

Code: Select all

``````% Normalized total reaction rates:

TOT_POWER                 (idx, [1:   2]) = [  1.21101E-11 0.00055 ];
TOT_GENRATE               (idx, [1:   2]) = [  9.18263E-01 0.00055 ];
TOT_FISSRATE              (idx, [1:   2]) = [  3.73240E-01 0.00056 ];
TOT_ABSRATE               (idx, [1:   2]) = [  1.00000E+00 0.00000 ];
TOT_LEAKRATE              (idx, [1:   2]) = [  0.00000E+00 0.00000 ];
TOT_LOSSRATE              (idx, [1:   2]) = [  1.00000E+00 0.00000 ];
TOT_FLUX                  (idx, [1:   2]) = [  1.76238E+02 0.00062 ];
TOT_RR                    (idx, [1:   2]) = [  5.89023E+01 0.00062 ];
``````
By default the normalization is set to total loss rate = 1 (for compatibility with MCNP, If I remember correctly). In a criticality source simulation, the total neutron generation rate (TOT_GENRATE) divided by the analog estimator of k-eff (ANA_KEFF) should give you the number of fission neutrons emitted per second.

Leakage rate calculated by the detector gives you the total number of neutrons escaped from the system per second. Whether or not this is equal to the leakage probability depends on what you mean by the term. Some small fraction of the escaped neutrons were born in (n,2n) reactions, and in theory, if you have completely non-absorbing system, your leakage rate can exceed the rate at which neutrons are generated in fissions. In any case, the value given by the detector should be equal to the total leakage rate (TOT_LEAKRATE) in the _res.m output (at least within statistics, there may be some small numerical difference because the estimator is different).
Hello,
I have a model of the small AmO2-Be neutron source. The problem of interest is the surface neutron current on the cylinder surface, I set the srcrate normalization (n/s), but surface detector scores more neutrons than neutron source generating.

Integral of the Cylinder Surface Detector is about 7.11E6 n/s that is equal to TOTAL_LEAKRATE and seems to be OK, but what is interesting - "Neutron balance" shows that N,XN reactions are missing in the system, capture rate is much more than the fission rate. So the main question is: Why the surface detector scores more neutrons than the neutron source generating and and TOTAL_LEAKRATE exceeds srcrate in this system?

Serpent 2.1.30 ENDF-B/VII.0

Code: Select all

``````% -----------------------------------------------------------------------------

set title "AmO2-576-2 Be Source"

% -----------------------------------------------------------------------------

% Geometry:

% -----------------------------------------------------------------------------

cell 1     0 AmBe      -1
cell 2     0 Steel       -2 1
cell 3     0 outside     2

% -----------------------------------------------------------------------------

surf 1      cylz 0.0 0.0 1.11  -2.415 2.415    % Source
surf 2      cylz 0.0 0.0 1.35  -2.900 2.900    % Container
surf 3      cylz 0.0 0.0 5.00  -10.00 10.00    % Boundary Air

% -----------------------------------------------------------------------------

% Material data:

% -----------------------------------------------------------------------------

mat AmBe       5.66E-02 rgb 166 0 0

4009.03c     9.95E-01
8016.03c     3.45E-03
8017.03c     1.31E-06
95241.03c     1.63E-03
95242.03c     6.21E-05
95243.03c     3.09E-05

mat Steel     -7.86 rgb 32 32 32
6000.03c    -0.000750
7014.03c    -0.001250
14000.03c    -0.005000
15031.03c    -0.000300
16000.03c    -0.000150
24000.03c    -0.180000
25055.03c    -0.087500
26000.03c    -0.675050
28000.03c    -0.050000

mat Air       -0.001205 rgb 255 255 255

6000.03c     0.000150
7014.03c     0.784431
8016.03c     0.210748
18040.03c    0.004671

% -----------------------------------------------------------------------------

% Calculation parameters:

% -----------------------------------------------------------------------------

% --- Geometry plot:

plot 3 500 500
plot 2 500 500

mesh 3 500 500
mesh 2 500 500

% --- Detectors
% --- Steel Container Detector
det SteelNetSideSurfaceCurrent
ds  2  0
de  EquiLet
ene EquiLet 1
0
1.00E-01
2.00E-01
3.00E-01
4.00E-01
5.00E-01
6.00E-01
7.00E-01
8.00E-01
9.00E-01
1.00E+00
1.10E+00
1.20E+00
1.30E+00
1.40E+00
1.50E+00
1.60E+00
1.70E+00
1.80E+00
1.90E+00
2.00E+00
2.10E+00
2.20E+00
2.30E+00
2.40E+00
2.50E+00
2.60E+00
2.70E+00
2.80E+00
2.90E+00
3.00E+00
3.10E+00
3.20E+00
3.30E+00
3.40E+00
3.50E+00
3.60E+00
3.70E+00
3.80E+00
3.90E+00
4.00E+00
4.10E+00
4.20E+00
4.30E+00
4.40E+00
4.50E+00
4.60E+00
4.70E+00
4.80E+00
4.90E+00
5.00E+00
5.10E+00
5.20E+00
5.30E+00
5.40E+00
5.50E+00
5.60E+00
5.70E+00
5.80E+00
5.90E+00
6.00E+00
6.10E+00
6.20E+00
6.30E+00
6.40E+00
6.50E+00
6.60E+00
6.70E+00
6.80E+00
6.90E+00
7.00E+00
7.10E+00
7.20E+00
7.30E+00
7.40E+00
7.50E+00
7.60E+00
7.70E+00
7.80E+00
7.90E+00
8.00E+00
8.10E+00
8.20E+00
8.30E+00
8.40E+00
8.50E+00
8.60E+00
8.70E+00
8.80E+00
8.90E+00
9.00E+00
9.10E+00
9.20E+00
9.30E+00
9.40E+00
9.50E+00
9.60E+00
9.70E+00
9.80E+00
9.90E+00
1.00E+01
1.01E+01
1.02E+01
1.03E+01
1.04E+01
1.05E+01
1.06E+01
1.07E+01
1.08E+01
1.09E+01
1.10E+01
1.11E+01
1.12E+01
1.13E+01
1.14E+01
1.15E+01
1.16E+01
1.17E+01
1.18E+01
1.19E+01
1.20E+01

% --- Source cylinder detector
det SourceNetSideSurfaceCurrent
ds  1  0
de  EquiLet

% --- Source

src AmBe576 n sc 1
sb 120 1.00E-01	4.55E+04
2.00E-01	7.15E+04
3.00E-01	7.16E+04
4.00E-01	7.34E+04
5.00E-01	7.69E+04
6.00E-01	8.25E+04
7.00E-01	8.61E+04
8.00E-01	8.21E+04
9.00E-01	7.40E+04
1.00E+00	6.26E+04
1.10E+00	5.68E+04
1.20E+00	5.34E+04
1.30E+00	5.01E+04
1.40E+00	4.56E+04
1.50E+00	4.09E+04
1.60E+00	3.58E+04
1.70E+00	3.33E+04
1.80E+00	3.67E+04
1.90E+00	4.70E+04
2.00E+00	5.31E+04
2.10E+00	5.65E+04
2.20E+00	5.97E+04
2.30E+00	6.14E+04
2.40E+00	6.15E+04
2.50E+00	6.12E+04
2.60E+00	6.12E+04
2.70E+00	6.69E+04
2.80E+00	7.34E+04
2.90E+00	8.40E+04
3.00E+00	1.02E+05
3.10E+00	1.22E+05
3.20E+00	1.38E+05
3.30E+00	1.43E+05
3.40E+00	1.32E+05
3.50E+00	1.23E+05
3.60E+00	1.17E+05
3.70E+00	1.13E+05
3.80E+00	1.10E+05
3.90E+00	1.08E+05
4.00E+00	1.06E+05
4.10E+00	1.05E+05
4.20E+00	1.04E+05
4.30E+00	1.03E+05
4.40E+00	1.03E+05
4.50E+00	1.04E+05
4.60E+00	1.05E+05
4.70E+00	1.07E+05
4.80E+00	1.11E+05
4.90E+00	1.15E+05
5.00E+00	1.19E+05
5.10E+00	1.16E+05
5.20E+00	1.15E+05
5.30E+00	1.13E+05
5.40E+00	1.10E+05
5.50E+00	1.04E+05
5.60E+00	1.01E+05
5.70E+00	1.03E+05
5.80E+00	1.01E+05
5.90E+00	9.24E+04
6.00E+00	8.16E+04
6.10E+00	7.08E+04
6.20E+00	6.46E+04
6.30E+00	6.07E+04
6.40E+00	5.64E+04
6.50E+00	5.94E+04
6.60E+00	5.95E+04
6.70E+00	5.65E+04
6.80E+00	5.02E+04
6.90E+00	4.45E+04
7.00E+00	4.18E+04
7.10E+00	4.17E+04
7.20E+00	4.34E+04
7.30E+00	4.65E+04
7.40E+00	5.06E+04
7.50E+00	5.51E+04
7.60E+00	5.88E+04
7.70E+00	6.13E+04
7.80E+00	6.25E+04
7.90E+00	6.25E+04
8.00E+00	6.12E+04
8.10E+00	5.87E+04
8.20E+00	5.37E+04
8.30E+00	4.61E+04
8.40E+00	3.87E+04
8.50E+00	3.31E+04
8.60E+00	2.81E+04
8.70E+00	2.35E+04
8.80E+00	1.94E+04
8.90E+00	1.60E+04
9.00E+00	1.37E+04
9.10E+00	1.25E+04
9.20E+00	1.28E+04
9.30E+00	1.46E+04
9.40E+00	1.75E+04
9.50E+00	2.09E+04
9.60E+00	2.31E+04
9.70E+00	2.37E+04
9.80E+00	2.44E+04
9.90E+00	2.44E+04
1.00E+01	2.42E+04
1.01E+01	2.36E+04
1.02E+01	2.24E+04
1.03E+01	2.06E+04
1.04E+01	1.80E+04
1.05E+01	1.52E+04
1.06E+01	1.28E+04
1.07E+01	1.08E+04
1.08E+01	8.07E+03
1.09E+01	4.98E+03
1.10E+01	1.96E+03
1.11E+01	1.13E+02
1.12E+01	0.00E+00
1.13E+01	0.00E+00
1.14E+01	0.00E+00
1.15E+01	0.00E+00
1.16E+01	0.00E+00
1.17E+01	0.00E+00
1.18E+01	0.00E+00
1.19E+01	0.00E+00
1.20E+01	0.00E+00

% --- Run Parameters

set bc 1
set srcrate 6.97E+06
set nps 10000000000

% --- Libraries:

set acelib "/home/iol4/serpent2/xs/endfb7/sss_endfb7u.xsdata"
set declib "/home/iol4/serpent2/xs/endfb7/sss_endfb7.dec"
set nfylib "/home/iol4/serpent2/xs/endfb7/sss_endfb7.nfy"``````
Output .res

Code: Select all

``````% Neutron balance (particles/weight):

BALA_SRC_NEUTRON_SRC     (idx, [1:  2])  = [ 1000000000 6.29264E+13 ];
BALA_SRC_NEUTRON_FISS    (idx, [1:  2])  = [ 876078 4.93156E+10 ];
BALA_SRC_NEUTRON_NXN     (idx, [1:  2])  = [ 0 1.56051E+12 ];
BALA_SRC_NEUTRON_VR      (idx, [1:  2])  = [ 0 0.00000E+00 ];
BALA_SRC_NEUTRON_TOT     (idx, [1:  2])  = [ 1000876078 6.45362E+13 ];

BALA_LOSS_NEUTRON_CAPT    (idx, [1:  2])  = [ 4858716 3.17041E+11 ];
BALA_LOSS_NEUTRON_FISS    (idx, [1:  2])  = [ 217222 1.28217E+10 ];
BALA_LOSS_NEUTRON_LEAK    (idx, [1:  2])  = [ 995800140 6.42063E+13 ];
BALA_LOSS_NEUTRON_CUT     (idx, [1:  2])  = [ 0 0.00000E+00 ];
BALA_LOSS_NEUTRON_ERR     (idx, [1:  2])  = [ 0 0.00000E+00 ];
BALA_LOSS_NEUTRON_TOT     (idx, [1:  2])  = [ 1000876078 6.45362E+13 ];

BALA_NEUTRON_DIFF         (idx, [1:  2])  = [ 0 0.00000E+00 ];

% Normalized total reaction rates (neutrons):

TOT_POWER                 (idx, [1:   2]) = [  4.77623E-08 0.00011 ];
TOT_POWDENS               (idx, [1:   2]) = [  0.00000E+00 0.0E+00 ];
TOT_GENRATE               (idx, [1:   2]) = [  5.43864E+03 0.00011 ];
TOT_FISSRATE              (idx, [1:   2]) = [  1.41407E+03 0.00011 ];
TOT_CAPTRATE              (idx, [1:   2]) = [  3.51070E+04 9.7E-05 ];
TOT_ABSRATE               (idx, [1:   2]) = [  3.65211E+04 9.7E-05 ];
TOT_SRCRATE               (idx, [1:   2]) = [  6.97000E+06 0.0E+00 ];
TOT_FLUX                  (idx, [1:   2]) = [  1.16941E+07 8.3E-05 ];
TOT_PHOTON_PRODRATE       (idx, [1:   4]) = [  0.00000E+00 0.0E+00  0.00000E+00 0.0E+00 ];
TOT_LEAKRATE              (idx, [1:   2]) = [  7.11177E+06 6.7E-06 ];
ALBEDO_LEAKRATE           (idx, [1:   2]) = [  0.00000E+00 0.0E+00 ];
TOT_LOSSRATE              (idx, [1:   2]) = [  7.14830E+06 6.8E-06 ];
TOT_CUTRATE               (idx, [1:   2]) = [  0.00000E+00 0.0E+00 ];
TOT_RR                    (idx, [1:   2]) = [  2.07222E+06 8.1E-05 ];
INI_FMASS                 (idx, 1)        =  0.00000E+00 ;
TOT_FMASS                 (idx, 1)        =  0.00000E+00 ;

% Criticality eigenvalues:

ANA_KEFF                  (idx, [1:   6]) = [  7.83088E-04 0.00253  0.00000E+00 0.0E+00  0.00000E+00 0.0E+00 ];
IMP_KEFF                  (idx, [1:   2]) = [  7.79692E-04 0.00011 ];
COL_KEFF                  (idx, [1:   2]) = [  7.79684E-04 0.00011 ];
ABS_KEFF                  (idx, [1:   2]) = [  7.79692E-04 0.00011 ];
ABS_KINF                  (idx, [1:   2]) = [ -3.98710E-02 6.5E-05 ];
ANA_EXT_K                 (idx, [1:  20]) = [  7.83163E-04 0.00253  6.93856E-04 0.09757  0.00000E+00 0.0E+00  0.00000E+00 0.0E+00  0.00000E+00 0.0E+00  0.00000E+00 0.0E+00  0.00000E+00 0.0E+00  0.00000E+00 0.0E+00  0.00000E+00 0.0E+00  0.00000E+00 0.0E+00 ];
SRC_MULT                  (idx, [1:   2]) = [  1.00078E+00 2.0E-06 ];
MEAN_NGEN                 (idx, [1:   2]) = [  0.00000E+00 0.0E+00 ];
PROMPT_CHAIN_LENGTH       (idx, [1:   2]) = [  1.00074E+00 5.9E-05 ];
GEOM_ALBEDO               (idx, [1:   6]) = [  1.00000E+00 0.0E+00  1.00000E+00 0.0E+00  1.00000E+00 0.0E+00 ];[code]
``````
Output Surface Detector

Code: Select all

``````DETSteelNetSideSurfaceCurrent = [
1    1    1    1    1    1    1    1    1    1 -3.22208E+04 0.00075
2    2    1    1    1    1    1    1    1    1 -1.10037E+05 0.00031
3    3    1    1    1    1    1    1    1    1 -1.10414E+05 0.00031
4    4    1    1    1    1    1    1    1    1 -1.15345E+05 0.00033
5    5    1    1    1    1    1    1    1    1 -1.17906E+05 0.00031
6    6    1    1    1    1    1    1    1    1 -1.18095E+05 0.00031
7    7    1    1    1    1    1    1    1    1 -1.09001E+05 0.00031
8    8    1    1    1    1    1    1    1    1 -1.08014E+05 0.00031
9    9    1    1    1    1    1    1    1    1 -9.39765E+04 0.00032
10   10    1    1    1    1    1    1    1    1 -8.31460E+04 0.00033
11   11    1    1    1    1    1    1    1    1 -7.55036E+04 0.00033
12   12    1    1    1    1    1    1    1    1 -7.16782E+04 0.00033
13   13    1    1    1    1    1    1    1    1 -6.71975E+04 0.00033
14   14    1    1    1    1    1    1    1    1 -6.24841E+04 0.00033
15   15    1    1    1    1    1    1    1    1 -5.78912E+04 0.00036
16   16    1    1    1    1    1    1    1    1 -5.35394E+04 0.00036
17   17    1    1    1    1    1    1    1    1 -5.20052E+04 0.00037
18   18    1    1    1    1    1    1    1    1 -5.51503E+04 0.00035
19   19    1    1    1    1    1    1    1    1 -6.32975E+04 0.00035
20   20    1    1    1    1    1    1    1    1 -6.76044E+04 0.00034
21   21    1    1    1    1    1    1    1    1 -6.96231E+04 0.00032
22   22    1    1    1    1    1    1    1    1 -7.18518E+04 0.00035
23   23    1    1    1    1    1    1    1    1 -7.26327E+04 0.00033
24   24    1    1    1    1    1    1    1    1 -7.18254E+04 0.00035
25   25    1    1    1    1    1    1    1    1 -7.12004E+04 0.00033
26   26    1    1    1    1    1    1    1    1 -7.12707E+04 0.00034
27   27    1    1    1    1    1    1    1    1 -7.52085E+04 0.00034
28   28    1    1    1    1    1    1    1    1 -7.74675E+04 0.00035
29   29    1    1    1    1    1    1    1    1 -8.66540E+04 0.00033
30   30    1    1    1    1    1    1    1    1 -1.02285E+05 0.00033
31   31    1    1    1    1    1    1    1    1 -1.17841E+05 0.00034
32   32    1    1    1    1    1    1    1    1 -1.30137E+05 0.00034
33   33    1    1    1    1    1    1    1    1 -1.33174E+05 0.00034
34   34    1    1    1    1    1    1    1    1 -1.23716E+05 0.00033
35   35    1    1    1    1    1    1    1    1 -1.15364E+05 0.00035
36   36    1    1    1    1    1    1    1    1 -1.10225E+05 0.00034
37   37    1    1    1    1    1    1    1    1 -1.06423E+05 0.00033
38   38    1    1    1    1    1    1    1    1 -1.03573E+05 0.00034
39   39    1    1    1    1    1    1    1    1 -1.01605E+05 0.00034
40   40    1    1    1    1    1    1    1    1 -9.99315E+04 0.00033
41   41    1    1    1    1    1    1    1    1 -9.91983E+04 0.00034
42   42    1    1    1    1    1    1    1    1 -9.83510E+04 0.00034
43   43    1    1    1    1    1    1    1    1 -9.76010E+04 0.00035
44   44    1    1    1    1    1    1    1    1 -9.68112E+04 0.00034
45   45    1    1    1    1    1    1    1    1 -9.79098E+04 0.00035
46   46    1    1    1    1    1    1    1    1 -9.85423E+04 0.00033
47   47    1    1    1    1    1    1    1    1 -1.00308E+05 0.00034
48   48    1    1    1    1    1    1    1    1 -1.03681E+05 0.00034
49   49    1    1    1    1    1    1    1    1 -1.06885E+05 0.00034
50   50    1    1    1    1    1    1    1    1 -1.09750E+05 0.00032
51   51    1    1    1    1    1    1    1    1 -1.06896E+05 0.00034
52   52    1    1    1    1    1    1    1    1 -1.05636E+05 0.00035
53   53    1    1    1    1    1    1    1    1 -1.03770E+05 0.00035
54   54    1    1    1    1    1    1    1    1 -1.00659E+05 0.00034
55   55    1    1    1    1    1    1    1    1 -9.54649E+04 0.00034
56   56    1    1    1    1    1    1    1    1 -9.26519E+04 0.00035
57   57    1    1    1    1    1    1    1    1 -9.37745E+04 0.00034
58   58    1    1    1    1    1    1    1    1 -9.15858E+04 0.00034
59   59    1    1    1    1    1    1    1    1 -8.36768E+04 0.00035
60   60    1    1    1    1    1    1    1    1 -7.40379E+04 0.00036
61   61    1    1    1    1    1    1    1    1 -6.46496E+04 0.00032
62   62    1    1    1    1    1    1    1    1 -5.92332E+04 0.00033
63   63    1    1    1    1    1    1    1    1 -5.57631E+04 0.00036
64   64    1    1    1    1    1    1    1    1 -5.20414E+04 0.00034
65   65    1    1    1    1    1    1    1    1 -5.44682E+04 0.00035
66   66    1    1    1    1    1    1    1    1 -5.42187E+04 0.00037
67   67    1    1    1    1    1    1    1    1 -5.14537E+04 0.00034
68   68    1    1    1    1    1    1    1    1 -4.59823E+04 0.00034
69   69    1    1    1    1    1    1    1    1 -4.11350E+04 0.00035
70   70    1    1    1    1    1    1    1    1 -3.89080E+04 0.00035
71   71    1    1    1    1    1    1    1    1 -3.90131E+04 0.00033
72   72    1    1    1    1    1    1    1    1 -4.05444E+04 0.00034
73   73    1    1    1    1    1    1    1    1 -4.33685E+04 0.00035
74   74    1    1    1    1    1    1    1    1 -4.69353E+04 0.00035
75   75    1    1    1    1    1    1    1    1 -5.08253E+04 0.00035
76   76    1    1    1    1    1    1    1    1 -5.39446E+04 0.00034
77   77    1    1    1    1    1    1    1    1 -5.58884E+04 0.00037
78   78    1    1    1    1    1    1    1    1 -5.67503E+04 0.00036
79   79    1    1    1    1    1    1    1    1 -5.64731E+04 0.00034
80   80    1    1    1    1    1    1    1    1 -5.50787E+04 0.00036
81   81    1    1    1    1    1    1    1    1 -5.25870E+04 0.00034
82   82    1    1    1    1    1    1    1    1 -4.79280E+04 0.00033
83   83    1    1    1    1    1    1    1    1 -4.11500E+04 0.00035
84   84    1    1    1    1    1    1    1    1 -3.45182E+04 0.00035
85   85    1    1    1    1    1    1    1    1 -2.95604E+04 0.00035
86   86    1    1    1    1    1    1    1    1 -2.51476E+04 0.00035
87   87    1    1    1    1    1    1    1    1 -2.11127E+04 0.00036
88   88    1    1    1    1    1    1    1    1 -1.75524E+04 0.00035
89   89    1    1    1    1    1    1    1    1 -1.46313E+04 0.00035
90   90    1    1    1    1    1    1    1    1 -1.27073E+04 0.00035
91   91    1    1    1    1    1    1    1    1 -1.17711E+04 0.00036
92   92    1    1    1    1    1    1    1    1 -1.21351E+04 0.00035
93   93    1    1    1    1    1    1    1    1 -1.37669E+04 0.00035
94   94    1    1    1    1    1    1    1    1 -1.63279E+04 0.00034
95   95    1    1    1    1    1    1    1    1 -1.92409E+04 0.00034
96   96    1    1    1    1    1    1    1    1 -2.10910E+04 0.00035
97   97    1    1    1    1    1    1    1    1 -2.15695E+04 0.00035
98   98    1    1    1    1    1    1    1    1 -2.21215E+04 0.00035
99   99    1    1    1    1    1    1    1    1 -2.20355E+04 0.00036
100  100    1    1    1    1    1    1    1    1 -2.17879E+04 0.00035
101  101    1    1    1    1    1    1    1    1 -2.11509E+04 0.00036
102  102    1    1    1    1    1    1    1    1 -2.00105E+04 0.00036
103  103    1    1    1    1    1    1    1    1 -1.82975E+04 0.00034
104  104    1    1    1    1    1    1    1    1 -1.59522E+04 0.00036
105  105    1    1    1    1    1    1    1    1 -1.34398E+04 0.00035
106  106    1    1    1    1    1    1    1    1 -1.12707E+04 0.00036
107  107    1    1    1    1    1    1    1    1 -9.44520E+03 0.00034
108  108    1    1    1    1    1    1    1    1 -7.01123E+03 0.00037
109  109    1    1    1    1    1    1    1    1 -4.29435E+03 0.00037
110  110    1    1    1    1    1    1    1    1 -1.67553E+03 0.00035
111  111    1    1    1    1    1    1    1    1 -9.70160E+01 0.00083
112  112    1    1    1    1    1    1    1    1 -4.67513E-01 0.13307
113  113    1    1    1    1    1    1    1    1 -4.27220E-01 0.14986
114  114    1    1    1    1    1    1    1    1 -4.97337E-01 0.12958
115  115    1    1    1    1    1    1    1    1 -4.85793E-01 0.15057
116  116    1    1    1    1    1    1    1    1 -3.12278E-01 0.16306
117  117    1    1    1    1    1    1    1    1 -3.18202E-01 0.16064
118  118    1    1    1    1    1    1    1    1 -2.70004E-01 0.16493
119  119    1    1    1    1    1    1    1    1 -3.34639E-01 0.16106
120  120    1    1    1    1    1    1    1    1 -2.72824E-01 0.18281
];

DETSteelNetSideSurfaceCurrentE = [
0.00000E+00  1.00000E-01  5.00000E-02
1.00000E-01  2.00000E-01  1.50000E-01
2.00000E-01  3.00000E-01  2.50000E-01
3.00000E-01  4.00000E-01  3.50000E-01
4.00000E-01  5.00000E-01  4.50000E-01
5.00000E-01  6.00000E-01  5.50000E-01
6.00000E-01  7.00000E-01  6.50000E-01
7.00000E-01  8.00000E-01  7.50000E-01
8.00000E-01  9.00000E-01  8.50000E-01
9.00000E-01  1.00000E+00  9.50000E-01
1.00000E+00  1.10000E+00  1.05000E+00
1.10000E+00  1.20000E+00  1.15000E+00
1.20000E+00  1.30000E+00  1.25000E+00
1.30000E+00  1.40000E+00  1.35000E+00
1.40000E+00  1.50000E+00  1.45000E+00
1.50000E+00  1.60000E+00  1.55000E+00
1.60000E+00  1.70000E+00  1.65000E+00
1.70000E+00  1.80000E+00  1.75000E+00
1.80000E+00  1.90000E+00  1.85000E+00
1.90000E+00  2.00000E+00  1.95000E+00
2.00000E+00  2.10000E+00  2.05000E+00
2.10000E+00  2.20000E+00  2.15000E+00
2.20000E+00  2.30000E+00  2.25000E+00
2.30000E+00  2.40000E+00  2.35000E+00
2.40000E+00  2.50000E+00  2.45000E+00
2.50000E+00  2.60000E+00  2.55000E+00
2.60000E+00  2.70000E+00  2.65000E+00
2.70000E+00  2.80000E+00  2.75000E+00
2.80000E+00  2.90000E+00  2.85000E+00
2.90000E+00  3.00000E+00  2.95000E+00
3.00000E+00  3.10000E+00  3.05000E+00
3.10000E+00  3.20000E+00  3.15000E+00
3.20000E+00  3.30000E+00  3.25000E+00
3.30000E+00  3.40000E+00  3.35000E+00
3.40000E+00  3.50000E+00  3.45000E+00
3.50000E+00  3.60000E+00  3.55000E+00
3.60000E+00  3.70000E+00  3.65000E+00
3.70000E+00  3.80000E+00  3.75000E+00
3.80000E+00  3.90000E+00  3.85000E+00
3.90000E+00  4.00000E+00  3.95000E+00
4.00000E+00  4.10000E+00  4.05000E+00
4.10000E+00  4.20000E+00  4.15000E+00
4.20000E+00  4.30000E+00  4.25000E+00
4.30000E+00  4.40000E+00  4.35000E+00
4.40000E+00  4.50000E+00  4.45000E+00
4.50000E+00  4.60000E+00  4.55000E+00
4.60000E+00  4.70000E+00  4.65000E+00
4.70000E+00  4.80000E+00  4.75000E+00
4.80000E+00  4.90000E+00  4.85000E+00
4.90000E+00  5.00000E+00  4.95000E+00
5.00000E+00  5.10000E+00  5.05000E+00
5.10000E+00  5.20000E+00  5.15000E+00
5.20000E+00  5.30000E+00  5.25000E+00
5.30000E+00  5.40000E+00  5.35000E+00
5.40000E+00  5.50000E+00  5.45000E+00
5.50000E+00  5.60000E+00  5.55000E+00
5.60000E+00  5.70000E+00  5.65000E+00
5.70000E+00  5.80000E+00  5.75000E+00
5.80000E+00  5.90000E+00  5.85000E+00
5.90000E+00  6.00000E+00  5.95000E+00
6.00000E+00  6.10000E+00  6.05000E+00
6.10000E+00  6.20000E+00  6.15000E+00
6.20000E+00  6.30000E+00  6.25000E+00
6.30000E+00  6.40000E+00  6.35000E+00
6.40000E+00  6.50000E+00  6.45000E+00
6.50000E+00  6.60000E+00  6.55000E+00
6.60000E+00  6.70000E+00  6.65000E+00
6.70000E+00  6.80000E+00  6.75000E+00
6.80000E+00  6.90000E+00  6.85000E+00
6.90000E+00  7.00000E+00  6.95000E+00
7.00000E+00  7.10000E+00  7.05000E+00
7.10000E+00  7.20000E+00  7.15000E+00
7.20000E+00  7.30000E+00  7.25000E+00
7.30000E+00  7.40000E+00  7.35000E+00
7.40000E+00  7.50000E+00  7.45000E+00
7.50000E+00  7.60000E+00  7.55000E+00
7.60000E+00  7.70000E+00  7.65000E+00
7.70000E+00  7.80000E+00  7.75000E+00
7.80000E+00  7.90000E+00  7.85000E+00
7.90000E+00  8.00000E+00  7.95000E+00
8.00000E+00  8.10000E+00  8.05000E+00
8.10000E+00  8.20000E+00  8.15000E+00
8.20000E+00  8.30000E+00  8.25000E+00
8.30000E+00  8.40000E+00  8.35000E+00
8.40000E+00  8.50000E+00  8.45000E+00
8.50000E+00  8.60000E+00  8.55000E+00
8.60000E+00  8.70000E+00  8.65000E+00
8.70000E+00  8.80000E+00  8.75000E+00
8.80000E+00  8.90000E+00  8.85000E+00
8.90000E+00  9.00000E+00  8.95000E+00
9.00000E+00  9.10000E+00  9.05000E+00
9.10000E+00  9.20000E+00  9.15000E+00
9.20000E+00  9.30000E+00  9.25000E+00
9.30000E+00  9.40000E+00  9.35000E+00
9.40000E+00  9.50000E+00  9.45000E+00
9.50000E+00  9.60000E+00  9.55000E+00
9.60000E+00  9.70000E+00  9.65000E+00
9.70000E+00  9.80000E+00  9.75000E+00
9.80000E+00  9.90000E+00  9.85000E+00
9.90000E+00  1.00000E+01  9.95000E+00
1.00000E+01  1.01000E+01  1.00500E+01
1.01000E+01  1.02000E+01  1.01500E+01
1.02000E+01  1.03000E+01  1.02500E+01
1.03000E+01  1.04000E+01  1.03500E+01
1.04000E+01  1.05000E+01  1.04500E+01
1.05000E+01  1.06000E+01  1.05500E+01
1.06000E+01  1.07000E+01  1.06500E+01
1.07000E+01  1.08000E+01  1.07500E+01
1.08000E+01  1.09000E+01  1.08500E+01
1.09000E+01  1.10000E+01  1.09500E+01
1.10000E+01  1.11000E+01  1.10500E+01
1.11000E+01  1.12000E+01  1.11500E+01
1.12000E+01  1.13000E+01  1.12500E+01
1.13000E+01  1.14000E+01  1.13500E+01
1.14000E+01  1.15000E+01  1.14500E+01
1.15000E+01  1.16000E+01  1.15500E+01
1.16000E+01  1.17000E+01  1.16500E+01
1.17000E+01  1.18000E+01  1.17500E+01
1.18000E+01  1.19000E+01  1.18500E+01
1.19000E+01  1.20000E+01  1.19500E+01
];

DETSourceNetSideSurfaceCurrent = [
1    1    1    1    1    1    1    1    1    1 -2.81874E+04 0.00081
2    2    1    1    1    1    1    1    1    1 -1.03049E+05 0.00032
3    3    1    1    1    1    1    1    1    1 -1.02877E+05 0.00032
4    4    1    1    1    1    1    1    1    1 -1.04796E+05 0.00035
5    5    1    1    1    1    1    1    1    1 -1.09505E+05 0.00033
6    6    1    1    1    1    1    1    1    1 -1.09683E+05 0.00033
7    7    1    1    1    1    1    1    1    1 -9.86126E+04 0.00034
8    8    1    1    1    1    1    1    1    1 -9.81687E+04 0.00033
9    9    1    1    1    1    1    1    1    1 -8.55252E+04 0.00033
10   10    1    1    1    1    1    1    1    1 -7.40710E+04 0.00036
11   11    1    1    1    1    1    1    1    1 -6.67774E+04 0.00036
12   12    1    1    1    1    1    1    1    1 -6.22431E+04 0.00036
13   13    1    1    1    1    1    1    1    1 -5.84065E+04 0.00036
14   14    1    1    1    1    1    1    1    1 -5.40244E+04 0.00036
15   15    1    1    1    1    1    1    1    1 -4.99934E+04 0.00037
16   16    1    1    1    1    1    1    1    1 -4.59656E+04 0.00038
17   17    1    1    1    1    1    1    1    1 -4.46066E+04 0.00037
18   18    1    1    1    1    1    1    1    1 -4.79235E+04 0.00037
19   19    1    1    1    1    1    1    1    1 -5.66807E+04 0.00036
20   20    1    1    1    1    1    1    1    1 -6.18322E+04 0.00035
21   21    1    1    1    1    1    1    1    1 -6.45830E+04 0.00034
22   22    1    1    1    1    1    1    1    1 -6.72067E+04 0.00035
23   23    1    1    1    1    1    1    1    1 -6.82631E+04 0.00034
24   24    1    1    1    1    1    1    1    1 -6.78951E+04 0.00035
25   25    1    1    1    1    1    1    1    1 -6.78672E+04 0.00035
26   26    1    1    1    1    1    1    1    1 -6.87823E+04 0.00035
27   27    1    1    1    1    1    1    1    1 -7.30563E+04 0.00035
28   28    1    1    1    1    1    1    1    1 -7.58129E+04 0.00036
29   29    1    1    1    1    1    1    1    1 -8.55931E+04 0.00034
30   30    1    1    1    1    1    1    1    1 -1.01786E+05 0.00033
31   31    1    1    1    1    1    1    1    1 -1.19220E+05 0.00034
32   32    1    1    1    1    1    1    1    1 -1.33002E+05 0.00034
33   33    1    1    1    1    1    1    1    1 -1.36924E+05 0.00034
34   34    1    1    1    1    1    1    1    1 -1.26990E+05 0.00032
35   35    1    1    1    1    1    1    1    1 -1.18873E+05 0.00034
36   36    1    1    1    1    1    1    1    1 -1.13253E+05 0.00033
37   37    1    1    1    1    1    1    1    1 -1.09579E+05 0.00032
38   38    1    1    1    1    1    1    1    1 -1.06794E+05 0.00034
39   39    1    1    1    1    1    1    1    1 -1.04898E+05 0.00033
40   40    1    1    1    1    1    1    1    1 -1.03174E+05 0.00033
41   41    1    1    1    1    1    1    1    1 -1.02543E+05 0.00034
42   42    1    1    1    1    1    1    1    1 -1.02058E+05 0.00034
43   43    1    1    1    1    1    1    1    1 -1.01046E+05 0.00034
44   44    1    1    1    1    1    1    1    1 -1.00522E+05 0.00033
45   45    1    1    1    1    1    1    1    1 -1.01365E+05 0.00034
46   46    1    1    1    1    1    1    1    1 -1.02385E+05 0.00033
47   47    1    1    1    1    1    1    1    1 -1.04214E+05 0.00033
48   48    1    1    1    1    1    1    1    1 -1.07893E+05 0.00033
49   49    1    1    1    1    1    1    1    1 -1.11560E+05 0.00033
50   50    1    1    1    1    1    1    1    1 -1.14766E+05 0.00031
51   51    1    1    1    1    1    1    1    1 -1.11665E+05 0.00033
52   52    1    1    1    1    1    1    1    1 -1.10580E+05 0.00035
53   53    1    1    1    1    1    1    1    1 -1.08595E+05 0.00034
54   54    1    1    1    1    1    1    1    1 -1.05608E+05 0.00033
55   55    1    1    1    1    1    1    1    1 -9.99786E+04 0.00033
56   56    1    1    1    1    1    1    1    1 -9.70295E+04 0.00034
57   57    1    1    1    1    1    1    1    1 -9.85468E+04 0.00034
58   58    1    1    1    1    1    1    1    1 -9.63359E+04 0.00034
59   59    1    1    1    1    1    1    1    1 -8.81193E+04 0.00035
60   60    1    1    1    1    1    1    1    1 -7.79629E+04 0.00035
61   61    1    1    1    1    1    1    1    1 -6.79729E+04 0.00031
62   62    1    1    1    1    1    1    1    1 -6.21762E+04 0.00032
63   63    1    1    1    1    1    1    1    1 -5.85064E+04 0.00034
64   64    1    1    1    1    1    1    1    1 -5.45144E+04 0.00033
65   65    1    1    1    1    1    1    1    1 -5.71726E+04 0.00034
66   66    1    1    1    1    1    1    1    1 -5.70384E+04 0.00036
67   67    1    1    1    1    1    1    1    1 -5.41728E+04 0.00034
68   68    1    1    1    1    1    1    1    1 -4.83677E+04 0.00033
69   69    1    1    1    1    1    1    1    1 -4.31596E+04 0.00034
70   70    1    1    1    1    1    1    1    1 -4.07809E+04 0.00034
71   71    1    1    1    1    1    1    1    1 -4.08124E+04 0.00032
72   72    1    1    1    1    1    1    1    1 -4.24391E+04 0.00033
73   73    1    1    1    1    1    1    1    1 -4.54094E+04 0.00034
74   74    1    1    1    1    1    1    1    1 -4.92148E+04 0.00034
75   75    1    1    1    1    1    1    1    1 -5.33573E+04 0.00034
76   76    1    1    1    1    1    1    1    1 -5.67065E+04 0.00033
77   77    1    1    1    1    1    1    1    1 -5.88242E+04 0.00037
78   78    1    1    1    1    1    1    1    1 -5.97701E+04 0.00035
79   79    1    1    1    1    1    1    1    1 -5.95517E+04 0.00033
80   80    1    1    1    1    1    1    1    1 -5.81254E+04 0.00035
81   81    1    1    1    1    1    1    1    1 -5.55628E+04 0.00034
82   82    1    1    1    1    1    1    1    1 -5.06863E+04 0.00032
83   83    1    1    1    1    1    1    1    1 -4.34948E+04 0.00033
84   84    1    1    1    1    1    1    1    1 -3.65034E+04 0.00035
85   85    1    1    1    1    1    1    1    1 -3.12451E+04 0.00034
86   86    1    1    1    1    1    1    1    1 -2.65611E+04 0.00035
87   87    1    1    1    1    1    1    1    1 -2.22881E+04 0.00035
88   88    1    1    1    1    1    1    1    1 -1.85060E+04 0.00035
89   89    1    1    1    1    1    1    1    1 -1.53914E+04 0.00034
90   90    1    1    1    1    1    1    1    1 -1.33331E+04 0.00035
91   91    1    1    1    1    1    1    1    1 -1.23102E+04 0.00035
92   92    1    1    1    1    1    1    1    1 -1.26710E+04 0.00035
93   93    1    1    1    1    1    1    1    1 -1.43830E+04 0.00034
94   94    1    1    1    1    1    1    1    1 -1.70839E+04 0.00034
95   95    1    1    1    1    1    1    1    1 -2.01940E+04 0.00033
96   96    1    1    1    1    1    1    1    1 -2.21802E+04 0.00034
97   97    1    1    1    1    1    1    1    1 -2.26907E+04 0.00034
98   98    1    1    1    1    1    1    1    1 -2.32926E+04 0.00035
99   99    1    1    1    1    1    1    1    1 -2.32157E+04 0.00035
100  100    1    1    1    1    1    1    1    1 -2.29703E+04 0.00035
101  101    1    1    1    1    1    1    1    1 -2.23148E+04 0.00035
102  102    1    1    1    1    1    1    1    1 -2.11325E+04 0.00035
103  103    1    1    1    1    1    1    1    1 -1.93447E+04 0.00033
104  104    1    1    1    1    1    1    1    1 -1.68763E+04 0.00035
105  105    1    1    1    1    1    1    1    1 -1.42199E+04 0.00034
106  106    1    1    1    1    1    1    1    1 -1.19373E+04 0.00035
107  107    1    1    1    1    1    1    1    1 -1.00204E+04 0.00033
108  108    1    1    1    1    1    1    1    1 -7.45178E+03 0.00035
109  109    1    1    1    1    1    1    1    1 -4.57977E+03 0.00036
110  110    1    1    1    1    1    1    1    1 -1.79471E+03 0.00034
111  111    1    1    1    1    1    1    1    1 -1.03989E+02 0.00080
112  112    1    1    1    1    1    1    1    1 -4.62140E-01 0.13298
113  113    1    1    1    1    1    1    1    1 -4.36415E-01 0.14769
114  114    1    1    1    1    1    1    1    1 -5.44984E-01 0.12398
115  115    1    1    1    1    1    1    1    1 -5.31346E-01 0.14238
116  116    1    1    1    1    1    1    1    1 -3.39211E-01 0.15867
117  117    1    1    1    1    1    1    1    1 -3.44551E-01 0.15368
118  118    1    1    1    1    1    1    1    1 -2.80356E-01 0.16069
119  119    1    1    1    1    1    1    1    1 -3.48695E-01 0.15687
120  120    1    1    1    1    1    1    1    1 -2.78995E-01 0.17952
];``````
Same task calculated with the MCNPX 2.7e ENDF-B/VII.0 showed the same neutron current on the surfaces. But the neutron balance is pretty different.

Code: Select all

``````      MCNPX Visual Editor Version X_25                                                     probid =   06/16/19 13:50:03

neutron creation       tracks   weight        energy            neutron loss           tracks   weight        energy
(per source particle)                                           (per source particle)

source            10000000    1.0000E+00    4.5067E+00          escape            10260482    1.0200E+00    4.2327E+00
nucl. interaction        0    0.            0.                  energy cutoff            0    0.            0.
particle decay           0    0.            0.                  time cutoff              0    0.            0.
weight window            0    0.            0.                  weight window            0    0.            0.
cell importance          0    0.            0.                  cell importance          0    0.            0.
weight cutoff            0    0.            0.                  weight cutoff            7    9.4951E-08    8.5605E-11
energy importance        0    0.            0.                  energy importance        0    0.            0.
dxtran                   0    0.            0.                  dxtran                   0    0.            0.
forced collisions        0    0.            0.                  forced collisions        0    0.            0.
exp. transform           0    0.            0.                  exp. transform           0    0.            0.
upscattering             0    0.            2.3521E-15          downscattering           0    0.            1.7298E-01
*photonuclear            0    0.            0.                  capture                  0    5.1782E-03    2.5749E-02
(n,xn)              509564    4.9317E-02    6.2698E-02          loss to (n,xn)      254782    2.4659E-02    1.3874E-01
prompt fission        7705    7.6441E-04    1.7659E-03          loss to fission       2005    1.9873E-04    9.8102E-04
delayed fission          7    6.9223E-07    5.7866E-07          nucl. interaction        0    0.            0.
*prompt photofiss        0    0.            0.                  particle decay           0    0.            0.
tabular boundary         0    0.            0.                  tabular boundary         0    0.            0.
tabular sampling         0    0.            0.
total         10517276    1.0501E+00    4.5711E+00              total         10517276    1.0501E+00    4.5711E+00

number of neutrons banked                  260489        average time of (shakes)              cutoffs
neutron tracks per source particle     1.0517E+00          escape            1.2576E+02          tco   1.0000E+34
neutron collisions per source particle 3.0433E-01          capture           5.8692E+01          eco   0.0000E+00
total neutron collisions                  3043343          capture or escape 1.2542E+02          wc1  -5.0000E-01
net multiplication              1.0252E+00 0.0000          any termination   1.2245E+02          wc2  -2.5000E-01

computer time so far in this run     0.67 minutes            maximum number ever in bank         5
computer time in mcrun               0.64 minutes            bank overflows to backup file       0
source particles per minute            1.5688E+07            dynamic storage         0 words,          0 bytes.
random numbers generated                127480902            most random numbers used was      178 in history  8966337

range of sampled source weights = 1.0000E+00 to 1.0000E+00

source efficiency = 0.6668  in cell     1
1source distribution frequency tables                                                                   print table 170

the expected values below do not include the effect of the rejection loop which samples position.``````
- Igor

Jaakko Leppänen
Posts: 2442
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

### Re: set power normalization

Serpent uses implicit (n,xn) by default, so even though these reactions are not shown in neutron balance, they appear in the weight balance. The treatment can be changed to analog using the set impl card.

Also note that the balance data refers to simulated particles and weight. Total source rate, fission rate, capture rate, etc. are normalized values.
- Jaakko