From ffe8db531252f8bca347a24c9d9c08acd7123b3a Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Tue, 16 Jul 2019 16:26:40 +0200 Subject: [PATCH] Test performance with full set of tallies --- e2-estia.i | 21 ++++++------- tallies/detector_f5.i | 2 +- tallies/selene1_f5.i | 70 +++++++++++++++++++++++++++++++++++++++++++ 3 files changed, 82 insertions(+), 11 deletions(-) create mode 100644 tallies/selene1_f5.i diff --git a/e2-estia.i b/e2-estia.i index e599772..bc2eff1 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -94,7 +94,7 @@ c ------------------------------------------------------- c --------------- WEIGHT CARDS -------------------------- c ------------------------------------------------------- c weight window generator, can be run with previous weight window files -WWGE:n 1e-3 1e-1 1 10 100 3e3 +c WWGE:n 1e-3 1e-1 1 10 100 3e3 c mesh based weight windows, can be reused after change of cell structure c WWG 15 0 0 c MESH GEOM=CYL ORIGIN=0 0 -500 REF=133.7 -176.7 13.5 AXS=0 0 1 @@ -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 255 61003 0 +c WWG 255 61003 0 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ c ----------------------------------------------------------- c ------------------- TALLY CARDS --------------------------- @@ -121,22 +121,23 @@ c read file=tallies/full_mesh_np.i echo c read file=tallies/NBPI_f4.i echo c read file=tallies/BBG_f4.i echo c read file=tallies/BBG_f5.i echo -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 +read file=tallies/inbunker_f4.i echo +read file=tallies/chopper_pit_doses.i echo +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/pre_selene1_mesh.i +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 -c read file=tallies/shutter_wall_f4.i echo +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 c read file=tallies/experimental_cave_f5.i echo c read file=tallies/cave_mesh.i echo +read file=tallies/selene1_f5.i echo c read file=tallies/selene1_f6.i echo c c =============================== @@ -168,7 +169,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 f255 0.05 +stop nps 1e9 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 @@ -207,7 +208,7 @@ c read file=sources/mcstas_source_mf.i echo c read file=sources/mcstas_gammas.i echo c read file=sources/mcstas_gammas_2.i echo c read file=sources/mcstas_source_shutter.i echo -prdmp 1e8 -30 1 2 1e6 +prdmp 1e8 -30 1 2 1e7 print 30 32 40 c 120 $O importance function quality c 128 $O universe map diff --git a/tallies/detector_f5.i b/tallies/detector_f5.i index fc3308d..7b90cdd 100644 --- a/tallies/detector_f5.i +++ b/tallies/detector_f5.i @@ -1,6 +1,6 @@ c f5 detector tally to estimate background at the detector fc255 neutron flux at detector position -f255:n 2280 -3180 -30 0.0 +f255:n 2280 -3180 -30 0.0 ND 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 diff --git a/tallies/selene1_f5.i b/tallies/selene1_f5.i new file mode 100644 index 0000000..4b6b0c6 --- /dev/null +++ b/tallies/selene1_f5.i @@ -0,0 +1,70 @@ +c tallies for Selene 1 component fluxes (without extra cells) +fc905 metrology cart center and top +f905:n 1209 -1665 -10 0 1209 -1665 24 0 ND +e905 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 +fq905 e f +fc915 robot center +f915:n 1162 -1698 -10 0 ND +e915 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 +fq915 e f +fc925 c-axis before and after Cu blocks +f925:n 963 -1341 -7 0 1038 -1444 -7 0 ND +e925 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 +fq925 e f +pd905 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 +pd915 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 +pd925 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