From 8f55b02ced76d56691254dcf140d422e0dbb400e Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Wed, 16 Feb 2022 15:18:14 +0100 Subject: [PATCH] Try fixing artifacts by merging mirror cells and add large mesh tallies --- e2-estia.i | 6 +++--- geometry/estia_mf.i | 36 ++++++++++++++++++------------------ geometry/estia_surfaces.i | 4 ++-- geometry/sm_reflectivity.i | 8 ++++---- tallies/middle_focus_mesh.i | 28 +++++++++++++++++++--------- 5 files changed, 46 insertions(+), 36 deletions(-) diff --git a/e2-estia.i b/e2-estia.i index 61e08eb..e166243 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -197,7 +197,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 291r + 0 279r 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 @@ -208,7 +208,7 @@ fcl:n 0 276r $ in monolith 0 13r -1.0 6r $ NFGA 0 206r -1 $ Glass in Selene 1 0 53r -1 7r $ S1 back plate and mineral cast in universe - 0 211r + 0 199r c c ----------------------------------------------------------- c ------------------------------------------------------- @@ -225,7 +225,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 842r +imp:n,p,h,/,z 1 276r 0 0 1 830r 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 diff --git a/geometry/estia_mf.i b/geometry/estia_mf.i index e4890fe..d816695 100644 --- a/geometry/estia_mf.i +++ b/geometry/estia_mf.i @@ -73,27 +73,27 @@ c Magnets right/left (230553 -230554): (230555 -230556): (230557 -230558): (230559 -230560)) c Mirrors (silicon) -230601 90 -2.33 230562 -230563 230610 -230609 -230603 230606 -230602 90 -2.33 230562 -230563 230609 -230608 -230602 230605 -230603 90 -2.33 230562 -230563 230608 -230607 -230601 230604 -230604 90 -2.33 230562 -230563 230620 -230619 -230613 230616 -230605 90 -2.33 230562 -230563 230619 -230618 -230612 230615 -230606 90 -2.33 230562 -230563 230618 -230617 -230611 230614 +230601 90 -2.33 230562 -230563 ((230610 -230609 -230603 230606): + (230609 -230608 -230602 230605): + (230608 -230607 -230601 230604)) +230604 90 -2.33 230562 -230563 ((230620 -230619 -230613 230616): + (230619 -230618 -230612 230615): + (230618 -230617 -230611 230614)) 230611 90 -2.33 230562 -230563 230622 -230623 230621 -230624 -230656 230612 90 -2.33 230562 -230563 230626 -230627 230625 -230628 -230666 c Mirror coatings -230651 89 -5.102 230562 -230563 230610 -230609 230603 -230653 -230652 89 -5.102 230562 -230563 230609 -230608 230602 -230652 -230653 89 -5.102 230562 -230563 230608 -230607 230601 -230651 -230654 89 -5.102 230562 -230563 230610 -230609 230656 -230606 -230655 89 -5.102 230562 -230563 230609 -230608 230655 -230605 -230656 89 -5.102 230562 -230563 230608 -230607 230654 -230604 -230661 89 -5.102 230562 -230563 230620 -230619 230613 -230663 -230662 89 -5.102 230562 -230563 230619 -230618 230612 -230662 -230663 89 -5.102 230562 -230563 230618 -230617 230611 -230661 -230664 89 -5.102 230562 -230563 230620 -230619 230666 -230616 -230665 89 -5.102 230562 -230563 230619 -230618 230665 -230615 -230666 89 -5.102 230562 -230563 230618 -230617 230664 -230614 +230651 89 -5.102 230562 -230563 ((230610 -230609 230603 -230653): + (230609 -230608 230602 -230652): + (230608 -230607 230601 -230651)) +230654 89 -5.102 230562 -230563 ((230610 -230609 230656 -230606): + (230609 -230608 230655 -230605): + (230608 -230607 230654 -230604)) +230661 89 -5.102 230562 -230563 ((230620 -230619 230613 -230663): + (230619 -230618 230612 -230662): + (230618 -230617 230611 -230661)) +230664 89 -5.102 230562 -230563 ((230620 -230619 230666 -230616): + (230619 -230618 230665 -230615): + (230618 -230617 230664 -230614)) 230671 89 -5.102 230562 -230563 230622 -230623 230671 -230621 -230656 230672 89 -5.102 230562 -230563 230622 -230623 230624 -230672 -230656 230673 89 -5.102 230562 -230563 230626 -230627 230673 -230625 -230666 diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index e54ca35..4996695 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -1338,8 +1338,8 @@ c Absorber Plates 230712 13 px -15.000 230713 3 pz -3.500 230714 3 pz 3.500 -230715 3 py 2.200 -230716 3 py -0.700 +230715 3 py 1.500 +230716 3 py -1.400 c Flight Tube 230801 13 cx 12.500 230802 13 cx 12.800 diff --git a/geometry/sm_reflectivity.i b/geometry/sm_reflectivity.i index f9be059..6914b5e 100644 --- a/geometry/sm_reflectivity.i +++ b/geometry/sm_reflectivity.i @@ -16,11 +16,11 @@ c Selene 2 REFLE41 208911 36 -206113 REFLE42 208912 36 -206113 c Polarizer -REFLE51 230611 55 -230663 -REFLE52 230612 55 -230662 +REFLE51 230611 55 -230661 +REFLE52 230612 55 -230661 REFLE53 230613 55 -230661 -REFLE54 230601 56 -230653 -REFLE55 230602 56 -230652 +REFLE54 230601 56 -230651 +REFLE55 230602 56 -230651 REFLE56 230603 56 -230651 REFLE57 230620 55 -230674 REFLE58 230624 56 -230672 diff --git a/tallies/middle_focus_mesh.i b/tallies/middle_focus_mesh.i index 93cc8a5..7181f64 100644 --- a/tallies/middle_focus_mesh.i +++ b/tallies/middle_focus_mesh.i @@ -1,13 +1,13 @@ c fine mesh around middle focus area to analyzer beam paths tmesh rmesh11:n dose 1 - cora11 2135 42i 2350 - corb11 -150 9i -50 22i -4 59i 2 23i 50 9i 150 - corc11 -200 6i -60 59i 60 11i 300 + cora11 1460 158i 4640 + corb11 -400 59i 800 + corc11 -200 39i 600 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 + 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 @@ -41,10 +41,11 @@ c Is described in ESS-0019931 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 - rmesh21:p dose 2 - cora21 2135 42i 2350 - corb21 -150 9i -50 22i -4 59i 2 23i 50 9i 150 - corc21 -200 6i -60 59i 60 11i 300 +c + rmesh31:p dose 2 + cora31 1460 158i 4640 + corb31 -400 59i 800 + corc31 -200 39i 600 c Official ESS dose conversion factors from gamma/cm²/s to Sv/h c Is described in ESS-0019931 mshmf2 0.000000 1.00000e-04 @@ -76,4 +77,13 @@ c Is described in ESS-0019931 4.00E+03 8.75E-07 5.00E+03 9.04E-07 6.00E+03 9.29E-07 8.00E+03 9.65E-07 1.00E+04 9.94E-07 +c Finer meshes around the middle focus components + rmesh51:n dose 1 + cora51 2135 42i 2350 + corb51 -150 9i -50 22i -4 59i 2 23i 50 9i 150 + corc51 -200 6i -60 59i 60 11i 300 + rmesh61:p dose 2 + cora61 2135 42i 2350 + corb61 -150 9i -50 22i -4 59i 2 23i 50 9i 150 + corc61 -200 6i -60 59i 60 11i 300 endmd