Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
古田 琢哉; 佐藤 達彦
Radiological Physics and Technology, 14(3), p.215 - 225, 2021/09
PHITSユーザー数は、集計を開始した2010年以降に順調な増加傾向を示している。その中で、近年は医学物理分野における新規ユーザー数の増加が顕著となっている。そのため、PHITSを利用した数多くの研究が発表されており、その内容も放射線治療応用、施設の遮蔽計算、放射線生物学応用、医療機器研究開発等、多岐に渡っている。本稿では、これら応用研究の具体例を示しつつレビューを行うとともに、医学物理応用研究で有益なPHITSの機能を紹介する。
藤本 望*; 多田 健一; Ho, H. Q.; 濱本 真平; 長住 達; 石塚 悦男
Annals of Nuclear Energy, 158, p.108270_1 - 108270_8, 2021/08
被引用回数：0 パーセンタイル：0.02(Nuclear Science & Technology)Japan Atomic Energy Agency has developed a new nuclear data processing code, namely FRENDY, to generate the ACE files from various nuclear libraries. A code-to-experiment verification of FRENDY processing was carried out in this study with criticality benchmark assessments of the high temperature engineering test reactor. The ACE files of the JENDL-4.0 and ENDF-B-VII.1 was generated successfully by FRENDY. These ACE files have been used in MCNP6 transportation calculation for various benchmark problems of the high temperature engineering test reactor. As a result, the k
and reaction rate obtained by MCNP6 calculation presented a good agreement compared to the experimental data. The proper ACE files generation by FRENDY was confirmed for the HTTR criticality calculations.
植木 太郎
Nuclear Science and Engineering, 195(2), p.214 - 226, 2021/02
被引用回数：1 パーセンタイル：42.23(Nuclear Science & Technology)極端な無秩序の下では、物理的な力学系は、逆冪乗則パワースペクトルで特徴づけられる状態に落ち着く。本論文は、このような無秩序状態の乱雑系に関して、ホワイトノイズを含む形へのワイエルシュトラス関数モデルの拡張および多種物質混合方式の開発を報告する。具体的な逆冪乗則パワースペクトルの式および多種物質系でも計算機上で実現可能な混合手順を提示する。SOLOMONモンテカルロソルバーおよびJENDL-4核データーライブラリを使用した数値計算例についても報告する。
Cs and 
Cs radioactivity measurementsMalins, A.; 今村 直広*; 新里 忠史; 高橋 純子*; Kim, M.; 佐久間 一幸; 篠宮 佳樹*; 三浦 覚*; 町田 昌彦
Journal of Environmental Radioactivity, 226, p.106456_1 - 106456_12, 2021/01
被引用回数：0 パーセンタイル：0(Environmental Sciences)Understanding the relationship between the distribution of radioactive 
Cs and 
Cs in forests and ambient dose equivalent rates (
*(10)) in the air is important for researching forests in eastern Japan affected by the Fukushima Dai-ichi Nuclear Power Plant (FDNPP) accident. This study used a large number of measurements from forest samples, including Cs and Cs radioactivity concentrations, densities and moisture contents, to perform Monte Carlo radiation transport simulations for *(10) between 2011 and 2017. Calculated *(10) at 0.1 and 1 m above the ground had mean residual errors of 19% and 16%, respectively, from measurements taken with handheld NaI(Tl) scintillator survey meters. Setting aside the contributions from natural background radiation, Cs and Cs in the organic layer and the top 5 cm of forest soil generally made the largest contributions to calculated *(10). The contributions from Cs and Cs in the forest canopy were calculated to be largest in the first two years following the accident. Uncertainties were evaluated in the simulation results due to the measurement uncertainties in the model inputs by assuming Gaussian measurement errors. The mean uncertainty (relative standard deviation) of the simulated *(10) at 1 m height was 11%. The main contributors to the total uncertainty in the simulation results were the accuracies to which the Cs and Cs radioactivities of the organic layer and top 5 cm of soil, and the vertical distribution of Cs and Cs within the 5 cm soil layers, were known. Radioactive cesium located in the top 5 cm of soil was the main contributor to *(10) at 1 m by 2016 or 2017 in the calculation results for all sites.
Malins, A.; 越智 康太郎; 町田 昌彦; 眞田 幸尚
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.147 - 154, 2020/10
Compton-to-peak analysis is a method for selecting suitable coefficients to convert count rates measured with in situ gamma ray spectrometry to radioactivity concentrations of Cs & Cs in the environment. The Compton-to-peak method is based on the count rate ratio of the spectral regions containing Compton scattered gamma rays to that with the primary Cs & Cs photopeaks. This is known as the Compton-to-peak ratio (RCP). RCP changes as a function of the vertical distribution of Cs & Cs within the ground. Inferring this distribution enables the selection of appropriate count rate to activity concentration conversion coefficients. In this study, the PHITS Monte Carlo radiation transport code was used to simulate the dependency of RCP on different vertical distributions of Cs & Cs within the ground. A model was created of a LaBr
(Ce) detector used in drone helicopter aerial surveys in Fukushima Prefecture. The model was verified by comparing simulated gamma ray spectra to measurements from test sources. Simulations were performed for the infinite half-space geometry to calculate the dependency of RCP on the mass depth distribution (exponential or uniform) of Cs & Cs within the ground, and on the altitude of the detector above the ground. The calculations suggest that the sensitivity of the Compton-to-peak method is greatest for the initial period following nuclear fallout when Cs & Cs are located close to the ground surface, and for aerial surveys conducted at low altitudes. This is because the relative differences calculated between RCP with respect to changes in the mass depth distribution were largest for these two cases. Data on the measurement height above and on the Cs & Cs activity ratio is necessary for applying the Compton-to-peak method to determine the distribution and radioactivity concentration of Cs & Cs in the ground.
久保 光太郎; Zheng, X.; 田中 洋一; 玉置 等史; 杉山 智之; Jang, S.*; 高田 孝*; 山口 彰*
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.308 - 315, 2020/10
動的確率論的リスク評価(PRA)PRAは従来のPRA手法の現実性と網羅性を向上させる手法の一つである、しかしながら、それらの向上と引き換えに膨大な計算コストが発生する。本稿では、複数のサンプリング手法を簡易的な事故シーケンスに対する動的PRAに対して適用した。具体的には、モンテカルロ法,ラテン超方格法,格子点サンプリング及び準モンテカルロ法を比較した。その結果、今回の検討の範囲においては、準モンテカルロ法が最も効率的であった。
原田 秀郎; 高山 直毅; 米田 政夫
Journal of Physics Communications (Internet), 4(8), p.085004_1 - 085004_17, 2020/08
原子炉を用いた放射化分析などで重要な中性子共鳴積分値を高精度化するため、熱外中性子スペクトルの新しい近似を定式化した。近似式の導出に当たっては、はじめにモンテカルロ計算コードMVP-3を用いて参照解となる中性子スペクトルを計算し、これから 型の関数型を導出した。従来の近似式に比較し、導出した関数型は、中性子共鳴積分値を高精度に決定できることを示した。この検討は、過去にJRR-3で行われた
Csの中性子共鳴積分値の測定データに基づき行われた。また、提唱した近似式に導入したパラメータ
及び
を実験的に決定するため、3種類のフラックスモニター(
Au, 
Co及び
Zr)を用いる手法を提唱するとともに、解析手法を定式化した。
Riyana, E. S.; 奥村 啓介; 寺島 顕一; 松村 太伊知; 坂本 雅洋
Mechanical Engineering Journal (Internet), 7(3), p.19-00543_1 - 19-00543_8, 2020/06
Prediction of the fuel debris location and distribution inside the primary containment vessel (PCV) of the Fukushima Daiichi Nuclear Power Plant is important to decide further decommissioning step and strategy. The radiation measurements in the past internal investigations have not yet provided enough information to predict fuel debris location and its distribution inside PCV. To support further measurement efforts, we simulate the detector response inside the PCV. The calculation result could provide a base on deciding suitable detector systems to assist the efforts on searching, localizing and defining distributions of the fuel debris.
長住 達; 松中 一朗*; 藤本 望*; 石井 俊晃; 石塚 悦男
JAEA-Technology 2020-003, 13 Pages, 2020/05
JAEA-Technology-2020-003.pdf:1.5MB
MVPコードを用いて、制御棒の幾何形状を実機の制御棒構造に近づけた詳細モデルを作成し、制御棒モデルの詳細化によるHTTRの核特性への影響について検討した。この結果、臨界制御棒位置は、制御棒モデルの詳細化により、従来のモデルと比べて11mm低くなり、実測値である1775mmに近づいた。また、制御棒先端のショックアブソーバーにより吸収された反応度は0.2%
k/kとなり、臨界制御棒位置にして14mmの差となることが分かった。さらに、制御棒の形状効果によるSRACコードの解析値に対する補正量は、制御棒モデルの詳細化とショックアブソーバーによる影響を考慮して、低温臨界時において反応度で-0.05%k/k、臨界制御棒位置にして-3mmとなった。
Malins, A.; 町田 昌彦; Vu, TheDang; 相澤 一也; 石田 武和*
Nuclear Instruments and Methods in Physics Research A, 953, p.163130_1 - 163130_7, 2020/02
被引用回数：5 パーセンタイル：91.19(Instruments & Instrumentation)Radiation transport simulations were used to analyse neutron imaging with the current-biased kinetic inductance detector (CB-KID). The PHITS Monte Carlo code was applied for simulating neutron, 
He, 
Li, photon and electron transport, 
B(n,)Li reactions, and energy deposition by particles within CB-KID. Slight blurring in simulated CB-KID images originated He and Li ions spreading out in random directions from the B conversion layer in the detector prior to causing signals in the 
and 
superconducting Nb nanowire meander lines. 478 keV prompt gamma rays emitted by Li nuclei from neutron-B reactions had negligible contribution to the simulated CB-KID images. Simulated neutron images of B dot arrays indicate that sub 10 
m resolution imaging should be feasible with the current CB-KID design. The effect of the geometrical structure of CB-KID on the intrinsic detection efficiency was calculated from the simulations. An analytical equation was then developed to approximate this contribution to the detection efficiency. Detection efficiencies calculated in this study are upper bounds for the reality as the effects of detector temperature, the bias current, signal processing and dead-time losses were not taken into account. The modelling strategies employed in this study could be used to evaluate modifications to the CB-KID design prior to actual fabrication and testing, conveying a time and cost saving.
Kim, M.; Malins, A.; 吉村 和也; 佐久間 一幸; 操上 広志; 北村 哲浩; 町田 昌彦; 長谷川 幸弘*; 柳 秀明*
Journal of Environmental Radioactivity, 210, p.105803_1 - 105803_10, 2019/12
被引用回数：2 パーセンタイル：21.19(Environmental Sciences)放射能被害地域において空間線量率のシミュレーション精度を向上させるには、環境内の放射性核種の異なる分布、例えば、農地, 都市, 森林の放射能レベル差を考慮して現実的にモデル化する必要がある。さらには建物, 樹木, 地形による
線の遮蔽効果をモデルに考慮すべきである。以下に、福島県の市街地及び農地の3次元モデルの作成システムの概要を述べる。線源設定はCs及びCsの放射能分布をモデルのさまざまな環境要素に異なる分布設定が可能である。構造物については、現地の建物モデルにおいては日本の典型的な9種類の建物モデルを用いて作成される。また、樹木については広葉樹と針葉樹モデル、地形モデルは、地形を考慮した地表面モデルを取り込んだ。計算対象のモデルの作成時は、数値標高モデル(DEM),数値表面モデル(DSM)及びユーザー編集の際にサポートする対象領域のオルソ画像で作られる。計算対象のモデルが作成されたら、放射線輸送解析計算コードであるPHITSに適したフォーマットでシステムから出力される。上記のシステムを用いて、福島第一原子力発電所から4km離れた地域でかつ、まだ除染作業が行われてない郊外を計算対象としてモデルを作成した。モデル作成後、PHITSによる空間線量率の計算結果は走行サーベイの実測値との相関があった。
Carter, L. M.*; Crawford, T. M.*; 佐藤 達彦; 古田 琢哉; Choi, C.*; Kim, C. H.*; Brown, J. L.*; Bolch, W. E.*; Zanzonico, P. B.*; Lewis, J. S.*
Journal of Nuclear Medicine, 60(12), p.1802 - 1811, 2019/12
被引用回数：9 パーセンタイル：80.16(Radiology, Nuclear Medicine & Medical Imaging)従来の内部被ばく研究では、CT画像などに基づき構築したボクセル人体ファントムを適用し、線量評価データを解析することが一般的であった。しかし、皮膚や膜状組織などをより精密に表現できるポリゴンメッシュ人体ファントムの開発が進められており、内部被ばく研究でもこのファントムを用いた評価が要求されている。しかし、モンテカルロ放射線輸送計算コードにこのファントムを取り込んで計算を行うことは、コードに精通していない場合は難しい現状となっている。そこで、PHITSの使い方に関する詳しい知識を持たないユーザーでも、このファントムを用いた内部被ばく線量計算が簡単にできるツールPARaDIMを開発した。このツールでは、グラフィカルインターフェースにより、使用する四面体メッシュファントムの選択や臓器毎の放射性核種の設定、そしてPHITSの実行が操作できる。このツールを用いた実行例をいくつか示し、先行研究との比較を行うことで、本ツールの有用性を確認した。
-rays elastic scattering by Monte Carlo simulationOmer, M.; 羽島 良一*
New Journal of Physics (Internet), 21(11), p.113006_1 - 113006_10, 2019/11
被引用回数：4 パーセンタイル：56.93(Physics, Multidisciplinary)Nuclear resonance fluorescence (NRF) is a promising technique for nondestructive assay (NDA) of nuclear materials. Its strength is improved when polarized -ray beams are used as probes because the polarized -rays provide an asymmetry in the intensity of the scattered photons. Nonetheless, NRF interaction takes place with other interactions such as elastic scattering of -rays. These interactions are unavoidable background which may limit the sensitivity of the NRF technique. Therefore, estimating polarization effects of elastic scattering is a significant factor in assessing the NRF method. Here, we implement a new Monte Carlo simulation to account for polarization effects of the elastic scattering. The simulation is based on explicit expressions driven in the Stokes parameters formalism. The scattering amplitudes of Rayleigh, nuclear Thomson, and Delbr
ck scattering processes have been superimposed into a two-orthogonal set of complex amplitudes. This set is then exploited to construct the core of the simulation in a way such that the simulation could handle arbitrary polarization states of the incoming beam and correspondingly generate polarization states for the outgoing beam. This work was supported by the subsidiary for promotion of strengthening nuclear security or the like of the Ministry of Education, Culture, Sports, Science, and Technology (MEXT), Japan.
古田 琢哉; El Basha, D.*; Iyer, S. S. R.*; Correa Alfonso, C. M.*; Bolch, W. E.*
Journal of Radiological Protection, 39(3), p.825 - 837, 2019/09
被引用回数：1 パーセンタイル：24.18(Environmental Sciences)人の眼球には多様性があるにもかかわらず、これまでの眼球に対する線量計算シミュレーションでは、ほぼ全ての研究で標準的な一つのモデルが採用されてきた。そこで、本研究では新たに開発した大きさ及び変形が可能な数値眼球モデルを利用し、モンテカルロ放射線輸送計算コードPHITSを用いて、標準的な照射場(AP, PA, RLAT, ROT等)での電子線,光子線および中性子線による眼球内組織の線量を計算した。ここでは、5種類(標準,大型,小型,近視形状型および遠視形状型)の極端な眼球モデルの計算結果の比較に基づき、組織の吸収線量への眼球の大きさや変形が与える影響を解析した。電子線では線量の集中性が高く、吸収線量は眼球表面からの組織の深さに大きく依存して変化する。このため、眼球の大きさや変形に伴い、組織全体の深さ位置が変化することによる吸収線量への顕著な影響が見られた。これに対し、光子線や中性子線では電子線に比べて線量の集中性が弱いため、眼球の大きさや変形による影響が小さいことがわかった。
長家 康展; 植木 太郎; 外池 幸太郎
Proceedings of 11th International Conference on Nuclear Criticality Safety (ICNC 2019) (Internet), 9 Pages, 2019/09
燃料デブリ体系に対する臨界安全解析のために新規モンテカルロ計算ソルバーSolomonを開発した。Solomonは、通常の臨界安全解析だけでなく、燃料デブリを含む損傷した原子炉の炉心の臨界安全解析もできるように設計されている。本論文では、Solomon開発の現状について述べ、乱雑化ワイエルシュトラス関数モデル, ボクセル形状を重ね合わせた乱雑化ワイエルシュトラス関数モデルの応用について紹介する。
Yeom, Y. S.*; Han, M. C.*; Choi, C.*; Han, H.*; Shin, B.*; 古田 琢哉; Kim, C. H.*
Health Physics, 116(5), p.664 - 676, 2019/05
被引用回数：6 パーセンタイル：77.19(Environmental Sciences)国際放射線防護委員会(ICRP)のタスクグループ103により、メッシュ形状の線量評価用人体ファントム(MRCPs)の開発が進められている。この人体ファントムは、将来的には線量評価で用いる標準人体モデルとして採用される予定である。そこで、このMRCPファントムに対するベンチマーク計算を主なモンテカルロ粒子輸送計算コード(Geant4, MCNP6およびPHITS)で行った。様々な粒子およびエネルギーで外部および内部被ばくの計算を実施し、計算時間やメモリ使用量をコード間で比較した。また、ボクセルファントムに対する計算も行い、コード毎の異なるメッシュ表現による性能の違いについて調べた。MRCPのメモリ使用量はGeant4およびMCNP6で10GB程度であったのに対し、PHITSでは1.2GBと顕著に少なかった。また、計算時間に関してもGeant4およびMCNP6ではボクセルファントムに比べてMRCPの計算時間は長くなる傾向を示したが、PHITSでは同程度もしくは短縮する傾向を示した。
Riyana, E. S.; 奥村 啓介; 寺島 顕一
Proceedings of 27th International Conference on Nuclear Engineering (ICONE-27) (Internet), 4 Pages, 2019/05
Determination of fuel debris location and distribution inside primary containment vessel (PCV) of Fukushima Daiichi Nuclear Power Station is important to decide further decommissioning step and strategy. Past measurements not yet provide enough information to determine fuel debris location and its distribution inside PCV. To support further measurements effort, we simulate detector response inside PCV based on calculated photon and neutron emission spectrum produced from fuel debris. The Calculation result could be use as basis for determination of suitable spectrometer system or detector for search, localized, define fuel debris distributions and treatment of fuel debris.
橋本 慎太郎; 佐藤 達彦
Journal of Nuclear Science and Technology, 56(4), p.345 - 354, 2019/04
被引用回数：1 パーセンタイル：24.18(Nuclear Science & Technology)モンテカルロ法に基づいた粒子輸送シミュレーションは、加速器施設の遮蔽計算等の様々な目的で利用されている。モンテカルロ法による計算結果の信頼性を定量的に評価するためには、試行回数によって決まる統計的不確かさに加えて、計算で使われる反応断面積や計算の入力情報となる遮蔽材の密度がもつ誤差の影響を系統的不確かさとして求める必要がある。本研究の評価方法は分散分析に基づいており、これを鉄やコンクリートを遮蔽材とする中性子の遮蔽計算の解析に適用し、統計的及び系統的不確かさの両方を評価可能であることを示した。その際、入力情報の値を誤差の範囲で変動させるランダム条件法の他、三条件(中央値,上限値及び下限値)のみ変動させる三条件法を提案した。ランダム条件法は計算体系に関わらず適切に系統的不確かさが評価できるものの長い計算時間を必要とするのに対し、入力情報の誤差の影響が複雑な場合を除いて、三条件法は計算時間を抑えてランダム条件法と同じ評価結果を与えることがわかった。さらに、試行回数を増加させた場合の収束状況を判断できる新しい基準値を示し、必要最小限の計算時間で収束した評価結果が得られることを明らかにした。
植木 太郎
Proceedings of International Conference on Mathematics and Computational Methods applied to Nuclear Science and Engineering (M&C 2019) (CD-ROM), p.151 - 160, 2019/00
物理・力学的な系は、極端な無秩序の下で、逆冪乗則パワースペクトルの状態に進展していくことが知られている。本著者は、確率的乱雑化ワイエルシュトラス関数による広範な逆冪乗則パワースペクトルの模擬を検討し、乱雑化無秩序媒質のモデル機能を、連続エネルギーモンテカルロ法ソルバーのデルタ追跡法を用いて検証した。数値計算例として、国際臨界安全ベンチマーク評価プロジェクトのTopsy炉心に基づく乱雑化体系の臨界性評価の揺らぎについて報告する。
木名瀬 栄; Mohammadi, A.*; G
mez-Ros, J.-M.*
Computational Anatomical Animal Models; Methodological Developments and Research Applications, p.5_1 - 5_9, 2018/12
被引用回数：0 パーセンタイル：0.04There are limited investigations on the computational frog models and the organ dose evaluations for frogs in environmental protection. In this article, computational frog models and their applications are reviewed to share some perspectives of frog model development in the near future. The authors hope that 3D printing frog phantoms with adequate tissue substitutes should be developed for the validation of the dosimetric quantities by the Monte Carlo simulations.
