diff --git a/e2-estia.i b/e2-estia.i index 2bfa138..a32290e 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -161,12 +161,12 @@ c read file=tallies/BBG_f5.i echo read file=tallies/inbunker_f4.i echo c read file=tallies/CPC1_f5.i echo c read file=tallies/CPC2_f5.i echo -c read file=tallies/chopper_pit_doses.i echo +read file=tallies/chopper_pit_doses.i echo c read file=tallies/bunker_feedthrough_f1.i echo read file=tallies/bunker_wall_f4.i echo c read file=tallies/bunker_wall_f5.i echo c read file=tallies/pre_selene1_mesh.i -c read file=tallies/e1_e2_f4.i echo +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 @@ -210,6 +210,9 @@ c dbcn 375642321 j j j j j j j j j j j 15291711 c stop calculation either at XeY histories, 6.5 days or tally 15 error <0.001 c stop nps 39242646 CTME 4e5 $f31 0.01 stop nps 1e9 CTME 4e5 $f31 0.01 +c write particle tracks +ssw 206505 pty=n $ record all neutron tracks crossing outward +c random number generator options rand gen=2 seed=29976355590901 $ default 19073486328125 c histp mode n p h / z $ neutron, photon, proton, pi+, pi0 @@ -249,6 +252,7 @@ c read file=sources/proton.i echo c * Source Term from Valentina for 2m adapted to size c read file=sources/e2_2m.i echo read file=sources/e2_2m_mod.i echo +c read file=sources/rssa.i echo c * Estimated source within bunker wall to optimize collimator efficiency c read file=sources/bunker_wall.i echo $ approximate, not for official results c * McStas source terms for gamma doeses diff --git a/geometry/ess_outside_cells.i b/geometry/ess_outside_cells.i index b413b92..fc7c3f4 100644 --- a/geometry/ess_outside_cells.i +++ b/geometry/ess_outside_cells.i @@ -644,8 +644,8 @@ c Main ESS plant geometry without the instrument 100264 200 -3.9 -100067 100066 -100111 100110 -100072 100068 100003 100265 200 -3.9 -100067 100066 -100112 100111 -100072 100068 100003 c ESS Magnadense Heavy Concrete -100266 200 -3.9 -100067 100066 -100113 100112 -100072 100068 100003 -100268 200 -3.9 -100067 100066 -100101 100114 -100072 100068 100003 +100266 200 -3.9 -100067 100066 -100113 100112 -206505 100068 100003 +100268 200 -3.9 -100067 100066 -100101 100114 -206505 100068 100003 100269 200 -3.9 (-100121 : -100065 : -100122 : 100120 : 100101 : 100123) (-100003 : -100069 : -100002 : -100070 : -100064 : 100071 : 100072) 100064 -100115 -100107 100103 100067 -100072 100003 (-100065 : diff --git a/geometry/estia_guides.i b/geometry/estia_guides.i index 4dd4675..c047b19 100644 --- a/geometry/estia_guides.i +++ b/geometry/estia_guides.i @@ -8,8 +8,8 @@ c -----Bunker Wall---- #(-206401 -206402 -206403 -206404) 205004 54 0.0833854 -100067 100066 -100114 100113 -206505 206405 #(206405 -206505 -206501 -206502 -206503 -206504) -205005 200 -3.9 -100067 100066 -100114 100113 -100072 206505 #(-206601 - -206602 -206603 -206604) +205005 200 -3.9 -100067 100066 -100101 100112 -100072 206505 #(-206601 + -206602 -206603 -206604 -100114) c -----Wall Cut-Out--- 206001 4 -1.2E-3 (((100068 -206205 -206201 -206202 -206203 -206204) :(206205 -206305 -206301 -206302 -206303 -206304)) 100003 226106 diff --git a/sbatch.tsh b/sbatch.tsh index 31a21bb..84e2fdc 100644 --- a/sbatch.tsh +++ b/sbatch.tsh @@ -24,6 +24,7 @@ mpirun -np $SLURM_NPROCS /home/ll_glavic/shielding/mcnp6/MCNP_SM/bin/mcnp6.mpi \ runtpe=runs/e2-estia_$CURRENT_VERSION.run \ mctal=results/e2-estia_$CURRENT_VERSION.mct \ mdata=results/e2-estia_$CURRENT_VERSION.mdt \ + wssa=results/e2-estia_$CURRENT_VERSION.ssa \ wwout=weight_windows/e2-estia_$CURRENT_VERSION.wwg \ wwone=weight_windows/e2-estia_$CURRENT_VERSION.wwo \ wwinp=weight_windows/current.wwg diff --git a/sources/rssa.i b/sources/rssa.i new file mode 100644 index 0000000..7b7b6c6 --- /dev/null +++ b/sources/rssa.i @@ -0,0 +1,8 @@ +c source definition based on surface source generated from 2m source term +cut:n j 0 +c cut all photons that can't produce neutrons +cut:p 1e+08 1.667 +c but adapted to model geometry and actual size of Estia beam extraction +c make sure to supply the correct WWINP file +WWP:n j j 4 j -1 j 1.0 1j 2.0E+13 1j +ssr diff --git a/tallies/bunker_wall_f4.i b/tallies/bunker_wall_f4.i index 513d340..8e9a436 100644 --- a/tallies/bunker_wall_f4.i +++ b/tallies/bunker_wall_f4.i @@ -17,31 +17,31 @@ e34 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 1e+3 1.58e+3 2.51e+3 fq34 e f c convert to dose rate Sv/h -de34 log 1.00000e-18 1.00000e-09 1.00000e-08 2.50000e-08 1.00000e-07 - 2.00000e-07 5.00000e-07 1.00000e-06 2.00000e-06 5.00000e-06 - 1.00000e-05 2.00000e-05 5.00000e-05 1.00000e-04 2.00000e-04 - 5.00000e-04 1.00000e-03 2.00000e-03 5.00000e-03 1.00000e-02 - 2.00000e-02 3.00000e-02 5.00000e-02 7.00000e-02 1.00000e-01 - 1.50000e-01 2.00000e-01 3.00000e-01 5.00000e-01 7.00000e-01 - 9.00000e-01 1.00000e+00 1.20000e+00 1.50000e+00 2.00000e+00 - 3.00000e+00 4.00000e+00 5.00000e+00 6.00000e+00 7.00000e+00 - 8.00000e+00 9.00000e+00 1.00000e+01 1.20000e+01 1.40000e+01 - 1.50000e+01 1.60000e+01 1.80000e+01 2.00000e+01 2.10000e+01 - 3.00000e+01 5.00000e+01 7.50000e+01 1.00000e+02 1.30000e+02 - 1.50000e+02 1.80000e+02 2.00000e+02 3.00000e+02 4.00000e+02 - 5.00000e+02 6.00000e+02 7.00000e+02 8.00000e+02 9.00000e+02 - 1.00000e+03 2.00000e+03 5.00000e+03 1.00000e+04 -df34 log 1.00000e-08 1.11000e-08 1.28000e-08 1.44000e-08 1.87000e-08 - 2.11000e-08 2.37000e-08 2.53000e-08 2.66000e-08 2.78000e-08 - 2.82000e-08 2.82000e-08 2.82000e-08 2.80000e-08 2.78000e-08 - 2.71000e-08 2.71000e-08 2.74000e-08 2.87000e-08 3.28000e-08 - 4.39000e-08 5.65000e-08 8.28000e-08 1.10000e-07 1.51000e-07 - 2.18000e-07 2.84000e-07 4.10000e-07 6.37000e-07 8.35000e-07 - 1.00000e-06 1.08000e-06 1.19000e-06 1.31000e-06 1.47000e-06 - 1.65000e-06 1.74000e-06 1.78000e-06 1.79000e-06 1.80000e-06 - 1.80000e-06 1.80000e-06 1.80000e-06 1.80000e-06 1.78000e-06 - 1.77000e-06 1.76000e-06 1.74000e-06 1.72000e-06 1.71000e-06 - 1.63000e-06 1.56000e-06 1.58000e-06 1.60000e-06 1.61000e-06 - 1.61000e-06 1.61000e-06 1.61000e-06 1.67000e-06 1.79000e-06 - 1.92000e-06 2.05000e-06 2.16000e-06 2.24000e-06 2.30000e-06 - 2.35000e-06 2.76000e-06 3.64000e-06 4.75000e-06 +c de34 log 1.00000e-18 1.00000e-09 1.00000e-08 2.50000e-08 1.00000e-07 +c 2.00000e-07 5.00000e-07 1.00000e-06 2.00000e-06 5.00000e-06 +c 1.00000e-05 2.00000e-05 5.00000e-05 1.00000e-04 2.00000e-04 +c 5.00000e-04 1.00000e-03 2.00000e-03 5.00000e-03 1.00000e-02 +c 2.00000e-02 3.00000e-02 5.00000e-02 7.00000e-02 1.00000e-01 +c 1.50000e-01 2.00000e-01 3.00000e-01 5.00000e-01 7.00000e-01 +c 9.00000e-01 1.00000e+00 1.20000e+00 1.50000e+00 2.00000e+00 +c 3.00000e+00 4.00000e+00 5.00000e+00 6.00000e+00 7.00000e+00 +c 8.00000e+00 9.00000e+00 1.00000e+01 1.20000e+01 1.40000e+01 +c 1.50000e+01 1.60000e+01 1.80000e+01 2.00000e+01 2.10000e+01 +c 3.00000e+01 5.00000e+01 7.50000e+01 1.00000e+02 1.30000e+02 +c 1.50000e+02 1.80000e+02 2.00000e+02 3.00000e+02 4.00000e+02 +c 5.00000e+02 6.00000e+02 7.00000e+02 8.00000e+02 9.00000e+02 +c 1.00000e+03 2.00000e+03 5.00000e+03 1.00000e+04 +c df34 log 1.00000e-08 1.11000e-08 1.28000e-08 1.44000e-08 1.87000e-08 +c 2.11000e-08 2.37000e-08 2.53000e-08 2.66000e-08 2.78000e-08 +c 2.82000e-08 2.82000e-08 2.82000e-08 2.80000e-08 2.78000e-08 +c 2.71000e-08 2.71000e-08 2.74000e-08 2.87000e-08 3.28000e-08 +c 4.39000e-08 5.65000e-08 8.28000e-08 1.10000e-07 1.51000e-07 +c 2.18000e-07 2.84000e-07 4.10000e-07 6.37000e-07 8.35000e-07 +c 1.00000e-06 1.08000e-06 1.19000e-06 1.31000e-06 1.47000e-06 +c 1.65000e-06 1.74000e-06 1.78000e-06 1.79000e-06 1.80000e-06 +c 1.80000e-06 1.80000e-06 1.80000e-06 1.80000e-06 1.78000e-06 +c 1.77000e-06 1.76000e-06 1.74000e-06 1.72000e-06 1.71000e-06 +c 1.63000e-06 1.56000e-06 1.58000e-06 1.60000e-06 1.61000e-06 +c 1.61000e-06 1.61000e-06 1.61000e-06 1.67000e-06 1.79000e-06 +c 1.92000e-06 2.05000e-06 2.16000e-06 2.24000e-06 2.30000e-06 +c 2.35000e-06 2.76000e-06 3.64000e-06 4.75000e-06 diff --git a/tallies/inbunker_f4.i b/tallies/inbunker_f4.i index ef7f3c6..e767e09 100644 --- a/tallies/inbunker_f4.i +++ b/tallies/inbunker_f4.i @@ -34,31 +34,31 @@ e24 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 1e+3 1.58e+3 2.51e+3 fq24 e f c convert to dose rate Sv/h -de24 log 1.00000e-18 1.00000e-09 1.00000e-08 2.50000e-08 1.00000e-07 - 2.00000e-07 5.00000e-07 1.00000e-06 2.00000e-06 5.00000e-06 - 1.00000e-05 2.00000e-05 5.00000e-05 1.00000e-04 2.00000e-04 - 5.00000e-04 1.00000e-03 2.00000e-03 5.00000e-03 1.00000e-02 - 2.00000e-02 3.00000e-02 5.00000e-02 7.00000e-02 1.00000e-01 - 1.50000e-01 2.00000e-01 3.00000e-01 5.00000e-01 7.00000e-01 - 9.00000e-01 1.00000e+00 1.20000e+00 1.50000e+00 2.00000e+00 - 3.00000e+00 4.00000e+00 5.00000e+00 6.00000e+00 7.00000e+00 - 8.00000e+00 9.00000e+00 1.00000e+01 1.20000e+01 1.40000e+01 - 1.50000e+01 1.60000e+01 1.80000e+01 2.00000e+01 2.10000e+01 - 3.00000e+01 5.00000e+01 7.50000e+01 1.00000e+02 1.30000e+02 - 1.50000e+02 1.80000e+02 2.00000e+02 3.00000e+02 4.00000e+02 - 5.00000e+02 6.00000e+02 7.00000e+02 8.00000e+02 9.00000e+02 - 1.00000e+03 2.00000e+03 5.00000e+03 1.00000e+04 -df24 log 1.00000e-08 1.11000e-08 1.28000e-08 1.44000e-08 1.87000e-08 - 2.11000e-08 2.37000e-08 2.53000e-08 2.66000e-08 2.78000e-08 - 2.82000e-08 2.82000e-08 2.82000e-08 2.80000e-08 2.78000e-08 - 2.71000e-08 2.71000e-08 2.74000e-08 2.87000e-08 3.28000e-08 - 4.39000e-08 5.65000e-08 8.28000e-08 1.10000e-07 1.51000e-07 - 2.18000e-07 2.84000e-07 4.10000e-07 6.37000e-07 8.35000e-07 - 1.00000e-06 1.08000e-06 1.19000e-06 1.31000e-06 1.47000e-06 - 1.65000e-06 1.74000e-06 1.78000e-06 1.79000e-06 1.80000e-06 - 1.80000e-06 1.80000e-06 1.80000e-06 1.80000e-06 1.78000e-06 - 1.77000e-06 1.76000e-06 1.74000e-06 1.72000e-06 1.71000e-06 - 1.63000e-06 1.56000e-06 1.58000e-06 1.60000e-06 1.61000e-06 - 1.61000e-06 1.61000e-06 1.61000e-06 1.67000e-06 1.79000e-06 - 1.92000e-06 2.05000e-06 2.16000e-06 2.24000e-06 2.30000e-06 - 2.35000e-06 2.76000e-06 3.64000e-06 4.75000e-06 +c de24 log 1.00000e-18 1.00000e-09 1.00000e-08 2.50000e-08 1.00000e-07 +c 2.00000e-07 5.00000e-07 1.00000e-06 2.00000e-06 5.00000e-06 +c 1.00000e-05 2.00000e-05 5.00000e-05 1.00000e-04 2.00000e-04 +c 5.00000e-04 1.00000e-03 2.00000e-03 5.00000e-03 1.00000e-02 +c 2.00000e-02 3.00000e-02 5.00000e-02 7.00000e-02 1.00000e-01 +c 1.50000e-01 2.00000e-01 3.00000e-01 5.00000e-01 7.00000e-01 +c 9.00000e-01 1.00000e+00 1.20000e+00 1.50000e+00 2.00000e+00 +c 3.00000e+00 4.00000e+00 5.00000e+00 6.00000e+00 7.00000e+00 +c 8.00000e+00 9.00000e+00 1.00000e+01 1.20000e+01 1.40000e+01 +c 1.50000e+01 1.60000e+01 1.80000e+01 2.00000e+01 2.10000e+01 +c 3.00000e+01 5.00000e+01 7.50000e+01 1.00000e+02 1.30000e+02 +c 1.50000e+02 1.80000e+02 2.00000e+02 3.00000e+02 4.00000e+02 +c 5.00000e+02 6.00000e+02 7.00000e+02 8.00000e+02 9.00000e+02 +c 1.00000e+03 2.00000e+03 5.00000e+03 1.00000e+04 +c df24 log 1.00000e-08 1.11000e-08 1.28000e-08 1.44000e-08 1.87000e-08 +c 2.11000e-08 2.37000e-08 2.53000e-08 2.66000e-08 2.78000e-08 +c 2.82000e-08 2.82000e-08 2.82000e-08 2.80000e-08 2.78000e-08 +c 2.71000e-08 2.71000e-08 2.74000e-08 2.87000e-08 3.28000e-08 +c 4.39000e-08 5.65000e-08 8.28000e-08 1.10000e-07 1.51000e-07 +c 2.18000e-07 2.84000e-07 4.10000e-07 6.37000e-07 8.35000e-07 +c 1.00000e-06 1.08000e-06 1.19000e-06 1.31000e-06 1.47000e-06 +c 1.65000e-06 1.74000e-06 1.78000e-06 1.79000e-06 1.80000e-06 +c 1.80000e-06 1.80000e-06 1.80000e-06 1.80000e-06 1.78000e-06 +c 1.77000e-06 1.76000e-06 1.74000e-06 1.72000e-06 1.71000e-06 +c 1.63000e-06 1.56000e-06 1.58000e-06 1.60000e-06 1.61000e-06 +c 1.61000e-06 1.61000e-06 1.61000e-06 1.67000e-06 1.79000e-06 +c 1.92000e-06 2.05000e-06 2.16000e-06 2.24000e-06 2.30000e-06 +c 2.35000e-06 2.76000e-06 3.64000e-06 4.75000e-06