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

Fix VS opening height run McStas rss with fine mesh around VS to check axis

Browse Source
This commit is contained in:
glavic_a
2020-03-03 14:20:54 +01:00
parent b1788d345c
commit 86b637444f
3 changed files with 63 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
e2-estia.i
View File

@@ -149,7 +149,8 @@ c read file=tallies/bunker_mesh_np.i echo
c read file=tallies/bunker_mesh_act.i echo
c read file=tallies/full_mesh.i echo
c read file=tallies/full_mesh_p.i echo
read file=tallies/full_mesh_np.i echo
c read file=tallies/full_mesh_np.i echo
read file=tallies/vs_fine_mesh.i echo
c read file=tallies/boundary_meshes.i echo
c read file=tallies/boundary_wall_meshes.i echo
c read file=tallies/selene_source_mesh_p.i echo
@@ -204,8 +205,8 @@ 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
c stop nps 392426460 CTME 4e5 $f31 0.01
stop nps 1e8 CTME 4e5 $f31 0.01
stop nps 39242646 CTME 4e5 $f31 0.01
c stop nps 1e8 CTME 4e5 $f31 0.01
rand gen=2 seed=29976355590901 $ default 19073486328125
c histp
mode n p h / z $ neutron, photon, proton, pi+, pi0
@@ -237,18 +238,18 @@ 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
c read file=sources/bunker_wall.i echo $ approximate, not for official results
c * McStas source terms for gamma doeses
c read file=sources/mcstas_rssa.i echo
read file=sources/mcstas_rssa.i echo
c read file=sources/mcstas_source_vs.i echo
c read file=sources/mcstas_source_sample.i echo
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 5e6 -30 1 2 5e5
prdmp 5e5 -30 1 2 1e5
print 30 32 40
c 120 $O importance function quality
c 128 $O universe map

6
geometry/estia_surfaces.i
View File

@@ -771,8 +771,8 @@ c VS Blades (open)
224502 p -0.2588190294771241 0.9659264084517616 0.0 -290.42
224503 p -0.2588190294771241 0.9659264084517616 0.0 -290.52
224504 p 0.2588190294771241 -0.9659264084517616 0.0 292.92
221501 pz 5.061
221502 pz 1.061
221501 pz 5.011
221502 pz 1.011
221503 pz -4.739
221504 pz -0.739
224505 p 0.9659264084517616 0.2588190294771241 0.0 1073.54
@@ -782,7 +782,7 @@ c VS Blades (open)
224513 p -0.2588190294771241 0.9659264084517616 0.0 -292.856
224514 p -0.2588190294771241 0.9659264084517616 0.0 -290.906
221511 pz 5.261
221512 pz 0.539
221512 pz -0.489
221513 pz 1.761
221514 pz -4.539
224515 p -0.9659264084517616 -0.2588190294771241 0.0 -1078.921

53
tallies/vs_fine_mesh.i Normal file
View File

@@ -0,0 +1,53 @@
c large and detailed mesh tally covering the beamline from monolith to full cave
c plus one very detailed tally to get the source details around VS
tmesh
rmesh1:n dose 1 trans 101
cora1 500 9i 1000 22i 1460
corb1 -700 11i -100 9i 100 3i 300
corc1 -200 1i -100 9i 100 5i 400
rmesh11:n dose 1 trans 101
cora11 1460 151i 4500
corb11 -400 59i 800
corc11 -200 39i 600
rmesh21:n dose 1 trans 101
cora21 1000 399i 1400
corb21 -50 99i 50
corc21 -50 99i 50
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
1.00000e-09 1.11000e-08 1.00000e-08 1.28000e-08
2.50000e-08 1.44000e-08 1.00000e-07 1.87000e-08
2.00000e-07 2.11000e-08 5.00000e-07 2.37000e-08
1.00000e-06 2.53000e-08 2.00000e-06 2.66000e-08
5.00000e-06 2.78000e-08 1.00000e-05 2.82000e-08
2.00000e-05 2.82000e-08 5.00000e-05 2.82000e-08
1.00000e-04 2.80000e-08 2.00000e-04 2.78000e-08
5.00000e-04 2.71000e-08 1.00000e-03 2.71000e-08
2.00000e-03 2.74000e-08 5.00000e-03 2.87000e-08
1.00000e-02 3.28000e-08 2.00000e-02 4.39000e-08
3.00000e-02 5.65000e-08 5.00000e-02 8.28000e-08
7.00000e-02 1.10000e-07 1.00000e-01 1.51000e-07
1.50000e-01 2.18000e-07 2.00000e-01 2.84000e-07
3.00000e-01 4.10000e-07 5.00000e-01 6.37000e-07
7.00000e-01 8.35000e-07 9.00000e-01 1.00000e-06
1.00000e+00 1.08000e-06 1.20000e+00 1.19000e-06
1.50000e+00 1.31000e-06 2.00000e+00 1.47000e-06
3.00000e+00 1.65000e-06 4.00000e+00 1.74000e-06
5.00000e+00 1.78000e-06 6.00000e+00 1.79000e-06
7.00000e+00 1.80000e-06 8.00000e+00 1.80000e-06
9.00000e+00 1.80000e-06 1.00000e+01 1.80000e-06
1.20000e+01 1.80000e-06 1.40000e+01 1.78000e-06
1.50000e+01 1.77000e-06 1.60000e+01 1.76000e-06
1.80000e+01 1.74000e-06 2.00000e+01 1.72000e-06
2.10000e+01 1.71000e-06 3.00000e+01 1.63000e-06
5.00000e+01 1.56000e-06 7.50000e+01 1.58000e-06
1.00000e+02 1.60000e-06 1.30000e+02 1.61000e-06
1.50000e+02 1.61000e-06 1.80000e+02 1.61000e-06
2.00000e+02 1.61000e-06 3.00000e+02 1.67000e-06
4.00000e+02 1.79000e-06 5.00000e+02 1.92000e-06
6.00000e+02 2.05000e-06 7.00000e+02 2.16000e-06
8.00000e+02 2.24000e-06 9.00000e+02 2.30000e-06
1.00000e+03 2.35000e-06 2.00000e+03 2.76000e-06
5.00000e+03 3.64000e-06 1.00000e+04 4.75000e-06
endmd