Automated burnup sequence example 1

From Serpent Wiki
Jump to: navigation, search

Input

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

set title "Serpent GC generation example"

% - 17x17 PWR fuel assembly, adopted with some modifications from the
%   MIT BEAVRS Benchmark: http://crpg.mit.edu/pub/beavrs
% - Single fuel type, no burnable absorber
% - Pin lattice modeled applying universe symmetry (45 degree sector)
%   --> number of unique pin positions reduced to 39
% - Automated division into depletion zones, each pin handled separately,
%   material volumes provided using the mvol card
% - Single assembly with reflective boundary conditions --> homogeneous
%   diffusion flux solver will not be run
% - B1 leakage correction run (WIMS 172-group structure)
% - Branches to : - low and high boron concentration
%                 - low and high coolant temperature (and density)
%                 - low and high fuel temperature (and density)
%                 - control rod insertion
% - Adjustment of boron concentrations using the repm variation
% - Adjustment of temperatures using the stp variation
% - Thermal expansion of fuel pins using the repu variation (pin and assembly
%   pitch not expanded)
% - Control rod insertion using the repu variation
% - All branch combinations included (3x3x3x2 = 54 combinations in matrix)

% NOTE: The values of temperature-adjusted densities and dimensions are
%       for demonstration only, may not reflect actual physics.

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

% ----- Geometry definitions --------------------------------------------------

% NOTE: Some of the definitions are for the branches

% --- Fuel pin at nominal conditions: (fuel at 900K, cladding 592K)

pin FUE
fuel     0.39480
helium   0.40137
zirc4    0.45870
cool

% --- Fuel pin at cold conditions: (fuel at 600K, cladding 556K)

pin FUEC
fuel     0.39350
helium   0.40121
zirc4    0.45852
cool

% --- Fuel pin at hot conditions: (fuel at 1200K, cladding 628K)

pin FUEH
fuel     0.39609
helium   0.40152
zirc4    0.45888
cool

% --- Empty guide tube:

pin GTU
cool     0.56134
zirc4    0.60198
cool

% --- Guide tube with control rod inserted: 

pin ROD
AIC      0.43310
helium   0.43688
steel    0.48387
cool     0.56134
zirc4    0.60198
cool

% --- Bare instrumentation tube:

pin ITU
cool     0.56134
zirc4    0.60198
cool

% --- Empty lattice cell:

pin XXX
void

% --- Assembly lattice:

% NOTE: Empty positions needed because of 45 deg symmetry

lat LAT1 1 5.0394 5.0394 9 9 1.25984
                   
ITU XXX XXX XXX XXX XXX XXX XXX XXX  
FUE FUE XXX XXX XXX XXX XXX XXX XXX  
FUE FUE FUE XXX XXX XXX XXX XXX XXX  
GTU FUE FUE GTU XXX XXX XXX XXX XXX  
FUE FUE FUE FUE FUE XXX XXX XXX XXX  
FUE FUE FUE FUE FUE GTU XXX XXX XXX  
GTU FUE FUE GTU FUE FUE FUE XXX XXX  
FUE FUE FUE FUE FUE FUE FUE FUE XXX  
FUE FUE FUE FUE FUE FUE FUE FUE FUE  

% --- Apply 45 degree symmetry:

set usym LAT1 z 2 0.0 0.0 45 45

% --- Assembly cell:

surf S01 sqc 0.0 0.0 10.70864

cell C01 A01 fill LAT1  -S01
cell C02 A01 cool        S01

% --- Global geometry:

surf B01 sqc 0.0 0.0  10.75182

cell G01 0 fill A01  -B01
cell G02 0 outside    B01

% --- Reflective boundary conditions:

set bc 2  

% --- Material subdivision: (handle each pin separately)

div fuel sep 1

% --- Volumes of sub-zones: (needed only for burnup calculation)

set mvol

fuel   1  1.95868
fuel   2  1.95868
fuel   3  1.95868
fuel   4  3.91737
fuel   5  1.95868
fuel   6  3.91737
fuel   7  3.91737
fuel   8  1.95868
fuel   9  3.91737
fuel  10  3.91737
fuel  11  3.91737
fuel  12  1.95868
fuel  13  1.95868
fuel  14  3.91737
fuel  15  3.91737
fuel  16  3.91737
fuel  17  3.91737
fuel  18  3.91737
fuel  19  3.91737
fuel  20  3.91737
fuel  21  3.91737
fuel  22  1.95868
fuel  23  1.95868
fuel  24  3.91737
fuel  25  3.91737
fuel  26  3.91737
fuel  27  3.91737
fuel  28  3.91737
fuel  29  3.91737
fuel  30  1.95868
fuel  31  1.95868
fuel  32  3.91737
fuel  33  3.91737
fuel  34  3.91737
fuel  35  3.91737
fuel  36  3.91737
fuel  37  3.91737
fuel  38  3.91737
fuel  39  1.95868

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

% ----- Material definitions --------------------------------------------------

% --- Fuel with 3.4 wt-% U-235: (at nominal 900K temperature) 

mat fuel  -10.2013  burn 1
 8016.09c   6.65955E-01
 8017.09c   2.53687E-04
92234.09c   8.67553E-05
92235.09c   1.07981E-02
92238.09c   3.22907E-01

% --- Cladding: (at nominal 592K temperature approximated as 600K)

mat zirc4  -6.5071
 8016.06c   7.08020E-03
24000.06c   1.74646E-03
26000.06c   3.41481E-03
40000.06c   9.76666E-01
50000.06c   1.10918E-02

% --- Steel:

mat steel  -7.7936
 6000.06c   1.15167E-02
14000.06c   5.41767E-03
25055.06c   1.33445E-02
15031.06c   4.46581E-04
42000.06c   2.54632E-03
26000.06c   9.61544E-01
28000.06c   5.18479E-03

% --- Gas gap:

mat helium -1.5981E-03
 2004.06c   1.00000

% --- Control rod absorber

mat AIC   -10.1608
47107.06c   4.19366E-01
47109.06c   3.89611E-01
49113.06c   6.11337E-03
49115.06c   1.36391E-01
48106.06c   6.06488E-04
48108.06c   4.31810E-04
48110.06c   6.05989E-03
48111.06c   6.21036E-03
48112.06c   1.17074E-02
48113.06c   5.92903E-03
48114.06c   1.39393E-02
48116.06c   3.63404E-03

% --- Water: (at nominal 550 K temperature)

mat water  -0.76973  tmp 550 moder lwtr 1001
 1001.03c   0.66667
 8016.03c   0.33333

therm lwtr lwe7.10t

% --- Natural boron:

mat boron   1.00000  tmp 550
 5010.03c   0.19900
 5011.03c   0.80100

% --- Coolant: (at nominal state with 500 ppm boron)

mix cool
water      -0.99950
boron      -500E-6

% --- Coolant with without boron:

mix cool_0B
water      -1.00000

% --- Coolant with 1000 ppm boron:

mix cool_1000B
water      -0.99900
boron      -1000E-6

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

% ----- Group constant generation ---------------------------------------------

% --- Universe for group constant generation:

set gcu 0

% --- Few-group structure: (two groups, boundary et 0.625 eV)

set nfg 2 0.625E-6

% --- Run B1 leakage correction: (WIMS 172-group structure, solution at end)

set fum wms172 2

% --- ADF calculation: (super-imposed surface and full symmetry)

surf ADF sqc 0.0 0.0  10.75182
set adf 0 ADF full 

% --- Calculate poison cross sections:

set poi 1

% --- Parameters included in coefficient output:

set coefpara 0

B1_RABSXS            % Reduced absorption cross section
B1_NSF               % Fission neutron production cross section
B1_DIFFCOEF          % Diffusion coefficient
B1_S0                % Scattering matrix
B1_KAPPA             % Energy deposited per fission (MeV)
B1_FISS              % Fission cross section
B1_INVV              % Inverse neutron speed
B1_I135_MICRO_ABS    % Microscopic absorption cross section of I-135
B1_XE135_MICRO_ABS   % Microscopic absorption cross section of Xe-135
B1_XE135_MACRO_ABS   % Macroscopic absorption cross section of Xe-135
B1_PM149_MICRO_ABS   % Microscopic absorption cross section of Pm-149
B1_I135_YIELD        % Fission yield of I-135
B1_XE135_YIELD       % Fission yield of Xe-135
B1_PM149_YIELD       % Fission yield of Pm-149
B1_SM149_YIELD       % Fission yield of Sm-149
B1_SM149_MICRO_ABS   % Microscopic absorption cross section of Sm-149
B1_SM149_MACRO_ABS   % Macroscopic absorption cross section of Sm-149
BETA_EFF             % Effective delayed neutron fraction
LAMBDA               % Delayed neutron decay constants
DF_SURF_DF           % Surface discontinuity factors
DF_CORN_DF           % Corner discontinuity factors

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

% ----- Branches -------------------------------------------------------------

% NOTE: Variables are used to pass additional information to output:
%       LIB - Cross section library used in the calculation
%       BHI - boron concentration during burnup calculation
%       VHI - coolant temperature during burnup calculation
%       BOR - boron concentration at branch
%       TMO - coolant temperature at branch
%       TFU - fuel temperature at branch
%       CR  - insertion of control rods (0 = out, 1 = in)

% --- History variables:

branch HIS
var LIB ENDFB7
var BHI 500
var VHI 550

% --- Nominal branch: (do nothing)

branch nom  

% --- Branch to low boron concentration:

branch Blo   
repm cool cool_0B
var BOR 0  

% --- Branch to high boron concentration:

branch Bhi   
repm cool cool_1000B
var BOR 1000  

% --- Branch to low coolant temperature:

branch Clo
stp water -0.84279 500 lwtr lwe7.08t lwe7.10t
stp boron  1.00000 500
var TMO 500  

% --- Branch to high coolant temperature:

branch Chi
stp water -0.66118 600 lwtr lwe7.10t lwe7.12t
stp boron  1.00000 600
var TMO 600  

% --- Branch to low fuel temperature:

branch Flo
stp fuel -10.2685 600
stp zirc4 -6.5123 556
repu FUE FUEC
var TFU 600  

% --- Branch to high fuel temperature:

branch Fhi
stp fuel -10.1348 1200
stp zirc4 -6.5020 628
repu FUE FUEH
var TFU 1200 

% --- Insertion of control rods:

branch CR 
repu GTU ROD
var CR 1

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

% ----- Burnup calculation ---------------------------------------------------

% --- Power density:

set powdens 38.6E-3

% --- Burnup algorithm: (linear extrapolation for predictor,
%     linear interpolation for corrector)

set pcc LELI

% --- Depletion history:

dep butot

 0.1
 0.5
 1.0
 2.0
 3.0
 4.0
 5.0
 6.0
 7.0
 8.0
 9.0
10.0
11.0
12.0
13.0
14.0
15.0
17.5
20.0
22.5
25.0
27.5
30.0
32.5
35.0
37.5
40.0
45.0
50.0

% --- Coefficient matrix: (full set of variations: 3x3x3x2)

coef 11 0 0.1 1 3 5 10 15 20 30 40 50

3 nom Blo Bhi
3 nom Clo Chi
3 nom Flo Fhi
2 nom CR
1 HIS

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

% ----- Misc. options ---------------------------------------------------------

% --- Criticality cycles: (5M histories in total)

set pop 25000 200 10

% --- Libraries:

set acelib "sss_endfb7u.xsdata"
set declib "sss_endfb7.dec"
set nfylib "sss_endfb7.nfy"

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

Output

In total the group constant output (.coe file) is more than 15000 lines long. Below are the first and the last two blocks.

1 594 1 54 1
5 nom nom nom nom HIS
6 VERSION 2.1.25 DATE 16/03/04 TIME 16:09:38 LIB ENDFB7 BHI 500 VHI 550
0 1 11
0 21
B1_RABSXS 2  8.87042E-03  8.38754E-02
B1_NSF 2  6.34881E-03  1.30902E-01
B1_DIFFCOEF 2  1.35858E+00  3.74262E-01
B1_S0 4  5.36648E-01  1.86897E-02  1.39199E-03  1.37185E+00
B1_KAPPA 2  2.03484E+02  2.02023E+02
B1_FISS 2  2.48523E-03  5.37209E-02
B1_INVV 2  5.63090E-08  2.50429E-06
B1_I135_MICRO_ABS 2  8.95379E-01  4.27630E+01
B1_XE135_MICRO_ABS 2  1.25147E+02  1.61930E+06
B1_XE135_MACRO_ABS 2  0.00000E+00  0.00000E+00
B1_PM149_MICRO_ABS 2  4.54345E+01  7.49439E+02
B1_I135_YIELD 2  6.48409E-02  6.28187E-02
B1_XE135_YIELD 2  1.86087E-03  2.56635E-03
B1_PM149_YIELD 2  1.24892E-02  1.08163E-02
B1_SM149_YIELD 2  1.30223E-12  1.71016E-12
B1_SM149_MICRO_ABS 2  9.84279E+01  4.84305E+04
B1_SM149_MACRO_ABS 2  0.00000E+00  0.00000E+00
BETA_EFF 9  6.86411E-03  1.90624E-04  1.09301E-03  1.05376E-03  3.18840E-03  1.00515E-03  3.33166E-04  0.00000E+00  0.00000E+00
LAMBDA 9  7.96989E-01  1.24907E-02  3.17022E-02  1.10018E-01  3.19883E-01  1.34694E+00  8.82720E+00  0.00000E+00  0.00000E+00
DF_SURF_DF 8  9.95351E-01  1.00532E+00  9.95351E-01  1.00532E+00  9.95351E-01  1.00532E+00  9.95351E-01  1.00532E+00
DF_CORN_DF 8  9.97021E-01  9.89742E-01  9.97021E-01  9.89742E-01  9.97021E-01  9.89742E-01  9.97021E-01  9.89742E-01
2 594 1 54 1
5 nom nom nom nom HIS
6 VERSION 2.1.25 DATE 16/03/04 TIME 16:10:37 LIB ENDFB7 BHI 500 VHI 550
0.1 2 11
0 21
B1_RABSXS 2  8.90522E-03  8.64748E-02
B1_NSF 2  6.35434E-03  1.29526E-01
B1_DIFFCOEF 2  1.36230E+00  3.74963E-01
B1_S0 4  5.37555E-01  1.88519E-02  1.43431E-03  1.37101E+00
B1_KAPPA 2  2.03480E+02  2.02030E+02
B1_FISS 2  2.48741E-03  5.31443E-02
B1_INVV 2  5.67012E-08  2.49361E-06
B1_I135_MICRO_ABS 2  9.01326E-01  4.25210E+01
B1_XE135_MICRO_ABS 2  1.25990E+02  1.60119E+06
B1_XE135_MACRO_ABS 2  2.80596E-07  3.19300E-03
B1_PM149_MICRO_ABS 2  4.56829E+01  7.45208E+02
B1_I135_YIELD 2  6.48345E-02  6.28220E-02
B1_XE135_YIELD 2  1.86624E-03  2.57662E-03
B1_PM149_YIELD 2  1.24838E-02  1.08180E-02
B1_SM149_YIELD 2  1.59412E-12  2.74691E-12
B1_SM149_MICRO_ABS 2  9.95010E+01  4.78917E+04
B1_SM149_MACRO_ABS 2  1.89051E-07  8.15555E-05
BETA_EFF 9  6.90190E-03  1.90933E-04  1.11859E-03  1.05900E-03  3.14686E-03  1.03962E-03  3.46911E-04  0.00000E+00  0.00000E+00
LAMBDA 9  8.15110E-01  1.24906E-02  3.16752E-02  1.10113E-01  3.19938E-01  1.34713E+00  8.84946E+00  0.00000E+00  0.00000E+00
DF_SURF_DF 8  9.96041E-01  1.00281E+00  9.96041E-01  1.00281E+00  9.96041E-01  1.00281E+00  9.96041E-01  1.00281E+00
DF_CORN_DF 8  9.99124E-01  9.84276E-01  9.99124E-01  9.84276E-01  9.99124E-01  9.84276E-01  9.99124E-01  9.84276E-01

(...)

593 594 54 54 1
5 Bhi Chi Fhi CR HIS
10 VERSION 2.1.25 DATE 16/03/05 TIME 03:35:02 BOR 1000 TMO 600 TFU 1200 CR 1 LIB ENDFB7 BHI 500 VHI 550
40 10 11
0 21
B1_RABSXS 2  1.45840E-02  1.20816E-01
B1_NSF 2  4.01023E-03  1.04049E-01
B1_DIFFCOEF 2  1.55011E+00  4.24660E-01
B1_S0 4  4.92427E-01  1.30168E-02  2.39991E-03  1.13277E+00
B1_KAPPA 2  2.06392E+02  2.06835E+02
B1_FISS 2  1.44318E-03  3.74459E-02
B1_INVV 2  5.11004E-08  2.34745E-06
B1_I135_MICRO_ABS 2  8.04398E-01  4.09391E+01
B1_XE135_MICRO_ABS 2  9.57297E+01  1.49080E+06
B1_XE135_MACRO_ABS 2  1.57715E-07  2.33607E-03
B1_PM149_MICRO_ABS 2  4.40999E+01  7.17546E+02
B1_I135_YIELD 2  6.72309E-02  6.54456E-02
B1_XE135_YIELD 2  3.00875E-03  7.41172E-03
B1_PM149_YIELD 2  1.42081E-02  1.22700E-02
B1_SM149_YIELD 2  2.28052E-10  5.03353E-10
B1_SM149_MICRO_ABS 2  8.46711E+01  4.48777E+04
B1_SM149_MACRO_ABS 2  1.52142E-06  7.67388E-04
BETA_EFF 9  4.84684E-03  1.01061E-04  8.58352E-04  7.31839E-04  2.03392E-03  8.75822E-04  2.45841E-04  0.00000E+00  0.00000E+00
LAMBDA 9  8.09453E-01  1.26344E-02  3.04614E-02  1.12210E-01  3.26038E-01  1.25631E+00  8.25938E+00  0.00000E+00  0.00000E+00
DF_SURF_DF 8  1.04465E+00  1.34650E+00  1.04465E+00  1.34650E+00  1.04465E+00  1.34650E+00  1.04465E+00  1.34650E+00
DF_CORN_DF 8  1.08144E+00  1.51008E+00  1.08144E+00  1.51008E+00  1.08144E+00  1.51008E+00  1.08144E+00  1.51008E+00
594 594 54 54 1
5 Bhi Chi Fhi CR HIS
10 VERSION 2.1.25 DATE 16/03/05 TIME 03:36:17 BOR 1000 TMO 600 TFU 1200 CR 1 LIB ENDFB7 BHI 500 VHI 550
50 11 11
0 21
B1_RABSXS 2  1.48326E-02  1.18601E-01
B1_NSF 2  3.76891E-03  9.62319E-02
B1_DIFFCOEF 2  1.56501E+00  4.23305E-01
B1_S0 4  4.92845E-01  1.29018E-02  2.34547E-03  1.13476E+00
B1_KAPPA 2  2.06883E+02  2.07488E+02
B1_FISS 2  1.33539E-03  3.40908E-02
B1_INVV 2  5.08111E-08  2.35886E-06
B1_I135_MICRO_ABS 2  7.97360E-01  4.11893E+01
B1_XE135_MICRO_ABS 2  9.40543E+01  1.50528E+06
B1_XE135_MACRO_ABS 2  1.41984E-07  2.16652E-03
B1_PM149_MICRO_ABS 2  4.39222E+01  7.21921E+02
B1_I135_YIELD 2  6.76783E-02  6.58339E-02
B1_XE135_YIELD 2  3.03661E-03  7.85162E-03
B1_PM149_YIELD 2  1.45667E-02  1.25037E-02
B1_SM149_YIELD 2  2.88049E-10  5.53854E-10
B1_SM149_MICRO_ABS 2  8.36402E+01  4.52185E+04
B1_SM149_MACRO_ABS 2  1.42463E-06  7.34828E-04
BETA_EFF 9  4.63295E-03  1.00960E-04  8.53696E-04  7.01447E-04  1.92259E-03  8.07486E-04  2.46765E-04  0.00000E+00  0.00000E+00
LAMBDA 9  8.12082E-01  1.26269E-02  3.02904E-02  1.12647E-01  3.26856E-01  1.23635E+00  8.14923E+00  0.00000E+00  0.00000E+00
DF_SURF_DF 8  1.04451E+00  1.33799E+00  1.04451E+00  1.33799E+00  1.04451E+00  1.33799E+00  1.04451E+00  1.33799E+00
DF_CORN_DF 8  1.08628E+00  1.49629E+00  1.08628E+00  1.49629E+00  1.08628E+00  1.49629E+00  1.08628E+00  1.49629E+00