From 75bbb7cda094a0067a321dff90b8c72739fbd8d7 Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Fri, 28 Jun 2019 15:50:19 +0200 Subject: [PATCH] Generate WW for detector position and test fine meshes behind collimation and bunker --- e2-estia.i | 9 ++++--- tallies/detector_f5.i | 25 +++++++++++++++++++ tallies/pre_selene1_mesh.i | 49 +++++++++++++++++++++++++++++++++++++ tallies/shutter_sample_f5.i | 39 ----------------------------- 4 files changed, 79 insertions(+), 43 deletions(-) create mode 100644 tallies/detector_f5.i create mode 100644 tallies/pre_selene1_mesh.i delete mode 100644 tallies/shutter_sample_f5.i diff --git a/e2-estia.i b/e2-estia.i index edc890e..e599772 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -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 diff --git a/tallies/detector_f5.i b/tallies/detector_f5.i new file mode 100644 index 0000000..fc3308d --- /dev/null +++ b/tallies/detector_f5.i @@ -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 diff --git a/tallies/pre_selene1_mesh.i b/tallies/pre_selene1_mesh.i new file mode 100644 index 0000000..3d99256 --- /dev/null +++ b/tallies/pre_selene1_mesh.i @@ -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 diff --git a/tallies/shutter_sample_f5.i b/tallies/shutter_sample_f5.i deleted file mode 100644 index f2d4366..0000000 --- a/tallies/shutter_sample_f5.i +++ /dev/null @@ -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