放射化の計算 ( dchain タリー ) について

DCHAIN について

• DCHAINは、PHITSに同梱されている 放射化計算用のコードパッケージ ．

• PHITS→DCHAINのスムーズな実行のために、DCHAINの インプットファイル生成用タリー がある．

特徴

• デフォルトのデータベースに含まれていない( γ, n ) 光核反応などは、取り扱うことができない

• 見積精度をだすためには、 データベースをきっちりと指定する必要あり ．

• 数値計算で簡便に評価できるため、指標を得る上で有用．

実行について

1. t-dchainを使用する準備

2. [t-dchain]タリーの準備

3. PHTISの実行

4. DCHAINの実行

t-dchain タリーの設定

設定項目

• 【忘れがち】 [Parameter] セクションに jmout=1 を設定する．

• file(21) : DCHAIN用データファイルの保存先を設定する．

• 領域 reg を使用して設定する場合、 領域の体積情報を与えることが必要 ．以下２つの方法がある．

  • 直接、記載する ([Volume] セクション) ．

  • モンテカルロ計算を利用して計算する ( [t-volume]セクション ) ．

• ビーム強度の指定 ( source / 秒 )

  • e.g.) sourceセクションの強度をビーム電流に設定、ビームon/off設定で 1時間など、照射時間を設定．

• ビーム on/off, 照射時間 の指定

• 評価タイミングの指定

設定例

$$ --------------------------------------------------------- $$
$$ --- volume settings                                   --- $$
$$ --------------------------------------------------------- $$

[volume]

  reg           vol
  301		1.0
  302		2.0
  303		3.0



$$ --------------------------------------------------------- $$
$$ --- dchain tally settings                             --- $$
$$ --------------------------------------------------------- $$

[T-dchain]

    $$ -- mesh / output settings -- $$
    mesh    = reg
    reg     = 301 302 303
    file    = out/tdchain.out

    $$ -- beam on/off settings   -- $$
    timeevo = 2       $$ #. of 
    100  h     1.0    $$ (value) (unit) (beam intensity)
    365  d     0.0

    $$ -- output time interval   -- $$
    outtime = 12
    1.0    h
    2.0    h
    5.0    h
    10.0   h
    1.0    d
    2.0    d
    5.0    d
    10.0   d
    20.0   d
    50.0   d
    100.0  d
    365.0  d

    $$ -- amplitude -- $$
    amp      = 1.0     $$ [ source / sec ]
                       $$ ==>  if [source's factor == uA etc., => 1.0]
    

実行例

• out/tdchain.outを指定してDCAHINを起動する．

  • linuxでは dchain.sh, windowsは右クリック

dchain.sh out/tdchain.out

参考

• 下記ディレクトリに詳細な情報があるとのこと．

  • 実行サンプル ( phitsrecommendationDCHAIN )

  • DCHAIN ( phitsdchain-spmanual )