Implement realistic shutter model with simple wall separator

This commit is contained in:
2021-03-11 13:48:33 +01:00
parent 34e0bb49fd
commit edb3007b90
5 changed files with 152 additions and 31 deletions

View File

@@ -18,6 +18,7 @@ c tallies
read file=geometry/estia_inbunker.i echo
read file=geometry/estia_guides.i echo
read file=geometry/estia_cave.i echo
read file=geometry/estia_shutter.i echo
read file=geometry/estia_universes.i echo
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
@@ -80,6 +81,10 @@ c Flange axis for Selene 1
TR21 1334.9578 0.36581 -2.4264 0.99992537 0 -0.012217
0 1 0
0.012217 0 0.99992537
c Shutter coordinate system centered in Y,Z flange and X to center
TR22 2149.5000 0.00000 -15.0000 1 0 0
0 1 0
0 0 1
c McStas coordinates with origin at virtual source
TR54 1114.9743 -2.1342 0.261 0 1 0
0.012217 0 0.99992537
@@ -168,9 +173,9 @@ c read file=tallies/bunker_wall_f5.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
read file=tallies/shutter_f4.i echo
c read file=tallies/shutter_f5.i echo
read file=tallies/shutter_wall_f4.i echo
c read file=tallies/shutter_wall_f4.i echo
read file=tallies/shutter_wall_mesh.i echo
c read file=tallies/detector_f5.i echo
c read file=tallies/selene1_walls_f4.i echo
@@ -190,7 +195,7 @@ ext:n 0 276r $ in monolith
0 146r 0.75v2 8r 0.75v3 4r $CPC1+2
0 28r 0.75v3 4r $ bunker wall
0 37r 0.75v3 7r $ BWC
0 230r
0 246r
c exponential transform reference vectors: origin, VS, shutter
vect v1 0 0 0 v2 1114.8 -1.75 0.5 v3 2137 -1.25 -12.5
c forced collisions in neutron guides
@@ -201,7 +206,7 @@ fcl:n 0 276r $ in monolith
0 13r -1.0 6r $ NFGA
0 208r -1 $ Glass in Selene 1
0 51r -1 7r $ S1 back plate and mineral cast in universe
0 150r
0 166r
c
c -----------------------------------------------------------
c -------------------------------------------------------
@@ -218,7 +223,7 @@ c random number generator options
rand gen=2 seed=29976355590901 $ default 19073486328125
c histp
mode n p h / z $ neutron, photon, proton, pi+, pi0
imp:n,p,h,/,z 1 276r 0 0 1 781r
imp:n,p,h,/,z 1 276r 0 0 1 797r
c Reduce model size for WWG runs until earlier stages of the model:
c void 201255 201257 203123 $ void cells downstream of BBG
c imp:n,p,h,/,z 1 276r 0 0 1 323r 0 1 0 1 34r 0 1 414r $ zero for above

View File

@@ -35,11 +35,10 @@ c Vacuum before and inside Selene Guide 1
#(208401 208402 208403 208404 -208405 208406)
#(208414 208415 208416 -208404 -208417 208418)
c Vacuum after Selene 1
206104 0 208436 -208100 -208922 100003
206105 5 0.0582256 -208426 208922 -208921 100003 -208100
206104 0 208436 -208103 -208922 100003
206105 5 0.0582256 -208426 208922 -208921 100003 -208103
c Vacuum from shutter to Selene 2
206111 0 ((208100 -208102 -208501 -208502 -208503 -208504)
:(208101 208635 -208932)) 100003
206111 0 208101 208635 -208932 100003
#100707 #208450
206112 5 0.0582256 -208426 -208625 208932 -208931 100003 208101
c Vacuum inside Selene 2
@@ -193,14 +192,12 @@ c area around tank with structrual beams (front, back, top, bottom, left, right)
0.012217 0 0.99992537)
208328 3 0.0876276 208441 208442 208445 208446 208424 -208447
#(208421 208422 208425 208426)
c Steel block in shutter wall
208401 54 0.0833854 (208100 -208101 -208221 -208222 -208223 -208224 100003)
#(208100 -208101 -208501 -208502 -208503 -208504)
c Shutter B4C and Lead
208402 20 -2.3 208102 -208103 -208501 -208502 -208503 -208504 100003
208403 13 -11.34 208103 -208101 -208501 -208502 -208503 -208504 100003
c 208402 0 208102 -208103 -208501 -208502 -208503 -208504 100003
c 208403 0 208103 -208101 -208501 -208502 -208503 -208504 100003
c Area feeding through shutter wall air, steel, concrete
208401 4 -1.2E-3 (208100 -208103 208921 -208923)
208402 54 0.0833854 (208100 -208101 -208221 -208222 -220206 -208224 208923)
#(208103 220002 -220003 220004 -220005)
208403 49 0.074981 (208100 -208101 -208221 -208222 220206 -208223)
#(208103 220002 -220003 220004 -220005)
c Copper in Selene 1
208411 14 -8.92 208111 -208112 100003 208402 ((208520 -208903)
:(208511 -208901)

38
geometry/estia_shutter.i Executable file
View File

@@ -0,0 +1,38 @@
c instrument shutter and surrounding elements
c outer boundary of shutter: X: 208103 -208101 Y: 220002 -220003 Z: 220004 -220005
c air around shutter
220001 4 -1.2E-3 208103 -208101 220002 -220003 220004 -220005
#220011 #220013
#(220000 -220001 220207 -220208 220204 -220205)
c Aluminum vacuum housing
220011 5 0.0582256 220000 -220001 220002 -220003 220004 -220204
220012 5 0.0582256 220000 -220001 220207 -220208 220204 -220205
#(220200 -220201 220202 -220203)
#(-220100)
220013 5 0.0582256 220000 -220001 220002 -220003 220205 -220206
c Copper inserts
220021 14 -8.92 220200 -220304 220202 -220203 220204 -220205
#(220300 -220301 220302 -220303)
220022 14 -8.92 220304 -220305 220202 -220306 220204 -220205
220023 14 -8.92 220304 -220305 220307 -220203 220204 -220205
220024 14 -8.92 220305 -220201 220202 -220203 220204 -220205
#(220300 -220301 220302 -220303)
c Shutter block B4C and Lead in Al casing
c closed shutter
220031 5 0.0582256 220400 -220401 220402 -220403 220404 -220405
#(220406 -220408 220410 -220411 220412 -220413)
220032 20 -2.3 220406 -220407 220410 -220411 220412 -220413
220033 13 -11.34 220407 -220408 220410 -220411 220412 -220413
c open shutter
c 220031 5 0.0582256 220400 -220401 220402 -220403 220420 -220423
c #(220406 -220408 220410 -220411 220421 -220422)
c 220032 20 -2.3 220406 -220407 220410 -220411 220421 -220422
c 220033 13 -11.34 220407 -220408 220410 -220411 220421 -220422
c Vacuum elements
220051 0 220000 -220200 -220100
220052 0 220200 -220304 220300 -220301 220302 -220303
220053 0 220304 -220305 220306 -220307 220204 -220205
#(220400 -220401 220402 -220403 220404 -220405) $ closed
c #(220400 -220401 220402 -220403 220420 -220423) $ opened
220054 0 220305 -220201 220300 -220301 220302 -220303
220055 0 220201 -220001 -220100

View File

@@ -429,18 +429,18 @@ c Vacuum Vessel
c Pre copper collor
c ----------- Ouside Bunker --------------
c Wall radii
208100 p 1.0000002116649709 1.148348833801549e-07 0.0 2112.01
208101 p 1.0000002116649709 1.148348833801549e-07 0.0 2157.17
208102 cz 2149.17
208103 cz 2150.17
208111 cz 1662.0
208112 cz 1672.0
208113 cz 1710.0
208114 cz 1720.0
208115 cz 1758.0
208116 cz 1768.0
208121 cz 2555.54
208122 cz 3275.49
208100 px 2112.0
208101 px 2164.0
208102 px 2132.0
208103 px 2135.0
208111 px 1662.0
208112 px 1672.0
208113 px 1710.0
208114 px 1720.0
208115 px 1758.0
208116 px 1768.0
208121 px 2550.0
208122 px 3275.0
c Shielding wall geometry
c first section outside
208201 p 1.148348833801549e-07 -1.0000002116649709 0.0 -148.964
@@ -462,8 +462,8 @@ c B4C concrete/steel interface
c wall breakthrough
208221 p -1.148348833801549e-07 1.0000002116649709 0.0 32.866
208222 p 1.148348833801549e-07 -1.0000002116649709 0.0 43.134
208223 pz 18.252
208224 p 0.0 0.0 -1.0 41.748
208223 pz 105.0
208224 p 0.0 0.0 -1.0 50.0
c Optical Cave Part
c Outside
208302 p -0.1201368660871055 -0.9927569472110325 0.0 -80.725
@@ -625,6 +625,7 @@ c Selene 2 Inside Mirrors Ellipses Top-Bottom/Right
c Vacuum vessels before S1/S1-S2
208921 21 cx 13.35
208922 21 cx 13.05
208923 21 cx 16.00
208931 3 c/x 0 -2 16.0
208932 3 c/x 0 -2 15.0
208933 3 c/x 0 -2 16.5
@@ -819,6 +820,60 @@ c Hatch SEE Chicane
c
c Shared Wall
210500 px 3629.979
c Shutter geometry
c outer boundary
220000 22 px -12.5
220001 22 px 12.5
220002 22 py -31.0
220003 22 py 31.0
220004 22 pz -17.0
220005 22 pz 118.0
c inner and outer vacuum flange cylinders
220100 22 cx 12.5
220101 22 cx 13.0
c Vacuum vessel inside (Al 2 cm thickness) and top
220200 22 px -10.5
220201 22 px 10.5
220202 22 py -26.0
220203 22 py 26.0
220204 22 pz -15.0
220205 22 pz 42.0
220206 22 pz 44.0
220207 22 py -28.0
220208 22 py 28.0
c Copper inserts
c window
220300 22 py -3.75
220301 22 py 3.75
220302 22 pz -8.0
220303 22 pz 8.0
c up-stream/center/downstream
220304 22 px -2.5
220305 22 px 5.5
c opening for shutter
220306 22 py -11.25
220307 22 py 11.25
c Shutter block outside
220400 22 px -1.46
220401 22 px 5.5
220402 22 py -11.0
220403 22 py 11.0
220404 22 pz -14.9
220405 22 pz 16.9
c Inside B-Al and Lead
220406 22 px -0.96
220407 22 px 0.0
220408 22 px 5.2
220410 22 py -6.0
220411 22 py 6.0
220412 22 pz -12.0
220413 22 pz 12.0
c vertical location when opened
220420 22 pz 10.1
220421 22 pz 13.0
220422 22 pz 37.0
220423 22 pz 41.9
c
c UNIVERSES
c Selene 1
c Chopper Pit final design geometry

26
pshutter.com Normal file
View File

@@ -0,0 +1,26 @@
origin 2150 0 -14.5 &
basis 1 0 0 0 1 0 &
extent 50 50 &
label 0 0 &
shade 4=LightYellow 9=LightGoldenrod 300=LightGoldenrod &
shade 5=green 10=DarkSeaGreen &
shade 49=SlateGray 16=DarkOliveGreen 200=DarkSlateGray &
shade 71=RosyBrown 3=red 54=RosyBrown 121=brown &
shade 14=orange &
shade 13=purple 120=DarkOliveGreen &
shade 18=magenta &
shade 12=cyan &
shade 2=PaleGreen 19=DarkOliveGreen 17=DarkOliveGreen &
shade 101=blue 115=blue 400=DodgerBlue &
shade 650=LightCyan 815=SlateGray 826=DarkSlateGray &
shade 870=LightYellow 1360=green 1361=green &
shade 1400=DarkOrchid 2629=VioletRosyBrown 2634=RosyBrown &
shade 2636=RosyBrown 2642=RosyBrown 2644=RosyBrown 2645=RosyBrown 2650=RosyBrown &
shade 2900=orange 6101=RosyBrown 6102=RosyBrown 6104=RosyBrown 6105=RosyBrown &
shade 6108=RosyBrown 6110=magenta 6113=magenta 6114=magenta 6115=magenta 6116=magenta &
shade 6117=magenta 6119=magenta 6120=coral 6121=coral 6122=coral 6123=coral &
shade 6125=coral 6129=coral 6132=purple 6133=purple 6136=purple 6138=purple &
shade 6140=purple 6141=purple 7410=DarkOliveGreen &
shade 31=pink 32=LightPink &
shade 601=SteelBlue 602=green 603=green 604=red 605=orange 606=purple &
shade 608=turquoise &