XGPT example

From Serpent Wiki
Jump to: navigation, search

Example on calculation of reaction-wise reactivity temperature coefficients using xGPT.


% --- Infinite homogeneous critical U235, U238, H1, B10 mixture

surf 1  inf

cell 1  0  fuel   -1
cell 2  0  outside 1

mat fuel         sum
 92235.03c  1.12633E-03
 92238.03c  2.18178E-02
  1001.03c  4.58880E-02
  5010.03c  8.61741e-5

% --- Cross section library

set acelib "/home/tpvtuomas/xsdata/xsss2-beta/sss_endfb71.xsdata"

% --- Neutron population

set pop 10000 100000 20 1 50 50

% --- Neutron and event buffer sizes

set nbuf 3 3

% --- The response we are interested in is k-effective

sens resp keff

% --- Use custom perturbations

%     For example:
%     CPT.txt contains relative difference of fuel material macroscopic capture XSs
%     in 300 K and 301 K (Sig_301 - Sig_300)/Sig_300
%     The value is scored when interactions happen in material "fuel"
%     and the sampled reaction is a capture reaction "capt", i.e. any reaction with
%     total neutron yield of 0.0

% --- This should yield the temperature reactivity effect of these three
%     (sum) reaction types

sens pert custom DopCPT  efunc "CPT.txt"  matlist fuel realist capt
sens pert custom DopFISS efunc "FISS.txt" matlist fuel realist fiss
sens pert custom DopELA  efunc "ELA.txt"  matlist fuel realist ela

% --- Number of latent generations to use

sens opt latgen 14

% --- Score perturbations directly to a result array (do not store events)

sens opt direct 0.20

% --- Define energy grid to be used (one that covers whole neutron energy range)

sens opt egrid myGrid

ene myGrid 1 1e-11 2e1

Each of the three files CPT.txt, FISS.txt and CPT.txt have a similar structure

E1 f1
E2 f2
E96219 f96219

where the E and f are energy-value pairs. To calculate the values for the energy dependent functions at each point, the user needs to obtain the material macroscopic cross sections at temperatures of 300 K and 301 K, based on which the relative difference between the two can be then calculated. The cross sections can be obtained, e.g. via the use of the set xsplot.