Kubo, Kotaro; Zheng, X.; Tanaka, Yoichi; Tamaki, Hitoshi; Sugiyama, Tomoyuki; Jang, S.*; Takata, Takashi*; Yamaguchi, Akira*
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.308 - 315, 2020/10
Dynamic probabilistic risk assessment (PRA) is a method for improving the realism and completeness of conventional PRA. However, enormous calculation costs are incurred by these improvements. One solution is to select an appropriate sampling method. In this paper, we applied the Monte Carlo, Latin hypercube, grid-point, and quasi-Monte Carlo sampling methods to the dynamic PRA of a simplified accident sequence and compared the results for each method. Quasi-Monte Carlo sampling was found to be the most effective method in this case.
Harada, Hideo; Takayama, Naoki; Komeda, Masao
Journal of Physics Communications (Internet), 4(8), p.085004_1 - 085004_17, 2020/08
A new convention of epithermal neutron spectrum is formulated for improving accuracy of resonance integrals. The new type function is proposed as an approximating function of epithermal neutron spectrum based on calculations by the state-of-art Monte Carlo code MVP-3. Bias effects on determination of resonance integrals due to utilizing approximating functions of the traditional types and the new type are compared. The other bias effect is also investigated, which is caused by neglecting position dependence of a neutron spectrum inside an irradiation capsule. For demonstrating the bias effects due to these assumptions on neutron spectrum quantitatively in a practical case, the thermal neutron-capture cross section and resonance integral of Cs measured at a research reactor JRR-3 are re-evaluated. A superior property of the proposed new convention is discussed. The experimental method is proposed to determine the new shape factor introduced in the convention by a combinational use of triple flux monitors (Au, Co and Zr), and its analytical methodology is formulated.
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
Nagasumi, Satoru; Matsunaka, Kazuaki*; Fujimoto, Nozomu*; Ishii, Toshiaki; Ishitsuka, Etsuo
JAEA-Technology 2020-003, 13 Pages, 2020/05
The influence of the control rod model on the nuclear characteristics of the HTTR has been evaluated, by creating detailed control rod model, in which geometric shape was close to that of the actual control rod structure, in MVP code. According to refinement of the control rod model, the critical control rod position was 11 mm lower than that of the conventional model, and this was close to the measured value of 1775 mm. The reactivity absorbed by the shock absorber located at the tip of the control rod was 0.2%k/k, and this was 14 mm difference at the critical control rod position. Considering the effect of refinement of the control rod and the effect of the shock absorber, the correction amount for the analysis value in SRAC code due to the shape effect of the control rod, is -0.05%k/k in reactivity, and -3 mm in the critical control rod position at low temperature criticality.
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.
Omer, M.; Hajima, Ryoichi*
New Journal of Physics (Internet), 21(11), p.113006_1 - 113006_10, 2019/11
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.
Nagaya, Yasunobu; Ueki, Taro; Tonoike, Kotaro
Proceedings of 11th International Conference on Nuclear Criticality Safety (ICNC 2019) (Internet), 9 Pages, 2019/09
A new Monte Carlo solver Solomon has been developed for the application to fuel-debris systems. It is designed not only for usual criticality safety analysis but also for criticality calculations of damaged reactor core including fuel debris. This paper describes the current status of Solomon and demonstrates the applications of the randomized Weierstrass function (RWF) model and the RWF model superimposed voxel geometry.
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.
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
A dynamical system under extreme physical disorder has the tendency of evolving toward the equilibrium state characterized by an inverse power law spectrum. In this paper, the author proposes a practically implementable modeling of random media under such a spectrum using a randomized form of the Weierstrass function. The proposed modeling is demonstrated by the continuous energy Monte Carlo particle transport with delta tracking for the criticality calculation of a randomized version of the Topsy spherical core in International Criticality Safety Benchmark Evaluation Project.
Rossi, F.; Rodriguez, D.; Takahashi, Toon; Seya, Michio; Koizumi, Mitsuo
Nihon Kaku Busshitsu Kanri Gakkai Dai-39-Kai Nenji Taikai Rombunshu (CD-ROM), 3 Pages, 2018/11
This paper (presentation) reports on the design study for seeking compactness of moderator (and reflector) structure for the Delayed Gamma-ray Spectroscopy Non Destructive Assay system to be applied for verification of fissile isotopic (Pu-239, Pu-241 and U-235) composition ratio in High Radioactive Nuclear Material sample such as dissolved solution of reprocessing facility. We first investigated advantages and disadvantages of using different compact neutron sources as DT and DD neutron generators. This led us to the conclusion of using DD sources to achieve a compact (and practical) NDA system. Based on MCNP simulations, with the use of a DD neutron generator, a preliminary optimization study was performed and will be shown here.
Rodriguez, D.; Rossi, F.; Seya, Michio; Koizumi, Mitsuo
Proceedings of 2017 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC 2017) (Internet), 3 Pages, 2018/11
El Basha, D.*; Furuta, Takuya; Iyer, S. S. R.*; Bolch, W. E.*
Physics in Medicine & 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 & 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.
Takada, Shusuke*; Okudaira, Takuya*; Goto, Fumiya*; Hirota, Katsuya*; Kimura, Atsushi; Kitaguchi, Masaaki*; Koga, Jun*; Nakao, Taro*; Sakai, Kenji; Shimizu, Hirohiko*; et al.
Journal of Instrumentation (Internet), 13(2), p.P02018_1 - P02018_21, 2018/02