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

Add mesh tallies for instrument boundaries with 20x20x20 grid

Browse Source
This commit is contained in:
glavic_a
2019-10-03 17:33:14 +02:00
parent 757fb0837d
commit 9fcff6c01a
4 changed files with 85 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
e2-estia.i
View File

@@ -92,6 +92,23 @@ 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
c Coordingates for instrument boundary meshes, y-direction points outside
c SKADI side wall starting at bunker boundary
TR121 962.0 -1155.0 -150.0 0.675590 -0.737277 0.000000
0.737277 0.675590 0.000000
0 0 1
c Fence on E01 side and left cave side
TR122 1089.0 -2723.0 -180.0 0.582123 -0.813101 0.000000
-0.813101 -0.582123 0.000000
0 0 1
c Downstream cave side
TR123 2841.0 -3360.0 -180.0 -0.813101 -0.582123 0.000000
0.582123 -0.813101 0.000000
0 0 1
c Dillitation joint
TR124 1104.0 -1820.0 -180.0 0.000000 -1.000000 0.000000
-1.000000 0.000000 0.000000
0 0 1
read file=geometry/ess_monolith_transformations.i echo
c -------------------------------------------------------
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
@@ -123,6 +140,7 @@ 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
c read file=tallies/full_mesh_np.i echo
read file=tallies/boundary_meshes.i echo
c read file=tallies/NBPI_f4.i echo
c read file=tallies/BBG_f4.i echo
c read file=tallies/BBG_f5.i echo
@@ -175,7 +193,7 @@ 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 1e6 CTME 4e5 f31 0.01
stop nps 1e7 CTME 4e5 $f31 0.01
rand gen=2 seed=29976355590901 $ default 19073486328125
c histp
mode n p h / z $ neutron, photon, proton, pi+, pi0

5
geometry/estia_surfaces.i
View File

@@ -844,6 +844,11 @@ c 250001 1 rpp 1339.25 2060.75 2 8 2 8
c 250002 1 rpp 2539.25 3260.75 -8 -2 -8 -2
c visualize Shutter Source
c 250003 3 rpp -170.5 -169.5 0 4 -6 6
c visualize boundary mesh geometry
c 250121 121 rpp -10 2910 -10 10 -10 310
c 250122 122 rpp -10 1550 -10 10 -10 310
c 250123 123 rpp -10 1050 -10 10 -10 310
c 250124 124 rpp -10 930 -10 10 -10 310
c Remote Handling Flange in local coordinate system
500000 so 50
c planes in beam direction

5
geometry/estia_universes.i
View File

@@ -16,6 +16,11 @@ c 18323 0 -18323 U=2
218325 0 -218325 U=2
c 18326 0 -18326 U=2
c
c just for plotting mesh geometry as transformation has bug in plot
c 300121 0 -250121
c 300122 0 -250122
c 300123 0 -250123
c 300124 0 -250124
c Remote handling flange with local coordinate system,
c has to be moved and rotated when filling universe and includes
c surrounding air and internal vacuum tube.

56
tallies/boundary_meshes.i Normal file
View File

@@ -0,0 +1,56 @@
c meshes around the instrument boundary using 20x20x20 grid size
tmesh
rmesh1:n dose 1 trans 121
cora1 -10 145i 2910
corb1 -10 10
corc1 -10 15i 310
rmesh11:n dose 1 trans 122
cora11 -10 77i 1550
corb11 -10 10
corc11 -10 15i 310
rmesh21:n dose 1 trans 123
cora21 -10 52i 1050
corb21 -10 10
corc21 -10 15i 310
rmesh31:n dose 1 trans 124
cora31 -10 46i 930
corb31 -10 10
corc31 -10 15i 310
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