electron-tantalum-ErD3 c c electron beam incident on tantalum converter--> ErD3 --> checking thickness c c -------------------------- Input Information ----------------------------- 80 c c c c ------------------------ End Input Information --------------------------- 80 c c ----------------------------- Cell Cards --------------------------------- 80 1 8 -7.85 -102 $SS Table Top 2 6 -16.65 -301 $Tantalum target 3 1 -3.8 -401 $Erbium Deuteride (ErD3) cylinder (50% packing density) 4 0 -201 $Vacuum region for tallying particles entering in front of Ta 5 0 -202 $Vacuum region for tallying particles entering behind Ta + ErD3 6 0 (-203 204 600 -205) $Vacuum region for tallying particles entering cylindrical shell around all 7 0 -602 -600 601 $Vacuum region showing location of source 8 7 -0.001225 -101 #1 #2 #3 #4 #5 #6 #7 $Air inside everything, outside all other materials 9 0 101 $Particle Graveyard c --------------------------- End Cell Cards ------------------------------- 80 c ---------------------------- Surface Cards ------------------------------- 80 c ROOM AIR $Room air box: coordinates of corner followed by vectors 101 BOX -15 -25 -5 50 0 0 0 50 0 0 0 30 c Stainless Steel Table 102 BOX -10 -20 0 40 0 0 0 40 0 0 0 0.3175 $Steel table top (1/8" thick) c Vacuum zones for tallies 201 rcc -1.1 0 11 0.05 0 0 2 $region for tallying electrons entering in front of Ta 202 rcc 24.5 0 11 0.1 0 0 10.5 $region for tallying particles entering behind Ta + ErD3 203 c/x 0 11 10.5 $region for tallying particles in surrounding cylindrical shell 204 c/x 0 11 10.4 205 px 24.45 c also uses px from 600 and px above. small gap between this and back tally region c Tantalum Target A 301 rcc -0.9 0 11 0.15 0 0 2 $Tatalum target, th ~ 0.15 cm c 401 rcc -0.74 0 11 0.10 0 0 10 $All ErD3, th= cm, Full cylider c Cylindrical region with a Surface Source near the tantalum side 600 px -1.2 $Slightly larger radius as tantalum - volume just to show where source is located 601 px -1.5 602 c/x 0 11 2.01 c c -------------------------- End Surface Cards ----------------------------- 80 c ----------------------------- Data Cards --------------------------------- 80 c Materials c Material 1 "ErD3" c Constituents (atom %): m1 68162. 0.000346 plib=.12p 68164. 0.004003 plib=.12p 68166. 0.083757 plib=.12p 68167. 0.057174 plib=.12p 68168. 0.067444 plib=.12p 68170. 0.037276 plib=.12p 1002. 0.750000 $ plib=.12p $Rounded down so that sum = 1 even c 50% packing fraction ~ density around (-3.8 g/mL)--------------------------- c c Material 6: "Tantalum" (-16.65 g/cm^3) c Constituents (atom %): c Ta181-73181 (1) m6 73181 1.0 c c Material 7: "Air" (-0.001225 "g/cm^3") m7 6012. 0.00015 7014. 0.784431 8016. 0.210748 18040. 0.004671 plib=.12p c c Material 8: "Stainless Steel 304" (-7.85 "g/cm^3") c Constituents (atom %): c Fe(nat)-26000 (0.71) Cr(nat)-24000 (0.18) Ni(nat)-28000 (0.08) c Mn55-25055 (0.02) Si(nat)-14000 (0.01) m8 26000 0.71 24000 0.18 28000 0.08 25055 0.02 14000 0.01 plib=.12p c c c Mode mode n e d t h s p a c c Cell Importances (total cells in CELL CARD) imp:n 1 1 1 1 1 1 1 1 0 $9 per line imp:e 1 1 1 1 1 1 1 1 0 $ imp:d 1 1 1 1 1 1 1 1 0 $ imp:t 1 1 1 1 1 1 1 1 0 $ imp:h 1 1 1 1 1 1 1 1 0 $ imp:s 1 1 1 1 1 1 1 1 0 $ imp:p 1 1 1 1 1 1 1 1 0 $ imp:a 1 1 1 1 1 1 1 1 0 $ c c ---------------------------Phys Cards---------------------------- c format: phys:n emax emcnf iunr J J J coilf cutn J J J c i_int_model i_els_model c in phys:n here, emax = 100MeV, emcn & iunr = 0 coilf = 4 c PHYS:N 100 0 0 J J J 4 -1 J J J 0 0 PHYS:E 100 0 0 0 0 1 1 1 1 0 0 J J 0.917 0.001 $ default values right now PHYS:D 100 J J J 1 PHYS:t 100 J J J 1 phys:h 100 0 -1 j 0 j 1 j j j 0 0 0 0.917 j 0 0 phys:p 100 0 0 -1 0 j 1 c c ----------------------------------------------------------------- c ---------------------------- Sources ---------------------------- c ----------------------------------------------------------------- c Source Definition c Distribution of Energy neutrons: sdef SUR=600 POS=-1.2 0 11 VEC 1 0 0 RAD=d1 DIR=1 par=e erg=10 si1 H 0 2.0001 $Source radius - similar to the petri dish outer radius sp1 -21 1 $power law: proportional to radius for even areal dist. c c Histories (or Computer Time Cutoff) nps 1e8 c c ----------------------------------------------------------------- c ---------------------------- Tallies ---------------------------- c ----------------------------------------------------------------- c ----electrons fc104 tally "F104e-Vac1" f104:e ( 4 ) E104 0.1 198i 20 c Energy Bins for incoming/outgoing electrons in the vacuum space before the Tantalum zone c fc114 tally "F114e-Ta" f114:e ( 2 ) E114 0.1 198i 20 c Energy Bins for incoming/outgoing electrons in the tantalum c fc124 tally "F124e-ErD3" f124:e ( 3 ) E124 0.1 198i 20 c Energy Bins for incoming/outgoing electrons in the ErD3 zone c fc134 tally "F134e-Vac2" f134:e ( 5 ) E134 0.1 198i 20 c Energy Bins for incoming/outgoing electrons in the vacuum space before the Tantalum zone c fc144 tally "F144e-Vac3" f144:e ( 6 ) E144 0.1 198i 20 c Energy Bins for incoming/outgoing electrons in the vacuum space surrounding (cylindrical shell) c c c ----photons fc204 tally "F204p-Vac1" f204:p ( 4 ) E204 0.1 198i 20 c Energy Bins for incoming/outgoing photons in the vacuum space before the Tantalum zone c fc214 tally "F214p-Ta" f214:p ( 2 ) E214 0.1 198i 20 c Energy Bins for incoming/outgoing photons in the tantalum c fc224 tally "F224p-ErD3" f224:p ( 3 ) E224 0.1 198i 20 c Energy Bins for incoming/outgoing photons in the ErD3 zone c fc234 tally "F234p-Vac2" f234:p ( 5 ) E234 0.1 198i 20 c Energy Bins for incoming/outgoing photons in the vacuum space before the Tantalum zone c fc244 tally "F244p-Vac3" f244:p ( 6 ) E244 0.1 198i 20 c Energy Bins for incoming/outgoing photons in the vacuum space surrounding (cylindrical shell) c c c ----neutrons fc304 tally "F304n-Vac1" f304:n ( 4 ) E304 0.1 198i 20 c Energy Bins for incoming/outgoing neutrons in the vacuum space before the Tantalum zone c fc314 tally "F314n-Ta" f314:n ( 2 ) E314 0.1 198i 20 c Energy Bins for incoming/outgoing neutrons in the tantalum c fc324 tally "F324n-ErD3" f324:n ( 3 ) E324 0.1 198i 20 c Energy Bins for incoming/outgoing neutrons in the ErD3 zone c fc334 tally "F334n-Vac2" f334:n ( 5 ) E334 0.1 198i 20 c Energy Bins for incoming/outgoing neutrons in the vacuum space before the Tantalum zone c fc344 tally "F344e-Vac3" f344:n ( 6 ) E344 0.1 198i 20 c Energy Bins for incoming/outgoing photons in the vacuum space surrounding (cylindrical shell) c c c c c --------------------------- CUT Cards ------------------------------------ 80 c Below: need to specifically define lower energy cutoffs instead of defaults cut:n j 0 c cut:e j 1e-5 $ testing lost particle errors w/ new geometry cut:d j 1e-3 cut:t j 1e-3 cut:h j 1e-3 $ These are the lower limits of MCNP cut:s j 1e-3 $ This sets the lower energy cutoff to 1keV. cut:p j 1e-3 $ c c -------------------------- Trying DBCN Card------------------------------- 80 DBCN 27J 200000 $Debug card. x28 deals with bank size. c max value of DBCN is 8-byte integer, so up to >18E18? c -------------------------- Trying DBCN Card------------------------------- 80 c c L'Ecuyer 63-bit random number generator (period=9.2E18) rand gen=2 c print -85 -86 -87 -98 c c --------------------------- End Data Cards ------------------------------- 80 c End MCNP Input