ESS-Estia/Estia-MCNP
1
0
Fork 0
Code Actions Packages Releases Activity

Generate WW for detector position and test fine meshes behind collimation and bunker

Browse Source
This commit is contained in:
glavic_a
2019-06-28 15:50:19 +02:00
parent 480d93162a
commit 75bbb7cda0
4 changed files with 79 additions and 43 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
e2-estia.i
View File

@@ -106,7 +106,7 @@ c JINTS= 2 4 4 8 8 4 4 6
c KMESH=298 303 309 315 360
c KINTS= 14 1 20 1 2
c generate cell based weight windows, needs to be updated on geometry change
WWG 54 61003 0
WWG 255 61003 0
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
c -----------------------------------------------------------
c ------------------- TALLY CARDS ---------------------------
@@ -125,12 +125,13 @@ c read file=tallies/inbunker_f4.i echo
c read file=tallies/chopper_pit_doses.i echo
c read file=tallies/bunker_wall_f4.i echo
c read file=tallies/bunker_wall_f5.i echo
read file=tallies/pre_selene1_mesh.i
c read file=tallies/e1_e2_f4.i echo
c read file=tallies/e1_e2_f5.i echo
c read file=tallies/shutter_f4.i echo
c read file=tallies/shutter_f5.i echo
read file=tallies/shutter_wall_f4.i echo
c read file=tallies/shutter_sample_f5.i echo
c read file=tallies/shutter_wall_f4.i echo
read file=tallies/detector_f5.i echo
c read file=tallies/selene1_walls_f4.i echo
c read file=tallies/optical_cave_f5.i echo
c read file=tallies/optical_cave_gammas.i echo
@@ -167,7 +168,7 @@ c --------------- PHYSICS CARDS --------------------------
c -------------------------------------------------------
c dbcn 375642321 j j j j j j j j j j j 15291711
c stop calculation either at 2e8 histories, 6.5 days or tally 15 error <0.001
stop nps 1e10 ctme 4.5e6 f54 0.05
stop nps 1e10 ctme 4.5e6 f255 0.05
rand gen=2 seed=29976355590901 $ default 19073486328125
c histp
mode n p h / z $ neutron, photon, proton, pi+, pi0

25
tallies/detector_f5.i Normal file
View File

@@ -0,0 +1,25 @@
c f5 detector tally to estimate background at the detector
fc255 neutron flux at detector position
f255:n 2280 -3180 -30 0.0
e255 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
1e-0 1.58 2.51 3.98 6.31
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
1e+3 1.58e+3 2.51e+3
fq255 e f
c reduce probability for contribution inside of bunker
pd255 0.001 276r $ in monolith
0.005 297r $ ess outside
0.001 0.001 78r $ NBOA to NFGA
0.05 75r $ up to after CPC1
0.1 34r $ up to after CPC2
0.25 91r $ up to end

49
tallies/pre_selene1_mesh.i Normal file
View File

@@ -0,0 +1,49 @@
c Fine 2D mesh behind bunker wall collimator upstream of Selene 1
tmesh
rmesh1:n dose 1 trans 1
cora1 1335 0i 1336
corb1 -100 100i 100
corc1 -150 150i 150
c second mesh just outside bunker wall
rmesh11:n dose 1 trans 1
cora11 1500 0i 1501
corb11 -100 100i 100
corc11 -150 150i 150
c Official ESS dose conversion factors from n/cm²/s to Sv/h
c Is described in ESS-0019931
mshmf1 0.000000000 1.00000e-08
1.00000e-09 1.11000e-08 1.00000e-08 1.28000e-08
2.50000e-08 1.44000e-08 1.00000e-07 1.87000e-08
2.00000e-07 2.11000e-08 5.00000e-07 2.37000e-08
1.00000e-06 2.53000e-08 2.00000e-06 2.66000e-08
5.00000e-06 2.78000e-08 1.00000e-05 2.82000e-08
2.00000e-05 2.82000e-08 5.00000e-05 2.82000e-08
1.00000e-04 2.80000e-08 2.00000e-04 2.78000e-08
5.00000e-04 2.71000e-08 1.00000e-03 2.71000e-08
2.00000e-03 2.74000e-08 5.00000e-03 2.87000e-08
1.00000e-02 3.28000e-08 2.00000e-02 4.39000e-08
3.00000e-02 5.65000e-08 5.00000e-02 8.28000e-08
7.00000e-02 1.10000e-07 1.00000e-01 1.51000e-07
1.50000e-01 2.18000e-07 2.00000e-01 2.84000e-07
3.00000e-01 4.10000e-07 5.00000e-01 6.37000e-07
7.00000e-01 8.35000e-07 9.00000e-01 1.00000e-06
1.00000e+00 1.08000e-06 1.20000e+00 1.19000e-06
1.50000e+00 1.31000e-06 2.00000e+00 1.47000e-06
3.00000e+00 1.65000e-06 4.00000e+00 1.74000e-06
5.00000e+00 1.78000e-06 6.00000e+00 1.79000e-06
7.00000e+00 1.80000e-06 8.00000e+00 1.80000e-06
9.00000e+00 1.80000e-06 1.00000e+01 1.80000e-06
1.20000e+01 1.80000e-06 1.40000e+01 1.78000e-06
1.50000e+01 1.77000e-06 1.60000e+01 1.76000e-06
1.80000e+01 1.74000e-06 2.00000e+01 1.72000e-06
2.10000e+01 1.71000e-06 3.00000e+01 1.63000e-06
5.00000e+01 1.56000e-06 7.50000e+01 1.58000e-06
1.00000e+02 1.60000e-06 1.30000e+02 1.61000e-06
1.50000e+02 1.61000e-06 1.80000e+02 1.61000e-06
2.00000e+02 1.61000e-06 3.00000e+02 1.67000e-06
4.00000e+02 1.79000e-06 5.00000e+02 1.92000e-06
6.00000e+02 2.05000e-06 7.00000e+02 2.16000e-06
8.00000e+02 2.24000e-06 9.00000e+02 2.30000e-06
1.00000e+03 2.35000e-06 2.00000e+03 2.76000e-06
5.00000e+03 3.64000e-06 1.00000e+04 4.75000e-06
endmd

39
tallies/shutter_sample_f5.i
View File

@@ -1,39 +0,0 @@
c f5 detector tally within the beam just outside the selene 1/2 separator wall
fc55 neutron flux behind instrument shutter
f55:n 1255.9 -1756.9 -12.6 0.0
e55 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
1e-0 1.58 2.51 3.98 6.31
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
1e+3 1.58e+3 2.51e+3
fq55 e f
c reduce probability for contribution inside of bunker
pd55 0.0001 265r 1 0.0001 454r 0.1 1r 1 2r 0.1 10r 1 76r
fc255 neutron flux at detector position
f255:n 2280 -3180 -30 0.0
e255 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
1e-0 1.58 2.51 3.98 6.31
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
1e+3 1.58e+3 2.51e+3
fq255 e f
c reduce probability for contribution inside of bunker
pd255 0.0001 265r 1 0.0001 454r 0.1 1r 1 2r 0.1 10r 1 76r