.. _app-xyzmesh-example: :ref:`4.13節 ` で使用するPHITS入力ファイル ---------------------------------------------------------------------------- :: [Title] Dose from 250MeV protons activating W/H2O/Fe block surrounded in Be [Parameters] icntl = 0 # (D=0) 3:ECH 5:NOR 6:SRC 7,8:GSH 11:DSH 12:DUMP maxcas = 1000 # (D=10) number of particles per one batch maxbch = 10 # (D=10) number of batches file(1) = c:/phits # (D=c:/phits) PHITS install folder name $ required options for DCHAIN file(21) = c:/phits/dchain-sp/data # (D=c:/phits/dchain-sp/data) DCHAIN data folder jmout = 1 # (D=0) Density echo, 0:input, 1:number density e-mode = 0 # (D=0) Event generator mode is not recommended for DCHAIN igamma = 3 # (D=2) 0:No, 1:Old, 2:EBITEM, 3:EBITEM+Isomer $ Photon source from DCHAIN output $ infl:{W_reg_target_gamma_t3.pht} $ This needs to come before the "off" source section or it will be skipped. [Source] s-type = 1 # mono-energetic axial source proj = proton # kind of incident particle e0 = 250.00 # energy of beam [MeV/u] r0 = 0.5000 # radius [cm] x0 = 0.0000 # (D=0.0) center position of x-axis [cm] y0 = 0.0000 # (D=0.0) center position of y-axis [cm] z0 = -10.000 # minimum position of z-axis [cm] z1 = -10.000 # maximum position of z-axis [cm] dir = 1.0000 # z-direction of beam [cosine] [Material] $ Tungsten, Density (g/cm3) = 19.300000 M1 74000 1.000000 $ W $ Water, Density (g/cm3) = 0.998207 ~= 1.0 M2 1000 0.666657 $ H 8000 0.333343 $ O $ Iron, Density (g/cm3) = 7.874000 M3 26000 1.000000 $ Fe $ Beryllium, Density (g/cm3) = 1.848000 M4 4000 1.000000 $ Be $ Air (dry, sea level), Density (g/cm3) = 0.001205 M10 6000 0.000150 $ C 7000 0.784431 $ N 8000 0.210748 $ O 18000 0.004671 $ Ar $ Concrete, Ordinary (NIST), Density (g/cm3) = 2.300 M11 1000 0.305330 $ H 6000 0.002880 $ C 8000 0.500407 $ O 11000 0.009212 $ Na 12000 0.000725 $ Mg 13000 0.010298 $ Al 14000 0.151042 $ Si 19000 0.003578 $ K 20000 0.014924 $ Ca 26000 0.001605 $ Fe set:c1[3] $ thickness of W portion set:c2[10] $ thickness of H2O portion set:c3[5] $ thickness of Fe portion set:c4[5] $ thickness of Be back end set:c5[4] $ half side length of inner rectangular prism set:c6[8] $ half side length of outer rectangular prism [Surface] $ Target components 1 rpp -c5 c5 -c5 c5 0 c1 $ W portion 2 rpp -c5 c5 -c5 c5 c1 c1+c2 $ H2O portion 3 rpp -c5 c5 -c5 c5 c1+c2 c1+c2+c3 $ Fe portion 4 rpp -c5 c5 -c5 c5 0 c1+c2+c3 $ Be inner border 5 rpp -c6 c6 -c6 c6 0 c1+c2+c3+c4 $ Be outer border $ Room structure 15 rpp -350 -300 -100 200 -100 150 $ inside shield 16 rpp -550 -500 -100 200 -100 50 $ doorway 20 rpp -500 100 -100 200 -100 300 $ inner wall 21 rpp -550 150 -150 250 -150 350 $ outer wall [Cell] $ Target 1 1 -19.3 -1 $ W 2 2 -1.0 -2 $ H20 3 3 -7.874 -3 $ Fe 4 4 -1.848 4 -5 $ Be $ Room (void) 10 0 -16:(5 -20 15) $ room + doorway air 11 0 -15:(20 -21 #10) $ outer + maze concrete walls $ Room (with materials) $ 10 10 -0.0012 -16:(5 -20 15) $ room + doorway air $ 11 11 -2.3 -15:(20 -21 #10) $ outer + maze concrete walls 99 -1 21 $ outer void / particle graveyard [Volume] $ required section for DCHAIN reg vol 1 c1*(2*c5)**2 2 c2*(2*c5)**2 3 c3*(2*c5)**2 4 ((c1+c2+c3+c4)*(2*c6)**2)-((c1+c2+c3)*(2*c5)**2) [Mat Name Color] mat name size color 0 void 1.00 lightgray 1 Tungsten 1.00 darkgreen 2 Water 1.00 pastelcyan 3 Iron 1.00 orange 4 Beryllium 1.00 pastelviolet 10 Air 1.00 {-0.9} 11 Concrete 1.00 gray $ Beam settings to be used in [T-Dchain] tallies set:c11[100.0] $ beam current (nA) set:c12[ c11 * 1.0e-9 / (1.602177e-19) ] $ beam power (source/sec) [T-DCHAIN] title = W target (reg) mesh = reg reg = 1 file = W_reg_target.out # file name of dchain input file timeevo = 2 # time evolution / irradiation schedule 10.0 m 1.0 50.0 m 0.0 outtime = 4 # output times 10.0 m -1 m -5 m -50 m amp = c12 # (D=1.0) Source Intensity(source/sec) idivs = 50 # (D=50) number of calculation substebs used in irradiation iphtout = 2 # (D=1) 2 = separate [Source] output to files for each time [T-DCHAIN] title = W target (xyz) mesh = xyz x-type = 2 nx = 8 xmin = -c5 xmax = c5 y-type = 2 ny = 8 ymin = -c5 ymax = c5 z-type = 2 nz = 3 zmin = 0 zmax = c1 file = W_xyz_target.out # file name of dchain input file timeevo = 2 # time evolution / irradiation schedule 10.0 m 1.0 50.0 m 0.0 outtime = 4 # output times 10.0 m -1 m -5 m -50 m amp = c12 # (D=1.0) Source Intensity(source/sec) idivs = 4 # (D=50) number of calculation substebs used in irradiation iphtout = 2 # (D=1) 2 = separate [Source] output to files for each time ipltmode = 1 # (D=0) enable DCHAIN 2D xy activity plot [T-DCHAIN] title = Whole target structure mesh = xyz x-type = 2 nx = 1 xmin = -c6 xmax = c6 y-type = 2 ny = 8 ymin = -c6 ymax = c6 z-type = 2 nz = 12 zmin = 0 zmax = c1+c2+c3+c4 file = whole_target_yz-view.out # file name of dchain input file timeevo = 2 # time evolution / irradiation schedule 10.0 m 1.0 50.0 m 0.0 outtime = 4 # output times 10.0 m -10 s -1 m -50 m amp = c12 # (D=1.0) Source Intensity(source/sec) idivs = 4 # (D=50) number of calculation substebs used in irradiation ipltmode = 4 # enable DCHAIN 2D yz activity plot [T-Track] title = proton, neutron, and photon distributions mesh = xyz x-type = 2 xmin = -10 xmax = 10 nx = 1 z-type = 2 zmin = -10 zmax = 30 nz = 200 y-type = 2 ymin = -20 ymax = 20 ny = 200 e-type = 2 emin = 0 emax = 2000 ne = 1 axis = yz file = yz-track.out part = proton neutron photon unit = 1 2D-type = 7 # 1:Cont, 2:Clust, 3:Color, 4:xyz, 5:mat, 6:Clust+Cont, 7:Col+Cont gshow = 3 # 0: no 1:bnd, 2:bnd+mat, 3:bnd+reg 4:bnd+lat epsout = 1 # automatically generate eps plot [T-Track] off title = Ambient dose equivalent H*(10) [mSv/hr] room map mesh = xyz x-type = 2 nx = 140 xmin = -550 xmax = 150 y-type = 2 ny = 1 ymin = -100 ymax = 100 z-type = 2 nz = 100 zmin = -150 zmax = 350 unit = 1 axis = xz e-type = 3 ne = 1 emin = 0.001 emax = 10 file = room-dose_reg-src.out set:c20[3600/1.0E+09] factor = c20 $ convert pSv/sec to mSv/hr multiplier = 6 $ number of regions using multiplier mat mset1 1 ( 0 -200 ) $ Zero out regions where we don't 2 ( 0 -200 ) $ care about the dose (inside walls 3 ( 0 -200 ) $ and the target). 4 ( 0 -200 ) 10 ( 1 -200 ) $ We only want to tally dose rate in air. 11 ( 0 -200 ) 2D-type = 7 # 1:Cont, 2:Clust, 3:Color, 4:xyz, 5:mat, 6:Clust+Cont, 7:Col+Cont gshow = 1 # 0: no 1:bnd, 2:bnd+mat, 3:bnd+reg 4:bnd+lat epsout = 1 $ automatically generate eps plot $ Set bounds of ANGEL plot color bar and make axis labels bigger angel = cmin(1.0E-5) cmax(1.0E+1) sangel = 2 x: {\Large z [cm]} y: {\Large x [cm]} [T-Gshow] mesh = xyz x-type = 2 nx = 280 xmin = -550 xmax = 150 y-type = 2 ny = 1 ymin = -5 ymax = 5 z-type = 2 nz = 200 zmin = -150 zmax = 350 axis = zx file = room_geometry_cross-section.out output = 2 epsout = 1 [T-Gshow] mesh = xyz x-type = 2 xmin = -20 xmax = 20 nx = 200 z-type = 2 zmin = -10 zmax = 30 nz = 200 y-type = 2 ymin = -10 ymax = 10 ny = 1 axis = zx file = target_geometry.out output = 2 epsout = 1 [End]