From dff3cdebbc1e22a651fc66307df0a3b322f54be7 Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Sat, 8 Jun 2019 13:12:56 +0200 Subject: [PATCH] Set physics, mode and model parameters after clarification with Alan Takibayevv --- e2-estia.i | 81 +++++++++++++++++++++--------------------------------- 1 file changed, 32 insertions(+), 49 deletions(-) diff --git a/e2-estia.i b/e2-estia.i index 8ccc8d2..d8fbeb1 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -94,7 +94,7 @@ c ------------------------------------------------------- c --------------- WEIGHT CARDS -------------------------- c ------------------------------------------------------- c weight window generator, can be run with previous weight window files -WWGE:n 1e-3 1e-1 1 10 100 3e3 +c WWGE:n 1e-3 1e-1 1 10 100 3e3 c mesh based weight windows, can be reused after change of cell structure c WWG 15 0 0 c MESH GEOM=CYL ORIGIN=0 0 -500 REF=133.7 -176.7 13.5 AXS=0 0 1 @@ -106,13 +106,7 @@ c JINTS= 2 4 4 8 8 4 4 6 c KMESH=298 303 309 315 360 c KINTS= 14 1 20 1 2 c generate cell based weight windows, needs to be updated on geometry change -WWG 35 61003 0 -c simple weight window game for other particles, high inside monolith -c source weight is 1e16, give large window with low split/roulette factor -c check only on collisions -c WWE:/,d,t,s,a 3e3 -c WWP:/,d,t,s,a 20 10 2 -1 0 -c WWN1:/,d,t,s,a 1e15 276r -1 -1 1e10 578r +c WWG 35 61003 0 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ c ----------------------------------------------------------- c ------------------- TALLY CARDS --------------------------- @@ -144,21 +138,18 @@ c read file=tallies/experimental_cave_f5.i echo c read file=tallies/cave_mesh.i echo c read file=tallies/selene1_f6.i echo c -c ========================================== -c -c ========================================== -c -c 1 2 3 4 5 6-14 15-22 23 24 25 26-34 35 36 37 38 39 -c 40-42 43-266 267 268-488 489-494 // -c 495 496 497 498 499 500 501 502 503 504 505 506 507 508-537 // -c 538 539 540-543 544-554 555-585 // -c 586 587-591 592 593 594-601 602 603-end// +c =============================== +c Global Variance Reduction Cards +c =============================== +c Exponential transform in collimators ext:n 0 276r $ in monolith 0 0 0 309r $ ess outside 0 132r 0.25v3 13r $CPC1+2 0 28r 0.25v3 4r $ bunker wall 0 16r 0.25v3 6r $ BWC 0 63r +c exponential transform reference vectors: origin, selene 1, sample +vect v1 0 0 0 v2 1000 -1394 -6.5 v3 2048 -2856 -30 c forced collisions in neutron guides fcl:n 0 276r $ in monolith 0 309r $ ess outside @@ -169,51 +160,43 @@ fcl:n 0 276r $ in monolith 0 10r -1 -1 $ S2 glass + carrier 0 22r -1 $ S1 carrier (universe) 0 8r -c exponential transform reference vectors: origin, selene 1, sample -vect v1 0 0 0 v2 1000 -1394 -6.5 v3 2048 -2856 -30 c c ----------------------------------------------------------- c ------------------------------------------------------- c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 -nps 1e9 +nps 1e4 c histp -mode n p h / - d $ deuteron - t $ triton - s $ helion - a $ alpha -c in ESS model from target: -c MODE N $ neutron -c P $ photon -c H $ proton -c / $ pion-+ -c | $ mu-- -c U $ nu-e -c V $ nu-mu -c Z $ pion-0 -c volume calculated for tally cells (spheres) -c and set for chopper disk/spindle and stages x1/x2/z1/z2/rot -c vol 1 506r 11j 1 183r 3848.5 10104 1j 1 7r 479.5 479.5 913 913 1734 -c 1 96r $ wall segments -c -234 -imp:h,/,d,t,s,a 1 276r 0 0 1 578r -imp:n,p 1 276r 0 0 1 578r -cut:h,/,d,t,s,a 1e+08 1e-3 +mode n p h / z $ neutron, photon, proton, pi+, pi0 +imp:n,p,h,/,z 1 276r 0 0 1 578r +cut:h,/,z 1e+08 1e-3 phys:n 3000 0 -phys:p 100 -phys:/,d,t,s,a 3000 -phys:h 3000 -lca 2 1 1 0023 1 1 0 1 1 0 -lea 1 4 1 0 1 0 0 1 +phys:p 3000 +phys:h,/,z 3000 +c Selection of used physics models, default from MCNP6.2 is CEM03+LAQGSM03 +LCA 2J +C 00 $ No ISABEL (default with INCL4) + 01 $ BERTINI for nucleons/pions and ISABEL for other particles +C 02 $ ISABEL for all particles + 5J +C 00 $ BERTINI or ISABEL + 01 $ CEM03 +C 02 $ INCL4 +C 00 $ ISABEL <940MeV/nucleon then LAQGSM03 + 01 $ LAQGSM03 for all heavy ion interactions +C ------------------------------------------------------------------------------ +LEA 6J + 00 $ DRESNER +C 02 $ ABLA + 1J c * For testing for geometry errors c read file=sources/geometry_test.i echo c * Source for simulation from proton beam (needs different WW and ESplit) -c read file=sources/proton.i echo +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 * McStas source terms for gamma doeses c read file=sources/mcstas_source_vs.i echo c read file=sources/mcstas_source_sample.i echo