Weight windows in Gamma transport

Separate section for discussion on gamma transport
Post Reply
Zsolti Elter
Posts: 2
Joined: Tue Oct 08, 2019 11:59 am
Security question 1: No
Security question 2: 19

Weight windows in Gamma transport

Post by Zsolti Elter » Tue Oct 08, 2019 12:16 pm

Hello,

I have been trying the weight window capabilities in 2.1.31. The example at http://serpent.vtt.fi/mediawiki/index.p ... _reduction works fine, however when I adapt it for a gamma problem I encounter the following error:

Fatal error in function GSconfigDataRead:

Unable to open file /home/elter/md7/gsconfig.dat for reading

Simulation aborted.

I am a bit uncertain what is this file. Could you please guide me?

Previously I have tried weight windows with 2.1.28, and I managed to get it work for gammas, and no such file was required.

Thank you!

zsolti

Zsolti Elter
Posts: 2
Joined: Tue Oct 08, 2019 11:59 am
Security question 1: No
Security question 2: 19

Re: Weight windows in Gamma transport

Post by Zsolti Elter » Wed Oct 09, 2019 8:18 pm

Hello,

sorry, I have overlooked the file while migrating from one version to an other. Now things work:)

best regards,

zsolti

kaartinen
Posts: 8
Joined: Wed Aug 01, 2018 7:53 am
Security question 1: No
Security question 2: 93
Location: Loviisa NPP

Re: Weight windows in Gamma transport

Post by kaartinen » Mon Oct 14, 2019 3:30 pm

I've had somewhat good experience with Serpent's weight window capabilites, especially the adaptive ww mesh.
Here is an example input that pretty much shows the power of adaptive ww mesh.
The case is a simple concrete labyrinth with a point Co-60 source. Running with only 1 core this input is able to get you a ww mesh over the whole geometry in less than 3 minutes.

If you try this make sure you change material id to correspond with your foton libraries

Code: Select all

% --------- SURFACES ------------ %
% a simple labyrinth
surf 1 cuboid 0 780 0 620 0 200
surf 2 cuboid 420 720 60 560 0 200
surf 3 cuboid 240 420 440 560 0 200
surf 4 cuboid 240 360 60 440 0 200
surf 5 cuboid 60 240 60 180 0 200
surf 6 suboid 60 180 180 440 0 200
surf 7 cuboid 0 180 440 560 0 200

% ---------- CELLS -----------------%
cell 999 0 outside 1 % outside
cell 1 0 M1 -1 2 3 4 5 6 7 % concrete
cell 2 0 M5 -2 % the rest is air
cell 3 0 M5 -3
cell 4 0 M5 -4
cell 5 0 M5 -5
cell 6 0 M6 -6
cell 7 0 M5 -7

% ------------- SOURCE ------------ %
% Co-60 point source with activity of 1E12 Bq
src 1 p sp 700 80 100
sb 7 0
0	0
3.471E-1	7.500E+7
8.261E-1	7.600E+7
1.173E+0	9.985E+11
1.332E+0	9.985E+11
2.159E+0	1.200E+7
2.506E+0	2.000E+4

set srcrate 1.998E+12

%---------- DETECTORS -----------------%
det 1 p	% a simple detector for flux mesh
det 2 p dr -200 void dv 1000 % dose rate mesh (Gy/h)
dx 0 780 78  % 10cm x 10cm x 10cm mesh
dy 620 20
dz 0 200 20
det 3 p dm M5  % a simple detector for flux mesh in air

% ---------- PLOT -------------------------%
plot 35 1E-10 1E10 -1 780 620 100 % plot weight window mesh
mesh 8 -13 1 3 780 620 % plot photon flux mesh
mesh 8 -13 3 3 780 620 % same but only in air

% ------------ SIMULATION CONTROL --------------%
set nps 100000 200	% yes, only 100k
set dt 0.0 0.0		% taking delta tracking away, its not use here
set cfe 1E18 -1 1E18 -1 % take virtual collisions away, they only slow us down
set bala 1		% doesnt really affect us if we use only 1 core

% --------- VARIANCE REDUCTION --------------%
wwgen 1 1E-10 1000 3 -1 1 % cartesian mesh
-10 790 10 % cover whole geometry
-10 630 10
-10 210 1

wwin 1
wi 3 6 1 % 6 iteration, tracking based adaptive
2 2 1 % only split in X and Y directions
1 1 % take adaptivity through material density off
10 50 % split cells if importance gradient bigger than 10
-10 10 10 200 % at smallest cells can be 10x10x200

%----- MATERIALS ----%
mat M1 -2.35 rgb 140 140 140 % CONCRETE
26000.84p 3.10E-4
1000.84p 7.86E-3
8000.84p 4.38E-2
11000.84p 1.05E-3
12000.84p 1.4E-4
13000.84p 2.39E-3
14000.84p 1.58E-2
19000.84p 6.9E-4
20000.84p 2.92E-2

mat M5 -0.001225 rgb 230 232 247 % AIR
7000.84p -0.755636
8000.84p -0.231475
18000.84p -0.012845

Post Reply