Riyana, E. S.; Okumura, Keisuke; Terashima, Kenichi; Matsumura, Taichi; Sakamoto, Masahiro
Mechanical Engineering Journal (Internet), 7(3), p.19-00543_1 - 19-00543_8, 2020/06
Iwamoto, Hiroki; Meigo, Shinichiro
Journal of Nuclear Science and Technology, 57(3), p.276 - 290, 2020/03
This paper presents an approach to uncertainty estimation of spallation particle multiplicity of lead (Pb), primarily focusing on proton-induced spallation neutron multiplicity () and its propagation to a neutron energy spectrum. The uncertainty is estimated from experimental proton-induced neutron-production double-differential cross sections (DDXs) and model calculations with the Particle and Heavy Ion Transport code System (PHITS). Uncertainties in multiplicities for , , and reactions are then inferred from the estimated uncertainty and the PHITS calculation. Using these uncertainties, uncertainty in a neutron energy spectrum produced from a thick Pb target bombarded with 500 MeV proton beams, measured in a previous experiment, is quantified by a random sampling technique, and propagation to the neutron energy spectrum is examined. Relatively large uncertainty intervals (UIs) were observed outside the lower limit of the measurement range, which is prominent in the backward directions. Our findings suggest that a reliable assessment of spallation neutron energy spectra requires systematic DDX experiments for detector angles and incident energies below 100 MeV as well as neutron energy spectrum measurements at lower energies below 1.4 MeV with an accuracy below the quantified UIs.
Meigo, Shinichiro; Matsuda, Hiroki; Iwamoto, Yosuke; Yoshida, Makoto*; Hasegawa, Shoichi; Maekawa, Fujio; Iwamoto, Hiroki; Nakamoto, Tatsushi*; Ishida, Taku*; Makimura, Shunsuke*
JPS Conference Proceedings (Internet), 28, p.061004_1 - 061004_6, 2020/02
no abstracts in English
Malins, A.; Machida, Masahiko; Vu, TheDang; Aizawa, Kazuya; Ishida, Takekazu*
Nuclear Instruments and Methods in Physics Research A, 953, p.163130_1 - 163130_7, 2020/02
Kim, M.; Malins, A.; Yoshimura, Kazuya; Sakuma, Kazuyuki; Kurikami, Hiroshi; Kitamura, Akihiro; Machida, Masahiko; Hasegawa, Yukihiro*; Yanagi, Hideaki*
Journal of Environmental Radioactivity, 210, p.105803_1 - 105803_10, 2019/12
To improve the accuracy of simulations for air dose rates over fallout contaminated areas, the distribution of the radionuclides within the environment should be modelled realistically, e.g. considering differences in radioactivity levels between agricultural land, urban surfaces, and forest compartments. Moreover simulations should model the shielding of rays by buildings, trees and land topography. Here we outline a system for generating three dimensional models of urban and rural areas in Fukushima Prefecture. The Cs and Cs radioactivity distribution can be set flexibly across the different components of the model. The models incorporate realistic representations of local buildings, based on nine common Japanese designs, individual conifer and broadleaf trees, and the topography of the land surface. Models are generated from Digital Elevation Model (DEM) and Digital Surface Model (DSM) datasets, and refined by users assisted with ortho-photographs of target sites. Completed models are exported from the system in a format suitable for the Particle and Heavy Ion Transport code System (PHITS) for the calculation of air dose rates and other radiological quantities. The system is demonstrated by modelling a suburban area 4 km from the Fukushima Daiichi Nuclear Power Plant that has yet to be decontaminated. Air dose rates calculated in PHITS were correlated with measurements taken across the site in a car-borne survey.
Carter, L. M.*; Crawford, T. M.*; Sato, Tatsuhiko; Furuta, Takuya; 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
Voxel human phantoms have been used for internal dose assessment. More anatomically accurate representation become possible for skins or layer tissues owing to recent developments of advanced polygonal mesh-type phantoms and thus internal dose assessment using those advanced phantoms are desired. However, the Monte Carlo transport calculation by implementing those phantoms require an advanced knowledge for the Monte Carlo transport codes and it is only limited to experts. We therefore developed a tool, PARaDIM, which enables users to conduct internal dose calculation with PHITS easily by themselves. With this tool, a user can select tetrahedral-mesh phantoms, set radionuclides in organs, and execute radiation transport calculation with PHITS. Several test cases of internal dosimetry calculations were presented and usefulness of this tool was demonstrated.
Tsai, P.-E.; Iwamoto, Yosuke; Hagiwara, Masayuki*
Transactions of the American Nuclear Society, 121(1), p.13 - 16, 2019/11
The importance of PKA characteristics to radiation damage assessment for proton accelerators has been demonstrated in the paper by the PHITS simulations. In order to obtain accurate results of displacement cross section calculated by PHITS, it is therefore critical to have the physics models implemented in PHITS be properly validated by experimental data. With the new measurement system, which has been being developed here at JAEA, it is expected to provide experimental PKA data with low measurement thresholds and good mass resolutions. The test result we obtained at CYRIC, Tohoku University, for the - gas ionization chamber as a part of the new measurement system was very promising. We were able to successfully distinguish C, B, Be, Li, and He elements with threshold energies lower than 1 MeV/nucleon in the configuration of 70-MeV proton beam and 920-nm thick C target. More experiments and tests shall continue in the future with heavier targets and higher proton beam energies.
Kim, M.; Malins, A.; Sakuma, Kazuyuki; Kitamura, Akihiro; Machida, Masahiko; Hasegawa, Yukihiro*; Yanagi, Hideaki*
Isotope News, (765), p.30 - 33, 2019/10
Here we outline a system for generating three dimensional models of urban and rural areas in Fukushima Prefecture. The Cs and Cs radioactivity distribution can be set flexibly across the different components of the model. The models incorporate realistic representations of local buildings, individual conifer and broadleaf trees, and the topography of the land surface. The system is demonstrated by modelling a suburban area 4 km from the Fukushima Daiichi Nuclear Power Plant that has yet to be decontaminated. Air dose rates calculated in PHITS were correlated with measurements taken across the site in a car-borne survey.
Miyahara, Shinya*; Ohdaira, Naoya*; Arita, Yuji*; Maekawa, Fujio; Matsuda, Hiroki; Sasa, Toshinobu; Meigo, Shinichiro
Nuclear Engineering and Design, 352, p.110192_1 - 110192_8, 2019/10
Lead-Bismuth Eutectic (LBE) is used as a spallation neutron target and coolant materials of Accelerator Driven System (ADS), and many kinds of elements are produced as spallation products. It is important to evaluate the release and transport behavior of the spallation products in the LBE. The inventories and the physicochemical composition of the spallation products produced in LBE have been investigated for an LBE loop in the ADS Target Test Facility (TEF-T) in J-PARC. The inventories of the spallation products in the LBE were estimated using the PHITS code. The physicochemical composition of the spallation products in the LBE was calculated using the Thermo-Calc code under the conditions of the operation temperatures of LBE from 350C to 500C and the oxygen concentrations in LBE from 10 ppb to 1 ppm. The calculation showed that the 5 elements of Rb, Tl, Tc, Os, Ir, Pt, Au and Hg were soluble in LBE under the all given conditions and any kinds of compound were not formed in LBE. It was suggested that the oxides of Ce, Sr, Zr and Y were stable as CeO, SrO, ZrO and YO in the LBE.
Furuta, Takuya; 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
Despite large variation of human eye, only one computational eye model has been adopted in almost all the radiation transport simulation studies. We thus adopted a new scalable and deformable eye model and studied the radiation exposure by electrons, photons, and neutrons in the standard radiation fields such as AP, PA, RLAT, ROT, by using Monte Carlo radiation transport code PHITS. We computed the radiation exposure for 5 eye models (standard, large, small, myopic, hyperopic) and analyzed influence of absorbed dose in ocular structures on eye size and shape. Dose distribution of electrons is conformal and therefore the absorbed doses in ocular structures depend on the depth location of each ocular structure. We thus found a significant variation of the absorbed doses for each ocular structure for electron exposure due to change of the depth location affected by eye size and shape. On the other hand only small variation was observed for photons and neutrons exposures because of less conformal dose distribution of those particles.
Abe, Shinichiro; Liao, W.*; Manabe, Seiya*; Sato, Tatsuhiko; Hashimoto, Masanori*; Watanabe, Yukinobu*
IEEE Transactions on Nuclear Science, 66(7, Part2), p.1374 - 1380, 2019/07
Single event upsets (SEUs) caused by secondary cosmic-ray neutrons have recognized as a serious reliability problem for microelectronic devices. Acceleration tests at neutron facilities are convenient to validate soft error rates (SERs) quickly, but some corrections caused from measurement conditions are required to derive realistic SERs at actual environment or to compare other measured data. In this study, the effect of irradiation side on neutron-induced SEU cross sections was investigated by performing neutron transport simulation using PHITS. SERs for 65-nm bulk CMOS SRAMs are estimated using the sensitive volume model. It was found from simulation that SERs for the sealant side irradiation are 30-50% larger than those for the board side irradiation. This difference comes from the difference of production yield and angular distribution of secondary H and He ions, which are the main cause of SEUs. Thus the direction of neutron irradiation should be reported when the result of acceleration tests are published. This result also indicates that SERs can be reduced by equipping device with sealant side facing downward.
Maeda, Makoto; Furutaka, Kazuyoshi; Kureta, Masatoshi; Ozu, Akira; Komeda, Masao; Toh, Yosuke
Journal of Nuclear Science and Technology, 56(7), p.617 - 628, 2019/07
Tsai, P.-E.; Heilbronn, L. H.*; Lai, B.-L.*; Iwata, Yoshiyuki*; Murakami, Takeshi*; Sheu, R.-J.*
Nuclear Instruments and Methods in Physics Research B, 449, p.62 - 70, 2019/06
The secondary neutrons produced from 100- and 230-MeV/nucleon He ions, respectively, stopping in the thick iron, PMMA and water targets are measured by the time of flight method combined with the pulse shape discrimination of liquid scintillators. The experimental data show that the secondary neutrons were contributed from breakup of projectiles, emission from overlap regions of projectile nuclei and target nuclei, and evaporation of projectiles and target nuclei. The measured double-differential thick target neutron yields, angular distributions, and total neutron yields per ion, were benchmarked by the PHITS, FLUKA, and MCNP model calculations. The default models in these codes agree the experimental data well at intermediate-to-large angles in the low-to-intermediate energy range. However, the physics models implemented in PHITS need further improvement for some particular nuclear interaction mechanisms, and this work can be an importance reference for future model development.
Yeom, Y. S.*; Han, M. C.*; Choi, C.*; Han, H.*; Shin, B.*; Furuta, Takuya; Kim, C. H.*
Health Physics, 116(5), p.664 - 676, 2019/05
Recently, Task Group 103 of the ICRP developed the mesh-type reference computational phantoms (MCRPs), which are planned for use in future ICRP dose coefficient calculation. Performance of major Monte Carlo particle transport codes (Geant4, MCNP6, and PHITS) were tested with MCRP. External and internal exposure of various particles and energies were calculated and the computational times and required memories were compared. Additionally calculation for voxel-mesh phantom was also conducted so that the influence of different mesh-representation in each code was studied. Memory usage of MRCP was as large as 10 GB with Geant4 and MCNP6 while it is much less with PHITS (1.2 GB). In addition, the computational time required for MRCP tends to increase compared to voxel-mesh phantoms with Geant4 and MCNP6 while it is equal or tends to decrease with PHITS.
Riyana, E. S.; Okumura, Keisuke; Terashima, Kenichi
Proceedings of 27th International Conference on Nuclear Engineering (ICONE-27) (Internet), 4 Pages, 2019/05
Hashimoto, Shintaro; Sato, Tatsuhiko
Journal of Nuclear Science and Technology, 56(4), p.345 - 354, 2019/04
Particle transport simulations based on the Monte Carlo method have been applied to shielding calculations. Estimation of not only statistical uncertainty related to the number of trials but also systematic one induced by unclear physical quantities is required to confirm the reliability of calculated results. In this study, we applied a method based on analysis of variance to shielding calculations. We proposed random- and three-condition methods. The first one determines randomly the value of the unclear quantity, while the second one uses only three values: the default value, upper and lower limits. The systematic uncertainty can be estimated adequately by the random-condition method, though it needs the large computational cost. The three-condition method can provide almost the same estimate as the random-condition method when the effect of the variation is monotonic. We found criterion to confirm convergence of the systematic uncertainty as the number of trials increases.
Iwamoto, Hiroki; Meigo, Shinichiro
Journal of Nuclear Science and Technology, 56(2), p.160 - 171, 2019/02
We present a new model to describe the fission probability of the high-energy fission model, as deduced from the intranuclear cascade calculation with the Intra-Nuclear Cascade model of Lige (INCL) version 4.6 and Prokofiev's phenomenological systematics of the proton-induced fission cross sections. This model is implemented in the de-excitation model of the Generalized Evaporation Model (GEM), and applied to Monte Carlo spallation reaction simulation using the Particle and Heavy Ion Transport code System (PHITS). Comparing with experimental data for subactinide nuclei shows that this model can provide a unified prediction of the proton-, neutron-, and deuteron-induced fission cross sections with markedly improved accuracy. The calculated fission fragments tend to shift to higher mass numbers. To account for the isotopic distributions of fission fragments within the framework of a coupled INCL/GEM, modification of INCL is required, especially for description of the highly-excited states of residual nuclei.
Ishitsuka, Etsuo; Kenzhina, I.*; Okumura, Keisuke; Ho, H. Q.; Takemoto, Noriyuki; Chikhray, Y.*
JAEA-Technology 2018-010, 33 Pages, 2018/11
As a part of study on the mechanism of tritium release to the primary coolant in research and testing reactors, tritium recoil release rate from Li and U impurities in the neutron reflector made by beryllium, aluminum and graphite were calculated by PHITS code. On the other hand, the tritium production from Li and U impurities in beryllium neutron reflectors for JMTR and JRR-3M were calculated by MCNP6 and ORIGEN2 code. By using both results, the amount of recoiled tritium from beryllium neutron reflectors were estimated. It is clear that the amount of recoiled tritium from Li and U impurities in beryllium neutron reflectors are negligible, and 2 and 5 orders smaller than that from beryllium itself, respectively.
Tsai, P.-E.; Iwamoto, Yosuke; Hagiwara, Masayuki*; Sato, Tatsuhiko; Ogawa, Tatsuhiko; Satoh, Daiki; Abe, Shinichiro; Ito, Masatoshi*; Watabe, Hiroshi*
Proceedings of 2017 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC 2017) (Internet), 3 Pages, 2018/11
The energy spectra of primary knock-on atoms (PKAs) are essential for radiation damage assessment in design of accelerator facilities. However up to date the experimental data are still limited, due to the poor mass resolution and the high measurement threshold energies in the conventional setup of nuclear physics experiments using solid state detectors, which are typically above a few MeV/nucleon. In this study, a novel detection system consisting of two time detectors and one dE-E energy detector is proposed and being constructed to measure the PKA spectra. The system and detector design was based on Monte Carlo simulations by using the PHITS code. The PHITS simulations show that the system is able to distinguish the PKA isotopes above 0.2-0.3 MeV/nucleon for A=2030 amu; the PKA mass identification thresholds decrease to 0.1 MeV/nucleon for PKAs lighter than 20 amu. The detection system will be tested in the summer of 2017, and the test results will be presented at the conference.