253 lines
11 KiB
OpenEdge ABL
253 lines
11 KiB
OpenEdge ABL
Input File ESTIA
|
|
c -------------------------------------------------------
|
|
c --------------- CELL CARDS --------------------------
|
|
c -------------------------------------------------------
|
|
read file=geometry/ess_monolith_cells.i echo
|
|
read file=geometry/ess_outside_cells.i echo
|
|
c tallies
|
|
100701 0 -100701:-101701
|
|
100702 0 -100702
|
|
100703 0 -100703
|
|
100704 0 -100704
|
|
100705 0 -100705
|
|
100706 0 -100706
|
|
100707 0 -100707
|
|
100708 0 -100708
|
|
100709 0 -100709
|
|
100710 0 -100710
|
|
read file=geometry/estia_inbunker.i echo
|
|
read file=geometry/estia_guides.i echo
|
|
read file=geometry/estia_cave.i echo
|
|
read file=geometry/estia_universes.i echo
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
|
|
c -------------------------------------------------------
|
|
c --------------- SURFACE CARDS -------------------------
|
|
c -------------------------------------------------------
|
|
c modified ESS geometry
|
|
read file=geometry/ess_monolith_surfaces.i echo
|
|
read file=geometry/ess_outside_surfaces.i echo
|
|
c tallies
|
|
c 700 s 117.0 -155.00 13.5 3.0
|
|
100701 1 s 542.0 -6.3 6.3 3.0
|
|
101701 1 s 542.0 -6.3 -6.3 3.0
|
|
100702 1 s 1010 -0.8 0 5.0
|
|
100703 1 s 1082 -25.0 0 5.0
|
|
100704 1 s 1082 -0.25 0 5.0
|
|
100705 1 s 1165 0.7 0 5.0
|
|
100706 21 sx -22 12.0
|
|
100707 s 2177 -2.5 -13 12.0
|
|
100708 s 1548 191 0 30.0
|
|
100709 s 1539 -165 0 30.0
|
|
100710 s 2177 -91 -13 12.0
|
|
c Estia geometry
|
|
read file=geometry/estia_surfaces.i echo
|
|
|
|
read file=geometry/materials.i echo
|
|
read file=geometry/materials_monolith.i echo
|
|
c
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
c -------------------------------------------------------
|
|
c --------------- TRANSFORM CARDS -----------------------
|
|
c Estia ISCS with origin at first focal point of feeder
|
|
TR1 15.05639956 -2.13420532 13.7 0.99992537 0 -0.012217
|
|
0 1 0
|
|
0.012217 0 0.99992537
|
|
1
|
|
C center of feeder vessel with x-axis parallel to cylinder
|
|
TR2 775.66787865 -7.64605191 4.3862 0.99992537 0 -0.012217
|
|
0 1 0
|
|
0.012217 0 0.99992537
|
|
1
|
|
c Beam parallel x-axis and origin at middle focus
|
|
TR3 2314.8848 -2.1342 -14.3991 0.99992537 0 -0.012217
|
|
0 1 0
|
|
0.012217 0 0.99992537
|
|
1
|
|
c Beam parallel x-axis and origin at virtual source
|
|
TR4 1114.9743 -2.1342 0.261 0.99992537 0 -0.012217
|
|
0 1 0
|
|
0.012217 0 0.99992537
|
|
1
|
|
c Parallel to c-axis but horizontal and centered at MF
|
|
TR13 2314.8848 -2.1342 -14.3991 1 0 0 0 1 0 0 0 1 1
|
|
c Parallel to c-axis but horizontal and centered at VS
|
|
TR14 1114.9743 -2.1342 0.261 1 0 0 0 1 0 0 0 1 1
|
|
c Parallel to c-axis but horizontal and centered at sample
|
|
TR15 3514.7959 -2.13570 -29.06 1 0 0 0 1 0 0 0 1 1
|
|
c Flange axis for Selene 1
|
|
TR21 1334.9578 0.36581 -2.4264 0.99992537 0 -0.012217
|
|
0 1 0
|
|
0.012217 0 0.99992537
|
|
c Coordinate system of the Estia beamport focal point
|
|
TR100 9.57163854 -4.09313273 13.7 0.99997563 -0.00698126 0
|
|
0.00698126 0.99997563 0
|
|
0 0 1 1
|
|
c Same coordinates but with origin at TCS
|
|
TR101 0 0 0 0.99997563 -0.00698126 0
|
|
0.00698126 0.99997563 0
|
|
0 0 1 1
|
|
c Transformation from SKADI model coordinate system here
|
|
TR102 -4.06550381 -2.91061485 1.0 0.99376792 -0.11146893 0
|
|
0.11146893 0.99376792 0
|
|
0 0 1 1
|
|
c Coordinate for surface source within bunker wall
|
|
c close to plane 206305 with same surface normal and center around beam
|
|
TR103 1202.6185 0.05148 -0.86 0.9999754884138854 -0.00698105188428727 0.0
|
|
0.00698105188428727 0.9999754884138854 0.0
|
|
0 0 1 1
|
|
c Coordingates for instrument boundary meshes, y-direction points outside
|
|
c SKADI side wall starting at bunker boundary
|
|
c TR121 962.0 -1155.0 -150.0 0.675590 -0.737277 0.000000
|
|
c 0.737277 0.675590 0.000000
|
|
c 0 0 1
|
|
c c Fence on E01 side and left cave side
|
|
c TR122 1089.0 -2723.0 -180.0 0.582123 -0.813101 0.000000
|
|
c -0.813101 -0.582123 0.000000
|
|
c 0 0 1
|
|
c c Downstream cave side
|
|
c TR123 2841.0 -3360.0 -180.0 -0.813101 -0.582123 0.000000
|
|
c 0.582123 -0.813101 0.000000
|
|
c 0 0 1
|
|
c c Dillitation joint
|
|
c TR124 1104.0 -1820.0 -180.0 0.000000 -1.000000 0.000000
|
|
c -1.000000 0.000000 0.000000
|
|
c 0 0 1
|
|
read file=geometry/ess_monolith_transformations.i echo
|
|
c -------------------------------------------------------
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
c -------------------------------------------------------
|
|
c --------------- WEIGHT CARDS --------------------------
|
|
c -------------------------------------------------------
|
|
c weight window generator, can be run with previous weight window files
|
|
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
|
|
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 255 61003 0
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
c -----------------------------------------------------------
|
|
c ------------------- TALLY CARDS ---------------------------
|
|
c
|
|
c read file=tallies/bunker_mesh.i echo
|
|
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_np.i echo
|
|
c read file=tallies/boundary_meshes.i echo
|
|
c read file=tallies/boundary_wall_meshes.i echo
|
|
c read file=tallies/selene_source_mesh_p.i echo
|
|
c read file=tallies/NBPI_f4.i echo
|
|
c 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_feedthrough_f1.i echo
|
|
c read file=tallies/bunker_wall_f4.i echo
|
|
c read file=tallies/bunker_wall_f5.i echo
|
|
c read file=tallies/pre_selene1_mesh.i
|
|
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
|
|
c read file=tallies/shutter_f5.i echo
|
|
c read file=tallies/shutter_wall_f4.i echo
|
|
c read file=tallies/detector_f5.i echo
|
|
c read file=tallies/selene1_walls_f4.i echo
|
|
c read file=tallies/optical_cave_f5.i echo
|
|
c read file=tallies/optical_cave_gammas.i echo
|
|
c read file=tallies/experimental_cave_f5.i echo
|
|
c read file=tallies/cave_mesh.i echo
|
|
c read file=tallies/selene1_f5.i echo
|
|
c read file=tallies/selene1_f6.i echo
|
|
c
|
|
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 17r 0.25v3 7r $ BWC
|
|
0 155r
|
|
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
|
|
0 -0.05 -0.1 4r $ NBOA
|
|
0 23r -0.1 $ BBG
|
|
0 13r -0.25 6r $ NFGA
|
|
0 173r -1 $ Glass in Selene 1
|
|
0 51r -1 -1 $ S1 back plate and mineral cast in universe
|
|
0 82r
|
|
c
|
|
c -----------------------------------------------------------
|
|
c -------------------------------------------------------
|
|
c --------------- PHYSICS CARDS --------------------------
|
|
c -------------------------------------------------------
|
|
c dbcn 375642321 j j j j j j j j j j j 15291711
|
|
c stop calculation either at XeY histories, 6.5 days or tally 15 error <0.001
|
|
stop nps 1e6 CTME 4e5 $f31 0.01
|
|
rand gen=2 seed=29976355590901 $ default 19073486328125
|
|
c histp
|
|
mode n p h / z $ neutron, photon, proton, pi+, pi0
|
|
imp:n,p,h,/,z 1 276r 0 0 1 672r
|
|
c imp:n,p,h,/,z 1 276r 0 0 1 487r 0 0 0 1 0 87r 1 9r
|
|
cut:h,/,z 1e+08 1e-3
|
|
phys:n 3000 0
|
|
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
|
|
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 * Estimated source within bunker wall to optimize collimator efficiency
|
|
c read file=sources/bunker_wall.i echo $ approximate, not for official results
|
|
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
|
|
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
|
|
c read file=sources/mcstas_source_shutter.i echo
|
|
prdmp 5e6 -30 1 2 5e5
|
|
print 30 32 40
|
|
c 120 $O importance function quality
|
|
c 128 $O universe map
|
|
c 150 $O DXTRAN diagnostics
|
|
c 160 $D TFC bin tally analysis
|
|
c 161 $D f(x) tally density plot
|
|
c 162 $D cumulative f(x) and tally density plot
|
|
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
|
|