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

Include scaled WW and run test against 5f at BBG, improve WWG ref

Browse Source
This commit is contained in:
glavic_a
2019-06-02 14:21:12 +02:00
parent 4cefa58ab3
commit b9ee1b5210
7 changed files with 37 additions and 26 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
e2-estia.i
View File

@@ -94,18 +94,25 @@ 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
WWGE:p 1.667 3e3
c mesh based weight windows, can be reused after change of cell structure
WWG 514 0 0
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 0 13.4 AXS=0 0 1
IMESH=200 550 1150 1500 2112.01 2157.17 4600
IINTS= 2 5 1 10 1 1 4
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 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
c source weight is 1e16, give large window with low split/roulette factor
c check only on collisions
c WWE:h,/,d,t,s,a 3e3
c WWP:h,/,d,t,s,a 100 10 2 -1 0
c WWN1:h,/,d,t,s,a 1e15 276r -1 -1 1e10 578r
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
c -----------------------------------------------------------
c ------------------- TALLY CARDS ---------------------------
@@ -117,8 +124,9 @@ 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/NBPI_f4.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
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
@@ -161,7 +169,7 @@ c -------------------------------------------------------
c --------------- PHYSICS CARDS --------------------------
c -------------------------------------------------------
c dbcn 375642321 j j j j j j j j j j j 15291711
nps 1e7
nps 1e4
c histp
mode n p h /
d $ deuteron

27
geometry/materials.i
View File

@@ -13,7 +13,7 @@ c Material : Stainless304 rho=0.0876276
m3 6000 0.00031864 14028 0.00170336 15031 6.95038e-05
16032 4.4752e-05 24000.50c 0.0174813 25055 0.00174159
26054 0.00338 26056 0.053455 26057 0.00128218
28000.50c 0.00815128 pnlib=70u
28000.50c 0.00815128
mx3:h model j j model 2r j j j model
c
C NIST STANDARD DRY AIR rho=1.2E-3
@@ -26,7 +26,7 @@ m5 13027 0.054381 12024 0.00207 12025 0.0002621
12026 0.0002885 14028 0.000217 22048 8.7e-05
24052.24c 7.8e-05 24053.24c 1e-05 25055 0.000408
26056 0.000214 26057 5e-06 29063 4e-05
29065 2e-05 30000 0.000145 pnlib=70u
29065 2e-05 30000 0.000145
mt5 al27.12t
mx5:h j model 2r j model j j model 4j model
c
@@ -46,7 +46,8 @@ c 1001 -0.0284
c mx8:h j j model j j
c mt8 lwtr.01t
c Helium
m9 2004 2.45e-05 pnlib=70u
m9 2004 2.45e-05
mx9:h 2004.71h
c
c BorAL 1100 alloy (2.7 g/cc)
m10 5010 -0.04 5011 -0.01 13027 -0.9288
@@ -56,7 +57,7 @@ mx10:h j model j j model 3r
c
c Material : H2O rho=0.100283
c m11 1001 0.0668456 1002 1.003e-05 8016 0.0334278
c pnlib=70u
c
c mt11 lwtr.01t
c
c bor-float glass 13% B2O3 2.51 g/cm3
@@ -121,14 +122,12 @@ m49 1001 0.00776555 1002 1.16501e-06 8016 0.0438499
14028 0.0158026 16032 5.63433e-05 19000.60c 0.000693104
20000.60c 0.00291501 26054 1.84504e-05 26056 0.000286826
26057 6.5671e-06 26058 8.75613e-07
pnlib=70u
mt49 lwtr.10t
mx49:h 3j model 1r j j model 2r 3j model
c
c Material : CastIron rho=0.0833854
m54 6000 0.00636 14028 0.0051404 25055 0.000386
26054 0.00419 26056 0.065789 26057 0.00152
pnlib=70u
mx54:h model j model j j j
c
c Material : Steel71 rho=0.0843223
@@ -137,34 +136,34 @@ m71 6000 0.000101187 14028 0.000295128 24000.50c 0.000168645
26057 0.00175516 26058 0.00023358 28058 0.000172218
28060 6.6335e-05 28061 2.88126e-06 28062 9.19272e-06
28064 2.33992e-06 29063
0.000116651 29065 5.19931e-05 pnlib=70u
0.000116651 29065 5.19931e-05
mx71:h model j model 1r j j j model j j j j j j j
c
c Material : Copper rho=0.0844385
c (rho=8.91001 g/cc)
m73 29063 0.0583892 29065 0.0260493 pnlib=70u
m73 29063 0.0583892 29065 0.0260493
c
c Material : Be5H2O rho=0.122325
c m118 1001 0.00334228 4009 0.117311 8016 0.00167114
c pnlib=70u
c
c mx118:h j model j
c mt118 be.60t lwtr.01t
c c
c c Material : Be10H2O rho=0.121164
c m119 1001 0.00668456 4009 0.111137 8016 0.00334278
c pnlib=70u
c
c c
c mx119:h j model j
c mt119 be.60t lwtr.01t
c
c Material : Tungsten_15.1g rho=0.0494621
c m120 74182 0.0131232 74183 0.00708653 74184 0.0151734
c 74186 0.014079 pnlib=70u
c 74186 0.014079
c
c Material : Iron_10H2O rho=0.0862739
c m121 1001 0.00668924 8016 0.00334462 26054 0.00445623
c 26056 0.0699533 26057 0.00161553 26058 0.000214997
c pnlib=70u
c
c mt121 lwtr.01t
c mx121:h 5j model
c
@@ -175,12 +174,12 @@ m122 1001 0.0229186 1002 1.23628e-06 5010 0.0001187
20040 0.0062902 20042 4.1981e-05 20043 8.759e-06
20044 0.0001353 20046 2.5956e-07 20048 1.21339e-05
26054 1.1341e-05 26056 0.00017803 26057 4.1115e-06
26058 5.47166e-07 pnlib=70u
26058 5.47166e-07
mx122:h j j j model model 5j model 4r 3j model
mt122 poly.10t
c
c Material : ParaOrtho%0.5 rho=0.041975
c m123 1001 0.000209875 1005.60c 0.0417651 pnlib=70u
c m123 1001 0.000209875 1005.60c 0.0417651
c mx123:h j 1001.70h
c mt123 hpara.10t hortho.10t
c

5
geometry/materials_monolith.i
View File

@@ -1286,7 +1286,7 @@ M6141 01001 0.302150460 $ SS316L + 41%VSMOW
42098 0.001921387
42100 0.000768078
C ==============================================================================
M7410 74180.50C 0.001200000 $ W
M7410 74180.80C 0.001200000 $ W
74182 0.265000000
74183 0.143100000
74184 0.306400000
@@ -1406,7 +1406,8 @@ C ======================================== REPLACEMENT NUCLIDES AND/OR MODELS ==
C ==============================================================================
MX0400:H MODEL
C ==============================================================================
MX0650:H 4J
MX0650:H 3J
MODEL
06012
C ==============================================================================
MX0870:H 06012

BIN
results/e2-estia_newmono_v02-0-g4cefa58.out.gz Normal file
View File

Binary file not shown.

8
sources/proton.i
View File

@@ -1,12 +1,14 @@
c Original definition of source from proton beam
c cut:n j 0.99e-06
cut:n j 0
cut:p 1e+08 0.5
cut:p 1e+08 1.667
c weight windows from external file generated by WWG
c make sure to supply the correct WWINP file
c WWP:n,p j j 10 j -1 j 8.7374e15 $ WW was generated for source scaling 1
WWP:n,p j j 2 j -1 j 1.0e2 1j 1e17 1j
c energy splitting for higher neutron energies and roulette below 1keV
esplt:n 0.05 1e-6 0.2 1e-3 1 10 3 100 10 3000
c esplt:n 0.05 1e-6 0.2 1e-3 1 10 3 100 10 3000
c energy splitting for higher neutron energies no roulette
c esplt:n 1 10 3 100 10 3000
c the source definition with protons/s
sdef dir=1 erg=2000 par=9 vec=0.0 1 0.0 x=d1 y=-30 z=d2 wgt=1.560377E16
si1 A -8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0

1
tallies/BBG_f5.i
View File

@@ -16,3 +16,4 @@ e15 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
1e+3 1.58e+3 2.51e+3
fq15 e f
pd15 0.01 276r 0 0 0.1 578r

2
weight_windows/current.wwg
View File

@@ -1 +1 @@
none.wwg
e2-estia_newmono_v02-0-g4cefa58.wwg