diff --git a/e2-estia.i b/e2-estia.i index 3cea4d4..90dc415 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -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 diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index 7924d23..35d90ee 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -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 diff --git a/geometry/estia_universes.i b/geometry/estia_universes.i index 0354245..83a5faa 100644 --- a/geometry/estia_universes.i +++ b/geometry/estia_universes.i @@ -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. diff --git a/tallies/boundary_meshes.i b/tallies/boundary_meshes.i new file mode 100644 index 0000000..b081a82 --- /dev/null +++ b/tallies/boundary_meshes.i @@ -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