Automated burnup sequence example 2
Revision as of 08:13, 1 June 2022 by Antti Rintala (talk | contribs)
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)
% - Historical variations: - low, nominal and high fuel temperature
% - 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
% - Control rod insertion using the repu variation
% - All branch combinations included (3x3x3x2 = 54 combinations in matrix)
% NOTE: The values of temperature-adjusted densities 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
% --- 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 (volume (10.75182*2)^2):
set poi 1 4.624065332496000e+02
% --- 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
% TFHI - fueltemperature 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
% --- History with low fuel temperature:
branch TFHIlo
stp fuel -10.2685 600
stp zirc4 -6.5123 556
var TFUHI 600
% --- History with nominal fuel temperature:
branch TFHInom
stp fuel -10.2013 900
stp zirc4 -6.5071 600
var TFUHI 900
% --- History with high fuel temperature:
branch TFHIhi
stp fuel -10.1348 1200
stp zirc4 -6.5020 628
var TFUHI 1200
% --- 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
var TFU 600
% --- Branch to nominal fuel temperature:
branch Fnom
stp fuel -10.2013 900
stp zirc4 -6.5071 600
var TFU 900
% --- Branch to high fuel temperature:
branch Fhi
stp fuel -10.1348 1200
stp zirc4 -6.5020 628
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
% --- Casematrix: (full set of histories: 3, full set of variations: 3x3x3x2)
casematrix cm_example
3 TFHIlo TFHInom TFHIhi
11 0 0.1 1 3 5 10 15 20 30 40 50
3 nom Blo Bhi
3 nom Clo Chi
3 Flo Fnom 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"
% -----------------------------------------------------------------------------
Run commands
Example 1
sss2 -casematrix cm_example 1 0 <input> sss2 -casematrix cm_example 2 0 <input> sss2 -casematrix cm_example 3 0 <input>
Example 2
First
sss2 -casematrix cm_example 1 -1 <input> sss2 -casematrix cm_example 2 -1 <input> sss2 -casematrix cm_example 3 -1 <input>
then
sss2 -casematrix cm_example 1 0 -coe <input> sss2 -casematrix cm_example 2 0 -coe <input> sss2 -casematrix cm_example 3 0 -coe <input>
Example 3
First
sss2 -casematrix cm_example 1 -1 <input> sss2 -casematrix cm_example 2 -1 <input> sss2 -casematrix cm_example 3 -1 <input>
then
for i in {1..54}
do
sss2 -casematrix cm_example 1 i <input>
done
for i in {1..54}
do
sss2 -casematrix cm_example 2 i <input>
done
for i in {1..54}
do
sss2 -casematrix cm_example 3 i <input>
done
Output
In total the group constant output (<input>_h<j>_r<i>.coe files, j=1, 2, 3; i=1, 2, ..., 54) is 46332 lines long. Below are the first and last two blocks.
1 594 1 54 1 5 nom nom Flo nom HIS 8 VERSION 2.2.1 DATE 22/05/13 TIME 15:35:29 TFU 600 LIB ENDFB7 BHI 500 VHI 550 TFUHI 600 0 1 11 0 21 B1_RABSXS 2 8.79100E-03 8.49362E-02 B1_NSF 2 6.43147E-03 1.32840E-01 B1_DIFFCOEF 2 1.35056E+00 3.73482E-01 B1_S0 4 5.37977E-01 1.87011E-02 1.27267E-03 1.37374E+00 B1_KAPPA 2 2.03732E+02 2.02270E+02 B1_FISS 2 2.51782E-03 5.45162E-02 B1_INVV 2 5.63800E-08 2.50615E-06 B1_I135_MICRO_ABS 2 9.06767E-01 3.88556E+01 B1_XE135_MICRO_ABS 2 1.26146E+02 1.47436E+06 B1_XE135_MACRO_ABS 2 0.00000E+00 0.00000E+00 B1_PM149_MICRO_ABS 2 4.60597E+01 6.80956E+02 B1_I135_YIELD 2 6.48392E-02 6.28187E-02 B1_XE135_YIELD 2 1.86145E-03 2.56635E-03 B1_PM149_YIELD 2 1.24878E-02 1.08163E-02 B1_SM149_YIELD 2 1.30260E-12 1.71016E-12 B1_SM149_MICRO_ABS 2 9.94630E+01 4.40187E+04 B1_SM149_MACRO_ABS 2 0.00000E+00 0.00000E+00 BETA_EFF 9 6.93533E-03 1.89233E-04 1.10077E-03 1.10507E-03 3.20762E-03 9.88703E-04 3.43933E-04 0.00000E+00 0.00000E+00 LAMBDA 9 8.00386E-01 1.24907E-02 3.16805E-02 1.10020E-01 3.20110E-01 1.34689E+00 8.84551E+00 0.00000E+00 0.00000E+00 DF_SURF_DF 8 9.95261E-01 1.00248E+00 9.95261E-01 1.00248E+00 9.95261E-01 1.00248E+00 9.95261E-01 1.00248E+00 DF_CORN_DF 8 9.97033E-01 9.90216E-01 9.97033E-01 9.90216E-01 9.97033E-01 9.90216E-01 9.97033E-01 9.90216E-01 2 594 1 54 1 5 nom nom Flo nom HIS 8 VERSION 2.2.1 DATE 22/05/13 TIME 15:36:41 TFU 600 LIB ENDFB7 BHI 500 VHI 550 TFUHI 600 0.1 2 11 0 21 B1_RABSXS 2 8.82667E-03 8.75524E-02 B1_NSF 2 6.42774E-03 1.31411E-01 B1_DIFFCOEF 2 1.35329E+00 3.74285E-01 B1_S0 4 5.39063E-01 1.88696E-02 1.30821E-03 1.37261E+00 B1_KAPPA 2 2.03725E+02 2.02277E+02 B1_FISS 2 2.51666E-03 5.39182E-02 B1_INVV 2 5.68287E-08 2.49521E-06 B1_I135_MICRO_ABS 2 9.14575E-01 3.85976E+01 B1_XE135_MICRO_ABS 2 1.27152E+02 1.45581E+06 B1_XE135_MACRO_ABS 2 2.83169E-07 3.24211E-03 B1_PM149_MICRO_ABS 2 4.64459E+01 6.76423E+02 B1_I135_YIELD 2 6.48286E-02 6.28219E-02 B1_XE135_YIELD 2 1.86825E-03 2.57627E-03 B1_PM149_YIELD 2 1.24789E-02 1.08180E-02 B1_SM149_YIELD 2 1.58945E-12 2.71187E-12 B1_SM149_MICRO_ABS 2 1.00392E+02 4.35421E+04 B1_SM149_MACRO_ABS 2 1.90871E-07 8.27847E-05 BETA_EFF 9 6.75361E-03 1.93081E-04 1.10606E-03 1.04009E-03 3.11166E-03 9.75198E-04 3.27524E-04 0.00000E+00 0.00000E+00 LAMBDA 9 7.95254E-01 1.24908E-02 3.16819E-02 1.10044E-01 3.20163E-01 1.34677E+00 8.86406E+00 0.00000E+00 0.00000E+00 DF_SURF_DF 8 9.94028E-01 1.00462E+00 9.94028E-01 1.00462E+00 9.94028E-01 1.00462E+00 9.94028E-01 1.00462E+00 DF_CORN_DF 8 9.91825E-01 9.94763E-01 9.91825E-01 9.94763E-01 9.91825E-01 9.94763E-01 9.91825E-01 9.94763E-01 (...) 593 594 54 54 1 5 Bhi Chi Fhi CR HIS 11 VERSION 2.2.1 DATE 22/05/14 TIME 05:19:31 BOR 1000 TMO 600 TFU 1200 CR 1 LIB ENDFB7 BHI 500 VHI 550 TFUHI 1200 40 10 11 0 21 B1_RABSXS 2 1.44616E-02 1.21087E-01 B1_NSF 2 4.00648E-03 1.05102E-01 B1_DIFFCOEF 2 1.56234E+00 4.25928E-01 B1_S0 4 4.89801E-01 1.30273E-02 2.39490E-03 1.13089E+00 B1_KAPPA 2 2.06637E+02 2.07081E+02 B1_FISS 2 1.44096E-03 3.78320E-02 B1_INVV 2 5.10865E-08 2.34671E-06 B1_I135_MICRO_ABS 2 8.09791E-01 3.90191E+01 B1_XE135_MICRO_ABS 2 9.68998E+01 1.43152E+06 B1_XE135_MACRO_ABS 2 1.60334E-07 2.36865E-03 B1_PM149_MICRO_ABS 2 4.43600E+01 6.82905E+02 B1_I135_YIELD 2 6.72070E-02 6.54430E-02 B1_XE135_YIELD 2 3.03223E-03 7.40090E-03 B1_PM149_YIELD 2 1.41847E-02 1.22685E-02 B1_SM149_YIELD 2 2.29430E-10 5.02156E-10 B1_SM149_MICRO_ABS 2 8.52364E+01 4.30738E+04 B1_SM149_MACRO_ABS 2 1.53671E-06 7.76567E-04 BETA_EFF 9 4.84854E-03 1.06453E-04 8.29470E-04 7.43615E-04 2.04635E-03 8.75171E-04 2.47486E-04 0.00000E+00 0.00000E+00 LAMBDA 9 8.23979E-01 1.26144E-02 3.04573E-02 1.12166E-01 3.25791E-01 1.25979E+00 8.55371E+00 0.00000E+00 0.00000E+00 DF_SURF_DF 8 1.04513E+00 1.33739E+00 1.04513E+00 1.33739E+00 1.04513E+00 1.33739E+00 1.04513E+00 1.33739E+00 DF_CORN_DF 8 1.08696E+00 1.48737E+00 1.08696E+00 1.48737E+00 1.08696E+00 1.48737E+00 1.08696E+00 1.48737E+00 594 594 54 54 1 5 Bhi Chi Fhi CR HIS 11 VERSION 2.2.1 DATE 22/05/14 TIME 05:21:17 BOR 1000 TMO 600 TFU 1200 CR 1 LIB ENDFB7 BHI 500 VHI 550 TFUHI 1200 50 11 11 0 21 B1_RABSXS 2 1.46975E-02 1.18910E-01 B1_NSF 2 3.77185E-03 9.74587E-02 B1_DIFFCOEF 2 1.57941E+00 4.24562E-01 B1_S0 4 4.90146E-01 1.29054E-02 2.33506E-03 1.13310E+00 B1_KAPPA 2 2.07123E+02 2.07720E+02 B1_FISS 2 1.33591E-03 3.45451E-02 B1_INVV 2 5.07503E-08 2.35873E-06 B1_I135_MICRO_ABS 2 8.01019E-01 3.93402E+01 B1_XE135_MICRO_ABS 2 9.49332E+01 1.44921E+06 B1_XE135_MACRO_ABS 2 1.44557E-07 2.20675E-03 B1_PM149_MICRO_ABS 2 4.40908E+01 6.88660E+02 B1_I135_YIELD 2 6.76457E-02 6.58237E-02 B1_XE135_YIELD 2 3.06545E-03 7.82440E-03 B1_PM149_YIELD 2 1.45349E-02 1.24977E-02 B1_SM149_YIELD 2 2.89658E-10 5.50926E-10 B1_SM149_MICRO_ABS 2 8.43392E+01 4.35188E+04 B1_SM149_MACRO_ABS 2 1.44546E-06 7.45853E-04 BETA_EFF 9 4.59341E-03 9.79931E-05 8.40226E-04 6.54586E-04 1.92507E-03 8.37270E-04 2.38265E-04 0.00000E+00 0.00000E+00 LAMBDA 9 8.10056E-01 1.26920E-02 3.02907E-02 1.12531E-01 3.27376E-01 1.22249E+00 8.16703E+00 0.00000E+00 0.00000E+00 DF_SURF_DF 8 1.04486E+00 1.34464E+00 1.04486E+00 1.34464E+00 1.04486E+00 1.34464E+00 1.04486E+00 1.34464E+00 DF_CORN_DF 8 1.08300E+00 1.50033E+00 1.08300E+00 1.50033E+00 1.08300E+00 1.50033E+00 1.08300E+00 1.50033E+00
