Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
小川 達彦; 平田 悠歩; 松谷 悠佑; 甲斐 健師; 佐藤 達彦; 岩元 洋介; 橋本 慎太郎; 古田 琢哉; 安部 晋一郎; 松田 規宏; et al.
EPJ Nuclear Sciences & Technologies (Internet), 10, p.13_1 - 13_8, 2024/11
放射線挙動解析コードPHITSは、モンテカルロ法に基づいてほぼ全ての放射線の挙動を解析することができる放射線挙動解析計算コードである。その最新版であるPHITS version 3.34の、飛跡構造解析機能に焦点を置いて説明する。飛跡構造解析とは、荷電粒子が物質中を運動する挙動を計算する手法の一つで、個々の原子反応を識別することにより原子スケールでの追跡を可能にするものである。従来の飛跡構造解析モデルは生体を模擬する水だけにしか適用できず、遺伝子への放射線損傷を解析するツールとして使われてきた飛跡構造解析であるが、PHITSにおいてはPHITS-ETS、PHITS-ETS for Si、PHITS-KURBUC、ETSART、ITSARという飛跡構造解析モデルが補い合うことにより、生体の放射線影響だけでなく、半導体や材料物質など任意物質に対する適用が可能になっている。実際にこれらのモデルを使って、放射線によるDNA損傷予測、半導体のキャリア生成エネルギー計算、DPAの空間配置予測など、新しい解析研究も発表されており、飛跡構造解析を基礎とするボトムアップ型の放射線影響研究の推進に重要な役割を果たすことが期待できる。
名越 康人*; 深堀 拓也*; 岡田 裕*; 高橋 昭如*; 下平 昌樹; 上田 貴志*; 小川 琢矢*; 八代醍 健志*; 高橋 由紀夫*; 大畑 充*
Transactions of the 27th International Conference on Structural Mechanics in Reactor Technology (SMiRT 27) (Internet), 9 Pages, 2024/03
日本溶接協会CAF小委員会では、塑性拘束効果を考慮した破壊評価手法ガイドラインの策定を目指している。この評価手法では、脆性破壊を評価するためのBereminモデルと延性亀裂成長を評価するためのGTNモデルを用いる。そこで、これらの評価モデルの適用性を検証するため、CAF小委員会の参加機関によるベンチマーク解析が行われた。ベンチマーク解析は、各機関が有する有限要素解析コードを用い、2種類の低合金鋼(A及びB)の破壊試験に対して実施されてきた。本発表では、低合金鋼Bに対する解析結果を報告する。Bereminモデルにおいて、一般的なワイブル形状母数( = 10, 20, 30)を用いた場合、各機関で計算されたワイブル応力が概ね一致することを確認した。また、Toughness Scaling Modelに基づいて、塑性拘束度が異なる2種類の試験片を用いてワイブル形状母数を算出した。算出されたワイブル形状母数は解析機関によりばらつきはあったものの、最終的に算出されるワイブル応力は一致することを確認した。GTNモデルに関して、評価に用いるパラメータを1T-C(T)試験片の室温での荷重-変位関係に基づいて最適化した。最適化されたパラメータを用いてGTNモデルに基づき評価されたJ-R曲線が各機関で一致することを確認した。
佐藤 達彦; 岩元 洋介; 橋本 慎太郎; 小川 達彦; 古田 琢哉; 安部 晋一郎; 甲斐 健師; 松谷 悠佑; 松田 規宏; 平田 悠歩; et al.
Journal of Nuclear Science and Technology, 61(1), p.127 - 135, 2024/01
被引用回数:53 パーセンタイル:99.93(Nuclear Science & Technology)放射線挙動解析コードPHITSは、モンテカルロ法に基づいてほぼ全ての放射線の挙動を解析することができる。その最新版であるPHITS version 3.31を開発し公開した。最新版では、高エネルギー核データに対する親和性や飛跡構造解析アルゴリズムなどが改良されている。また、PHIG-3DやRT-PHITSなど、パッケージに組み込まれた外部ソフトウェアも充実している。本論文では、2017年にリリースされたPHITS3.02以降に導入された新しい機能について説明する。
廣田 貴俊*; 名越 康人*; 北条 公伸*; 岡田 裕*; 高橋 昭如*; 勝山 仁哉; 上田 貴志*; 小川 琢矢*; 八代醍 健志*; 大畑 充*; et al.
Proceedings of ASME 2021 Pressure Vessels and Piping Conference (PVP 2021) (Internet), 9 Pages, 2021/07
In order to establish a guideline for fracture evaluation by considering plastic constraint in the ductile-brittle transition temperature (DBTT) region, the CAF (Constraint-Based Assessment of Fracture in Ductile-Brittle Transition Temperature Region) subcommittee has been launched in 2018 in the Japan Welding Engineering Society. In the committee, fracture tests are conducted using C(T), SE(B), and 50mm-thick flat plate with a surface flaw subjected to bending load or tensile load to verify fracture evaluation methods. Since simulation results are easily affected by analysis conditions, benchmark analysis is essential for the potential users of the guideline. Therefore, benchmark analyses are executed on brittle and ductile damages by Beremin and Gurson-Tvergaard-Needleman (GTN) models implemented in the finite element (FE) codes. The benchmark analyses are carried out in four steps; Step 0 is to confirm the output of FE codes in each member with the same input data and the same FE model. Step 1 is to confirm the result of Weibull stress analysis for C(T) specimens tested at -125C. The Weibull parameter, m, was fixed in this step. At step 2, sensitivity analyses are conducted on Weibull stresses in different conditions. The outputs by the GTN model are also confirmed. At the final step, the fracture simulation will be run for flat plate specimens with less plastic constraint than the standard fracture toughness specimen. As the results of the benchmark analyses up to step 2, a significant difference is not observed in the Weibull stress computed by committee members with the same input data and FE model and it is confirmed that the effects of element type, nonlinear deformation theory employed in FE analysis. For the calculation of the Weibull parameter m by using the fracture toughness test results and the developed programs by committee members, the converged values of m show good agreement among them.
井澤 一彦; 石井 淳一; 大久保 卓哉; 小川 和彦; 外池 幸太郎
Proceedings of 11th International Conference on Nuclear Criticality Safety (ICNC 2019) (Internet), 9 Pages, 2019/09
日本原子力研究開発機構は臨界実験装置STACYを非均質体系に変更する更新計画を進めている。更新されたSTACYでは、福島第一原子力発電所事故によって生成された燃料デブリを想定した臨界解析結果の検証のための実験が計画されている。STACY更新炉の初臨界は2021年初めに予定されている。初臨界後は、STACYでは「基本炉心」を構成し、運転員の習熟と実験結果の不確かさを把握するための一連の運転が行われる。STACY更新炉の初臨界に先立ち、基本炉心の核特性を把握するための一連の核特性解析を行った。本発表では、基本炉心の炉心構成条件を示すとともに、新規制基準のもとで炉心に課される諸条件をあきらかにする。
佐藤 達彦; 岩元 洋介; 橋本 慎太郎; 小川 達彦; 古田 琢哉; 安部 晋一郎; 甲斐 健師; Tsai, P.-E.; 松田 規宏; 岩瀬 広*; et al.
Journal of Nuclear Science and Technology, 55(6), p.684 - 690, 2018/06
被引用回数:880 パーセンタイル:99.99(Nuclear Science & Technology)粒子・重イオン輸送計算コードPHITS Version 3.02を開発して公開した。核反応モデルや原子反応モデルを改良することにより、適応エネルギー範囲や精度を向上させた。また、計算効率を向上させる新しいタリーや、放射性同位元素を線源とする機能、医学物理分野にPHITSを応用する際に有用となるソフトウェアなど、様々なユーザーサポート機能を開発して、その利便性を大幅に向上させた。これらの改良により、PHITSは、加速器設計、放射線遮へい及び防護、医学物理、宇宙線研究など、幅広い分野で3000名以上のユーザーに利用されている。本論文では、2013年に公開したPHITS2.52以降に追加された新機能を中心に紹介する。
橋本 慎太郎; 佐藤 達彦; 岩元 洋介; 小川 達彦; 古田 琢哉; 安部 晋一郎; 仁井田 浩二*
核データニュース(インターネット), (120), p.26 - 34, 2018/06
放射線の挙動を模擬できる粒子・重イオン輸送計算コードPHITSは、加速器の遮へい設計や宇宙線・地球科学、放射線防護研究等の広範な分野で利用されている。PHITSでは、放射線が引き起こす物理現象を輸送過程と衝突過程の2つに分けて記述しており、さらに核反応が含まれる衝突過程を「核反応の発生」、「前平衡過程」、「複合核過程」の3段階的に分けて計算することで模擬している。これらの計算では断面積モデルKurotamaや核反応モデルINCL4.6やJQMDといった数多くの物理モデルが使用されており、PHITSの信頼性を高めるために、我々は各モデルの高度化を行ってきた。本稿は、日本原子力学会2018年春の年会の企画セッション「我が国における核データ計算コード開発の現状と将来ビジョン」における報告を踏まえ、PHITSにおける核反応モデルの役割を解説するとともに、近年我々が行ってきた様々なモデルの高度化や今後の展開についてまとめたものである。
佐藤 達彦; 仁井田 浩二*; 岩元 洋介; 橋本 慎太郎; 小川 達彦; 古田 琢哉; 安部 晋一郎; 甲斐 健師; 松田 規宏; 奥村 啓介; et al.
EPJ Web of Conferences, 153, p.06008_1 - 06008_6, 2017/09
被引用回数:6 パーセンタイル:94.69(Nuclear Science & Technology)粒子・重イオン輸送計算コードPHITSは、原子力機構を中心とする日本やヨーロッパの研究所が共同で開発しているモンテカルロ放射線輸送計算コードである。PHITSには、様々な物理モデルやデータライブラリが含まれており、ほぼ全ての放射線の挙動を1TeVまで解析可能である。現在、国内外2500名以上のユーザーが、加速器設計、放射線遮へい、放射線防護、医学物理、地球惑星科学など様々な分野で利用している。本発表では、PHITS2.52からPHITS2.82までの間に導入した新しい物理モデルや計算機能などを紹介する。
佐藤 達彦; 岩元 洋介; 橋本 慎太郎; 小川 達彦; 古田 琢哉; 安部 晋一郎; 甲斐 健師; 松田 規宏; 岩瀬 広*; 仁井田 浩二*
放射線, 43(2), p.55 - 58, 2017/05
粒子・重イオン輸送計算コードPHITSは、原子力機構が中心となって開発しているモンテカルロ放射線輸送計算コードである。PHITSには、様々な物理モデルやデータライブラリが含まれており、ほぼ全ての放射線の挙動を1TeVまで解析可能である。これまでに、放射線施設の設計、医学物理計算、放射線防護研究、宇宙線・地球惑星科学など、工学,医学,理学の様々な分野で国内外2,500名以上の研究者・技術者に利用され、現在もその応用範囲は拡がっている。本稿では、その最新版であるPHITS2.88の特徴をまとめる。
岩元 洋介; 佐藤 達彦; 橋本 慎太郎; 小川 達彦; 古田 琢哉; 安部 晋一郎; 甲斐 健師; 松田 規宏; 細山田 龍二*; 仁井田 浩二*
Journal of Nuclear Science and Technology, 54(5), p.617 - 635, 2017/05
被引用回数:107 パーセンタイル:99.74(Nuclear Science & Technology)粒子・重イオン輸送計算コードPHITSは、加速器遮蔽設計、医学物理計算など様々な目的で利用されている。本研究では、最新のPHITSバージョン2.88を用いて、核反応による粒子生成断面積、中性子輸送計算、電磁カスケード等に対する58ケースのベンチマーク計算を実施した。本稿では、このうち特徴的な結果を示した22のケースについて詳細に報告する。ベンチマーク計算の結果、100MeV以上の陽子、中性子、重イオン等の入射反応や電磁カスケードについては、PHITSによる計算結果は実験値を概ね再現した。一方、100MeV未満の粒子入射反応は、PHITSが設定を推奨している核内カスケードモデルINCL4の適用範囲外のエネルギー領域であるため、再現性が悪いことが確認された。また、100MeV以上の陽子入射反応のうち、核分裂成分の収率に対して、蒸発モデルGEMの高エネルギー核分裂過程の取り扱いの問題から実験値を再現しない場合があった。これらの課題は、評価済み核データJENDL4.0/HEをPHITSに組み込むこと等で改善していく予定である。以上のように、本研究でPHITSの様々な適用分野における計算精度を検証するとともに、今後の効率的な改善に重要な指針を得ることができた。
岩元 洋介; 佐藤 達彦; 仁井田 浩二*; 橋本 慎太郎; 小川 達彦; 古田 琢哉; 安部 晋一郎; 甲斐 健師; 松田 規宏; 岩瀬 広*; et al.
JAEA-Conf 2016-004, p.63 - 69, 2016/09
粒子・重イオン輸送計算コードシステムPHITSは、原子力機構を中心に、複数の国内外の研究機関の協力の下で開発が進められている。同コードは、様々な核反応モデルやデータライブラリにより、中性子・陽子・重イオン・電子・光子等のほとんどの種類の粒子の輸送を取り扱うことが可能で、ソースプログラム,実行ファイル,データライブラリといった全ての構成要素が一つのパッケージにまとめられて配布されている。現在、1800人を超える研究者等がPHITSユーザーとして登録されており、核技術、加速器施設の設計、医学物理といった様々な分野で利用している。本発表ではPHITSに組み込まれている物理モデルを簡単にまとめ、最新のミューオン核反応モデル、脱励起モデルEBITEMについて紹介する。また、IAEA-CRP「初期の放射線損傷断面積」の活動の下で行っている、PHITSを用いた材料のはじき出し断面積、はじき出し原子エネルギースペクトル、カーマの計算結果について紹介を行う。
中嶋 享*; 寺本 高啓*; 赤木 浩; 藤川 高志*; 間嶋 拓也*; 峰本 紳一郎*; 小川 奏*; 酒井 広文*; 富樫 格*; 登野 健介*; et al.
Scientific Reports (Internet), 5, p.14065_1 - 14065_11, 2015/09
被引用回数:38 パーセンタイル:82.92(Multidisciplinary Sciences)X線自由電子レーザーを利用することで、レーザーによって整列したI分子からのX線光電子回折(XPD)パターンを観測した。X線自由電子レーザーの偏光方向に整列したI分子のXPDは我々の理論計算とよく一致した。さらに、実験で得られるXPDを分子構造決定に利用する際の適応基準を提案した。
佐藤 達彦; 仁井田 浩二*; 松田 規宏; 橋本 慎太郎; 岩元 洋介; 古田 琢哉; 野田 秀作; 小川 達彦; 岩瀬 広*; 中島 宏; et al.
Annals of Nuclear Energy, 82, p.110 - 115, 2015/08
被引用回数:38 パーセンタイル:95.18(Nuclear Science & Technology)原子力機構が中心となり日欧の複数機関が協力して汎用モンテカルロ粒子・重イオン輸送計算コードPHITSを開発している。PHITSは、幅広いエネルギー範囲のほぼすべての放射線の挙動を扱うことができ、そのすべてのコンポーネントは1つのパッケージにまとめられ、RIST, OECD/NEA Databank, RSICCを通じて全世界に配布されている。その結果、PHITSユーザー数は国内外で総計1,000名を超え、工学・理学・医学のさまざまな分野で利用されている。本論文では、PHITSの概要について紹介するとともに、イベントジェネレータモードやビーム輸送機能などPHITSに組み込まれた幾つかの重要な機能について解説する。
山本 真人*; 鬼沢 邦雄; 吉本 賢太郎*; 小川 琢矢*; 馬渕 靖宏*; Valo, M.*; Lambrecht, M.*; Viehrig, H.-W.*; 三浦 直樹*; 曽根田 直樹*
Small Specimen Test Techniques; 6th Volume (ASTM STP 1576), p.53 - 69, 2015/05
4mm厚のミニチュア破壊靭性試験片(微小C(T))によるマスターカーブ法の適用性を検証するため、日本の原子炉圧力容器鋼から採取した破壊靭性試験片を用いて、国内外の研究機関や産業界の参加を得てラウンドロビン試験を実施した。ASTM E1921規格に従って試験を行い、得られた破壊靭性参照温度の比較から、微小C(T)から得られる参照温度のばらつきは規格に示されている不確実さの範囲とほぼ同等であり、微小試験片に対するマスターカーブ法の有効性を確認した。また、本報では各機関が行った試験結果を取りまとめて統計処理を行い、参照温度には大きな相違はなく、標準サイズの試験片との寸法効果も認められないことが確認できた。
古田 琢哉; 佐藤 達彦; 小川 達彦; 仁井田 浩二*; 石川 顕一*; 野田 茂穂*; 高木 周*; 前山 拓哉*; 福西 暢尚*; 深作 和明*; et al.
Proceedings of Joint International Conference on Mathematics and Computation, Supercomputing in Nuclear Applications and the Monte Carlo Method (M&C + SNA + MC 2015) (CD-ROM), 9 Pages, 2015/04
粒子・重イオン輸送計算コードPHITSには計算時間短縮のために、二種類の並列計算機能が組み込まれている。一つはメッセージパッシングインターフェイス(MPI)を利用した分散メモリ型並列計算機能であり、もう一つはOpenMP指示文を利用した共有メモリ型並列計算機能である。それぞれの機能には利点と欠点があり、PHITSでは両方の機能を組み込むことで、利用者のニーズに合わせた並列計算が可能である。また、最大並列数が8から16程度のノードを一つの単位として、数千から数万というノード数で構成されるスーパーコンピュータでは、同一ノード内ではOpenMP、ノード間ではMPIの並列機能を使用するハイブリッド型での並列計算も可能である。それぞれの並列機能の動作について解説するとともにワークステーションや京コンピュータを使用した適用例について示す。
岩元 洋介; 佐藤 達彦; 仁井田 浩二*; 松田 規宏; 橋本 慎太郎; 古田 琢哉; 野田 秀作; 小川 達彦; 岩瀬 広*; 中島 宏; et al.
JAEA-Conf 2014-002, p.69 - 74, 2015/02
粒子・重イオン輸送計算コードシステムPHITSは、原子力機構を中心に、複数の国内外の研究機関の協力の下で開発が進められている。PHITSは、様々な核反応モデルやデータライブラリにより、中性子・陽子・重イオン・電子・光子等のほとんどの種類の粒子の輸送を取り扱うことができる。PHITSはFortran言語で記述されており、ソースプログラム、実行ファイル、データライブラリといった全ての構成要素が一つのパッケージにまとめられている。このパッケージは、国内へは高度情報科学技術研究機構、国外へは経済協力開発機構原子力機関のデータバンク、または米国放射線安全情報計算センターを通して、利用希望者へ配布されている。現在、1000人を超える研究者がPHITSユーザーとして登録されており、核技術、加速器施設の設計、医学物理といった様々な分野で利用している。本発表ではPHITSに組み込まれている物理モデルを簡単にまとめ、重要な機能のイベントジェネレータモード、データライブラリを用いた計算、適用例として材料の放射線損傷を計算するための機能を紹介する。また、PHITSのさらなる高度化に必要な核データコミュニティへの要望を示す。
山本 真人*; 木村 晃彦*; 鬼沢 邦雄; 吉本 賢太郎*; 小川 琢矢*; 馬渕 靖宏*; Viehrig, H.-W.*; 三浦 直樹*; 曽根田 直樹*
Proceedings of 2014 ASME Pressure Vessels and Piping Conference (PVP 2014) (DVD-ROM), 7 Pages, 2014/07
破壊靭性評価のためのマスターカーブ法は、最近試験規格として標準化され、原子炉圧力容器の信頼性を確保するための強力なツールであると期待されている。現行の監視試験において、マスターカーブ法のためのデータを得るためには、シャルピー試験片の試験後の半片から採取可能な小型の試験片の活用が重要である。著者らは、4mm厚のミニチュア破壊靭性試験片(微小C(T))によるマスターカーブ法の適用性を検証するため、典型的な日本の原子炉圧力容器鋼を用いて、国内の学界、産業界や研究機関の参加を得てラウンドロビン試験を進め、微小C(T)の有効性を確認した。本報ではブラインド試験により試験温度の選択を各研究機関が独自に行った結果を取りまとめ、試験温度依存性を比較した。その結果、得られた参照温度には大きな相違はなく、破壊靭性評価にあたって試験温度の選択は大きな影響を及ぼさないことが確認できた。
橋本 慎太郎; 仁井田 浩二*; 松田 規宏; 岩元 洋介; 岩瀬 広*; 佐藤 達彦; 野田 秀作; 小川 達彦; 中島 宏; 深堀 智生; et al.
医学物理, 33(2), p.88 - 95, 2013/10
粒子重イオン輸送計算コードPHITSは、任意の3次元空間における粒子の挙動を解析することができる汎用のモンテカルロ計算コードである。重イオンを含むほぼすべての粒子を非常に幅広いエネルギー範囲で取り扱うことができ、タリーと呼ばれる仮想の検出器を設定することで粒子フラックスやエネルギー付与などさまざまな物理量を導出することができる。さらに、最近、生体内の微視的領域のエネルギー付与を計算できるマイクロドジメトリタリー機能を開発し、医学物理分野における利用方法の幅が広がった。本稿では、2012年12月に公開した最新のPHITS2.52にある機能のうち、取り扱うことができる物理現象と核反応モデル、及びマイクロドジメトリタリー機能について解説する。加えて、医学物理分野においてPHITSを用いて得られた近年の成果として、マイクロドジメトリタリー機能を用いた新しい生物学的線量評価法の開発について紹介し、粒子線治療やX線治療、CT診断に関する研究への利用例も紹介する。
佐藤 達彦; 仁井田 浩二*; 松田 規宏; 橋本 慎太郎; 岩元 洋介; 野田 秀作; 小川 達彦; 岩瀬 広*; 中島 宏; 深堀 智生; et al.
Journal of Nuclear Science and Technology, 50(9), p.913 - 923, 2013/09
被引用回数:582 パーセンタイル:99.97(Nuclear Science & Technology)粒子・重イオン挙動解析コードPHITSを改良し、その最新版(バージョン2.52)をリリースした。最新版の主な特徴は、(1)新しい核反応モデルの組込と改良、(2)メモリ共有型並列計算機能の導入、(3).統計誤差計算方法の改良と再開始計算機能の導入、(4)残留放射能計算機能の組込、(5)核データ及び原子データライブラリの改訂、(6)ユーザーサポートツールの強化などである。これらの成果により、PHITSの計算精度・速度・機能などが格段に向上され、そのさらなる利用拡大が期待できる。本稿では、その最新版の特徴と幾つかの実験データに対するベンチマーク計算結果について記載する。
山本 真人*; 鬼沢 邦雄; 吉本 賢太郎*; 小川 琢矢*; 馬渕 靖宏*; 三浦 直樹*
Proceedings of 2013 ASME Pressure Vessels and Piping Conference (PVP 2013) (DVD-ROM), 8 Pages, 2013/07
破壊靭性評価のためのマスターカーブ法は、最近試験規格として標準化され、原子炉圧力容器の信頼性を確保するための強力なツールであると期待されている。現行の監視試験において、マスターカーブ法のためのデータを得るためには、シャルピー試験片の試験後の半片から採取可能な小型の試験片の活用が重要である。著者らは、4mm厚のミニチュア破壊靭性試験片(ミニCT)によるマスターカーブ法の適用性を検証するため、典型的な日本の原子炉圧力容器鋼を用いて、国内の学界,産業界や研究機関の参加を得てラウンドロビン試験を開始した。試験機の影響及び試験手順を検討した第1報に引き続き、本報では試験データを拡充するとともに、負荷速度の影響を評価した。得られた結果に基づき、塑性変形を考慮した負荷速度の定義や結果の相違を検討した結果、ミニCTによる破壊靭性参照温度には負荷速度は大きな影響を及ぼさないことが確認された。