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

Test run with old 2m source term and no WW

Browse Source
This commit is contained in:
glavic_a
2019-06-03 19:09:08 +02:00
parent f25c87dcb0
commit bd6fa655eb
5 changed files with 28 additions and 24 deletions
Download Patch File Download Diff File Expand all files Collapse all files

28
e2-estia.i
View File

@@ -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

8
geometry/estia_cave.i
View File

@@ -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

8
geometry/estia_surfaces.i
View File

@@ -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 ++++++++++++++++++++++++++++

6
sources/e2_2m.i
View File

@@ -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

2
tallies/BBG_f4.i
View File

@@ -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