4.14. 例題：PHITSを使用しないDCHAIN単体計算

DCHAINは主としてPHITSの[T-Dchain]と組み合わせて使用するが、PHITSの[T-Track]・[T-Yield]出力を必要としない単純な崩壊計算であれば、DCHAIN単体でも実行できる。この場合は、DCHAINの既定値からいくつかのパラメータを変更し、DCHAINデータフォルダーのパスを1行で記載したdch_link.datを計算フォルダーに置く必要がある。既定の内容を以下に示す。

c:/phits/dchain-sp/data/

以下は、\(^{183}\)Hf、\(^{160m}\)Tm、\(^{169}\)Ta、\(^{164}\)Hfをそれぞれ初期放射能1 MBqとし、10日間の崩壊を複数の出力時刻で評価する例である。単体崩壊計算で通常必要となる主要なパラメータだけを示しており、必要に応じて他の出力パラメータも追加できる。変更時に重要な値は入力内のコメントで説明している。

原則として、jmode= -1、inmtcf= 0、hnfluxs= noの設定が必要である。ただし、jmode= -1を指定すれば、DCHAINが他の2つを強制的に上書きするため、実際にはjmodeだけでも単体崩壊計算を実行できる。

htitle = Standalone DCHAIN decay calculation example

  imode =       2        ! output inventory, decay heats, and photon spectra
  jmode =      -1        ! decay only, no external source
idcylib =       5        ! default decay library using data from JENDL and ENDF

  acmin =  1.0000E-20    ! threshold activity [Bq] to appear in *.act output file
 istabl =       0        ! toggle output of stable nuclides
iwrtchn =       1        ! decay chain scheme settings for *.dcs file
chrlvth =  1.0000E-04    !   threshold for determining if chains printed to *.dcs file
iwrchdt =       1        !   add extra information for contribution by reaction

! --- irradiation time ---
 itstep =       1        ! a decay calculation is just a single 'cooling' time step
   10.0   d  0.0         ! this time step must be >= to the final requested output time

! --- output time ---
  itout =        9       ! request times when inventory information is outputted
     1  m
     5  m
     10 m
     30 m
     3  h
     12 h
     1  d
     5  d
     10 d

! --- irradiation condition ---
 inmtcf =       0       ! set to 0 so DCHAIN won't look for [T-Yield] file

! --- calculation region details ---
 iregon =      1
!1)HRGCMM 2)IREGS 3)ITGNCLS 4)FLUXS 5)HNFLUXS 6)VOLUMES
DUMMY001  -1  4  0.0  no  1.0E+0 ! 4 nuclides to be tracked; set HNFLUXS to 'no' so
Hf-183       b1.00E+06           !  DCHAIN won't look for [T-Track] file; vol = 1 cm^3.
Tm-160m      b1.00E+06           ! Inventory is normally in 1E24 atoms/cm^3, but putting
Ta-169       b1.00E+06           ! 'b' before the value changes it to an activity in Bq.
Hf-164       b1.00E+06           !  Here each nuclide has a starting activity of 1 MBq.

end