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

Run closed shutter mesh with MF derived source

Browse Source
This commit is contained in:
glavic_a
2019-05-04 17:18:03 +02:00
parent d6b178d599
commit 4486b279d6
4 changed files with 141 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
e2-estia.i
View File

@@ -107,7 +107,7 @@ c read file=tallies/bunker_mesh_p.i echo
c read file=tallies/bunker_mesh_np.i echo
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
read file=tallies/full_mesh_p.i echo
c read file=tallies/full_mesh_np.i echo
c read file=tallies/inbunker_f4.i echo
c read file=tallies/chopper_pit_doses.i echo
@@ -115,7 +115,7 @@ c read file=tallies/bunker_wall_f4.i echo
c read file=tallies/e1_e2_f4.i echo
c read file=tallies/e1_e2_f5.i echo
c read file=tallies/shutter_f4.i echo
read file=tallies/shutter_f5.i echo
c read file=tallies/shutter_f5.i echo
c read file=tallies/shutter_sample_f5.i echo
c read file=tallies/selene1_walls_f4.i echo
c read file=tallies/optical_cave_f5.i echo
@@ -147,7 +147,7 @@ c -------------------------------------------------------
c --------------- PHYSICS CARDS --------------------------
c -------------------------------------------------------
c dbcn 375642321 j j j j j j j j j j j 15291711
nps 1e10
nps 2e8
c histp
mode n p n h / d t s a
c volume calculated for tally cells (spheres)
@@ -168,7 +168,7 @@ 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
c * Source Term from Valentina for 2m adapted to size
read file=sources/e2_2m.i echo
c 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
@@ -176,6 +176,7 @@ c read file=sources/mcstas_source_sample.i echo
c read file=sources/mcstas_source_mf.i echo
c read file=sources/mcstas_gammas.i echo
c read file=sources/mcstas_gammas_2.i echo
read file=sources/mcstas_source_shutter.i echo
prdmp 1e8 -60 1 2 1e7
print 30 32 40 $ 120 128 150 160 161 162
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++

1
geometry/estia_cave.i
View File

@@ -39,3 +39,4 @@ c area around cave
c visualize Selene gamma sources (comment for real runs)
c 50001 0 -50001
c 50002 0 -50002
c 50003 0 -50003

2
geometry/estia_surfaces.i
View File

@@ -957,4 +957,6 @@ c
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 visualize Shutter Source
c 50003 3 rpp -170.5 -169.5 0 4 -6 6
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++

133
sources/mcstas_source_shutter.i Normal file
View File

@@ -0,0 +1,133 @@
c cold neutron source from McStas simulation for middle focus position
c activation (delayed particles)
c act nonfiss=p,a,e dg=lines
cut:n j 1e-11 0 0
cut:p j 0.01
c as this is primary for gamma shielding purposes,
c no directional dependence has been considered
sdef dir=1 erg=d11 x=-170. y=d23 z=d24 vec=1 0 0 tr=3
par=n ara=72.0 wgt=3.33085e+09
c scaling wgt is to flux from McStas nps
si1 -1 0.99976203 1
sp1 0 0 1
si11 0.00000e+00
9.009970e-11 9.070338e-11 9.131316e-11 9.192910e-11 9.255130e-11
9.317984e-11 9.381480e-11 9.445628e-11 9.510435e-11 9.575912e-11
9.642068e-11 9.708911e-11 9.776452e-11 9.844700e-11 9.913665e-11
9.983357e-11 1.005379e-10 1.012497e-10 1.019690e-10 1.026961e-10
1.034309e-10 1.041737e-10 1.049245e-10 1.056835e-10 1.064507e-10
1.072263e-10 1.080104e-10 1.088032e-10 1.096047e-10 1.104151e-10
1.112345e-10 1.120630e-10 1.129009e-10 1.137482e-10 1.146051e-10
1.154716e-10 1.163481e-10 1.172346e-10 1.181312e-10 1.190382e-10
1.199556e-10 1.208837e-10 1.218226e-10 1.227725e-10 1.237336e-10
1.247059e-10 1.256898e-10 1.266854e-10 1.276929e-10 1.287124e-10
1.297441e-10 1.307884e-10 1.318453e-10 1.329150e-10 1.339978e-10
1.350939e-10 1.362035e-10 1.373269e-10 1.384642e-10 1.396157e-10
1.407816e-10 1.419621e-10 1.431576e-10 1.443683e-10 1.455943e-10
1.468361e-10 1.480938e-10 1.493677e-10 1.506582e-10 1.519654e-10
1.532897e-10 1.546315e-10 1.559909e-10 1.573683e-10 1.587640e-10
1.601784e-10 1.616118e-10 1.630645e-10 1.645369e-10 1.660293e-10
1.675421e-10 1.690757e-10 1.706304e-10 1.722067e-10 1.738049e-10
1.754255e-10 1.770688e-10 1.787354e-10 1.804256e-10 1.821399e-10
1.838787e-10 1.856425e-10 1.874319e-10 1.892472e-10 1.910890e-10
1.929579e-10 1.948543e-10 1.967788e-10 1.987320e-10 2.007144e-10
2.027266e-10 2.047692e-10 2.068428e-10 2.089481e-10 2.110857e-10
2.132563e-10 2.154605e-10 2.176991e-10 2.199727e-10 2.222822e-10
2.246282e-10 2.270115e-10 2.294330e-10 2.318934e-10 2.343937e-10
2.369345e-10 2.395170e-10 2.421418e-10 2.448101e-10 2.475227e-10
2.502806e-10 2.530849e-10 2.559366e-10 2.588368e-10 2.617865e-10
2.647870e-10 2.678393e-10 2.709447e-10 2.741045e-10 2.773198e-10
2.805921e-10 2.839226e-10 2.873127e-10 2.907640e-10 2.942778e-10
2.978556e-10 3.014992e-10 3.052100e-10 3.089897e-10 3.128401e-10
3.167629e-10 3.207600e-10 3.248331e-10 3.289844e-10 3.332158e-10
3.375293e-10 3.419271e-10 3.464115e-10 3.509846e-10 3.556489e-10
3.604068e-10 3.652608e-10 3.702136e-10 3.752677e-10 3.804261e-10
3.856916e-10 3.910671e-10 3.965558e-10 4.021609e-10 4.078857e-10
4.137336e-10 4.197082e-10 4.258131e-10 4.320522e-10 4.384294e-10
4.449489e-10 4.516148e-10 4.584317e-10 4.654042e-10 4.725369e-10
4.798348e-10 4.873031e-10 4.949472e-10 5.027725e-10 5.107849e-10
5.189904e-10 5.273952e-10 5.360058e-10 5.448290e-10 5.538719e-10
5.631419e-10 5.726465e-10 5.823937e-10 5.923920e-10 6.026500e-10
6.131767e-10 6.239817e-10 6.350749e-10 6.464664e-10 6.581673e-10
6.701887e-10 6.825426e-10 6.952412e-10 7.082975e-10 7.217250e-10
7.355381e-10 7.497515e-10 7.643809e-10 7.794428e-10 7.949542e-10
8.109333e-10 8.273991e-10 8.443715e-10 8.618716e-10 8.799214e-10
8.985442e-10 9.177645e-10 9.376082e-10 9.581026e-10 9.792762e-10
1.001160e-09 1.023785e-09 1.047186e-09 1.071398e-09 1.096460e-09
1.122412e-09 1.149296e-09 1.177158e-09 1.206046e-09 1.236010e-09
1.267104e-09 1.299387e-09 1.332920e-09 1.367768e-09 1.404000e-09
1.441691e-09 1.480921e-09 1.521774e-09 1.564341e-09 1.608719e-09
1.655013e-09 1.703335e-09 1.753804e-09 1.806549e-09 1.861711e-09
1.919438e-09 1.979892e-09 2.043248e-09 2.109694e-09 2.179435e-09
2.252693e-09 2.329707e-09 2.410739e-09 2.496073e-09 2.586020e-09
2.680917e-09 2.781136e-09 2.887082e-09 2.999198e-09 3.117975e-09
3.243949e-09 3.377716e-09 3.519930e-09 3.671319e-09 3.832688e-09
4.004935e-09 4.189061e-09 4.386182e-09 4.597551e-09 4.824575e-09
5.068840e-09 5.332136e-09 5.616492e-09 5.924217e-09 6.257943e-09
6.620683e-09 7.015902e-09 7.447595e-09 7.920393e-09 8.439684e-09
9.011772e-09 9.644063e-09 1.034531e-08 1.112592e-08 1.199833e-08
1.297753e-08 1.408168e-08 1.533297e-08 1.675878e-08 1.839312e-08
2.027874e-08 2.246991e-08 2.503641e-08 2.806911e-08
sp11 0.00000e+00
1.214840e+04 1.256345e+04 1.302044e+04 1.326152e+04 1.379966e+04
1.421472e+04 1.475672e+04 1.528101e+04 1.583387e+04 1.642523e+04
1.689077e+04 1.734300e+04 1.778793e+04 1.864549e+04 1.934852e+04
2.014583e+04 2.083637e+04 2.110292e+04 2.219619e+04 2.249020e+04
2.357387e+04 2.393510e+04 2.534934e+04 2.582712e+04 2.715238e+04
2.786449e+04 2.887503e+04 2.949026e+04 3.083067e+04 3.151190e+04
3.245292e+04 3.405610e+04 3.445021e+04 3.619378e+04 3.735710e+04
3.828247e+04 3.996155e+04 4.109754e+04 4.278032e+04 4.438702e+04
4.546359e+04 4.713878e+04 4.830520e+04 4.981591e+04 5.213248e+04
5.434103e+04 5.570910e+04 5.739441e+04 5.832051e+04 6.047995e+04
6.407619e+04 6.473340e+04 6.728596e+04 6.979223e+04 7.170640e+04
7.417202e+04 7.763219e+04 7.943879e+04 8.224429e+04 8.616590e+04
8.848352e+04 9.056017e+04 9.319461e+04 9.511508e+04 9.947364e+04
1.036847e+05 1.084932e+05 1.116687e+05 1.146604e+05 1.193212e+05
1.225996e+05 1.269869e+05 1.324391e+05 1.357706e+05 1.401152e+05
1.464706e+05 1.468557e+05 1.539999e+05 1.613889e+05 1.656908e+05
1.693471e+05 1.760142e+05 1.812553e+05 1.873262e+05 1.984413e+05
2.015337e+05 2.106013e+05 2.119435e+05 2.196635e+05 2.285301e+05
2.368854e+05 2.422217e+05 2.543786e+05 2.633250e+05 2.706326e+05
2.750603e+05 2.910239e+05 3.045884e+05 3.130142e+05 3.188584e+05
3.285789e+05 3.399626e+05 3.507469e+05 3.637172e+05 3.741806e+05
3.904184e+05 4.008069e+05 4.254911e+05 4.293372e+05 4.523024e+05
4.622288e+05 4.807719e+05 4.938701e+05 5.078877e+05 5.231608e+05
5.412427e+05 5.589387e+05 5.817409e+05 5.971860e+05 6.122883e+05
6.425329e+05 6.618982e+05 6.809655e+05 7.039152e+05 7.354101e+05
7.549776e+05 7.766671e+05 8.026105e+05 8.429900e+05 8.740523e+05
8.908045e+05 9.191396e+05 9.662563e+05 1.001041e+06 1.030041e+06
1.060946e+06 1.089684e+06 1.150674e+06 1.164706e+06 1.200898e+06
1.279452e+06 1.271447e+06 1.333980e+06 1.390647e+06 1.440428e+06
1.471303e+06 1.539177e+06 1.586709e+06 1.634066e+06 1.695469e+06
1.741327e+06 1.806084e+06 1.869876e+06 1.933938e+06 1.979604e+06
2.093888e+06 2.166930e+06 2.265165e+06 2.294742e+06 2.354508e+06
2.464106e+06 2.520649e+06 2.670378e+06 2.722620e+06 2.768233e+06
2.867300e+06 3.003317e+06 3.028797e+06 3.179778e+06 3.309525e+06
3.465320e+06 3.576304e+06 3.655310e+06 3.846793e+06 3.923541e+06
4.092392e+06 4.248305e+06 4.518226e+06 4.487231e+06 4.671564e+06
4.996738e+06 5.108565e+06 5.315782e+06 5.570328e+06 5.744651e+06
5.872509e+06 5.987757e+06 6.341752e+06 6.470877e+06 6.772215e+06
6.946503e+06 7.278125e+06 7.640013e+06 7.714839e+06 8.039751e+06
8.494502e+06 8.869999e+06 9.033183e+06 9.488911e+06 9.766928e+06
1.022764e+07 1.049437e+07 1.083413e+07 1.103215e+07 1.157572e+07
1.212642e+07 1.268534e+07 1.294905e+07 1.330881e+07 1.407757e+07
1.441326e+07 1.519522e+07 1.535392e+07 1.634973e+07 1.694721e+07
1.747633e+07 1.814262e+07 1.910661e+07 1.957293e+07 2.046833e+07
2.118353e+07 2.226533e+07 2.303698e+07 2.394613e+07 2.531498e+07
2.605995e+07 2.715415e+07 2.846069e+07 2.953247e+07 3.073770e+07
3.225625e+07 3.393882e+07 3.470952e+07 3.674492e+07 3.897454e+07
4.005168e+07 4.103969e+07 4.364481e+07 4.538396e+07 4.664388e+07
4.986418e+07 5.144487e+07 5.514519e+07 5.669231e+07 5.826391e+07
6.109010e+07 6.432554e+07 6.770208e+07 7.007861e+07 7.545150e+07
7.780215e+07 8.083321e+07 8.390743e+07 8.773106e+07 9.210111e+07
9.338437e+07 9.794151e+07 1.017470e+08 1.055494e+08 1.079075e+08
1.089422e+08 1.128917e+08 1.072069e+08 9.669178e+07 7.435091e+07
4.410439e+07 1.706532e+07 4.223000e+06 6.719128e+05 5.827995e+04
2.447921e+03 4.470446e+01 6.579453e-01 1.203861e-03 1.592778e-06
8.984822e-10 3.613184e-14 2.100864e-19 3.859123e-29 0.000000e+00
0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
c approximate beam size from x and y projections
si23 0 4
sp23 0 1
si24 -6 6
sp24 0 1