Test performance with full set of tallies

This commit is contained in:
2019-07-16 16:26:40 +02:00
parent b1adeaf714
commit ffe8db5312
3 changed files with 82 additions and 11 deletions

View File

@@ -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

View File

@@ -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

70
tallies/selene1_f5.i Normal file
View File

@@ -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