soluble absorber in serpent 2

Questions, comments and complaints about the User's Manual
wangxi
Posts: 8
Joined: Mon Jul 29, 2019 10:39 pm
Security question 1: No
Security question 2: 19

soluble absorber in serpent 2

Post by wangxi » Mon Jul 29, 2019 10:47 pm

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 on-line manual. Anyone know where to find "set abs" in Serpent 2? Thanks,

User avatar
Jaakko Leppänen
Site Admin
Posts: 2434
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: soluble absorber in serpent 2

Post by Jaakko Leppänen » Tue Jul 30, 2019 8:49 am

There is another feature that can be used for the same purpose:

http://serpent.vtt.fi/mediawiki/index.p ... t_iter_nuc
- Jaakko

wangxi
Posts: 8
Joined: Mon Jul 29, 2019 10:39 pm
Security question 1: No
Security question 2: 19

Re: soluble absorber in serpent 2

Post by wangxi » Tue Jul 30, 2019 3:07 pm

Great. Thanks.

wangxi
Posts: 8
Joined: Mon Jul 29, 2019 10:39 pm
Security question 1: No
Security question 2: 19

Re: soluble absorber in serpent 2

Post by wangxi » Tue Jul 30, 2019 6:39 pm

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,

wangxi
Posts: 8
Joined: Mon Jul 29, 2019 10:39 pm
Security question 1: No
Security question 2: 19

Re: soluble absorber in serpent 2

Post by wangxi » Tue Jul 30, 2019 9:39 pm

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?

User avatar
Jaakko Leppänen
Site Admin
Posts: 2434
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: soluble absorber in serpent 2

Post by Jaakko Leppänen » Wed Jul 31, 2019 8:49 am

Can you post your full input here?
- Jaakko

wangxi
Posts: 8
Joined: Mon Jul 29, 2019 10:39 pm
Security question 1: No
Security question 2: 19

Re: soluble absorber in serpent 2

Post by wangxi » Wed Jul 31, 2019 3:30 pm

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.

Ville Valtavirta
Posts: 514
Joined: Fri Sep 07, 2012 1:43 pm
Security question 1: No
Security question 2: 92

Re: soluble absorber in serpent 2

Post by Ville Valtavirta » Wed Jul 31, 2019 5:25 pm

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:

Code: Select all

set title "Infinite LWR-like system iterating B-10 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
-Ville

wangxi
Posts: 8
Joined: Mon Jul 29, 2019 10:39 pm
Security question 1: No
Security question 2: 19

Re: soluble absorber in serpent 2

Post by wangxi » Tue Aug 06, 2019 4:07 pm

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.351708E-07
nuc = 6000.06c mt = 2 Emin = 1.000000E-09 xs = 4.783051E+00 sum = 4.783051E+00
...
nuc = 6000.06c mt = 1 Emin = 1.000000E-09 xs = 4.784168E+00 sum = 4.784168E+00
diff = -0.000000E+00

nuc = 5010.06c mt = 2 Emin = 1.000000E-09 xs = 2.199080E+00 sum = 2.199080E+00
...
nuc = 5010.06c mt = 1 Emin = 1.000000E-09 xs = 1.263040E+03 sum = 1.263040E+03
diff = -0.000000E+00

nuc = 5011.06c mt = 2 Emin = 1.000000E-09 xs = 4.888831E+00 sum = 4.888831E+00
...
nuc = 5011.06c mt = 1 Emin = 1.000000E-09 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.094520E-02 mt = 2 Emin = 1.000000E-09 xs = 2.436735E-01 sum = 2.436735E-01
...
nuc = 5010.06c adens = 8.809502E-03 mt = 113 Emin = 1.000000E-09 xs = 2.023902E-05 sum = 1.154390E+01

Optimization mode = 4
Reconstruction tolerance = 1.000000E-04
totxs = 1.154390E+01 sum = 1.154390E+01 diff = 2.220446E-16
majorant = 7.119618E+00 diff = 6.214217E-01

***** 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.000000E-09 xs = 5.945720E-04 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.117159E-01
8016.06c -8.859261E-01
8017.06c -3.580223E-04
5010.06c -3.686189E-04
5011.06c -1.631381E-03
therm lwtr lwj3.13t
%
mat b4c -1.81
6000.06c -1.0 %-2.172177E-01
5010.06c -1.442742E-01
5011.06c -6.385081E-01
% options
set pop 100000 1000 50

% burn
set iter nuc 10 1.0 1 50100 1 water2
dep daytot
10
20
...

Riku Tuominen
Posts: 68
Joined: Tue Feb 27, 2018 10:52 am
Security question 1: No
Security question 2: 93

Re: soluble absorber in serpent 2

Post by Riku Tuominen » Fri Aug 09, 2019 1:11 pm

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?

Post Reply