From f950e5a639865555f2a5df1769c76fa6ed242a56 Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Wed, 1 May 2019 19:00:14 +0200 Subject: [PATCH] Separate cave into segments, air all sides and replace SKADI side wall with heavy concrete --- e2-estia.i | 10 +++++----- geometry/ess_cells.i | 1 + geometry/estia_cave.i | 32 +++++++++++++++++++++----------- geometry/estia_guides.i | 2 +- geometry/estia_surfaces.i | 9 ++++++--- tallies/full_mesh.i | 2 +- tallies/full_mesh_np.i | 4 ++-- tallies/full_mesh_p.i | 2 +- 8 files changed, 38 insertions(+), 24 deletions(-) diff --git a/e2-estia.i b/e2-estia.i index 6ecbeae..a236b82 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -137,7 +137,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 82r -0.5V1 49r 0 24r -0.1V1 4r 0 13r 0.1V3 9r $ chopper pit to bunker wall - 0 2r 0.1V3 1r 0 43r + 0 2r 0.1V3 1r 0 56r 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 @@ -146,15 +146,15 @@ c ------------------------------------------------------- c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 -nps 2e6 +nps 2e8 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 506r 11j 1 183r 3848.5 10104 1j 1 7r 479.5 479.5 913 913 1734 - 1 81r $ wall segments -imp:h,/,d,t,s,a 0 0 1 797r -imp:n,p 0 0 1 797r + 1 94r $ wall segments +imp:h,/,d,t,s,a 0 0 1 810r +imp:n,p 0 0 1 810r cut:h,/,d,t,s,a 1e+08 1e-3 phys:n 3000 0 phys:p 100 diff --git a/geometry/ess_cells.i b/geometry/ess_cells.i index a709dd6..d847db6 100644 --- a/geometry/ess_cells.i +++ b/geometry/ess_cells.i @@ -8,6 +8,7 @@ c Main ESS plant geometry without the instrument #(8709 8720 112 -104 8703 8701) #(-104 112 72 8710 8209 8203) #( 104 8708 8707 8710 8703 8701) + #( 104 610 8720 8721 8703 8203) #(68 -72 8209 -70 -104 -2) 3 118 0.122325 24 4 -19 4 0 159 -137 -21 153 -154 diff --git a/geometry/estia_cave.i b/geometry/estia_cave.i index 97d3132..51ffacd 100644 --- a/geometry/estia_cave.i +++ b/geometry/estia_cave.i @@ -1,16 +1,7 @@ c Experimental cave c First largest area: up-stream/down-stream/roof/floor/walls c then cut corner and add extended roof -8701 49 0.074981 ((8709 8710 8701 8703 8706 8707 8708 8709 8710 (-8302 : 8705)) - :(8709 8710 8702 8703 8711 8712) - :(8709 8701 8703 -8302 8805 8806)) - #((-8809 -8810 -8801 -8803 -8806 -8807 -8808 -8809 -8810) - :(-8809 -8810 -8802 -8803 -8811 -8812) - :(-8311 8312 8709 8809 -8313 -8803) - :(-8801 -8803 -8809 8815 8312 8806)) - #(( 72 -8101 8201 8202 8203 8204 3) - :(8101 -8709 8301 -8302 8303 8304 3)) -8702 4 -1.2E-3 ((-8809 -8810 -8801 -8803 -8806 -8807 -8808 -8809 -8810 +8701 4 -1.2E-3 ((-8809 -8810 -8801 -8803 -8806 -8807 -8808 -8809 -8810 (-8806 : -8815)) :(-8809 -8810 -8802 -8803 -8811 -8812) :(-8311 8312 8709 8809 -8313 -8803) @@ -18,11 +9,30 @@ c then cut corner and add extended roof :(-8312 -8302 8809 8709 -8701 -8313)) 3 #(8621 8622 8623 8624 8626) #6111 #8501 #8511 #9001 +8702 4 -1.2e-3 104 610 8721 8703 8203 -8708 +8703 4 -1.2e-3 -8707 8710 8721 8703 8203 8708 +8704 4 -1.2e-3 104 -8710 8720 8721 8703 8203 +c floor +8711 49 0.074981 8703 8803 8707 8708 8709 (8706:(8805 -8302)) 8710 +c walls +8712 49 0.074981 -8803 -8801 8709 -8705 -8312 -8302 +8713 49 0.074981 -8803 -8801 -8302 8806 8705 8805 +8714 200 -3.9 -8803 -8801 -8805 -8810 8806 8706 $ heavy concrete +8715 49 0.074981 -8803 -8801 8706 8707 8810 8710 +8716 49 0.074981 -8803 -8801 -8810 8707 8807 -8808 +8717 49 0.074981 -8803 -8801 8707 8808 8708 8709 +8718 49 0.074981 -8803 -8801 8311 8809 8709 -8808 +c ceiling +8719 49 0.074981 8801 8701 8710 8707 (8706:(8805 -8302)) 8708 8709 + #(8712 -8810 8711) #(-8311 8809 8312) +8720 49 0.074981 8801 8702 8709 8712 8711 8710 + #(-8810 -8812 -8811 -8802) #(-8311 8809 8312) +8721 49 0.074981 8313 8702 8709 8809 8711 8712 c sample position 9001 19 -8.65 -9002 9003 -9004 c 9001 4 -1.2E-3 -9001 c area around cave -9900 4 -1.2E-3 8709 8710 -104 (-8708:-8707) 8703 8701 +9900 4 -1.2E-3 8709 8710 -104 (-8708:-8707) 8703 8203 9901 4 -1.2E-3 8709 8710 112 (8302:-8805) 8703 8701 -8706 9902 4 -1.2E-3 8720 -8710 112 -104 8703 8701 c visualize Selene gamma sources (comment for real runs) diff --git a/geometry/estia_guides.i b/geometry/estia_guides.i index 37ca674..bc4476f 100644 --- a/geometry/estia_guides.i +++ b/geometry/estia_guides.i @@ -146,7 +146,7 @@ c Selene 2 Vacuum tank (stainless steel) c #(8626 8636 -8932) c Air outside the guide shielding 8601 4 -1.2E-3 ((72 -8709 -104 -8201 8203 8204) - :(-8701 8203 8710 8709 -104 -8302)) + :(-8701 8203 8710 8709 8708 8707 -8302)) #(101 -610 -67) #(-8305 -610) #(-8302 8702 ((8712 8711):(8809 -8311 -8313))) #708 diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index 9fcc9a8..3cddee8 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -681,9 +681,12 @@ c Inside 8803 p 0.000000 0.000000 -1.000000 166.960 8804 p 0.813101 0.582123 0.000000 102.134 8805 p -0.582123 0.813101 0.000000 -3414.795 -8806 p 0.736097 0.676876 0.000000 -236.513 -8807 p -0.813101 -0.582123 0.000000 637.866 -8808 p -0.981311 0.192430 0.000000 -1964.017 +c 8806 p 0.736097 0.676876 0.000000 -236.513 +c 8807 p -0.813101 -0.582123 0.000000 637.866 +c 8808 p -0.981311 0.192430 0.000000 -1964.017 +8806 p 0.736097 0.676876 0.000000 -246.513 +8807 p -0.813101 -0.582123 0.000000 627.866 +8808 p -0.981311 0.192430 0.000000 -1974.017 8809 p -0.582123 0.813101 0.000000 -3059.795 8810 p 0.582123 -0.813101 0.000000 4314.795 8811 p 0.813101 0.582123 0.000000 82.134 diff --git a/tallies/full_mesh.i b/tallies/full_mesh.i index 421cb5d..540a81b 100644 --- a/tallies/full_mesh.i +++ b/tallies/full_mesh.i @@ -2,7 +2,7 @@ c large and detailed mesh tally covering the beamline from monolith to full cave tmesh rmesh1:n dose 1 trans 1 cora1 500 10i 1000 40i 1500 60i 4500 - corb1 -715 13i -65 24i 55 13i 705 + corb1 -715 13i -65 24i 55 15i 805 corc1 -200 1i -100 40i 100 5i 600 c Official ESS dose conversion factors from n/cm²/s to Sv/h c Is described in ESS-0019931 diff --git a/tallies/full_mesh_np.i b/tallies/full_mesh_np.i index 4b1c11e..07445d4 100644 --- a/tallies/full_mesh_np.i +++ b/tallies/full_mesh_np.i @@ -2,7 +2,7 @@ c large and detailed mesh tally covering the beamline from monolith to full cave tmesh rmesh1:n dose 1 trans 1 cora1 500 10i 1000 40i 1500 60i 4500 - corb1 -715 13i -65 24i 55 13i 705 + corb1 -715 13i -65 24i 55 15i 805 corc1 -200 1i -100 40i 100 5i 600 c Official ESS dose conversion factors from n/cm²/s to Sv/h c Is described in ESS-0019931 @@ -44,7 +44,7 @@ c Is described in ESS-0019931 c rmesh11:p dose 2 trans 1 cora11 500 10i 1000 20i 1500 60i 4500 - corb11 -715 13i -65 12i 55 13i 705 + corb11 -715 13i -65 24i 55 15i 805 corc11 -200 2i -100 20i 100 5i 600 c Official ESS dose conversion factors from gamma/cm²/s to Sv/h c Is described in ESS-0019931 diff --git a/tallies/full_mesh_p.i b/tallies/full_mesh_p.i index 1151645..66841f6 100644 --- a/tallies/full_mesh_p.i +++ b/tallies/full_mesh_p.i @@ -2,7 +2,7 @@ c large and detailed mesh tally covering the beamline from monolith to full cave tmesh rmesh1:p dose 1 trans 1 cora1 500 10i 1000 40i 1500 60i 4500 - corb1 -715 13i -65 24i 55 13i 705 + corb1 -715 13i -65 24i 55 15i 805 corc1 -200 1i -100 40i 100 5i 600 c Official ESS dose conversion factors from gamma/cm²/s to Sv/h c Is described in ESS-0019931