From bd6fa655eb2bb9c43167bdb3b6df4f1c55b3930e Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Mon, 3 Jun 2019 19:09:08 +0200 Subject: [PATCH] Test run with old 2m source term and no WW --- e2-estia.i | 28 ++++++++++++++-------------- geometry/estia_cave.i | 8 +++++--- geometry/estia_surfaces.i | 8 +++++--- sources/e2_2m.i | 6 +++--- tallies/BBG_f4.i | 2 +- 5 files changed, 28 insertions(+), 24 deletions(-) diff --git a/e2-estia.i b/e2-estia.i index 1d99fb1..99ef251 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -94,17 +94,17 @@ 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 -WWG 15 0 0 +c WWG 15 0 0 c MESH GEOM=CYL ORIGIN=0 0 -500 REF=133.7 -176.7 13.5 AXS=0 0 1 -MESH GEOM=CYL ORIGIN=0 0 -500 REF=0 -1 2.5 AXS=0 0 1 - IMESH=20 200 550 1150 1500 2112.01 2157.17 4600 - IINTS=2 2 5 1 10 1 1 4 - JMESH=100 300 460 500 540 700 900 1200 - JINTS= 2 4 4 8 8 4 4 6 - KMESH=298 303 309 315 360 - KINTS= 14 1 20 1 2 +c MESH GEOM=CYL ORIGIN=0 0 -500 REF=0 -1 2.5 AXS=0 0 1 +c IMESH=20 200 550 1150 1500 2112.01 2157.17 4600 +c IINTS=2 2 5 1 10 1 1 4 +c JMESH=100 300 460 500 540 700 900 1200 +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 c WWG 55 1042 0 c simple weight window game for other particles, high inside monolith @@ -125,8 +125,8 @@ 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 c read file=tallies/NBPI_f4.i echo -c read file=tallies/BBG_f4.i echo -read file=tallies/BBG_f5.i echo +read file=tallies/BBG_f4.i echo +c read file=tallies/BBG_f5.i echo c read file=tallies/inbunker_f4.i echo c read file=tallies/chopper_pit_doses.i echo c read file=tallies/bunker_wall_f4.i echo @@ -169,7 +169,7 @@ c ------------------------------------------------------- c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 -nps 12e7 +nps 5e7 c histp mode n p h / d $ deuteron @@ -202,9 +202,9 @@ lea 1 4 1 0 1 0 0 1 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) -read file=sources/proton.i echo +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.i echo c read file=sources/e2_2m_act.i echo c * McStas source terms for gamma doeses c read file=sources/mcstas_source_vs.i echo diff --git a/geometry/estia_cave.i b/geometry/estia_cave.i index 92d6339..eaf4492 100644 --- a/geometry/estia_cave.i +++ b/geometry/estia_cave.i @@ -37,7 +37,9 @@ c area around cave 209901 4 -1.2E-3 208709 208710 100112 (208302:-208805) 208703 208701 -208706 209902 4 -1.2E-3 208720 -208710 100112 -100104 208703 208701 +c visualize ESS 2m source term (comment for real runs) +c 250000 0 -250000 c visualize Selene gamma sources (comment for real runs) -c 50001 0 -50001 -c 50002 0 -50002 -c 50003 0 -50003 +c 250001 0 -250001 +c 250002 0 -250002 +c 250003 0 -250003 diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index eb032f3..2813d04 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -970,9 +970,11 @@ c 221514 pz -4.539 224515 p -0.351842 0.936060 0.000000 -1078.921 224516 p -0.351842 0.936060 0.000000 -1079.221 +c visualize 2m source +c 250000 101 rpp 199.49 199.51 -9.0 1.0 3.25 19.75 c visualize Selene gamma sources -c 50001 1 rpp 1339.25 2060.75 2 8 2 8 -c 50002 1 rpp 2539.25 3260.75 -8 -2 -8 -2 +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 50003 3 rpp -170.5 -169.5 0 4 -6 6 +c 250003 3 rpp -170.5 -169.5 0 4 -6 6 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ diff --git a/sources/e2_2m.i b/sources/e2_2m.i index dce8925..f32c798 100644 --- a/sources/e2_2m.i +++ b/sources/e2_2m.i @@ -1,12 +1,12 @@ c source definition based on tally_181.txt -cut:n j 1e-11 +cut:n j 0 c cut all photons that can't produce neutrons cut:p 1e+08 1.667 c but adapted to model geometry and actual size of Estia beam extraction c make sure to supply the correct WWINP file -WWP:n j j 5 j -1 j 1.0 $1.79 +c WWP:n j j 2 j -1 j 1.0e-3 1j 1e17 1j $ weights are lower then from proton c the source definition with neutrons/s -sdef dir=d1 erg=fdir d2 x=199.98 y=d23 z=d24 vec=1 0 0 tr=101 +sdef dir=d1 erg=fdir d2 x=199.5 y=d23 z=d24 vec=1 0 0 tr=101 par=n ara=165.0 wgt=6.436556E13 si1 -1 0 0.5 0.866 0.9659 0.9848 0.99619 0.9975640 0.998629 0.999390 0.999847 1 diff --git a/tallies/BBG_f4.i b/tallies/BBG_f4.i index 9c58d41..3348d74 100644 --- a/tallies/BBG_f4.i +++ b/tallies/BBG_f4.i @@ -1,6 +1,6 @@ c f4 tally sphere at BBG fc14 neutron flux at BBG -f14:n 701 +f14:n 100701 e14 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8