7.10. [ T-DPA ] セクション

「原子あたりのはじき出し数」で評価される放射線損傷は、原子炉や加速器施設等で使用される材料の劣化を表す指標であり、大強度ビーム、特に重イオンに関して重要な課題の一つです。 原子のはじき出しは入射粒子と標的原子の衝突を起源として起こり、飛ばされた原子、一次はじき出し原子である PKA が再度別の原子に衝突することによって連鎖的に発生します。 このタリーではその頻度を DPA として出力します。 輸送荷電粒子のクーロン散乱断面積からの寄与を含んでおり、重イオン入射の場合にも適切に評価されます。 また、 e-mode>=1 を用いて、低エネルギーの中性子入射反応から生成する荷電粒子による DPA の寄与も計算できます。 なお、 [t-dpa] セクションでは各タリー結果の合計値は出力されません。

表 7.10.1 mesh

value

explanation

reg, r-z, xyz, tet

メッシュ型。 メッシュ型サブセクションが必要です。

表 7.10.2 volume

value

explanation

(省略可)

reg メッシュに対して各領域の体積を定義します。 このオプションの下に体積定義サブセクションが必要です。 ここで定義しない場合は、 [volume] で定義した値が用いられます。 reg=( ) を用いた特殊な領域指定を行い、この volume サブセクションを定義しない場合は、内部で定義された領域番号が input echo に出力されます。

reg vol

体積定義。 詳細は [volume] セクションを参照してください。

表 7.10.3 iechrl

value

explanation

72 (default)

体積入力 echo における1行あたりの最大カラム数。

表 7.10.4 part

value

explanation

all (default), 粒子名

タリーする粒子。

表 7.10.5 material

value

explanation

(省略可)

タリーする物質を指定します。

all

all が default で、未定義と同じです。

物質数

物質数を指定し、次の行で物質番号を定義します。 この値を負値にすると、指定した物質をタリー対象から除外します。

(次の行)

2 5 8 のように物質番号を指定します。

表 7.10.6 unit

value

explanation

1, 2

1: DPA [x1.E-24/source] 2: DPA [1/source]

表 7.10.7 idpa

value

explanation

0 (default), 1

0: 欠陥の非熱的再結合補正を考慮しない dpa、NRT-dpa です。 1: 欠陥の非熱的再結合補正を考慮する dpa、arc-dpa です。

表 7.10.8 axis

value

explanation

reg, x, y, z, r

出力データの \(x\) 軸。

tet

出力データの \(x\) 軸。 mesh=tet でのみ有効です。

xy, yz, xz, rz

2次元表示。

表 7.10.9 samepage

value

explanation

(省略可, D=part

画像出力ファイルで同じページに表示するデータの種類。 axis で定義できるパラメータを指定できます。

表 7.10.10 file

value

explanation

file name

出力ファイル名を定義します。 これは axis の各設定に対して必要です。

表 7.10.11 resfile

value

explanation

(省略可, D=file

restart 計算で過去の tally を読み込むファイル名を定義します。 axis を複数定義していても、指定する resfile は1つだけです。

表 7.10.12 output

value

explanation

dpa

total DPA を出力します。

表 7.10.13 mother

value

explanation

(省略可)

DPA に寄与する母核を限定します。 複数定義可能です。

all, 数

all はデフォルトです。 数を指定した場合、その数だけの核種を次の行に記述します。 負の場合は、それらの母核を対象から外すことを意味します。

表 7.10.14 item

value

explanation

208Pb Pb

質量数を指定すれば、その核です。 質量数を指定しなければ、その核の同位体全てを指定します。 複数の母核群を指定したいときは、複数の [t-dpa] セクションを定義します。

表 7.10.15 factor

value

explanation

(省略可, D=1.0

規格化定数。

表 7.10.16 title

value

explanation

(省略可)

タイトル。

表 7.10.17 angel

value

explanation

(省略可)

ANGEL パラメータ。

表 7.10.18 sangel

value

explanation

(省略可)

ANGEL パラメータの特別な書式。

表 7.10.19 2d-type

value

explanation

1, 2, 3, 4, 5, 6, 7

2次元表示のオプション。

(省略可, D=3

表 7.10.20 gshow

value

explanation

0 (default), 1, 2, 3, 4, 5

mesh=xyz かつ axis=xy,yz,xz の場合に、このオプションで領域境界(1)、物質名(2)、領域名(3)、LAT番号(4)を表示します。 gshow=5 は、 icntl=8 の場合にピクセル形式で物質色のみを表示します。

表 7.10.21 rshow

value

explanation

0 (default), 1, 2, 3

mesh=reg,tet かつ axis=xy,yz,xz の場合に、このオプションで領域境界(1)、物質名(2)、領域名(3)を表示します。 このオプションの下に xyz メッシュサブセクションが必要です。

表 7.10.22 x-txt

value

explanation

(省略可)

\(x\) 軸タイトル。

表 7.10.23 y-txt

value

explanation

(省略可)

\(y\) 軸タイトル。

表 7.10.24 z-txt

value

explanation

(省略可)

\(z\) 軸タイトル。

表 7.10.25 resol

value

explanation

1 (default)

gshow または rshow オプションの時、領域境界を求める分解能を各辺 resol 倍します。

表 7.10.26 width

value

explanation

0.5 (default)

線の太さを定義します。

表 7.10.27 volmat

value

explanation

(省略可, D=9

xyz メッシュで material を定義した場合に、各メッシュの体積値を補正するオプションです。 0 は補正なしを意味します。 volmat の値は、1つの xyz メッシュ辺を何分割して走査するかを表します。

表 7.10.28 epsout

value

explanation

0 (default), 1, 2

epsout=1 では、結果を eps ファイルとして出力します。 ファイル名は拡張子を .eps に置き換えたものになります。 epsout=2 では、2次元型の axis=xy, yz, xz, rz を除き、eps ファイルにエラーバーも表示します。

表 7.10.29 bmpout

value

explanation

0 (default), 1

2次元タリー出力の Bitmap 図を生成します。 このファイル名は拡張子を .bmp に置き換えたものになります。 mesh=xyz かつ axis=xy, yz, xz の場合に利用できます。

表 7.10.30 vtkout

value

explanation

0 (default), 1

xyz メッシュのタリー結果を ParaView の入力形式で出力します。 このファイル名は拡張子を .vtk に置き換えたものになります。 mesh=xyz かつ axis=xy, yz, xz の場合に利用できます。

表 7.10.31 vtkfmt

value

explanation

0 (default), 1

ParaView 用出力ファイルの形式。 0: ascii、 1: binary。

表 7.10.32 foamout

value

explanation

0 (default), 1, 2

foamout=1 では、数値データを OpenFOAM field data 形式で出力します。 このファイル名は拡張子を .foam に置き換えたものになります。 foamout=2 では、要素番号、重心の x 座標、y 座標、z 座標、体積、出力データ、相対誤差を CSV 形式で出力します。 このファイル名は拡張子を .csv に置き換えたものになります。 mesh=tet かつ axis=tet の場合に利用できます。

表 7.10.33 maxangel

value

explanation

part の数 (default)

eps ファイルに表示する part の数を指定します。 このパラメータはタリー eps ファイルに表示する粒子数だけを制限し、数値データファイルには影響しません。

表 7.10.34 ctmin(i)

value

explanation

(省略可, D=-9999

i 番目の counter の最小値。

表 7.10.35 ctmax(i)

value

explanation

(省略可, D=9999

i 番目の counter の最大値。

表 7.10.36 chmin(i)

value

explanation

(省略可, D=-9999

i 番目の history-counter の最小値。 このパラメータは batch variance mode の istdev=1 では指定できません。

表 7.10.37 chmax(i)

value

explanation

(省略可, D=9999

i 番目の history-counter の最大値。 このパラメータは batch variance mode の istdev=1 では指定できません。

表 7.10.38 trcl

value

explanation

(省略可)

r-z または xyz メッシュに対する座標変換番号、またはその定義。

表 7.10.39 dump

value

explanation

データ数

Dump機能により出力する粒子情報の数を定義します。 mesh = reg のみ有効です。

次行: データ定義

出力する粒子情報の順番を 表 6.7.1 に示す粒子情報番号で定義します。

表 7.10.40 gslat

value

explanation

2 (default)

gshow または rshow を指定した際の lattice および連続四面体の境界線に関するオプション。

0

描画しない。

1

描画する。 ただし、ボクセルファントムの線は明瞭に描画されないことがあります。

2

同じセル内の線は描画しない。

3

同じ物質内の線は描画しない。 この場合、同じ物質で満たされた隣接セルの境界は、lattice や連続四面体でなくても描画されません。 また、 gshow>=3 の場合でもセル番号や LAT 番号は表示されなくなります。

表 7.10.41 stdcut

value

explanation

(省略可, D=-1

STD cut off のしきい値。

stdcut を指定すると、PHITS は STD、すなわち標準偏差の値に応じて自動的 に計算を停止します。 この機能は stdcut が正で、 [parameters] セクションにおいて itall=0,1 を指定した場合に利用できます。 1バッチの最後に、タリー結果の STD の相対値がすべて 0 より大きく stdcut より小さい場合、計算は停止します。 2つ以上のタリーセクションで stdcut を設定した場合は、それらすべてが 条件を満たした時にこの機能が動作します。