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

Remove all shielding but SKADI around Selene 1, add vacuum tank and simulate gamma dose

Browse Source
This commit is contained in:
glavic_a
2018-11-14 12:28:22 +01:00
parent 6fbf2b5009
commit f43517c8d4
5 changed files with 70 additions and 43 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
e2-estia.i
View File

@@ -31,10 +31,10 @@ c 700 s 117.0 -155.00 13.5 3.0
703 1 s 1082 -25.0 0 5.0
704 1 s 1082 -0.25 0 5.0
705 1 s 1165 0.7 0 5.0
706 1 s 1328 5 0 12.0
706 1 s 1320 5 0 12.0
707 1 s 2190 3.5 0 12.0
708 s 730 -1380 0 30.0
709 s 1020 -1170 0 30.0
709 s 1030 -1155 0 30.0
710 s 1328 -1756 -13 12.0
c Estia geometry
read file=geometry/estia_surfaces.i echo
@@ -94,7 +94,7 @@ c -----------------------------------------------------------
c ------------------- TALLY CARDS ---------------------------
c
c read file=tallies/bunker_mesh.i echo
c read file=tallies/bunker_mesh_p.i echo
read file=tallies/bunker_mesh_p.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
@@ -122,7 +122,7 @@ ext:n 0 0 0.5V2 0 0.5V2 0 8r 0.5V2 7r 0 0.5V2 0 0.5V2 8r 0 0.5V2 0 0.5V2 0
0 0.5V2 0 0.5V2 0 0 0.5V2 0 0.5V2 0 0.5V2 0 0.5V2 0 22r
0.5V2 10r 0 4r 0.5V2 3r 0 0.5V2 13r 0 48r
-0.5V1 2r 0 4r -0.1V1 9r 0 0.1V3 3r 0 11r 0.1V3 2r
0 0.1V3 0 0.1V3 0 0.1V3 5r 0 30r
0 0.1V3 0 0 0.1V3 0 0.1V3 5r 0 30r
c exponential transform reference vectors: origin, selene 1, inst. shutter
vect v1 0 0 0 v2 1000 -1394 -6.5 v3 1282 -1793 -13
c
@@ -131,15 +131,15 @@ c -------------------------------------------------------
c --------------- PHYSICS CARDS --------------------------
c -------------------------------------------------------
c dbcn 375642321 j j j j j j j j j j j 15291711
nps 2e9
nps 1e9
c histp
mode n p n h / d t s a
c volume calculated for tally cells (spheres)
c and set for chopper disk/spindle and stages x1/x2/z1/z2/rot
vol 1 507r 10j 1 107r 3848.5 10104 1 13r 479.5 479.5 913 913 1734
1 18r 3.34e6 1 2r 3.5e6 1 3.5e6 1 6r 2.05e5 1 12r $ wall segments
imp:h,/,d,t,s,a 0 0 1 691r
imp:n,p 0 0 1 691r
1 18r 3.34e6 1 2r 3.5e6 1 3.5e6 1 6r 2.05e5 1 13r $ wall segments
imp:h,/,d,t,s,a 0 0 1 692r
imp:n,p 0 1 1 692r
cut:h,/,d,t,s,a 1e+08 1e-3
c cut:n j 0.99E-6
cut:n j 1e-11
@@ -153,11 +153,11 @@ lea 1 4 1 0 1 0 0 1
c * Source for simulation from proton beam (needs different WW and ESplit)
c read file=sources/proton.i echo
c * Source Term from Valentina for 2m adapted to size
read file=sources/e2_2m.i echo
c read file=sources/e2_2m.i echo
c * McStas source terms for gamma doeses
c read file=sources/mcstas_source_vs.i echo
c read file=sources/mcstas_source_sample.i echo
c read file=sources/mcstas_gammas.i echo
read file=sources/mcstas_gammas.i echo
prdmp 5e7 -60 1 2 5e6
print 30 32 $ 120 128 150 160 161 162
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++

51
geometry/estia_guides.i
View File

@@ -13,19 +13,20 @@ c -----Wall Cut-Out---
6003 4 -1.2E-3 (6103 -72 -6601 -6602 -6603 -6604 -6605) 3
#(7103 8405 -8921)
#(7102 -7103 7205 7206 7207 7208)
#(-8405 8401 8402 8403 8404)
#(8421 8422 8423 8424 8425 8426)
c #(-8405 8401 8402 8403 8404)
c Vacuum before and inside Selene Guide 1
6101 0 ((7103 8405 -8922 3)
:(-8405 8406 -8411 8412 8413 8402 3))
#706 #8311 #8411 #8412 #8413
6102 5 0.0582256 (7103 8405 8922 -8921 3)
6101 0 (7103 8435 -8922 3) #706
6102 5 0.0582256 (7103 -8425 8922 -8921 3)
6103 0 (-8431 -8432 -8433 -8434 -8435 -8436 (-8441:-8442) (-8443:-8444))
#8301 #8311 #8411 #8412 #8413
c Vaccum after Selene 1, until and inside Selen 2
6111 0 ((-8406 -8100 -8932 3)
6111 0 ((8436 -8100 -8932 3)
:(8100 -8101 -8501 -8502 -8503 -8504 3)
:(8101 -8121 -8932 3)
:(8121 -8122 8601 -8611 -8612 -8613 3))
#707 #8511
6112 5 0.0582256 (-8406 -8121 8932 -8931 3 #(8100 -8101))
6112 5 0.0582256 (-8426 -8121 8932 -8931 3 #(8100 -8101))
c -----Bunker Wall Collimator (BWC)-----
7001 0 (7100 -7103 -7301 -7302 -7303 -7304 3)
#(7101 -7401 7402)
@@ -51,7 +52,7 @@ c concrete floor
:(72 -8101 8209 -59 8201 8202)
:(8101 -8709 8209 -59 8301 8302)) 3
c ceiling outside
8002 49 0.074981 ( 72 -8101 8207 8203 8201 8202)
8002 4 -1.2E-3 ( 72 -8101 8207 8203 8201 8202)
c left outer wall (E01)
8003 4 -1.2E-3 ( 72 -8100 8201 -8205 59 -8207 3)
c Selene 1/2 separator
@@ -62,31 +63,37 @@ c Optical cave
#(8101 -8709 -8311 -8312 -8313 -8314)
c air
8008 4 -1.2E-3 (72 -8100 -8211 -8212 -8213 -8214 3)
#(-8406 -8931)
#(-8405 8406 8401 8402 8403 8404)
#(8421 8422 8423 8424 8425 8426)
#(-8426 -8931)
c #(-8405 8406 8401 8402 8403 8404)
c #8370 #8371 #8372
8009 4 -1.2E-3 (8101 -8709 -8311 -8312 -8313 -8314 3)
#(8121 -8122 8601 8602 8603 8604)
#(-8121 -8931)
#707 #710
c B4C concrete inner layers
8011 4 -1.2E-3 (8211 8215 -8213 -8214 72 -8100 3)
8012 16 -2.4 (8205 8202 8213 -8218 72 -8100 3)
8013 4 -1.2E-3 ((-8216 8111 -8100):(72 -8111 18202))
8212 -8213 -8214 3
8014 16 -2.4 (8202 -8206 8217 -8213 -8214 8111 -8100 3)
c Steel middle layers
8020 71 0.0843223 (-18202 8202 -8213 -8214 72 -8111 3)
8021 71 0.0843223 (-8215 8205 -8213 -8214 72 -8100 3)
8022 71 0.0843223 (8205 8202 8218 -8207 72 -8100 3)
8023 71 0.0843223 (8216 8202 (-8217:8206) -8213 -8214 8111 -8100 3)
c Selene 1 shielding
8011 4 -1.2E-3 (8211 8205 -8213 -8214 72 -8100 3)
8012 4 -1.2E-3 (8205 8202 8213 -8218 72 -8100 3)
8013 49 0.074981 72 8202 8212 8217 -8213 -8214 3
c 8014 4 -1.2E-3 (8202 -8206 8217 -8213 -8214 8111 -8100 3)
c
8020 49 0.074981 (-8217 -8100 8215 8202 -8213 -8214 3)
8021 4 -1.2E-3 (-8217 -8100 -8215 8212 -8213 -8214 3)
8022 4 -1.2E-3 (8205 8202 8218 -8207 72 -8100 3)
c 8023 4 -1.2E-3 (8216 8202 (-8217:8206) -8213 -8214 8111 -8100 3)
c
c Selene 1 carrier
8301 18 -2.4 -8405 8406 ((8401 8402 8403 8404)
#(-8411 8402 8412 8413))
c Selene 1 glass (change material!!)
c Selene 1 glass
8311 12 -2.51 -8405 8406 -8903 -8904 -8518
#(-8901 -8902) #(-8515 -8517) 3
c Selene 1 Vacuum tank (stainless steel)
8321 3 0.0876276 ( 8421 8422 8423 8424 8425 8426)
#(-8431 -8432 -8433 -8434 -8435 -8436)
#(8425 8435 -8922)
#(8426 8436 -8932)
8422 13 -11.34 -8435 -8436 -8432 -8433 ((8441 8442):(8443 8444))
c Reference sensor material for Selene 1 dose simulations
c 8370 71 0.0843223 -8370
c 8371 71 0.0843223 -8371

34
geometry/estia_surfaces.i
View File

@@ -471,7 +471,7 @@ c Sides 4
c Sides 5
6601 p -0.809017 -0.587785 0.000000 85.500
6602 p 0.809017 0.587785 0.000000 41.800
6603 pz 76.000
6603 pz 85.0
6604 p 0.000000 0.000000 -1.000000 81.000
6605 p -0.987688 0.156434 0.000000 -895.551
c Cu-collimator
@@ -521,7 +521,7 @@ c 8102 cz 3023.41
c Shielding wall geometry
c first section outside
8201 p 0.813101 0.582123 0.000000 -148.964
8202 p -0.776046 -0.630676 0.000000 1.917
8202 p -0.771625 -0.636078 0.000000 -21.127
18202 p -0.776046 -0.630676 0.000000 16.917
8203 p 0.000000 0.000000 -1.000000 -166.400
8204 pz -200.000
@@ -533,14 +533,14 @@ c 8208 p -0.776046 -0.630676 0.000000 31.917
8209 pz -300.000
c inside
8211 p -0.813101 -0.582123 0.000000 68.794
8212 p 0.813101 0.582123 0.000000 71.037
8212 p 0.790690 0.612217 0.000000 0.000
8213 pz 85.900
8214 p 0.000000 0.000000 -1.000000 200.000
c B4C concrete/steel interface
8215 p 0.813101 0.582123 0.000000 -88.964
8216 p 0.813101 0.582123 0.000000 91.037
8217 p 0.813101 0.582123 0.000000 101.037
8218 pz 105.900
8215 p 0.739506 0.673150 0.000000 -102.437
8216 pz 50.000
8217 p -0.612217 0.790690 0.000000 -1958.974
8218 pz 105.900
c 8219 p -0.776046 -0.630676 0.000000 111.917
c wall breakthrough
8221 p -0.813101 -0.582123 0.000000 32.866
@@ -578,6 +578,26 @@ c Inside
8411 p -0.813101 -0.582123 0.000000 14.866
8412 p -0.007112 0.009934 -0.999925 -38.883
8413 p 0.007112 -0.009934 0.999925 -11.117
c Selene 1 Vacuum
c Outside
8421 p 0.813101 0.582123 0.000000 -66.866
8422 p -0.813101 -0.582123 0.000000 -35.134
8423 p -0.007112 0.009934 -0.999925 -96.883
8424 p 0.007112 -0.009934 0.999925 -52.617
8425 p 0.582080 -0.813040 -0.012217 1345.888
8426 p -0.582080 0.813040 0.012217 -2097.888
c Inside
8431 p -0.813101 -0.582123 0.000000 64.166
8432 p 0.813101 0.582123 0.000000 31.934
8433 p 0.007112 -0.009934 0.999925 94.183
8434 p -0.007112 0.009934 -0.999925 50.917
8435 p -0.582080 0.813040 0.012217 -1347.388
8436 p 0.582080 -0.813040 -0.012217 2096.388
c Anti-gamma lead shield concept
8441 p 0.813101 0.582123 0.000000 28.934
8442 p 0.007112 -0.009934 0.999925 44.183
8443 p 0.007112 -0.009934 0.999925 91.183
8444 p 0.813101 0.582123 0.000000 -9.934
c Beam through wall
c Inside
8501 p -0.792291 -0.610143 -0.000426 79.405

4
tallies/bunker_mesh.i
View File

@@ -2,8 +2,8 @@ c large mesh tally covering the beamline from monolith to S1/2 separator
tmesh
rmesh1:n dose 1 trans 1
cora1 500 10i 1000 20i 1500 16i 2300
corb1 -305 5i -55 4i 45 5i 295
corc1 -200 1i -100 5i 100 2i 300
corb1 -305 5i -55 4i 45 9i 495
corc1 -200 1i -100 5i 100 4i 500
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

4
tallies/bunker_mesh_p.i
View File

@@ -2,8 +2,8 @@ c large mesh tally covering the beamline from monolith to S1/2 separator
tmesh
rmesh1:p dose 1 trans 1
cora1 500 10i 1000 20i 1500 16i 2300
corb1 -305 5i -55 4i 45 5i 295
corc1 -200 1i -100 5i 100 2i 300
corb1 -305 5i -55 4i 45 9i 495
corc1 -200 1i -100 5i 100 5i 600
c Official ESS dose conversion factors from gamma/cm²/s to Sv/h
c Is described in ESS-0019931
mshmf1 0.000000 1.00000e-04