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.
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.
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.
Kim, M.; Malins, A.; Sakuma, Kazuyuki; Kitamura, Akihiro; Machida, Masahiko; Hasegawa, Yukihiro*; Yanagi, Hideaki*
RIST News, (64), p.3 - 16, 2018/09
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.
ANS RPSD 2018; 20th Topical Meeting of the Radiation Protection and Shielding Division of ANS (CD-ROM), 5 Pages, 2018/08
Recently we introduced a function in PHITS to treat tetrahedral-mesh geometry. Tetrahedral-mesh geometry is a structure composed of combination of tetrahedrons and able to construct complex objects. Tetrahedral-mesh objects can be obtained by converting polygon data using mesh generation software such as TetGen. We also implemented a function in PHITS to export tally results into the format of the three-dimensional visualization software ParaView. TetGen is able to convert the polygon data into ParaView format. Together with these tools, three-dimensional analysis can be realized for PHITS simulation using a polygon objects.
Matsuda, Hiroki; Meigo, Shinichiro; Iwamoto, Hiroki
Journal of Nuclear Science and Technology, 55(8), p.955 - 961, 2018/08
We have started an experimental program to measure activation cross sections systematically in the proton-induced spallation reaction in structural materials commonly used in high-intensity proton accelerator-based facilities, such as Japan Proton Accelerator Research Complex (J-PARC). As the first step of the program, aluminum (Al) was chosen to verify the adequacy of the measurement technique implemented in a J-PARC proton beam environment because data of Al have been relatively well studied both by experimental measurement and simulation. Activation cross sections of Be, Na, and Na in Al were measured at proton energy points from 0.4, 1.3, 2.2 to 3.0 GeV, which could be delivered smoothly from the synchrotron. The validity of experimental data has been verified by introducing an effective proton numbers determination procedure. We compared the measured data with existing experimental data, the evaluated data (JENDL-HE/2007), and the calculations with several intra-nuclear cascade models by the Particle and Heavy Ion Transport code System (PHITS) code. Although the experimental data agreed with JENDL-HE/2007, the calculations underestimated about 40%. This could come from the evaporation model (generalized evaporation model) being implemented in the PHITS code. We found that the calculations agreed with the experimental data by an upgraded PHITS code.
Matsuda, Hiroki; Meigo, Shinichiro; Iwamoto, Hiroki
Journal of Physics; Conference Series, 1021(1), p.012017_1 - 012017_4, 2018/06
Spallation neutron at 180 degrees is of important for an evaluation of radiation protection for ADS (Accelerator-Driven System) and the nuclear physics. It was, however, quite difficult to measure it. We measured the energy spectrum of spallation neutron at 180 degrees at the proton transport beam line (3NBT) to MLF (Materials and Life Science Experimental Facility) on J-PARC by the NE213 liquid scintillator. The irradiated proton energy was 3 GeV, and the intensity was 110 protons above. The neutron energy was determined by Time-Of-Flight method with n-gamma discrimination. We also simulated the energy spectrum by using PHITS code and compared with measured spectrum. In this paper, the overview of the experiment and the results are described.
Hashimoto, Shintaro; Sato, Tatsuhiko; Iwamoto, Yosuke; Ogawa, Tatsuhiko; Furuta, Takuya; Abe, Shinichiro; Niita, Koji*
Kaku Deta Nyusu (Internet), (120), p.26 - 34, 2018/06
Particle and heavy-ion transport code system PHITS has been used for calculations of radiation shielding in accelerator facilities. PHITS describes physical phenomena induced by radiation as combination of transport and collision processes. The collision process including nuclear reactions is simulated by the three-step calculation: a generation of a reaction, pre-equilibrium, and compound processes. In the simulation, many physics models are used. This report explains roles of the models in PHITS and shows their developments we recently performed.
Meigo, Shinichiro; Matsuda, Hiroki; Iwamoto, Yosuke; Iwamoto, Hiroki; Hasegawa, Shoichi; Maekawa, Fujio; Yoshida, Makoto*; Ishida, Taku*; Makimura, Shunsuke*; Nakamoto, Tatsushi*
Proceedings of 9th International Particle Accelerator Conference (IPAC '18) (Internet), p.499 - 501, 2018/06
no abstracts in English
El Basha, D.*; Furuta, Takuya; Iyer, S. S. R.*; Bolch, W. E.*
Physics in Medicine and Biology, 63(10), p.105017_1 - 105017_13, 2018/05
With recent changes in the recommended annual limit on eye lens exposures to ionizing radiation by International Commission on Radiological Protection, there is considerable interest in predictive computational dosimetry models of the human eye and its various ocular structures. Several computational eye models to date have been constructed for this purpose but they are typically constructed of nominal size and of a roughly spherical shape associated with the emmetropic eye. We therefore constructed a geometric eye model that is both scalable (allowing for changes in eye size) and deformable (allowing for changes in eye shape), and that is suitable for use in radiation transport studies of ocular exposures and radiation treatments of eye disease. As an example, electron and photon anterior-posterior radiation transport with the constructed eye model was conducted and analyzed resultant energy-dependent dose profiles. Due to anterior-posterior irradiation, the energy dose response was shifted to higher energy for a larger-size eye or an axially deformed eye in prolate shape because the structures were located in deeper depth compared to the normal eye.
Han, M. C.*; Yeom, Y. S.*; Lee, H. S.*; Shin, B.*; Kim, C. H.*; Furuta, Takuya
Physics in Medicine and Biology, 63(9), p.09NT02_1 - 09NT02_9, 2018/05
The multi-threading computation performances of the Geant4, MCNP6, and PHITS codes were evaluated using three tetrahedral-mesh phantoms with different complexity. Photon and neutron transport simulations were conducted and the initialization time, calculation time, and memory usage were measured as a function of the number of threads N used in the simulation. The initialization time significantly increases with the complexity of the phantom, but not much with the number of the threads. For the calculation time, Geant4 showed good parallelization efficiency with multi-thread computation (30 times speed-up factor for N = 40) adopting the private tallies while saturation of the speed-up factor were observed in MCNP6 and PHITS (10 and a few times for N = 40) due to the time delay for the sharing tallies. On the other hand, Geant4 requires larger memory specification and the memory usage rapidly increases with the number of threads compared to MCNP6 or PHITS. It is notable that when compared to the other codes, the memory usage of PHITS is much smaller, regardless of both the complexity of the phantom and the number of the threads.
Meigo, Shinichiro; Matsuda, Hiroki; Iwamoto, Hiroki
Proceedings of 13th International Topical Meeting on Nuclear Applications of Accelerators (AccApp '17) (Internet), p.396 - 402, 2018/05
no abstracts in English
Abe, Shinichiro; Sato, Tatsuhiko
Meson, (47), p.34 - 39, 2018/03
Some researchers actively investigate secondary cosmic-ray muon radiography for imaging of large structures, detecting nuclear matter, and so on. The simulation is one of effective methods to estimate the required time for muon radiography or to optimize detection system. In this report, we introduce particle transport code PHITS and analytical model for estimating the atmospheric cosmic-ray spectra PARMA. Benchmark results are also shown. PHITS describes intensities of terrestrial muons in the ground well. PHITS also shows good agreement with measured data for particle generations from several interactions between muon and material. These results indicate the practicability of PHITS for simulation of muon radiography. PARMA reproduces cosmic-ray spectra on several spots well. Excellent agreement is seen between the measurements and PARMA for muon fluxes with large zenith angles and high energies, which are important for muon radiography.
Tsai, P.-E.; Lai, B.-L.*; Heilbronn, L. H.*; Sheu, R.-J.*
Nuclear Instruments and Methods in Physics Research B, 416, p.16 - 29, 2018/02
Fifteen thin target experiments were selected for this benchmark study of the neutron production cross sections. The studied cases include a mix combination of C, Ne Ar, Kr, and Xe ions bombarding Li, C, Al, Cu, and Pb target with projectile energies between 135 and 600 MeV/nucleon. The experimental data were compared to the model calculations performed by (1) PHITS version 2.73 with JQMD and GEM model, (2) PHITS version 2.82 with revised JQMD 2.0 and GEM model, (3) FLUKA version 2011.2c with RQMD 2.4 and FLUKA's own de-excitation model, and (4) MCNP6 version 1.0 with LAQGSM 03.03 and GEM2 model. This study provides useful information not only for code users but also for model developers, and it will lead to future improvements of the PHITS - JQMD model for heavy-ion induced reactions, which is important for accelerator facilities, heavy-ion radiotherapy, and space radiation.
Tada, Kenichi; Nagaya, Yasunobu; Kunieda, Satoshi; Suyama, Kenya; Fukahori, Tokio
EPJ Web of Conferences (Internet), 146, p.02028_1 - 02028_5, 2017/09
JAEA has started to develop new nuclear data processing system FRENDY (FRom Evaluated Nuclear Data libralY to any application). In this presentation, the outline of the development of FRENDY is presented. And functions and performances of FRENDY are demonstrated by generation and validation of the continuous energy cross section data libraries for MVP, PHITS and MCNP codes.