Multi-Purpose Particle and Heavy Ion Transport code System

Slides:



Advertisements
Similar presentations
P HI T S 陽子ビームで雪だるまを溶かそう Multi-Purpose Particle and Heavy Ion Transport code System title 年 3 月改訂.
Advertisements

P HI T S α 線, β 線, γ 線,中性子線を止めるに は? Multi-Purpose Particle and Heavy Ion Transport code System title 年 5 月改訂.
P HI T S PSFC4PHITS の使い方 Multi-Purpose Particle and Heavy Ion Transport code System Title 年 2 月改訂.
シーベルトとは? 高エネルギー加速器研究機構 平山 英夫 日本原子力学会 2014 年秋の大会 保健物理・環境科学部会企画セッション.
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用
ボクセルファントムを用いた X線治療のシミュレーション
PHITS 講習会 基礎実習(III): 計算条件の設定
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
PHITS講習会 基礎実習(I): 体系及び線源の定義
PHITS講習会 基礎実習(II): Tally(タリー)の定義
強度変調回転放射線治療(VMAT)の シミュレーション
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
相対論的重イオン衝突実験 PHENIXにおける Aerogel Cherenkov Counterの シミュレーションによる評価
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
物理化学(メニュー) 0-1. 有効数字 0-2. 物理量と単位 0-3. 原子と原子量 0-4. 元素の周期表 0-5.
Multi-Purpose Particle and Heavy Ion Transport code System
放射能計算の基礎と DCHAIN-SPの応用実習
第9回 星間物質その2(星間塵) 東京大学教養学部前期課程 2012年冬学期 宇宙科学II 松原英雄(JAXA宇宙研)
Multi-Purpose Particle and Heavy Ion Transport code System
W e l c o m ! いい天気♪ W e l c o m ! 腹減った・・・ 暑い~ 夏だね Hey~!! 暇だ。 急げ~!!
PHITSによるX線治療シミュレーション 基本操作復習編
ParaViewを用いたPHITS 計算結果の3次元表示
物質(人体含む)が吸収した放射線のエネルギー
放射線(エックス線、γ線)とは? 高エネルギー加速器研究機構 平山 英夫.
In situ cosmogenic seminar
スクリプト言語を用いたPHITSの連続実行
埼玉大学大学院理工学研究科 物理機能系専攻 物理学コース 06MP111 吉竹 利織
PHITS 講習会 基礎実習(III): 計算条件の設定
Multi-Purpose Particle and Heavy Ion Transport code System
誘導放射能の評価方法: PHITSとDCHAIN-SPの接続計算
α線,β線,γ線,中性子線を止めるには?
In situ cosmogenic seminar
原子核物理学 第4講 原子核の液滴模型.
HLab meeting 6/03/08 K. Shirotori.
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用 A
X線CTにおけるファントム中の エネルギー変化についての検討
2.5MeV単色中性子に対する TGCの動作特性の研究
応用実習用資料 Environmental radioactivity
鉄核媒質中の閾値近傍における 中性π中間子生成実験
Multi-Purpose Particle and Heavy Ion Transport code System
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用 A
PHITS 講習会 基礎実習(III): 計算条件の設定
放射線防護に用いられる線量概念 ー実効線量ー 岩井 敏 原子力安全推進協会 保健物理・環境科学部会セッション
応用実習用資料 Neutron target
IAEA phase space fileを用いた X線治療シミュレーション
Azimuthal distribution (方位角分布)
研究背景 電荷移行反応とは・・・ 核融合(重水素 + 三重水素→ヘリウム原子核+中性子) ・・・しかし、
理研RIBFにおける 中性子過剰Ne同位体の核半径に関する研究
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
治療用フィルムによる線量分布測定の 基礎的検討Ⅱ
Multi-Purpose Particle and Heavy Ion Transport code System
In situ cosmogenic seminar
東北大学電子光理学研究センター(ELPH) 菊永英寿
Charmonium Production in Pb-Pb Interactions at 158 GeV/c per Nucleon
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
STOF check.
Multi-Purpose Particle and Heavy Ion Transport code System
SMILE-2用simulation 上野一樹.
卒業論文発表 中性子ハロー核14Beの分解反応 物理学科4年 中村研究室所属   小原雅子.
α線,β線,γ線,中性子線を止めるには?
計算と実測値の比較 高エネルギー加速器研究機構 平山 英夫.
Multi-Purpose Particle and Heavy Ion Transport code System
Simulation study for drift region
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用 B
荷電粒子の物質中でのエネルギー損失と飛程
Presentation transcript:

Multi-Purpose Particle and Heavy Ion Transport code System PHITS Multi-Purpose Particle and Heavy Ion Transport code System 遮蔽設計の演習 PHITS講習会 入門実習 2017年1月改訂 title 1

実習目的 ? 高エネルギー中性子を効率良く遮蔽する壁材の検討をPHITSを用いて行ってみよう 高エネルギー中性子 「実習」 前の基本的な話から 線量の評価は実効線量で行うこととする Contents 2

実効線量とは? 物理量 実用量 防護量 吸収線量(Gy) フルエンス 線量換算係数 周辺線量当量(Sv) 個人線量当量(Sv) シミュレーション 人体模型(ファントム) 放射線加重係数[WR] 組織加重係数[WT] シミュレーション ICRU球 線質係数[Q(L)] 線量換算係数 実用量 防護量 周辺線量当量(Sv) 個人線量当量(Sv) 実効線量(Sv) 校正 関連 「実習」 前の基本的な話から 放射線測定器の指示値 放射線による健康リスク サーベイメータ、個人線量計 がん罹患率、致死率 線量換算係数を用いた計算には[T-deposit]ではなく[T-track]を使う Effective dose 3

shield.inpの確認 基本計算条件 計算体系 trackXZ.eps cross.eps 入射粒子: 体系: タリー: 200MeV陽子(半径0.01cmのペンシルビーム) 半径50cmで10cm厚の円柱が10個並ぶ(中身は空気) それ以外の領域も空気で満たされている [t-track]によるフラックス空間分布(二次元xz、一次元z) [t-cross]により上記円柱面毎のエネルギー分布 10個 陽子フラックス (1ページ目) 領域100 ⇒領域1 陽子 中性子 … 200MeV 陽子 空気 10cm 計算体系 trackXZ.eps cross.eps Check Input File 4

ステップ1:中性子を発生させよう タングステン標的を設定し、陽子線照射によって中性子を発生させよう Z軸中心の半径5cmで厚み5cm(z=- -10 ~ -5)の円柱(領域20) タングステンは登録済み(マテリアル2番)、密度は19.25g/cm3 領域100から標的領域を除く 中性子フラックス (2ページ目) 陽子フラックス (1ページ目) trackXZ.eps 入射陽子は 標的内で停止 標的との衝突で 中性子が発生 Step 1 5

ステップ2:実効線量に換算しよう 2番目 Multiplierセクションの換算係数を使用し、 フラックスを実効線量に換算しよう [ T - T r a c k ] title = Track Z ... y-txt = Effective dose [pSv/source] multiplier = all part = neutron emax = 1000.0 mat mset1 all (1.0 -201) part = photon all (1.0 -202) … [ M u l t i p l i e r ] number = -201 interpolation = log ne = 68 1.0E-9 3.09 1.0E-8 3.55 2番目 Multiplierセクションの換算係数を使用し、 フラックスを実効線量に換算しよう グラフのy軸の 表記を変更 Multiplier サブセクション追加 使用するMultiplier番号 規格化定数 線量換算係数[ICRP116] (フラックス⇒実効線量) 1/cm2    pSv trackZ.eps 中性子の寄与がほとんど Step 2 6

ステップ3:陽子ビーム電流を調整する ビームを電流値1Aで連続照射している時の実効線量率(Sv/h)を計算してみよう ヒント 2番目 実効線量は換算係数を掛けた結果pSv/sourceの単位で出力されている 1A(アンペア)は,1秒当たり1C(クーロン)の電流が流れる状態を表す 素電荷(陽子1つ当たりの電荷)は1.6x10-19Cとする (マイクロ)は10-6でp(ピコ)は10-12 [ T - T r a c k ] ... y-txt = Effective dose [pSv/source] multiplier = all part = neutron emax = 1000.0 mat mset1 all (1.0 -201) part = photon all (1.0 -202) 2番目 1アンペアで1秒照射するときの発生陽子数は 1.0 / 1.6e-19 = 6.25e18(個) 1Aで1時間照射する際の発生陽子数は 6.25e18 x 3600 x 1.0e-6 = 2.25e16(個) Sv/hで出力するための規格化定数 2.25e16 x 1.0e-12 = 2.25e4 [Sv/h] 2.25e4 距離100~105cm ⇒ trackZ.outの84行目 6.8435E+01 Sv/h Step 3 7

ステップ4:壁で遮蔽しよう 領域1~2 (合計20cm)の材質を変更しよう コンクリート(物質3、密度2.2g/cm3) 縦軸を揃えるために二番目の[T-track]タリーにangel = ymin(1.0e-2) ymax(1.0e3) を追加 物質を確認しやすくするため一番目の[T-track]タリーのgshowオプションを2に変更 コンクリート(物質3、密度2.2g/cm3) 鉄(物質4、密度7.7g/cm3) 2.3065E+01 2.6781E+01 trackZ.eps Step 4 8

ステップ5:壁の厚みを増やしてみよう 領域1~10 (合計100cm)の材質を変更しよう コンクリート(物質3、密度2.2g/cm3) trackXZ.eps 中性子深層透過計算になるため統計精度が悪い Step 5 9

ステップ6:中性子を遠くまで飛ばそう インポータンス設定して、端まで中性子が飛ぶようにしよう コンクリート trackXZ.eps [importance] set: c1[1.0] part = neutron photon reg imp 100 c1**0 20 c1**0 2 c1**1 3 c1**2 4 c1**3 5 c1**4 6 c1**5 7 c1**6 8 c1**7 9 c1**8 10 c1**9 200 c1**9 1以上の値に設定してみる trackXZ.eps あまり大きな値を設定するとメッセージが表示され急激に計算が遅くなる jbnk = 0, ibnk = 1 ... **** Warning: Too many secondary particles created **** **** MAXBNK overflowed thus HDD is used 10 times **** 距離100~105cmの実効線量⇒ 2.2579E+00 Sv/h 1.6603E+00 Sv/h コンクリート 鉄 c1=2.0 Step 6 10

密度が高ければいいのか? 1.6603E+0 4.6125E+0 trackZ.eps ∝ A2/3 × 密度/A Step 6 11 鉄(密度7.7g/cm3)の代わりに鉛(密度11.34g/cm3)で試してみよう 鉛を物質5として追加し、領域1~10 に使用する [ M a t e r i a l ] MAT[5] Pb 1.0 鉄(物質4、密度7.7g/cm3) 1.6603E+0 鉛(物質5、密度11.34g/cm3) 4.6125E+0 鉄に比べて 遮蔽効果が小さい trackZ.eps 高エネルギー中性子の 反応断面積(遮蔽効果) 原子核当りの 反応断面積 単位体積当たり の原子核数 ∝ × ∝ A2/3 ×  密度/A 鉄 2.01 鉛 1.91 Step 6 11

ステップ7:素材を組み合わせてみよう 領域1~5 を鉄(物質4、密度7.7g/cm3) とした場合の実効線量を鉄、コンクリート単体の場合と比較しよう 鉄とコンクリートの場所を反対(コンクリートが上流で鉄が下流)にするとどう変わるか? 2.3258E-01 鉄 コンクリート 1.6048E+00 コンクリート 鉄 trackZ.eps Step 7 12

透過した中性子のスペクトルは? cross.eps (コンクリート⇒鉄) cross.eps (鉄⇒コンクリート) 空気⇒コンクリート コンクリート20cm コンクリート⇒鉄 鉄30cm 鉄⇒空気 cross.eps (コンクリート⇒鉄) 空気⇒鉄 鉄20cm 鉄⇒コンクリート コンクリート30cm コンクリート⇒空気 cross.eps (鉄⇒コンクリート) 鉄(高密度)でエネルギーを落としてから、コンクリート(水素元素を含む)で低エネルギー中性子を止めることでより遮蔽できる Tally 13

ステップ8:壁の放射化を調べよう [t-dchain] tallyを有効にし、1時間照射した後の50年の線量率の推移を10年毎に調べよう 但し、領域1~5 を鉄、領域6~10をコンクリートとする jmout = 1 file(21)= c:/phits/dchain-sp/data e-mode = 0 [parameters]セクション に追加 Offをとる PHITSを実行して得られるtdchain.outをインプットとしてDCHAINを実行 [ T - D C H A I N ] $ must section for DCHAIN title = Induced radiation mesh = reg reg = 5 6 file = tdchain.out timeevo = 2 1.0 h 1.0 50.0 y 0.0 outtime = 6 1.0 h 10.0 y 20.0 y 30.0 y 40.0 y 50.0 y $ beam current (nA) set:c21[1000.0] amp = c21*1.0e-9/1.602e-19 26Alが支配的 鉄 コンクリート tdchain.eps(6ページ目) Step 8 14

微量不純物の影響を調べよう tdchain.out (50行目付近) 極微少量(1ppm)の微量不純物(59Co)を鉄に加えて(tdchain.outを修正)、DCHAINの実行をやり直してみよう !1)HRGCMM 2)IREGS 3)ITGNCLS ... DUMMY001 5 2 2.4071E+08 ... Fe 8.3055E-02 Co-59 1.0000E-06 構成元素数 微量不純物(59Co)なし 微量不純物(59Co)あり コンクリート 鉄 数十年後の線量率寄与は微量元素から生成された60Coが支配的 鉄 コンクリート tdchain.eps(6ページ目) Step 8 15

まとめ [Multiplier]セクションと[T-track]タリーにより、 実効線量を計算できる 高エネルギー中性子の場合、高密度の物質 でエネルギーを落としてから、原子番号の軽 い原子を含む物質で遮蔽すると効果的 長期間の誘導放射能の評価には、長寿命 核種を生む微量不純物を含めた評価が重要 Summary 16

宿題(超難題!) ターゲット(タングステン)の放射化を計算する 現在のビーム設定で1時間照射し、1日後の値を調べる 放射化ターゲットから1mの距離での実効線量(SV/h)を計算 ヒント(手順) PHITS⇒DCHAIN⇒PHITSの順で計算 1st PHITS [t-dchain]の設定を変更する [volume]セクションでターゲットの体積を定義する 2nd PHITS DCHAINの出力(tdchain.pht)から[source]セクションをコピーする 壁を全て空気に変更し、[Impotance]セクションをoffにする [t-track]タリーのmultiplierサブセクションの規格化は3600x10-6=3.6E-3 カラーバーのテキストは”z-txt = ***”で変更可能 Homework 17

回答例(answer-step1.inp, answer-step2.inp) trackZ.eps trackXZ.eps DCHAINの簡易見積もり結果(tdchain.actの1121行目)の値よりも一桁程小さい total g-ray dose-rate 2.42797E+03 [uSv/h*m^2] ターゲット自身による自己遮蔽の効果 Homework 18