Detector time binning

This category replaces the missing input manual
Post Reply
User avatar
Jaakko Leppänen
Site Admin
Posts: 2354
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Detector time binning

Post by Jaakko Leppänen » Mon Oct 22, 2012 4:04 pm

Time dependence in external source calculation mode was revised in update 2.1.9 and detector time binning was added. Time bins are added in detectors using the "di" entry:

Code: Select all

det <name> di <tbin>
where <tbin> is defined by a new tme card:

Code: Select all

tme <tbin> 1 <t1> <t2> <t3> ... <tn>
tme <tbin> <type> <N> <tmin> <tmax>
The second paramer defines the binning type. For type 1, an arbitrary binning is used, and the type number is followed by the interval boundaries in ascending order. Linear uniform binning is used for type 2 and logarithmic uniform binning for type 3. For these two types the type number is followed by the number of bins and the minimum and maximum time.

When time binning is used, an additional column is added in detector outputs (the results are found in column 12 instead of 11). All time values are in seconds.
- Jaakko

EdsonHenrice
Posts: 10
Joined: Sun Mar 09, 2014 5:18 pm
Security question 1: No
Security question 2: 92

Re: Detector time binning

Post by EdsonHenrice » Wed May 07, 2014 9:50 pm

I' m using external neutron source. As cited in other post, when I use "time binning" the simulation is limited only millisecond. There are other way to obtain this in other time scale For example 50 seconds. If not, still under development it?

Follows part off my code.

set nps 20000 100 intervalos
set nbuf 5

det contador dm fuelin di intervalos
tme intervalos 2 500 0 0.1

When I use it, only first row show-me results.

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

Re: Detector time binning

Post by Jaakko Leppänen » Wed May 07, 2014 11:51 pm

I assume you are simulating a sub-critical system?

Unless you define time dependence for the source, all your source neutrons are born at t=0s. The result is that the neutron histories die out very fast, in the order of milliseconds in thermal systems and even faster in fast systems. So there are simply no neutrons left to make contributions in the time bins later on.

If you are looking for the contribution of delayed neutrons, you should note that delayed neutron emission is switched off by default in external source simulations, and can be switched on with:

Code: Select all

set delnu 1
Also note that delayed neutron emission is a relatively rare event, so you probably need to run a lot of histories to get good statistics. I've never tried this type of simulation myself. Let me know if it works.
- Jaakko

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

Re: Detector time binning

Post by Jaakko Leppänen » Thu May 08, 2014 10:03 pm

Actually, there is another problem... On line 409 of score.c, the neutron time is set to zero before the detectors are scored:

Code: Select all

     t = 0.0;
If you remove this line, time is correctly binned. I'll fix the problem for the next update.
- Jaakko

jesse.johns
Posts: 114
Joined: Wed Apr 20, 2011 5:19 am

Re: Detector time binning

Post by jesse.johns » Sun Aug 02, 2015 7:48 pm

Hello,

I am wondering if there's a simple was to sample the entire production and loss rate, or if I need to set up a detector for each reaction and a surface for leakage? I'm currently doing this, but it would be splendid if the _res.m indexed for time step for things like TOT_FISSRATE.

Thanks,

Jesse

subterfug3x
Posts: 32
Joined: Thu Oct 04, 2018 7:53 am
Security question 1: No
Security question 2: 96
Location: Tomsk, Russia

Re: Detector time binning

Post by subterfug3x » Sat Nov 23, 2019 6:15 pm

Jaakko Leppänen wrote:
Wed May 07, 2014 11:51 pm
I assume you are simulating a sub-critical system?

Unless you define time dependence for the source, all your source neutrons are born at t=0s. The result is that the neutron histories die out very fast, in the order of milliseconds in thermal systems and even faster in fast systems. So there are simply no neutrons left to make contributions in the time bins later on.

If you are looking for the contribution of delayed neutrons, you should note that delayed neutron emission is switched off by default in external source simulations, and can be switched on with:
The use of st card with src leads to error:
Input error:

Source sampling failed because of low efficiency

Code: Select all

set nps 1000000000 200 promt
src DT n sc 4001 se 14.1 sz 0 300 srad 0 17 st 0.0 1.16E-02 % Source in Plasma ZONE CELL 4001 14.1 MeV
set srcrate 4.5e17
set bc 1
set gcu -1

set tcut 3600

tme global 2 2000 0.0 2000 
tme promt  2 4000 0.0 1.16E-02
tme delay  2 500 1.16E-02 1000

set delnu 1

set bala 1
set nbuf 100000
set impl 0 0 0

set mcvol 100000000
System is subcritical. Calculation works normal without st card.
- Igor

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

Re: Detector time binning

Post by Jaakko Leppänen » Mon Nov 25, 2019 12:51 pm

Could you post the full input, please.
- Jaakko

subterfug3x
Posts: 32
Joined: Thu Oct 04, 2018 7:53 am
Security question 1: No
Security question 2: 96
Location: Tomsk, Russia

Re: Detector time binning

Post by subterfug3x » Mon Nov 25, 2019 2:52 pm

Jaakko Leppänen wrote:
Mon Nov 25, 2019 12:51 pm
Could you post the full input, please.
Full Input:

Code: Select all

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

set title "HGTRU Reactor Core with TNS"  

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

% Geometry:

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

%####################################
%##### Particles ####################
%####################################

% --- Coated particle:

particle 1

fuel    0.0175                   % --- Fuel Kernel
PyC     0.0265                   % --- PyC Buffer
Coating 0.0300                   % --- Ti3SiC2 Coating
matrix                           % --- Graphite graphite

% --- Microfuel Cell

cell 101  11  matrix  -100       % ---  graphite cell
pbed 101  11 "/home/iol4/calculations/particles/17pf-80cm-hgtru_core_part.inp"     % ---  Fill graphite with particles

%##### Particle 2 ####################

particle 2

fuel     0.0175                  % --- Fuel Kernel
PyC      0.0265                  % --- PyC Buffer
Coating  0.0300                  % --- Ti3SiC2 Coating
matrix                           % --- Graphite graphite

cell 102  12  matrix  -100        % ---  graphite cell
pbed 102  12 "/home/iol4/calculations/particles/13pf-80cm-hgtru_core_profile.inp"     % ---  Fill graphite with particles

%##### Particle 3 ####################

particle 3

fuel     0.0175                  % --- Fuel Kernel
PyC      0.0265                  % --- PyC Buffer
Coating  0.0300                  % --- Ti3SiC2 Coating
matrix                           % --- Graphite graphite

cell 103  13  matrix  -100        % ---  graphite cell
pbed 103  13 "/home/iol4/calculations/particles/19pf-80cm-hgtru_core_profile.inp"      % ---  Fill graphite with particles

%##### Particle 4 ####################

particle 4

fuel     0.0175                  % --- Fuel Kernel
PyC      0.0265                  % --- PyC Buffer
Coating  0.0300                  % --- Ti3SiC2 Coating
matrix                           % --- Graphite graphite

cell 104  14  matrix  -100        % ---  graphite cell
pbed 104  14 "/home/iol4/calculations/particles/21pf-80cm-hgtru_core_profile.inp"      % ---  Fill graphite with particles

%##### Particle 5 ####################

particle 5

fuel     0.0175                  % --- Fuel Kernel
PyC      0.0265                  % --- PyC Buffer
Coating  0.0300                  % --- Ti3SiC2 Coating
matrix                           % --- Graphite graphite

cell 105  15  matrix  -100        % ---  graphite cell
pbed 105  15 "/home/iol4/calculations/particles/23pf-80cm-hgtru_core_part.inp"     % ---  Fill graphite with particles

%########## 17%pf ###########

cell 211 31 fill 101  -201         % --- Fuel Compact 17%
cell 212 31 SiC       201 -202    % --- Compact Coating
cell 213 31 graphite        202    % --- Graphite outside compact

cell 214 40 graphite  -100
cell 215 40 graphite        100    % --- Empty lattice position:

%########## 15%pf ###########

cell 221 32 fill 102  -201         % --- Fuel Compact 17%
cell 222 32 SiC       201 -202    % --- Compact Coating
cell 223 32 graphite        202    % --- Graphite outside compact

cell 224 41 ingraphite      -100   % Not used now
cell 225 41 ingraphite       100   % --- Empty lattice position:

%########## 19%pf ###########

cell 231 33 fill 103  -201         % --- Fuel Compact 19%
cell 232 33 SiC       201 -202    % --- Compact Coating
cell 233 33 graphite        202    % --- Graphite outside compact

%########## 21%pf ###########

cell 241 34 fill 104  -201         % --- Fuel Compact 21%
cell 242 34 SiC       201 -202    % --- Compact Coating
cell 243 34 graphite        202    % --- Graphite outside compact

%########## 23%pf ###########

cell 251 35 fill 105  -201         % --- Fuel Compact 23%
cell 252 35 SiC       201 -202    % --- Compact Coating
cell 253 35 graphite        202    % --- Graphite outside compact

%########## Heterogeneous burnable poison ###########

cell 261 37 BPGd      -201         % --- Fuel Compact
cell 262 37 SiC        201 -202    % --- Compact Coating
cell 263 37 graphite        202    % --- Graphite outside compact

%####################################
%##### Fuel block and Reflector ##### LOW ENRICHMENT
%####################################

% --- Fuel Block 

cell 311 51 fill FB17  -301 311 312 313 314 315 316 317
cell 312 51 graphite  301

% --- Fuel Block Coolant Holes

cell 313 51 helium   -311
cell 314 51 helium   -312
cell 315 51 helium   -313
cell 316 51 helium   -314
cell 317 51 helium   -315
cell 318 51 helium   -316
cell 319 51 helium   -317

%####################################
%##### Inner Fuel block ############# HIGH ENRICHMENT
%####################################

cell 321 52 fill FB15  -301 311 312 313 314 315 316 317
cell 322 52 graphite 301

% --- Fuel Block Coolant Holes

cell 323 52 helium   -311
cell 324 52 helium   -312
cell 325 52 helium   -313
cell 326 52 helium   -314
cell 327 52 helium   -315
cell 328 52 helium   -316
cell 329 52 helium   -317

%####################################
%##### Inner Fuel block #############
%####################################

cell 331 53 fill FB19  -301 311 312 313 314 315 316 317
cell 332 53 ingraphite 301

% --- Fuel Block Coolant Holes

cell 333 53 helium   -311
cell 334 53 helium   -312
cell 335 53 helium   -313
cell 336 53 helium   -314
cell 337 53 helium   -315
cell 338 53 helium   -316
cell 339 53 helium   -317

%####################################
%##### Inner Fuel block #############
%####################################

cell 341 54 fill FB21  -301 311 312 313 314 315 316 317
cell 342 54 ingraphite 301

% --- Fuel Block Coolant Holes

cell 343 54 helium   -311
cell 344 54 helium   -312
cell 345 54 helium   -313
cell 346 54 helium   -314
cell 347 54 helium   -315
cell 348 54 helium   -316
cell 349 54 helium   -317

%####################################
%##### Inner Fuel block #############
%####################################

cell 351 55 fill FB23  -301 311 312 313 314 315 316 317
cell 352 55 ingraphite 301

% --- Fuel Block Coolant Holes

cell 353 55 helium   -311
cell 354 55 helium   -312
cell 355 55 helium   -313
cell 356 55 helium   -314
cell 357 55 helium   -315
cell 358 55 helium   -316
cell 359 55 helium   -317

%####################################
%##### Control Rod Block ############
%####################################

cell 3101 90 absorber_graphite  -301 311 312 313 314 315 316 317
cell 3102 90 absorber_graphite   301

% --- Coolant Holes for Reflector

cell 3103 90 helium             -311
cell 3104 90 absorber           -312
cell 3105 90 absorber           -313
cell 3106 90 absorber           -314
cell 3107 90 absorber           -315
cell 3108 90 absorber           -316
cell 3109 90 absorber           -317

%####################################
%##### Blocks in Core with 7 holes #####
%####################################

cell 3211 57 graphite  -301 311 312 313 314 315 316 317
cell 3212 57 graphite   301

% --- Coolant Holes for Reflector

cell 3213 57 helium    -311
cell 3214 57 helium    -312
cell 3215 57 helium    -313
cell 3216 57 helium    -314
cell 3217 57 helium    -315
cell 3218 57 helium    -316
cell 3219 57 helium    -317

%####################################
%##### Reflector with 7 holes UP AND BOT########
%####################################

cell 3311 58 reflector  -302 311 312 313 314 315 316 317
cell 3312 58 reflector   302

% --- Coolant Holes for Reflector

cell 3313 58 reflector_coolant    -311
cell 3314 58 reflector_coolant    -312
cell 3315 58 reflector_coolant    -313
cell 3316 58 reflector_coolant    -314
cell 3317 58 reflector_coolant    -315
cell 3318 58 reflector_coolant    -316
cell 3319 58 reflector_coolant    -317

%####################################
%####### Side Reflector #############
%####################################

% --- Reflector with 1 coolant hole CORE LAYER

cell 3411 60 reflector_coolant  -311
cell 3412 60 reflector          -301 311
cell 3413 60 reflector           301 

% --- Reflector with 1 coolant hole CORE LAYER ABSOBER

cell 3421 80 absorber           -320   % -311 Absorber
cell 3422 80 reflector          -301 320
cell 3423 80 reflector           301

% --- Upper and Bottom LAYER

cell 3514 61 reflector_coolant  -311
cell 3515 61 reflector          -302 311
cell 3516 61 reflector           302 

% --- Upper and Bottom LAYER ABSOBER

cell 3524 81 absorber           -320  % -311 Absorber
cell 3525 81 reflector          -302 320
cell 3526 81 reflector           302 

%####################################
%####### Outside World ##############
%####################################

% --- Helium outside zones in core layers

cell 3511 70 reflector_coolant -301
cell 3512 70 reflector_coolant  301

% --- Helium outside zones in reflector layers

cell 3521 71 reflector_coolant -302
cell 3522 71 reflector_coolant  302

%###################################
%############ Core Design ##########
%###################################
 
% --- Core Fill Cell

cell 401 20 fill CoreLayer    -401
cell 402 20 graphite           401

% --- Reflector Fill Cell

cell 411 21 fill TopRef       -402
cell 412 21 reflector          402

% --- Core Axial Lattice

lat Core 9 0 0 5 
0   21   % --- Bot Refl Level 0
30  20   % --- Core Level 1
110 20   % --- Core Level 2
190 20   % --- Core Level 3
270 21   % --- Upp Refl Level 5

%###################################
%########### CORE + TNS ############
%###################################

% --- TNS Cell


cell 4001 0 plasma -4000 4008 -4007
cell 4002 0 void    4000 -4001 4008 -4007
cell 4003 0 steel   4001 -4002 4008 -4007
cell 4004 0 air     4002 -4003 4008 -4007
cell 4005 0 CuH2O  -4004 4003 4010 -4009
cell 4006 0 void   -4004 4003 4012 -4013

% ---- Solenoid coils + Solenoid Coolant

cell	4007	0	CuH2O	-4004	4003	4015	-4014 
cell	4008	0	CuH2O	-4004	4003	4017	-4016
cell	4009	0	CuH2O	-4004	4003	4019	-4018
cell	4010	0	CuH2O	-4004	4003	4021	-4020
cell	4011	0	CuH2O	-4004	4003	4023	-4022
cell	4012	0	CuH2O	-4004	4003	4025	-4024
cell	4013	0	CuH2O	-4004	4003	4027	-4026
cell	4014	0	CuH2O	-4004	4003	4029	-4028
cell	4015	0	CuH2O	-4004	4003	4031	-4030
cell	4016	0	CuH2O	-4004	4003	4033	-4032
cell	4017	0	CuH2O	-4004	4003	4035	-4034
cell	4018	0	CuH2O	-4004	4003	4037	-4036
cell	4019	0	CuH2O	-4004	4003	4039	-4038
cell	4020	0	CuH2O	-4004	4003	4041	-4040
cell	4021	0	CuH2O	-4004	4003	4043	-4042
cell	4022	0	CuH2O	-4004	4003	4045	-4044
cell	4023	0	CuH2O	-4004	4003	4047	-4046
cell	4024	0	CuH2O	-4004	4003	4049	-4048
cell	4025	0	CuH2O	-4004	4003	4051	-4050
cell	4026	0	CuH2O	-4004	4003	4053	-4052
cell	4027	0	CuH2O	-4004	4003	4055	-4054
cell	4028	0	CuH2O	-4004	4003	4057	-4056
cell	4029	0	CuH2O	-4004	4003	4059	-4058
cell	4030	0	CuH2O	-4004	4003	4061	-4060
cell	4031	0	CuH2O	-4004	4003	4063	-4062
cell	4032	0	CuH2O	-4004	4003	4065	-4064
cell	4033	0	CuH2O	-4004	4003	4067	-4066
cell	4034	0	CuH2O	-4004	4003	4069	-4068
cell	4035	0	CuH2O	-4004	4003	4071	-4070
cell	4036	0	CuH2O	-4004	4003	4073	-4072
cell	4037	0	CuH2O	-4004	4003	4075	-4074
cell	4038	0	CuH2O	-4004	4003	4077	-4076
cell	4039	0	CuH2O	-4004	4003	4079	-4078
cell	4040	0	CuH2O	-4004	4003	4081	-4080
cell	4041	0	CuH2O	-4004	4003	4083	-4082
cell	4042	0	CuH2O	-4004	4003	4085	-4084

% ---- Vaccum holes 

cell	4043	0	void	-4004	4003	4087	-4088
cell	4044	0	void	-4004	4003	4089	-4090
cell	4045	0	void	-4004	4003	4091	-4092
cell	4046	0	void	-4004	4003	4093	-4094
cell	4047	0	void	-4004	4003	4095	-4096
cell	4048	0	void	-4004	4003	4097	-4098
cell	4049	0	void	-4004	4003	4099	-4100
cell	4050	0	void	-4004	4003	4101	-4102
cell	4051	0	void	-4004	4003	4103	-4104
cell	4052	0	void	-4004	4003	4105	-4106
cell	4053	0	void	-4004	4003	4107	-4108
cell	4054	0	void	-4004	4003	4109	-4110
cell	4055	0	void	-4004	4003	4111	-4112
cell	4056	0	void	-4004	4003	4113	-4114
cell	4057	0	void	-4004	4003	4115	-4116
cell	4058	0	void	-4004	4003	4117	-4118
cell	4059	0	void	-4004	4003	4119	-4120
cell	4060	0	void	-4004	4003	4121	-4122
cell	4061	0	void	-4004	4003	4123	-4124
cell	4062	0	void	-4004	4003	4125	-4126
cell	4063	0	void	-4004	4003	4127	-4128
cell	4064	0	void	-4004	4003	4129	-4130
cell	4065	0	void	-4004	4003	4131	-4132
cell	4066	0	void	-4004	4003	4133	-4134
cell	4067	0	void	-4004	4003	4135	-4136
cell	4068	0	void	-4004	4003	4137	-4138
cell	4069	0	void	-4004	4003	4139	-4140
cell	4070	0	void	-4004	4003	4141	-4142
cell	4071	0	void	-4004	4003	4143	-4144
cell	4072	0	void	-4004	4003	4145	-4146
cell	4073	0	void	-4004	4003	4147	-4148
cell	4074	0	void	-4004	4003	4149	-4150
cell	4075	0	void	-4004	4003	4151	-4152
cell	4076	0	void	-4004	4003	4153	-4154
cell	4077	0	void	-4004	4003	4155	-4156
cell	4078	0	void	-4004	4003	4157	-4158
cell	4079	0	void	-4004	4003	4159	-4160

cell 4080 0 steel   -4004 4003 4008 -4007 #4005 #4006 #4007 #4008 #4009 #4010 #4011 #4012 #4013 #4014 #4015 #4016 #4017 #4018 #4019 #4020 #4021 #4022 #4023 #4024 #4025 #4026 #4027 #4028 #4029 #4030 #4031 #4032 #4033 #4034 #4035 #4036 #4037 #4038 #4039 #4040 #4041 #4042 #4043 #4044 #4045 #4046 #4047 #4048 #4049 #4050 #4051 #4052 #4053 #4054 #4055 #4056 #4057 #4058 #4059 #4060 #4061 #4062 #4063 #4064 #4065 #4066 #4067 #4068 #4069 #4070 #4071 #4072 #4073 #4074 #4075 #4076 #4077 #4078 #4079

cell 4081 0 steel -4005 4004 -4007 4008
cell 4082 0 steel -4006 4005 -4007 4008 

% --- Core with TNS

cell 500 0 fill Core   -500 4006 -4007 4008
cell 501 0 outside          500 -4007 4008
cell 4084 0 outside    -4008
cell 4085 0 outside     4007

% --- Fuel Rod lattice

lat FB15 2 0.0 0.0 19 19 1.8
      40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 1
       40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 2
        40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 3
         40 40 40 40 40 40 40 40 40 32 32 32 32 32 32 32 40 40 40  %ROW 4
          40 40 40 40 40 40 40 40 32 40 40 32 32 40 40 32 40 40 40  %ROW 5
           40 40 40 40 40 40 40 32 40 40 40 32 40 40 40 32 40 40 40  %ROW 6
            40 40 40 40 40 40 32 32 40 40 32 32 40 40 32 32 40 40 40  %ROW 7
             40 40 40 40 40 32 32 32 32 32 32 32 32 32 32 32 40 40 40  %ROW 8
              40 40 40 40 32 40 40 32 32 40 40 40 32 40 40 32 40 40 40  %ROW 9
               40 40 40 32 40 40 40 32 40 40 40 32 40 40 40 32 40 40 40  %ROW 10
                40 40 40 32 40 40 32 40 40 40 32 32 40 40 32 40 40 40 40  %ROW 11
                 40 40 40 32 32 32 32 32 32 32 32 32 32 32 40 40 40 40 40  %ROW 12
                  40 40 40 32 32 40 40 32 32 40 40 32 32 40 40 40 40 40 40  %ROW 13
                   40 40 40 32 40 40 40 32 40 40 40 32 40 40 40 40 40 40 40  %ROW 14
                    40 40 40 32 40 40 32 32 40 40 32 40 40 40 40 40 40 40 40  %ROW 15
                     40 40 40 32 32 32 32 32 32 32 40 40 40 40 40 40 40 40 40  %ROW 16
                      40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 17
                       40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 18
                        40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 19

lat FB17 2 0.0 0.0 19 19 1.8
      40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 1
       40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 2
        40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 3
         40 40 40 40 40 40 40 40 40 31 31 31 31 31 31 31 40 40 40  %ROW 4
          40 40 40 40 40 40 40 40 31 40 40 31 31 40 40 31 40 40 40  %ROW 5
           40 40 40 40 40 40 40 31 40 40 40 31 40 40 40 31 40 40 40  %ROW 6
            40 40 40 40 40 40 31 31 40 40 31 31 40 40 31 31 40 40 40  %ROW 7
             40 40 40 40 40 31 31 31 31 31 31 31 31 31 31 31 40 40 40  %ROW 8
              40 40 40 40 31 40 40 31 31 40 40 40 31 40 40 31 40 40 40  %ROW 9
               40 40 40 31 40 40 40 31 40 40 40 31 40 40 40 31 40 40 40  %ROW 10
                40 40 40 31 40 40 31 40 40 40 31 31 40 40 31 40 40 40 40  %ROW 11
                 40 40 40 31 31 31 31 31 31 31 31 31 31 31 40 40 40 40 40  %ROW 12
                  40 40 40 31 31 40 40 31 31 40 40 31 31 40 40 40 40 40 40  %ROW 13
                   40 40 40 31 40 40 40 31 40 40 40 31 40 40 40 40 40 40 40  %ROW 14
                    40 40 40 31 40 40 31 31 40 40 31 40 40 40 40 40 40 40 40  %ROW 15
                     40 40 40 31 31 31 31 31 31 31 40 40 40 40 40 40 40 40 40  %ROW 16
                      40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 17
                       40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 18
                        40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 19

% --- Fuel Rod lattice

lat FB19 2 0.0 0.0 19 19 1.8
      40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 1
       40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 2
        40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 3
         40 40 40 40 40 40 40 40 40 33 33 33 33 33 33 33 40 40 40  %ROW 4
          40 40 40 40 40 40 40 40 33 40 40 33 33 40 40 33 40 40 40  %ROW 5
           40 40 40 40 40 40 40 33 40 40 40 33 40 40 40 33 40 40 40  %ROW 6
            40 40 40 40 40 40 33 33 40 40 33 33 40 40 33 33 40 40 40  %ROW 7
             40 40 40 40 40 33 33 33 33 33 33 33 33 33 33 33 40 40 40  %ROW 8
              40 40 40 40 33 40 40 33 33 40 40 40 33 40 40 33 40 40 40  %ROW 9
               40 40 40 33 40 40 40 33 40 40 40 33 40 40 40 33 40 40 40  %ROW 10
                40 40 40 33 40 40 33 40 40 40 33 33 40 40 33 40 40 40 40  %ROW 11
                 40 40 40 33 33 33 33 33 33 33 33 33 33 33 40 40 40 40 40  %ROW 12
                  40 40 40 33 33 40 40 33 33 40 40 33 33 40 40 40 40 40 40  %ROW 13
                   40 40 40 33 40 40 40 33 40 40 40 33 40 40 40 40 40 40 40  %ROW 14
                    40 40 40 33 40 40 33 33 40 40 33 40 40 40 40 40 40 40 40  %ROW 15
                     40 40 40 33 33 33 33 33 33 33 40 40 40 40 40 40 40 40 40  %ROW 16
                      40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 17
                       40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 18
                        40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 19

lat FB21 2 0.0 0.0 19 19 1.8
      40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 1
       40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 2
        40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 3
         40 40 40 40 40 40 40 40 40 34 34 34 34 34 34 34 40 40 40  %ROW 4
          40 40 40 40 40 40 40 40 34 40 40 34 34 40 40 34 40 40 40  %ROW 5
           40 40 40 40 40 40 40 34 40 40 40 34 40 40 40 34 40 40 40  %ROW 6
            40 40 40 40 40 40 34 34 40 40 34 34 40 40 34 34 40 40 40  %ROW 7
             40 40 40 40 40 34 34 34 34 34 34 34 34 34 34 34 40 40 40  %ROW 8
              40 40 40 40 34 40 40 34 34 40 40 40 34 40 40 34 40 40 40  %ROW 9
               40 40 40 34 40 40 40 34 40 40 40 34 40 40 40 34 40 40 40  %ROW 10
                40 40 40 34 40 40 34 40 40 40 34 34 40 40 34 40 40 40 40  %ROW 11
                 40 40 40 34 34 34 34 34 34 34 34 34 34 34 40 40 40 40 40  %ROW 12
                  40 40 40 34 34 40 40 34 34 40 40 34 34 40 40 40 40 40 40  %ROW 13
                   40 40 40 34 40 40 40 34 40 40 40 34 40 40 40 40 40 40 40  %ROW 14
                    40 40 40 34 40 40 34 34 40 40 34 40 40 40 40 40 40 40 40  %ROW 15
                     40 40 40 34 34 34 34 34 34 34 40 40 40 40 40 40 40 40 40  %ROW 16
                      40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 17
                       40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 18
                        40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 19

lat FB23 2 0.0 0.0 19 19 1.8
      40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 1
       40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 2
        40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 3
         40 40 40 40 40 40 40 40 40 35 35 35 35 35 35 35 40 40 40  %ROW 4
          40 40 40 40 40 40 40 40 35 40 40 35 35 40 40 35 40 40 40  %ROW 5
           40 40 40 40 40 40 40 35 40 40 40 35 40 40 40 35 40 40 40  %ROW 6
            40 40 40 40 40 40 35 35 40 40 35 35 40 40 35 35 40 40 40  %ROW 7
             40 40 40 40 40 35 35 35 35 35 35 35 35 35 35 35 40 40 40  %ROW 8
              40 40 40 40 35 40 40 35 35 40 40 40 35 40 40 35 40 40 40  %ROW 9
               40 40 40 35 40 40 40 35 40 40 40 35 40 40 40 35 40 40 40  %ROW 10
                40 40 40 35 40 40 35 40 40 40 35 35 40 40 35 40 40 40 40  %ROW 11
                 40 40 40 35 35 35 35 35 35 35 35 35 35 35 40 40 40 40 40  %ROW 12
                  40 40 40 35 35 40 40 35 35 40 40 35 35 40 40 40 40 40 40  %ROW 13
                   40 40 40 35 40 40 40 35 40 40 40 35 40 40 40 40 40 40 40  %ROW 14
                    40 40 40 35 40 40 35 35 40 40 35 40 40 40 40 40 40 40 40  %ROW 15
                     40 40 40 35 35 35 35 35 35 35 40 40 40 40 40 40 40 40 40  %ROW 16
                      40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 17
                       40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 18
                        40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40  %ROW 19

% --- Upper and boottom Reflector lattice

lat TopRef 3 0.0 0.0 19 19 20.8
      71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71  %ROW 1 
      71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71  %ROW 2 
      71 71 71 71 71 71 71 71 71 61 61 61 61 61 61 61 61 71 71  %ROW 3 
      71 71 71 71 71 71 71 71 61 61 61 61 61 61 61 61 61 71 71  %ROW 4 
      71 71 71 71 71 71 71 61 61 58 58 58 58 58 58 61 61 71 71  %ROW 5 
      71 71 71 71 71 71 61 61 58 58 58 58 58 58 58 61 61 71 71  %ROW 6 
      71 71 71 71 71 61 61 58 58 58 58 58 58 58 58 61 61 71 71  %ROW 7 
      71 71 71 71 61 61 58 58 58 58 58 58 58 58 58 61 61 71 71  %ROW 8 
      71 71 71 61 61 58 58 58 58 58 58 58 58 58 58 61 61 71 71  %ROW 9 
      71 71 61 61 58 58 58 58 58 58 58 58 58 58 58 61 61 71 71  %ROW 10 
      71 71 61 61 58 58 58 58 58 58 58 58 58 58 61 61 71 71 71  %ROW 11 
      71 71 61 61 58 58 58 58 58 58 58 58 58 61 61 71 71 71 71  %ROW 12 
      71 71 61 61 58 58 58 58 58 58 58 58 61 61 71 71 71 71 71  %ROW 13 
      71 71 61 61 58 58 58 58 58 58 58 61 61 71 71 71 71 71 71  %ROW 14 
      71 71 61 61 58 58 58 58 58 58 61 61 71 71 71 71 71 71 71  %ROW 15 
      71 71 61 61 61 61 61 61 61 61 61 71 71 71 71 71 71 71 71  %ROW 16 
      71 71 61 61 61 61 61 61 61 61 71 71 71 71 71 71 71 71 71  %ROW 17 
      71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71  %ROW 18 
      71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71  %ROW 19 

%lat TopRef 3 0.0 0.0 19 19 20.8
%      71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71  %ROW 1 
%      71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71  %ROW 2 
%      71 71 71 71 71 71 71 71 71 61 61 61 61 61 61 61 61 71 71  %ROW 3 
%      71 71 71 71 71 71 71 71 61 61 81 61 81 61 81 61 61 71 71  %ROW 4 
%      71 71 71 71 71 71 71 61 81 58 58 58 58 58 58 81 61 71 71  %ROW 5 
%      71 71 71 71 71 71 61 61 58 58 58 58 58 58 58 61 61 71 71  %ROW 6 
%      71 71 71 71 71 61 81 58 58 58 58 58 58 58 58 81 61 71 71  %ROW 7 
%      71 71 71 71 61 61 58 58 58 58 58 58 58 58 58 61 61 71 71  %ROW 8 
%      71 71 71 61 81 58 58 58 58 58 58 58 58 58 58 81 61 71 71  %ROW 9 
%      71 71 61 61 58 58 58 58 58 58 58 58 58 58 58 61 61 71 71  %ROW 10 
%      71 71 61 81 58 58 58 58 58 58 58 58 58 58 81 61 71 71 71  %ROW 11 
%      71 71 61 61 58 58 58 58 58 58 58 58 58 61 61 71 71 71 71  %ROW 12 
%      71 71 61 81 58 58 58 58 58 58 58 58 81 61 71 71 71 71 71  %ROW 13 
%      71 71 61 61 58 58 58 58 58 58 58 61 61 71 71 71 71 71 71  %ROW 14 
%      71 71 61 81 58 58 58 58 58 58 81 61 71 71 71 71 71 71 71  %ROW 15 
%      71 71 61 61 81 61 81 61 81 61 61 71 71 71 71 71 71 71 71  %ROW 16 
%      71 71 61 61 61 61 61 61 61 61 71 71 71 71 71 71 71 71 71  %ROW 17 
%      71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71  %ROW 18 
%      71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71 71  %ROW 19 
      
%% --- Reactor Core lattice w/o abs rods
%
lat CoreLayer 3 0.0 0.0 19 19 20.8
      70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70  %ROW 1 
      70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70  %ROW 2 
      70 70 70 70 70 70 70 70 70 60 60 60 60 60 60 60 60 70 70  %ROW 3 
      70 70 70 70 70 70 70 70 60 60 60 60 60 60 60 60 60 70 70  %ROW 4 
      70 70 70 70 70 70 70 60 60 51 51 51 51 51 51 60 60 70 70  %ROW 5 
      70 70 70 70 70 70 60 60 51 55 54 54 54 55 51 60 60 70 70  %ROW 6 
      70 70 70 70 70 60 60 51 54 51 52 52 51 54 51 60 60 70 70  %ROW 7 
      70 70 70 70 60 60 51 54 52 57 57 57 52 54 51 60 60 70 70  %ROW 8 
      70 70 70 60 60 51 54 52 57 57 57 57 52 54 51 60 60 70 70  %ROW 9 
      70 70 60 60 51 55 51 57 57 57 57 57 51 55 51 60 60 70 70  %ROW 10 
      70 70 60 60 51 54 52 57 57 57 57 52 54 51 60 60 70 70 70  %ROW 11 
      70 70 60 60 51 54 52 57 57 57 52 54 51 60 60 70 70 70 70  %ROW 12 
      70 70 60 60 51 54 51 52 52 51 54 51 60 60 70 70 70 70 70  %ROW 13 
      70 70 60 60 51 55 54 54 54 55 51 60 60 70 70 70 70 70 70  %ROW 14 
      70 70 60 60 51 51 51 51 51 51 60 60 70 70 70 70 70 70 70  %ROW 15 
      70 70 60 60 60 60 60 60 60 60 60 70 70 70 70 70 70 70 70  %ROW 16 
      70 70 60 60 60 60 60 60 60 60 70 70 70 70 70 70 70 70 70  %ROW 17 
      70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70  %ROW 18 
      70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70  %ROW 19 

% --- Reactor Core lattice with abs rods

%lat CoreLayer 3 0.0 0.0 19 19 20.8
%      70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70  %ROW 1 
%      70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70  %ROW 2 
%      70 70 70 70 70 70 70 70 70 60 60 60 60 60 60 60 60 70 70  %ROW 3 
%      70 70 70 70 70 70 70 70 60 60 80 60 80 60 80 60 60 70 70  %ROW 4 
%      70 70 70 70 70 70 70 60 80 51 51 51 51 51 51 80 60 70 70  %ROW 5 
%      70 70 70 70 70 70 60 60 51 55 54 54 54 55 51 60 60 70 70  %ROW 6 
%      70 70 70 70 70 60 80 51 54 51 52 52 51 54 51 80 60 70 70  %ROW 7 
%      70 70 70 70 60 60 51 54 52 90 90 90 52 54 51 60 60 70 70  %ROW 8 
%      70 70 70 60 80 51 54 52 90 90 90 90 52 54 51 80 60 70 70  %ROW 9 
%      70 70 60 60 51 55 51 90 90 90 90 90 51 55 51 60 60 70 70  %ROW 10 
%      70 70 60 80 51 54 52 90 90 90 90 52 54 51 80 60 70 70 70  %ROW 11 
%      70 70 60 60 51 54 52 90 90 90 52 54 51 60 60 70 70 70 70  %ROW 12 
%      70 70 60 80 51 54 51 52 52 51 54 51 80 60 70 70 70 70 70  %ROW 13 
%      70 70 60 60 51 55 54 54 54 55 51 60 60 70 70 70 70 70 70  %ROW 14 
%      70 70 60 80 51 51 51 51 51 51 80 60 70 70 70 70 70 70 70  %ROW 15 
%      70 70 60 60 80 60 80 60 80 60 60 70 70 70 70 70 70 70 70  %ROW 16 
%      70 70 60 60 60 60 60 60 60 60 70 70 70 70 70 70 70 70 70  %ROW 17 
%      70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70  %ROW 18 
%      70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70  %ROW 19 

%///////////////////////////////////////////////////////////////////////////////
% ############################# SURFACES #######################################
%///////////////////////////////////////////////////////////////////////////////

surf 100 inf

% ---- Fuel Pellets and Fuel Block surfaces

surf 201  cyl 0.0 0.0 0.47850  % Pellet 1050
surf 202  cyl 0.0 0.0 0.50850  % SiC 300 mkm ZrB2/SiC

% ---- Upper and Botoom Reflector surfaces

surf 301 hexyprism 0.0 0.0 10.4 0 80
surf 303 hexyprism 0.0 0.0 10.4 0 80
surf 302 hexyprism 0.0 0.0 10.4 0 30

% ---- Coolant Surfaces

surf 311    cyl       0.0     0.0   1.2      % Coolant Hole 1
surf 312    cyl      -7.2     0.0   1.2      % Coolant Hole 2
surf 313    cyl       7.2     0.0   1.2      % Coolant Hole 3
surf 314    cyl      -3.6    -6.3   1.2      % Coolant Hole 4
surf 315    cyl       3.6    -6.3   1.2      % Coolant Hole 5
surf 316    cyl      -3.6     6.3   1.2      % Coolant Hole 6
surf 317    cyl       3.6     6.3   1.2      % Coolant Hole 7

surf 320    cyl       0.0     0.0   1.2      % Absorber Hole

% ---- Core and Reflector surfaces

surf 401 hexxprism 0.0 0.0 160   0 80
surf 402 hexxprism 0.0 0.0 160   0 30

%%% - TIN %%%%

surf 411   cylz 0 0 17.0
surf 412   cylz 0 0 18.5
surf 413   cylz 0 0 19.5 
surf 414   cylz 0 0 20.0
surf 415   cylz 0 0 30.0 

%%% - TIN + Solenoid %%%%

% -- surf 4000   cone 0 0 0 17 2247 

surf 	4000   	cyl 	0 0 17.0
surf 	4001   	cyl 	0 0 18.5
surf 	4002   	cyl 	0 0 19.5 
surf 	4003   	cyl 	0 0 20.5
surf 	4004   	cyl 	0 0 28.5
surf 	4005   	cyl 	0 0 29.5
surf 	4006   	cyl 	0 0 30.0 

surf 	4007   	pz  	300
surf 	4008   	pz  	0

surf	4009	pz		152
surf	4010	pz		148
surf	4012	pz		153.14
surf	4013	pz		154.86

%%% - TIN + Solenoid Coils %%%%

surf	4014	pz	160
surf	4015	pz	156
surf	4016	pz	168
surf	4017	pz	164
surf	4018	pz	176
surf	4019	pz	172
surf	4020	pz	184
surf	4021	pz	180
surf	4022	pz	192
surf	4023	pz	188
surf	4024	pz	200
surf	4025	pz	196
surf	4026	pz	208
surf	4027	pz	204
surf	4028	pz	216
surf	4029	pz	212
surf	4030	pz	224
surf	4031	pz	220
surf	4032	pz	232
surf	4033	pz	228
surf	4034	pz	240
surf	4035	pz	236
surf	4036	pz	248
surf	4037	pz	244
surf	4038	pz	256
surf	4039	pz	252
surf	4040	pz	264
surf	4041	pz	260
surf	4042	pz	272
surf	4043	pz	268
surf	4044	pz	280
surf	4045	pz	276
surf	4046	pz	288
surf	4047	pz	284
surf	4048	pz	296
surf	4049	pz	292
surf	4050	pz	144
surf	4051	pz	140
surf	4052	pz	136
surf	4053	pz	132
surf	4054	pz	128
surf	4055	pz	124
surf	4056	pz	120
surf	4057	pz	116
surf	4058	pz	112
surf	4059	pz	108
surf	4060	pz	104
surf	4061	pz	100
surf	4062	pz	96
surf	4063	pz	92
surf	4064	pz	88
surf	4065	pz	84
surf	4066	pz	80
surf	4067	pz	76
surf	4068	pz	72
surf	4069	pz	68
surf	4070	pz	64
surf	4071	pz	60
surf	4072	pz	56
surf	4073	pz	52
surf	4074	pz	48
surf	4075	pz	44
surf	4076	pz	40
surf	4077	pz	36
surf	4078	pz	32
surf	4079	pz	28
surf	4080	pz	24
surf	4081	pz	20
surf	4082	pz	16
surf	4083	pz	12
surf	4084	pz	8
surf	4085	pz	4
surf	4086	pz	0

%%% - TIN + Vaccuum Coils %%%%

surf	4087	pz	161.14
surf	4088	pz	162.86
surf	4089	pz	169.14
surf	4090	pz	170.86
surf	4091	pz	177.14
surf	4092	pz	178.86
surf	4093	pz	185.14
surf	4094	pz	186.86
surf	4095	pz	193.14
surf	4096	pz	194.86
surf	4097	pz	201.14
surf	4098	pz	202.86
surf	4099	pz	209.14
surf	4100	pz	210.86
surf	4101	pz	217.14
surf	4102	pz	218.86
surf	4103	pz	225.14
surf	4104	pz	226.86
surf	4105	pz	233.14
surf	4106	pz	234.86
surf	4107	pz	241.14
surf	4108	pz	242.86
surf	4109	pz	249.14
surf	4110	pz	250.86
surf	4111	pz	257.14
surf	4112	pz	258.86
surf	4113	pz	265.14
surf	4114	pz	266.86
surf	4115	pz	273.14
surf	4116	pz	274.86
surf	4117	pz	281.14
surf	4118	pz	282.86
surf	4119	pz	289.14
surf	4120	pz	290.86
surf	4121	pz	297.14
surf	4122	pz	298.86
surf	4123	pz	145.14
surf	4124	pz	146.86
surf	4125	pz	137.14
surf	4126	pz	138.86
surf	4127	pz	129.14
surf	4128	pz	130.86
surf	4129	pz	121.14
surf	4130	pz	122.86
surf	4131	pz	113.14
surf	4132	pz	114.86
surf	4133	pz	105.14
surf	4134	pz	106.86
surf	4135	pz	97.14
surf	4136	pz	98.86
surf	4137	pz	89.14
surf	4138	pz	90.86
surf	4139	pz	81.14
surf	4140	pz	82.86
surf	4141	pz	73.14
surf	4142	pz	74.86
surf	4143	pz	65.14
surf	4144	pz	66.86
surf	4145	pz	57.14
surf	4146	pz	58.86
surf	4147	pz	49.14
surf	4148	pz	50.86
surf	4149	pz	41.14
surf	4150	pz	42.86
surf	4151	pz	33.14
surf	4152	pz	34.86
surf	4153	pz	25.14
surf	4154	pz	26.86
surf	4155	pz	17.14
surf	4156	pz	18.86
surf	4157	pz	9.14
surf	4158	pz	10.86
surf	4159	pz	1.14
surf	4160	pz	2.86

surf 500   hexxprism 0.0 0.0 160   0 300

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

% Material data:

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

% -- Fuel:
mat fuel      -10.4 tmp 1520 rgb 166 0 0 %burn 1

%wAtFr_crit(Gd2O3)=0.16 //// wAtFr_0.95(Gd2O3)=0.298
64154.15c	0.00274
64155.15c	0.01763
64156.15c	0.02462
64157.15c	0.01870
64158.15c	0.02959
64160.15c	0.02601
90232.15c	0.20200
94239.15c	0.47000
94240.15c	0.02500
94241.15c	0.00500
8016.15c	2.17870

mat fuel1      -10.4 tmp 1520 rgb 166 0 0 % burn 1 vol 2.17E+05

 90232.15c     0.500
 94239.15c     0.470
 94240.15c     0.025
 94241.15c     0.005
  8016.15c     2

mat fuel2      -10.4 tmp 1520 rgb 166 0 0 % burn 1 vol 2.17E+05 % Inner Fuel Block Fuel

 90232.15c     0.500
 94239.15c     0.470
 94240.15c     0.025
 94241.15c     0.005
  8016.15c     2

mat fuel3      -10.4 tmp 1520 rgb 166 0 0 % burn 1 vol 2.17E+05 % Inner Fuel Block Fuel

 90232.15c     0.500
 94239.15c     0.470
 94240.15c     0.025
 94241.15c     0.005
  8016.15c     2

mat fuel4      -10.4 tmp 1520 rgb 166 0 0 % burn 1 vol 2.17E+05 % Inner Fuel Block Fuel

 90232.15c     0.500
 94239.15c     0.470
 94240.15c     0.025
 94241.15c     0.005
  8016.15c     2

mat fuel5      -10.4 tmp 1520 rgb 166 0 0 % burn 1 vol 2.17E+05 % Inner Fuel Block Fuel

 90232.15c     0.500
 94239.15c     0.470
 94240.15c     0.025
 94241.15c     0.005
  8016.15c     2

mat BPGd -6.0 tmp 1500 burn 1  % p crit = 6.0 g/cm3

 64154.15c     0.0436   %2*0.0218
 64155.15c     0.2960   %2*0.148
 64156.15c     0.4094   %2*0.2047
 64157.15c     0.3130   %2*0.1565
 64158.15c     0.4968   %2*0.2484
 64160.15c     0.4372   %2*0.2186
  8016.15c     3

% --- Carbon buffer layer:

mat PyC       -1.2 tmp 1520 rgb 192 192 192

 6000.15c      1
 
% --- Ti3SiC2:

mat Coating   -4.2 tmp 1520 rgb 32 32 32

 22000.15c     3
 14000.15c     1
  6000.15c     2
 
% --- Graphite graphite:

mat graphite    -2.2 moder grph 6000 tmp 1370 rgb 158 158 158 % Fuel Pellet graphite

 6000.12c      1
 
mat matrix   -2.2 moder grph 6000 tmp 1370 rgb 166 0 0        % Kern graphite

 6000.12c      1

mat ingraphite   -2.2 moder grph 6000 tmp 1370 rgb 204 102 0 % Inner Fuel Block

 6000.12c      1

%///////////////////////////////////////////////////////////////////////////////
% ############################# REFLECTOR MAT ##################################
%///////////////////////////////////////////////////////////////////////////////
 
mat reflector    -2.2 moder grph 6000 tmp 1370 rgb 158 158 158 % reflector graphite

 6000.12c      1
 
mat reflector_coolant   2.65156E-05 rgb 0 0  169 % reflector helium

 2000.12c    2.65156E-05
 
% --- Silicon carbide layer:

mat SiC       -3.18 tmp 1370 rgb 32 32 32

 14000.12c    1
  6000.12c    1
 
% --- ZrB2 Layer:

mat zrb2      -6.08 tmp 1370 rgb 32 32 32 % burn 1

  5010.12c     1.33E-01 
  5011.12c     5.33E-01
 40000.12c     3.33E-01 

% --- Thermal scattering data:

therm grph gre7.22t

mat helium   2.65156E-05 rgb 0 0  169

 2000.12c    2.65156E-05
 5010.12c    1.33E-01
 5011.12c    5.33E-01
 
%% TIN MAT

% -- Plasma:

mat plasma     5E-11

 1002.18c      0.5
 1003.18c      0.5

% --- 08X10HT:

mat steel       8.633E-02

  6000.06c     -8.00E-04
 14000.06c     -8.00E-03
 25055.06c     -2.00E-02
 28000.06c     -9.00E-02
 16000.06c     -2.00E-04
 15031.06c     -3.50E-04
 24000.06c     -1.90E-01
 29000.06c     -3.00E-03
 26000.06c     -6.88E-01
 
% --- Air:

mat air        -0.0006157
             
  6000.06c     -0.000124     
  7014.06c     -0.755268     
  8016.06c     -0.231781     
 18040.06c     -0.012827     
 
% --- Cu + H2O:

mat CuH2O       0.08242

 29000.06c     -0.9860692
  1001.06c     -0.0015479
  8016.06c     -0.0123829
  
% ------ Control Rod (total incorrect in specs) -----

mat absorber   9.8481E-02 moder grph 6000 tmp 1370 rgb 255 0 0

  5010.12c      6.8220E-03
  5011.12c      2.4974E-02
  6000.12c      6.6685E-02

% ------ Graphite in Control Rod Block  -----

mat absorber_graphite    -2.2 moder grph 6000 tmp 1370 rgb 158 158 158 % Fuel Pellet graphite

 6000.12c      1
 
% -----------------------------------------------------------------------------

% Calculation parameters:

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

% --- Geometry plot:

plot 30 5000 5000
plot 20 1000 1000

mesh 3 1000 1000   %1
mesh 2 1000 1000   %2

det thermal_n   dr -1 void de thermal
det resonance_n dr -1 void de resonance
det fast_n      dr -1 void de fast
det ffast_n     dr -1 void de ffast

ene thermal   1 0.0 0.625e-06  
ene resonance 1 0.625e-06 0.1
ene fast      1 0.1 20
ene ffast     1 2 20

% --- Fuel Spectra

set cpd 3

det CoreFissRate
dl CoreLayer
dr -6 void

det PINFissRate15
dl FB15
dr -6 void

det PINFissRate17
dl FB17
dr -6 void

%det PINFissRate19
%dl FB19
%dr -6 void

det PINFissRate21
dl FB21
dr -6 void

det PINFissRate23
dl FB23
dr -6 void

det FissionEnergy                       
dr -8 void
dn 1 0.0 160 100 0 360 1 0 300 30

det FluxRadialDistr                         
dr -1 void
dn 1 0.0 160 100 0 360 1 0 300 30

det ThermalFluxRadialDistr                         
dr -1 void
de thermal
dn 1 0.0 160 100 0 360 1 0 300 30

det ResonanceFluxRadialDistr                         
dr -1 void
de resonance
dn 1 0.0 160 100 0 360 1 0 300 30

det FastFluxRadialDistr                         
dr -1 void
de fast
dn 1 0.0 160 100 0 360 1 0 300 30

det FuelSpectraIn172
dm  fuel1
de  wms172
ene wms172 4 wms172 

det FuelSpectra
dm  fuel
de  wms172

det FuelSpectraIn69
dm  fuel1
de  wms69
ene wms69 4 wms69 

det FuelSpectraOut69
dm  fuel
de  wms69

det FisRateIn
dm  fuel1                               
dr -6 void 

det FisRateOut
dm  fuel                               
dr -6 void

det CaptRateIn
dm  fuel1                              
dr -2 void

det CaptRateOut
dm  fuel                              
dr -2 void

det Flux
dr -15 void

det 1 di global du 0
det 2 di promt  du 0
det 3 di delay  du 0

% --- Parametres:

%set nps 1000000
%src DT n sc 4001 se 14.1 % Source in Plasma ZONE CELL 4000 14.1 MeV
%set srcrate 4.5e17

set nps 1000000000 200 promt
src DT n sc 4001 se 14.1 sz 0 300 srad 0 17 st 0.0 1.16E-02 % Source in Plasma ZONE CELL 4001 14.1 MeV
set srcrate 4.5e17
set bc 1
set gcu -1

set tcut 3600

tme global 2 2000 0.0 2000 
tme promt  2 4000 0.0 1.16E-02
tme delay  2 500 1.16E-02 1000

set delnu 1

%set keff 0.95

set bala 1
set nbuf 100000
set impl 0 0 0

%set pop 25000 250 30 1
%set gcu -1
%set srcrate 4.5e17
%set bc 1

set mcvol 100000000

% --- Libraries:

set acelib "/home/iol4/serpent2/xs/endfb7/sss_endfb7u.xsdata"
set declib "/home/iol4/serpent2/xs/endfb7/sss_endfb7.dec"
set nfylib "/home/iol4/serpent2/xs/endfb7/sss_endfb7.nfy"

%% --- Nuclide inventory:

%set inventory all
%dep butot
%       0.1
%       0.2
%       0.3
%       0.5
%       1
%       2
%       3
%       5
%       10
%       20
%       30
%       40
%       50
%       60
%       70
%       80
%       90
%       100
%       110
%       120
%       130
%       140
%       150
%       160
%       170
%       180
%       190
%       200
%       210
%       220
%       230
%       240
%%%       250
%%       260
%%       270
%%       280
%%       290
%%       300
%%       310
%%       320
%%       330
%%       340
%%       350
%%       360
%%       370
%%       380
%%       390
%%       400
Pbed Info:

Code: Select all

particle 1
pbed info: % Packing fraction = 0.17 Cylinder r=0.4785, z0=0, z1=80 cm particle universe=1 particle radius = 0.03 cm
particle 2
pbed info: % Packing fraction = 0.13 Cylinder r=0.4785, z0=0, z1=80 cm particle universe=2 particle radius = 0.03 cm
particle 3
pbed info: % Packing fraction = 0.19 Cylinder r=0.4785, z0=0, z1=80 cm particle universe=3 particle radius = 0.03 cm
particle 4
pbed info: % Packing fraction = 0.21 Cylinder r=0.4785, z0=0, z1=80 cm particle universe=4 particle radius = 0.03 cm
particle 5
pbed info: % Packing fraction = 0.23 Cylinder r=0.4785, z0=0, z1=80 cm particle universe=5 particle radius = 0.03 cm
- Igor

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

Re: Detector time binning

Post by Jaakko Leppänen » Mon Nov 25, 2019 7:20 pm

The problem is that your source particles are emitted between 0 and 1.16E-2 seconds, but the time bin width is only 2.9E-6 s. So the probability of sampling a source point on that time interval is a bit low.

You can fix this issue by setting the st entry in your source card into "st 0.0 2.9E-6", but the root cause is that Serpent expects that when the sequential population control is used in a time-dependent problem, all your source particles must be born on the first interval.
- Jaakko

Post Reply