Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Riyana, E. S.; Okumura, Keisuke; Sakamoto, Masahiro; Matsumura, Taichi; Terashima, Kenichi
Journal of Nuclear Science and Technology, 59(4), p.424 - 430, 2022/04
Kubo, Kotaro; Jang, S.*; Takata, Takashi*; Yamaguchi, Akira*
Journal of Nuclear Science and Technology, 59(3), p.357 - 367, 2022/03
Times Cited Count：0Dynamic probabilistic risk assessment (PRA), which handles epistemic and aleatory uncertainties by coupling the thermal-hydraulics simulation and probabilistic sampling, enables a more realistic and detailed analysis than 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 station blackout sequence in a boiling water reactor and compared each method. The result indicated that quasi-Monte Carlo sampling method handles the uncertainties most effectively in the assumed scenario.
Furuta, Takuya
Igaku Butsuri, 41(4), P. 194, 2021/12
Number of medical uses of Particle and Heavy Ion Transport code System (PHITS) has been increased due to the recent high demands of medical use of radiations. The summary of such research works was described in the review article on medical application of Particle and Heavy Ion Transport code System PHITS published in Radiological Physics and Technology in 2021. There was a request from the editorial board of Japan Society of Medical Physics (JSMP) for writing an introductory article of this article in their internal journal. The research works on medical applications described in the review article, useful functions for medical application in PHITS, and newly opened user forum of PHITS have been introduced.
Furuta, Takuya; Sato, Tatsuhiko
Radiological Physics and Technology, 14(3), p.215 - 225, 2021/09
Number of the PHITS users has steadily increased since 2010 from when it is officially counted. Among them, increase of new users in medical physics is outstanding. Many research works in medical physics using PHITS have been published and the applications are widely spread in different fields such as applications to different types of radiotherapy, shielding calculations of medical facilities, application to radiation biology, and research and development of medical tools. In this article, we will introduce useful functions for medical application in PHITS by referring to examples of various medical applications.
Fujimoto, Nozomu*; Tada, Kenichi; Ho, H. Q.; Hamamoto, Shimpei; Nagasumi, Satoru; Ishitsuka, Etsuo
Annals of Nuclear Energy, 158, p.108270_1 - 108270_8, 2021/08
Times Cited Count：0 Percentile：0.01(Nuclear Science & Technology)Li, C.; Uchibori, Akihiro; Takata, Takashi; Pellegrini, M.*; Erkan, N.*; Okamoto, Koji*
Dai-25-Kai Doryoku, Enerugi Gijutsu Shimpojiumu Koen Rombunshu (Internet), 4 Pages, 2021/07
The capability of stable cooling and avoiding re-criticality on the debris bed are the main issues for achieving IVR (In-Vessel Retention). In the actual situation, the debris bed is composed of mixed-density debris particles. Hence, when these mixed-density debris particles were launched to re-distribute, the debris bed would possibly form a density-stratified distribution. For the proper evaluation of this scenario, the multi-physics model of CFD-DEM-Monte-Carlo based neutronics is established to investigate the coolability and re-criticality on the heterogeneous density-stratified debris bed with considering the particle relocation. The CFD-DEM model has been verified by utilizing water injection experiments on the mixed-density particle bed in the first portion of this research. In the second portion, the coupled system of the CFD-DEM-Monte-Carlo based neutronics model is applied to reactor cases. Afterward, the debris particles' movement, debris particles' and coolant's temperature, and the k-eff eigenvalue are successfully tracked. Ultimately, the relocation and stratification effects on debris bed's coolability and re-criticality had been quantitatively confirmed.
Ueki, Taro
Nuclear Science and Engineering, 195(2), p.214 - 226, 2021/02
Times Cited Count：1 Percentile：39.17(Nuclear Science & Technology)A dynamical system under extreme physical disorder has the tendency of evolving toward the equilibrium state characterized by an inverse power law power spectrum. In this paper, a practically implementable three-dimensional model is proposed for the random media formed by multi-materials mixture under such a power spectrum using a randomized form of Weierstrass function, its extension covering the white noise, and partial volumes pairings of constituent materials. The proposed model is implemented in the SOLOMON Monte Carlo solver with delta tracking. Two sets of numerical results are shown using the JENDL-4 nuclear data libraries.
Cs and 
Cs radioactivity measurementsMalins, A.; Imamura, Naohiro*; Niizato, Tadafumi; Takahashi, Junko*; Kim, M.; Sakuma, Kazuyuki; Shinomiya, Yoshiki*; Miura, Satoru*; Machida, Masahiko
Journal of Environmental Radioactivity, 226, p.106456_1 - 106456_12, 2021/01
Times Cited Count：0 Percentile：0(Environmental Sciences)Malins, A.; Ochi, Kotaro; Machida, Masahiko; Sanada, Yukihisa
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.147 - 154, 2020/10
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
JAEA-Technology-2020-003.pdf:1.5MB
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
Times Cited Count：6 Percentile：89.87(Instruments & Instrumentation)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
Times Cited Count：3 Percentile：20.11(Environmental Sciences)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
Times Cited Count：9 Percentile：78.97(Radiology, Nuclear Medicine & Medical Imaging)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.
-rays elastic scattering by Monte Carlo simulationOmer, M.; Hajima, Ryoichi*
New Journal of Physics (Internet), 21(11), p.113006_1 - 113006_10, 2019/11
Times Cited Count：4 Percentile：55.7(Physics, Multidisciplinary)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
Times Cited Count：1 Percentile：23.13(Environmental Sciences)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
Times Cited Count：7 Percentile：81.33(Environmental Sciences)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.
