soluble absorber in serpent 2
soluble absorber in serpent 2
I'm trying to run depletion with soluble boron in the moderator using set abs , but get message saying Parameter is not supported in Serpent 2, and I couldn't find the new parameter setting for set abs in the online manual. Anyone know where to find "set abs" in Serpent 2? Thanks,
 Jaakko Leppänen
 Site Admin
 Posts: 2440
 Joined: Thu Mar 18, 2010 10:43 pm
 Security question 2: 0
 Location: Espoo, Finland
 Contact:
Re: soluble absorber in serpent 2
There is another feature that can be used for the same purpose:
http://serpent.vtt.fi/mediawiki/index.p ... t_iter_nuc
http://serpent.vtt.fi/mediawiki/index.p ... t_iter_nuc
 Jaakko
Re: soluble absorber in serpent 2
Great. Thanks.
Re: soluble absorber in serpent 2
however, when used set iter nuc, get error" Fatal error in function CheckRealistSum::
Run without set iter nuc, it works. I'm wondering if I could find any successful examples with set iter nuc in the manual or other places?
THanks,
Run without set iter nuc, it works. I'm wondering if I could find any successful examples with set iter nuc in the manual or other places?
THanks,
Re: soluble absorber in serpent 2
Tried some cases. Looks to me that if the absorber only exist in the material in the iteration, it works. If the absorber also exist in some material not in the iteration, error occurs.
E.g. set iter nuc 10 1.0 1 50100 1 water
dep daytot
10
20...
It works if B10 only exist in water. If there are absorber rods made with B4C, the burnup calculation failed.
Is my observation correct?
E.g. set iter nuc 10 1.0 1 50100 1 water
dep daytot
10
20...
It works if B10 only exist in water. If there are absorber rods made with B4C, the burnup calculation failed.
Is my observation correct?
 Jaakko Leppänen
 Site Admin
 Posts: 2440
 Joined: Thu Mar 18, 2010 10:43 pm
 Security question 2: 0
 Location: Espoo, Finland
 Contact:
Re: soluble absorber in serpent 2
Thanks, Jaakko. My work hasn't been published yet, I don't think it's a good idea for me to post the input file here yet. I'll keep working on it and let you know what's going on.

 Posts: 518
 Joined: Fri Sep 07, 2012 1:43 pm
 Security question 1: No
 Security question 2: 92
Re: soluble absorber in serpent 2
Hi,
you can specify which materials are included in the boron iteration (see the input manual). This will allow you to only iterate the boron atomic density in coolant/moderator.
A minimal working input for boron iteration in a single material case is here:
Ville
you can specify which materials are included in the boron iteration (see the input manual). This will allow you to only iterate the boron atomic density in coolant/moderator.
A minimal working input for boron iteration in a single material case is here:
Code: Select all
set title "Infinite LWRlike system iterating B10 density"
%  Geometry is just a cube
surf 2 cuboid 100.0 100.0 100.0 100.0 100.0 100.0
%  Cell definitions:
cell 3 0 fuel 2
cell 99 0 outside 2 % Outside world
%  Fuel materials:
mat fuel 5.0
8016.03c 0.316667
1001.03c 0.633333
92235.03c 0.0015
92238.03c 0.0485
5010.03c 0.0001
%  Iterate a multiplicative factor for 50100 atomic density
% to get a keff of 1.0
set iter nuc 20 1.0 1 50100 1 fuel
%  Reflective boundary condition:
set bc 2 2 2
%  Neutron population:
set pop 5000 100 10
Re: soluble absorber in serpent 2
I have material water2 and b4c, and both have B10. I tried a case with boron commented out in b4c, i.e. only have boron in water2, the case was successful. Otherwise , I got fatal error in CheckRealistSum, looks to me it's the summing over b4c instead of water2.
I'm puzzled. Please see the error message between , and the input at the end. Does anyone ever experienced similar situation? Thanks,

Sampling initial source...
CheckRealistSum: Fatal error encountered...
mat = b4c E = 2.351708E07
nuc = 6000.06c mt = 2 Emin = 1.000000E09 xs = 4.783051E+00 sum = 4.783051E+00
...
nuc = 6000.06c mt = 1 Emin = 1.000000E09 xs = 4.784168E+00 sum = 4.784168E+00
diff = 0.000000E+00
nuc = 5010.06c mt = 2 Emin = 1.000000E09 xs = 2.199080E+00 sum = 2.199080E+00
...
nuc = 5010.06c mt = 1 Emin = 1.000000E09 xs = 1.263040E+03 sum = 1.263040E+03
diff = 0.000000E+00
nuc = 5011.06c mt = 2 Emin = 1.000000E09 xs = 4.888831E+00 sum = 4.888831E+00
...
nuc = 5011.06c mt = 1 Emin = 1.000000E09 xs = 4.890638E+00 sum = 4.890638E+00
diff = 0.000000E+00
totxs = 1.154390E+01 (ures [1.000000E+37 1.000000E+37])
majorant = 7.119618E+00 (ures [0.000000E+00 0.000000E+00])
Calculating sum over nuclides and reactions:
nuc = 6000.06c adens = 5.094520E02 mt = 2 Emin = 1.000000E09 xs = 2.436735E01 sum = 2.436735E01
...
nuc = 5010.06c adens = 8.809502E03 mt = 113 Emin = 1.000000E09 xs = 2.023902E05 sum = 1.154390E+01
Optimization mode = 4
Reconstruction tolerance = 1.000000E04
totxs = 1.154390E+01 sum = 1.154390E+01 diff = 2.220446E16
majorant = 7.119618E+00 diff = 6.214217E01
***** Tue Aug 6 08:44:16 2019:
 MPI task = 8
 OpenMP thread = 9
 RNG parent seed = 1565095384
 RNG history seed = 16818534749929305560
 RNG history idx = 74974
... ...
Simulation aborted.
 RNG parent seed = 1565095384
 RNG history seed = 17757758591319832024
 RNG history idx = 99952
Emin = 1.000000E09 xs = 5.945720E04 sum = 4.752207E+00
nuc = 6000.06c mt = 107 Fatal error in function CheckRealistSum:
Error in total or majorant, or reaction sampling failedEmin = 6.181241E+00 xs = 0.000000E+00 sum = 4.752207E+00
nuc = 6000.06c mt = 1
Simulation aborted.

mat water2 0.7207 moder lwtr 1001
1001.06c 1.117159E01
8016.06c 8.859261E01
8017.06c 3.580223E04
5010.06c 3.686189E04
5011.06c 1.631381E03
therm lwtr lwj3.13t
%
mat b4c 1.81
6000.06c 1.0 %2.172177E01
5010.06c 1.442742E01
5011.06c 6.385081E01
% options
set pop 100000 1000 50
% burn
set iter nuc 10 1.0 1 50100 1 water2
dep daytot
10
20
...
I'm puzzled. Please see the error message between , and the input at the end. Does anyone ever experienced similar situation? Thanks,

Sampling initial source...
CheckRealistSum: Fatal error encountered...
mat = b4c E = 2.351708E07
nuc = 6000.06c mt = 2 Emin = 1.000000E09 xs = 4.783051E+00 sum = 4.783051E+00
...
nuc = 6000.06c mt = 1 Emin = 1.000000E09 xs = 4.784168E+00 sum = 4.784168E+00
diff = 0.000000E+00
nuc = 5010.06c mt = 2 Emin = 1.000000E09 xs = 2.199080E+00 sum = 2.199080E+00
...
nuc = 5010.06c mt = 1 Emin = 1.000000E09 xs = 1.263040E+03 sum = 1.263040E+03
diff = 0.000000E+00
nuc = 5011.06c mt = 2 Emin = 1.000000E09 xs = 4.888831E+00 sum = 4.888831E+00
...
nuc = 5011.06c mt = 1 Emin = 1.000000E09 xs = 4.890638E+00 sum = 4.890638E+00
diff = 0.000000E+00
totxs = 1.154390E+01 (ures [1.000000E+37 1.000000E+37])
majorant = 7.119618E+00 (ures [0.000000E+00 0.000000E+00])
Calculating sum over nuclides and reactions:
nuc = 6000.06c adens = 5.094520E02 mt = 2 Emin = 1.000000E09 xs = 2.436735E01 sum = 2.436735E01
...
nuc = 5010.06c adens = 8.809502E03 mt = 113 Emin = 1.000000E09 xs = 2.023902E05 sum = 1.154390E+01
Optimization mode = 4
Reconstruction tolerance = 1.000000E04
totxs = 1.154390E+01 sum = 1.154390E+01 diff = 2.220446E16
majorant = 7.119618E+00 diff = 6.214217E01
***** Tue Aug 6 08:44:16 2019:
 MPI task = 8
 OpenMP thread = 9
 RNG parent seed = 1565095384
 RNG history seed = 16818534749929305560
 RNG history idx = 74974
... ...
Simulation aborted.
 RNG parent seed = 1565095384
 RNG history seed = 17757758591319832024
 RNG history idx = 99952
Emin = 1.000000E09 xs = 5.945720E04 sum = 4.752207E+00
nuc = 6000.06c mt = 107 Fatal error in function CheckRealistSum:
Error in total or majorant, or reaction sampling failedEmin = 6.181241E+00 xs = 0.000000E+00 sum = 4.752207E+00
nuc = 6000.06c mt = 1
Simulation aborted.

mat water2 0.7207 moder lwtr 1001
1001.06c 1.117159E01
8016.06c 8.859261E01
8017.06c 3.580223E04
5010.06c 3.686189E04
5011.06c 1.631381E03
therm lwtr lwj3.13t
%
mat b4c 1.81
6000.06c 1.0 %2.172177E01
5010.06c 1.442742E01
5011.06c 6.385081E01
% options
set pop 100000 1000 50
% burn
set iter nuc 10 1.0 1 50100 1 water2
dep daytot
10
20
...

 Posts: 72
 Joined: Tue Feb 27, 2018 10:52 am
 Security question 1: No
 Security question 2: 93
Re: soluble absorber in serpent 2
I tried to reproduce the error with a very simple input but I did not succeed. Can you provide us some input that produces the error? Otherwise it is a bit difficult to help you. Also, what Serpent version are you using?