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

Add new surface source in bunker wall and test it

Browse Source
This commit is contained in:
glavic_a
2019-10-03 11:56:33 +02:00
parent f1dd1fd74d
commit 757fb0837d
4 changed files with 265 additions and 13 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
e2-estia.i
View File

@@ -87,6 +87,11 @@ TR101 0.0 0.0 0.0 0.58778525 -0.80901699 0.0 &
0.0 0.0 1.0 1
c Transformation from SKADI model coordinate system here
TR102 0 5.0 1.0 0.66913061 -0.74314483 0 0.74314483 0.66913061 0 0 0 1 1
c Coordinate for surface source within bunker wall
c close to plane 206305 with same surface normal and center around beam
TR103 700.03 -977.88 -0.86 0.587785 -0.809017 0.000000
0.809017 0.587785 0.000000
0 0 1
read file=geometry/ess_monolith_transformations.i echo
c -------------------------------------------------------
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
@@ -123,15 +128,15 @@ c read file=tallies/BBG_f4.i echo
c read file=tallies/BBG_f5.i echo
c read file=tallies/inbunker_f4.i echo
c read file=tallies/chopper_pit_doses.i echo
read file=tallies/bunker_feedthrough_f1.i echo
c read file=tallies/bunker_wall_f4.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
c read file=tallies/shutter_wall_f4.i echo
read file=tallies/shutter_wall_f4.i echo
c read file=tallies/detector_f5.i echo
c read file=tallies/selene1_walls_f4.i echo
c read file=tallies/optical_cave_f5.i echo
@@ -170,12 +175,12 @@ c --------------- PHYSICS CARDS --------------------------
c -------------------------------------------------------
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
stop nps 1e9 CTME 4e5 f31 0.01
stop nps 1e6 CTME 4e5 f31 0.01
rand gen=2 seed=29976355590901 $ default 19073486328125
c histp
mode n p h / z $ neutron, photon, proton, pi+, pi0
c imp:n,p,h,/,z 1 276r 0 0 1 589r
imp:n,p,h,/,z 1 276r 0 0 1 487r 0 0 0 1 0 87r 1 9r
imp:n,p,h,/,z 1 276r 0 0 1 589r
c imp:n,p,h,/,z 1 276r 0 0 1 487r 0 0 0 1 0 87r 1 9r
cut:h,/,z 1e+08 1e-3
phys:n 3000 0
phys:p 3000
@@ -202,7 +207,9 @@ 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
c read file=sources/e2_2m.i echo
read file=sources/e2_2m_mod.i echo
c read file=sources/e2_2m_mod.i echo
c * Estimated source within bunker wall to optimize collimator efficiency
read file=sources/bunker_wall.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
@@ -210,7 +217,7 @@ c read file=sources/mcstas_source_mf.i echo
c read file=sources/mcstas_gammas.i echo
c read file=sources/mcstas_gammas_2.i echo
c read file=sources/mcstas_source_shutter.i echo
prdmp 1e8 -30 1 2 1e7
prdmp 1e6 -30 1 2 1e5
print 30 32 40
c 120 $O importance function quality
c 128 $O universe map

8
geometry/estia_guides.i
View File

@@ -18,10 +18,10 @@ c #(-5102 5501 5502 -5505 5504)
c #(5102 -7100 -5531)
206002 4 -1.2E-3 206305 -207205 -206401 -206402 -206403 -206404
#207022 #207024 #207025 #207031 #207032 #207033
c fill=50 (703.47 -983.22 -0.89
c 0.58207953 -0.81304008 0.0
c -0.81310076 -0.58212297 0.0
c 0.0 0.0 1.0)
fill=50 (703.47 -983.22 -0.89
0.58207953 -0.81304008 0.0
-0.81310076 -0.58212297 0.0
0.0 0.0 1.0)
206003 4 -1.2E-3 207205 -207222 -206401 -206402 -206403 -206404
#(207201 -207202 -207203 207204)
206004 4 -1.2E-3 207206 -206405 -206401 -206402 -206403 -206404 208921

1
geometry/estia_surfaces.i
View File

@@ -381,6 +381,7 @@ c Segmentation surfaces for bunker feedthrough f1 tally
206702 4 c/x 3.0 0.0 7
206703 4 c/x 3.0 0.0 13
206704 4 c/x 3.0 0.0 30
c 206705 103 rpp 0 2 -10 10 -10 10 $ to show location of surface cell TR103
c Cu-collimator
c Radii
207102 cz 1271.50

244
sources/bunker_wall.i Normal file
View File

@@ -0,0 +1,244 @@
c source definition based on tally_181.txt, modified for high E flux
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 1.0E+8 1j
c the source definition with neutrons/s
sdef tr=103 pos=0 0 0 vec=1 0 0
dir=d1 rad=fdir d2 erg=fdir d3
par=n ara=11309.7 wgt=2.0e+8
c distribution of radial direction
si1 0
5.00000E-01 8.66025E-01 9.65926E-01
9.96195E-01 9.99391E-01 1.00000E+00
sp1 0.00000000e+00
4.61452024e+07 1.70198257e+08 1.08329950e+08
6.38706911e+07 9.10952408e+06 2.26490368e+06
c directional dependent distribution of position
ds2 q 0 21
5.00000E-01 22
8.66025E-01 23
9.65926E-01 24
9.96195E-01 25
9.99391E-01 26
1.00000E+00 27
c Energy dependent probability for each direction pair
ds3 q 0 31
5.00000E-01 32
8.66025E-01 33
9.65926E-01 34
9.96195E-01 35
9.99391E-01 36
1.00000E+00 37
c direction probabilities for each segment
si21 0 1
sp21 0 1.0000E-20
c
si22 0 2 7 13 30 122.51
sp22 0 2.30230E+04 2.07305E+05 5.25726E+05 1.23189E+06 2.10846E+07
c
si23 0 2 7 13 30 122.51
sp23 0 6.76418E+04 1.27453E+06 2.04084E+06 6.08084E+06 7.56352E+07
c
si24 0 2 7 13 30 122.51
sp24 0 8.63986E+04 7.45148E+05 1.24766E+06 2.35238E+06 4.97334E+07
c
si25 0 2 7 13 30 122.51
sp25 0 2.43288E+04 4.55438E+05 5.77578E+05 1.32580E+06 2.95522E+07
c
si26 0 2 7 13 30 122.51
sp26 0 1.63116E+05 1.04990E+05 1.76650E+05 4.00354E+05 3.70965E+06
c
si27 0 2 7 13 30 122.51
sp27 0 1.68399E+04 4.20802E+04 9.41964E+04 8.99954E+04 8.89340E+05
c energy dependent probability for each direction
si31 0 1
sp31 0 1.0000E-20
C ------------------------------------------------------------------------
si32 0.00000E+00
1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10
1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09
1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08
1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07
1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06
1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05
1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04
1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03
1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02
1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01
1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00
1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01
1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02
1.00000E+03 1.58000E+03 2.51000E+03
sp32 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
3.14844E+05 1.50026E+05 0.00000E+00 0.00000E+00 0.00000E+00
5.04168E+05 0.00000E+00 2.47386E+05 0.00000E+00 0.00000E+00
1.19420E+05 0.00000E+00 0.00000E+00 2.83700E+05 5.18008E+05
4.50826E+05 1.17780E+05 0.00000E+00 0.00000E+00 3.24702E+05
0.00000E+00 3.04862E+05 5.05966E+04 1.27270E+05 3.51030E+05
3.12144E+05 2.54320E+05 4.39610E+05 1.89490E+05 3.07168E+05
5.67617E+05 5.24070E+05 6.84092E+05 1.08998E+06 1.05406E+06
1.56784E+06 2.06740E+06 1.63649E+06 1.28594E+06 9.86652E+05
9.00414E+05 8.92366E+05 1.09287E+06 1.32177E+06 1.10393E+06
6.95039E+05 1.87933E+05 4.50658E+04 1.73907E+03 0.00000E+00
0.00000E+00 0.00000E+00 2.30726E+07
c
si33 0.00000E+00
1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10
1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09
1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08
1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07
1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06
1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05
1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04
1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03
1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02
1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01
1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00
1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01
1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02
1.00000E+03 1.58000E+03 2.51000E+03
sp33 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 8.96446E+04
1.24880E+05 2.83058E+05 2.51954E+05 6.87496E+05 7.15686E+05
6.94112E+05 9.17851E+05 1.41215E+05 6.94870E+04 3.30944E+05
2.24454E+05 1.15589E+06 3.71718E+05 6.70040E+05 6.24094E+05
9.71906E+05 7.78151E+05 3.52868E+05 2.61442E+05 2.02862E+05
3.12048E+05 4.24589E+05 4.77852E+05 1.01942E+06 3.71538E+05
5.74270E+05 5.68819E+05 7.53604E+05 9.31588E+05 5.57942E+05
8.05481E+05 8.61544E+05 8.70714E+05 8.55040E+05 9.96667E+05
1.90695E+06 1.83161E+06 2.15864E+06 2.40250E+06 4.03989E+06
3.85865E+06 4.51882E+06 4.32260E+06 3.50218E+06 2.75366E+06
2.17294E+06 2.62505E+06 4.26290E+06 6.18216E+06 7.27704E+06
6.33514E+06 3.45029E+06 1.01333E+06 1.69420E+05 1.41528E+04
3.21214E+02 0.00000E+00 8.50991E+07
c
si34 0.00000E+00
1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10
1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09
1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08
1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07
1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06
1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05
1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04
1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03
1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02
1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01
1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00
1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01
1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02
1.00000E+03 1.58000E+03 2.51000E+03
sp34 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 3.50044E+05 2.11036E+05 7.25122E+05 4.01714E+05
2.85638E+05 0.00000E+00 2.89634E+05 5.13630E+05 2.83914E+05
7.57076E+05 1.06405E+05 1.20970E+05 3.91716E+05 1.32458E+05
0.00000E+00 1.28368E+05 0.00000E+00 3.83210E+05 3.27540E+05
4.95046E+05 4.01546E+05 1.45736E+05 3.95288E+05 4.07088E+05
2.15020E+05 2.84168E+05 3.56171E+05 2.78774E+05 4.20790E+05
5.99099E+05 4.42646E+05 1.72771E+05 4.67428E+05 3.10148E+05
6.13976E+05 8.44794E+05 9.57095E+05 1.36761E+06 1.50384E+06
1.43597E+06 2.37344E+06 1.74013E+06 1.49268E+06 1.20460E+06
8.77910E+05 1.15478E+06 1.96657E+06 4.33218E+06 6.25133E+06
6.14042E+06 5.83465E+06 2.37386E+06 7.54202E+05 1.23595E+05
1.91279E+04 1.98462E+03 5.41650E+07
c
si35 0.00000E+00
1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10
1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09
1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08
1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07
1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06
1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05
1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04
1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03
1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02
1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01
1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00
1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01
1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02
1.00000E+03 1.58000E+03 2.51000E+03
sp35 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 5.80822E+04
7.82170E+04 0.00000E+00 0.00000E+00 0.00000E+00 1.18801E+05
2.03312E+05 9.17614E+04 0.00000E+00 0.00000E+00 3.13736E+05
0.00000E+00 2.85096E+04 1.01698E+05 0.00000E+00 0.00000E+00
4.15596E+04 1.82898E+05 0.00000E+00 1.31384E+05 4.61576E+04
1.54784E+05 3.45572E+04 1.52229E+05 0.00000E+00 1.33941E+04
1.05392E+05 1.80908E+05 2.09700E+05 8.77208E+04 3.44191E+05
1.95796E+05 2.88313E+05 3.23630E+05 4.55422E+05 6.71276E+05
5.39635E+05 8.32866E+05 2.83956E+05 5.40580E+05 4.28812E+05
3.37518E+05 3.94562E+05 7.58188E+05 1.78063E+06 3.63320E+06
4.66280E+06 7.41933E+06 3.86641E+06 1.47653E+06 2.91329E+05
5.72558E+04 1.82814E+04 3.19354E+07
c
si36 0.00000E+00
1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10
1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09
1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08
1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07
1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06
1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05
1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04
1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03
1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02
1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01
1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00
1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01
1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02
1.00000E+03 1.58000E+03 2.51000E+03
sp36 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
1.14870E+05 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
1.77094E+03 1.95123E+03 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 6.35549E+04 2.26664E+04 0.00000E+00 0.00000E+00
2.35992E+04 7.97576E+04 5.61072E+04 3.09262E+04 4.60783E+04
5.68676E+04 9.41971E+04 2.47496E+04 1.27729E+04 6.29456E+04
3.05540E+04 9.95486E+04 5.93783E+04 2.44680E+05 5.31452E+05
5.44316E+05 1.09089E+06 7.82100E+05 3.84272E+05 7.48612E+04
1.54286E+04 4.46045E+03 4.55476E+06
c
si37 0.00000E+00
1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10
1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09
1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08
1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07
1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06
1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05
1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04
1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03
1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02
1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01
1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00
1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01
1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02
1.00000E+03 1.58000E+03 2.51000E+03
sp37 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
5.63064E-10 0.00000E+00 0.00000E+00 0.00000E+00 4.99236E-02
0.00000E+00 0.00000E+00 1.17927E-04 0.00000E+00 1.87258E+05
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
0.00000E+00 1.51526E-02 0.00000E+00 0.00000E+00 1.82373E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
1.94199E+03 0.00000E+00 4.92952E+04 0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.85002E+04
0.00000E+00 2.37029E+04 0.00000E+00 1.65305E+04 0.00000E+00
1.11586E+04 7.41496E+04 8.46998E+03 5.06353E+03 5.09402E+01
9.70764E+01 1.31903E+04 2.21342E+04 3.09494E+04 5.67470E+04
9.29058E+04 2.40455E+05 1.69563E+05 7.28883E+04 2.05490E+04
6.57676E+03 2.72816E+02 1.13245E+06