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

Test generation of surface source in combination with f4 tallies

Browse Source
This commit is contained in:
glavic_a
2020-04-29 10:35:57 +02:00
parent 2272904a66
commit fab54aff69
7 changed files with 75 additions and 62 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
e2-estia.i
View File

@@ -161,12 +161,12 @@ c read file=tallies/BBG_f5.i echo
read file=tallies/inbunker_f4.i echo read file=tallies/inbunker_f4.i echo
c read file=tallies/CPC1_f5.i echo c read file=tallies/CPC1_f5.i echo
c read file=tallies/CPC2_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 c read file=tallies/bunker_feedthrough_f1.i echo
read file=tallies/bunker_wall_f4.i echo read file=tallies/bunker_wall_f4.i echo
c read file=tallies/bunker_wall_f5.i echo c read file=tallies/bunker_wall_f5.i echo
c read file=tallies/pre_selene1_mesh.i 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/e1_e2_f5.i echo
c read file=tallies/shutter_f4.i echo c read file=tallies/shutter_f4.i echo
c read file=tallies/shutter_f5.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 calculation either at XeY histories, 6.5 days or tally 15 error <0.001
c stop nps 39242646 CTME 4e5 $f31 0.01 c stop nps 39242646 CTME 4e5 $f31 0.01
stop nps 1e9 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 rand gen=2 seed=29976355590901 $ default 19073486328125
c histp c histp
mode n p h / z $ neutron, photon, proton, pi+, pi0 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 * Source Term from Valentina for 2m adapted to size
c read file=sources/e2_2m.i echo c read file=sources/e2_2m.i echo
read file=sources/e2_2m_mod.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 * Estimated source within bunker wall to optimize collimator efficiency
c read file=sources/bunker_wall.i echo $ approximate, not for official results c read file=sources/bunker_wall.i echo $ approximate, not for official results
c * McStas source terms for gamma doeses c * McStas source terms for gamma doeses

4
geometry/ess_outside_cells.i
View File

@@ -644,8 +644,8 @@ c Main ESS plant geometry without the instrument
100264 200 -3.9 -100067 100066 -100111 100110 -100072 100068 100003 100264 200 -3.9 -100067 100066 -100111 100110 -100072 100068 100003
100265 200 -3.9 -100067 100066 -100112 100111 -100072 100068 100003 100265 200 -3.9 -100067 100066 -100112 100111 -100072 100068 100003
c ESS Magnadense Heavy Concrete c ESS Magnadense Heavy Concrete
100266 200 -3.9 -100067 100066 -100113 100112 -100072 100068 100003 100266 200 -3.9 -100067 100066 -100113 100112 -206505 100068 100003
100268 200 -3.9 -100067 100066 -100101 100114 -100072 100068 100003 100268 200 -3.9 -100067 100066 -100101 100114 -206505 100068 100003
100269 200 -3.9 (-100121 : -100065 : -100122 : 100120 : 100101 : 100123) 100269 200 -3.9 (-100121 : -100065 : -100122 : 100120 : 100101 : 100123)
(-100003 : -100069 : -100002 : -100070 : -100064 : 100071 : 100072) (-100003 : -100069 : -100002 : -100070 : -100064 : 100071 : 100072)
100064 -100115 -100107 100103 100067 -100072 100003 (-100065 : 100064 -100115 -100107 100103 100067 -100072 100003 (-100065 :

4
geometry/estia_guides.i
View File

@@ -8,8 +8,8 @@ c -----Bunker Wall----
#(-206401 -206402 -206403 -206404) #(-206401 -206402 -206403 -206404)
205004 54 0.0833854 -100067 100066 -100114 100113 -206505 206405 205004 54 0.0833854 -100067 100066 -100114 100113 -206505 206405
#(206405 -206505 -206501 -206502 -206503 -206504) #(206405 -206505 -206501 -206502 -206503 -206504)
205005 200 -3.9 -100067 100066 -100114 100113 -100072 206505 #(-206601 205005 200 -3.9 -100067 100066 -100101 100112 -100072 206505 #(-206601
-206602 -206603 -206604) -206602 -206603 -206604 -100114)
c -----Wall Cut-Out--- c -----Wall Cut-Out---
206001 4 -1.2E-3 (((100068 -206205 -206201 -206202 -206203 -206204) 206001 4 -1.2E-3 (((100068 -206205 -206201 -206202 -206203 -206204)
:(206205 -206305 -206301 -206302 -206303 -206304)) 100003 226106 :(206205 -206305 -206301 -206302 -206303 -206304)) 100003 226106

1
sbatch.tsh
View File

@@ -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 \ runtpe=runs/e2-estia_$CURRENT_VERSION.run \
mctal=results/e2-estia_$CURRENT_VERSION.mct \ mctal=results/e2-estia_$CURRENT_VERSION.mct \
mdata=results/e2-estia_$CURRENT_VERSION.mdt \ mdata=results/e2-estia_$CURRENT_VERSION.mdt \
wssa=results/e2-estia_$CURRENT_VERSION.ssa \
wwout=weight_windows/e2-estia_$CURRENT_VERSION.wwg \ wwout=weight_windows/e2-estia_$CURRENT_VERSION.wwg \
wwone=weight_windows/e2-estia_$CURRENT_VERSION.wwo \ wwone=weight_windows/e2-estia_$CURRENT_VERSION.wwo \
wwinp=weight_windows/current.wwg wwinp=weight_windows/current.wwg

8
sources/rssa.i Normal file
View File

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

56
tallies/bunker_wall_f4.i
View File

@@ -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 1e+3 1.58e+3 2.51e+3
fq34 e f fq34 e f
c convert to dose rate Sv/h c convert to dose rate Sv/h
de34 log 1.00000e-18 1.00000e-09 1.00000e-08 2.50000e-08 1.00000e-07 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 2.35000e-06 2.76000e-06 3.64000e-06 4.75000e-06

56
tallies/inbunker_f4.i
View File

@@ -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 1e+3 1.58e+3 2.51e+3
fq24 e f fq24 e f
c convert to dose rate Sv/h c convert to dose rate Sv/h
de24 log 1.00000e-18 1.00000e-09 1.00000e-08 2.50000e-08 1.00000e-07 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 c 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 2.35000e-06 2.76000e-06 3.64000e-06 4.75000e-06