Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Tanaka, Minori*; Watabane, Masashi*; Machida, Masahiko; Yamada, Susumu; Enomoto, Yota*; Gunji, Kota*; Arikawa, Taro*
Doboku Gakkai Rombunshu, B2 (Kaigan Kogaku) (Internet), 76(2), p.I_103 - I_108, 2020/11
no abstracts in English
Nagao, Fumiya; Niizato, Tadafumi; Sasaki, Yoshito; Ito, Satomi; Watanabe, Takayoshi; Dohi, Terumi; Nakanishi, Takahiro; Sakuma, Kazuyuki; Hagiwara, Hiroki; Funaki, Hironori; et al.
JAEA-Research 2020-007, 249 Pages, 2020/10
JAEA-Research-2020-007.pdf:15.83MB
The accident of the Fukushima Daiichi Nuclear Power Station, Tokyo Electric Power Company Holdings, Inc. occurred due to the Great East Japan Earthquake, Sanriku offshore earthquake, of 9.0 magnitude and the accompanying tsunami. As a result, large amount of radioactive materials was released into the environment. Under these circumstances, Japan Atomic Energy Agency (JAEA) has been conducting "Long-term Assessment of Transport of Radioactive Contaminants in the Environment of Fukushima" concerning radioactive materials released in environment, especially migration behavior of radioactive cesium since November 2012. This report is a summary of the research results that have been obtained in environmental dynamics research conducted by JAEA in Fukushima Prefecture.
Yamada, Susumu; Machida, Masahiko; Kurikami, Hiroshi
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.140 - 146, 2020/10
no abstracts in English
Idomura, Yasuhiro; Onodera, Naoyuki; Yamada, Susumu; Yamashita, Susumu; Ina, Takuya*; Imamura, Toshiyuki*
Supa Kompyuthingu Nyusu, 22(5), p.18 - 29, 2020/09
A communication avoiding multigrid preconditioned conjugate gradient method (CAMGCG) is applied to the pressure Poisson equation in a multiphase CFD code JUPITER, and its computational performance and convergence property are compared against the conventional Krylov methods. The CAMGCG solver has robust convergence properties regardless of the problem size, and shows both communication reduction and convergence improvement, leading to higher performance gain than CA Krylov solvers, which achieve only the former. The CAMGCG solver is applied to extreme scale multiphase CFD simulations with 90 billion DOFs, and its performance is compared against the preconditioned CG solver. In this benchmark, the number of iterations is reduced to 
, and 
speedup is achieved with keeping excellent strong scaling up to 8,000 nodes on the Oakforest-PACS.
Machida, Masahiko; Yamada, Susumu; Iwata, Ayako; Otosaka, Shigeyoshi; Kobayashi, Takuya; Watanabe, Masahisa; Funasaka, Hideyuki; Morita, Takami*
Journal of Nuclear Science and Technology, 57(8), p.939 - 950, 2020/08
Times Cited Count：0 Percentile：100(Nuclear Science & Technology)After direct discharges of highly contaminated water of the Fukushima Daiichi Nuclear Power Plant (1F) from April to May 2011, Kanda suggested that relatively small amounts of run-off of radionuclides from the 1F port into the Fukushima coastal region subsequently continued by his estimation method. However, the estimation period was limited to up to September 2012. Therefore, this paper estimates the discharge inventory up to June 2018. In the missing period, the Japanese government and Tokyo Electric Power Company Holdings have continued efforts to stop the discharge, and consequently, the radionuclide concentration in seawater inside the 1F port has gradually diminished. We show the monthly discharge inventory of 
Cs up to June 2018 by two methods, i.e., Kanda method partially improved by the authors and a more sophisticated method using Voronoi tessellation reflecting the increase in the number of monitoring points inside the 1 F port. The results show that the former always yields overestimated results compared with the latter, but the ratio of the former to the latter is less than one order of magnitude. Using these results, we evaluate the impact of the discharge inventory from the 1F port into the coastal area and radiation dose upon fish ingestion.
Yamada, Susumu; Machida, Masahiko; Imamura, Toshiyuki*
Parallel Computing; Technology Trends, p.105 - 113, 2020/00
no abstracts in English
Matsumoto, Kazuya*; Idomura, Yasuhiro; Ina, Takuya*; Mayumi, Akie; Yamada, Susumu
Journal of Supercomputing, 75(12), p.8115 - 8146, 2019/12
Times Cited Count：0 Percentile：100(Computer Science, Hardware & Architecture)A communication-avoiding generalized minimum residual method (CA-GMRES) is implemented on a hybrid CPU-GPU cluster, targeted for the performance acceleration of iterative linear system solver in the gyrokinetic toroidal five-dimensional Eulerian code GT5D. In addition to the CA-GMRES, we implement and evaluate a modified variant of CA-GMRES (M-CA-GMRES) proposed in our previous study to reduce the amount of floating-point calculations. This study demonstrates that beneficial features of the CA-GMRES are in its minimum number of collective communications and its highly efficient calculations based on dense matrix-matrix operations. The performance evaluation is conducted on the Reedbush-L GPU cluster, which contains four NVIDIA Tesla P100 GPUs per compute node. The evaluation results show that the M-CA-GMRES is 1.09x, 1.22x and 1.50x faster than the CA-GMRES, the generalized conjugate residual method (GCR), and the GMRES, respectively, when 64 GPUs are used.
Yamada, Susumu
Koyuchi Keisan To Tokuichi Keisan; Keisan Rikigaku Rekucha Kosu, p.77 - 104, 2019/12
no abstracts in English
Machida, Masahiko; Yamada, Susumu; Iwata, Ayako; Otosaka, Shigeyoshi; Kobayashi, Takuya; Watanabe, Masahisa; Funasaka, Hideyuki; Morita, Takami*
Nippon Genshiryoku Gakkai Wabun Rombunshi, 18(4), p.226 - 236, 2019/12
After direct discharges of highly-contaminated water from Unit 2 and 3 in Fukushima Daiichi Nuclear Power Plant (1F) in April to May 2011, Kanda suggested that relatively small run-off of radionuclides from 1F port into Fukushima coastal region has subsequently continued by using his estimation scheme. However, the estimation period was limited until September 2012, and there has been no report on the issue since the work. Therefore, this paper focuses on discharge inventory from 1F port until June 2018. In the missing period, the central government and Tokyo Electric Power Company Holdings have done continuous efforts to stop the discharge, and consequently sea water concentration inside 1F port has diminished gradually. We show monthly discharge inventory of Cs-137 until June 2018 by two schemes, i.e., Kanda's scheme partially improved by authors and more sophisticated one using Voronoi tessellation reflecting the increment of the number of monitoring points inside 1F port. The results show that the former always presents overestimated results compared to the latter but the ratio of former to latter is less than one order. Based on these results, we evaluate impact of discharge inventory from 1F port into the coastal area and radiation does via fish digestion.
Nagao, Fumiya; Niizato, Tadafumi; Sasaki, Yoshito; Ito, Satomi; Watanabe, Takayoshi; Dohi, Terumi; Nakanishi, Takahiro; Sakuma, Kazuyuki; Hagiwara, Hiroki; Funaki, Hironori; et al.
JAEA-Research 2019-002, 235 Pages, 2019/08
JAEA-Research-2019-002.pdf:21.04MB
The accident of the Fukushima Daiichi Nuclear Power Station (hereinafter referred to 1F), Tokyo Electric Power Company Holdings, Inc. occurred due to the Great East Japan Earthquake, Sanriku offshore earthquake, of 9.0 magnitude and the accompanying tsunami. As a result, large amount of radioactive materials was released into the environment. Under these circumstances, JAEA has been conducting Long-term Environmental Dynamics Research concerning radioactive materials released in environment, especially migration behavior of radioactive cesium since November 2012. This report is a summary of the research results that have been obtained in environmental dynamics research conducted by JAEA in Fukushima Prefecture.
Wang, H.*; Otsu, Hideaki*; Chiga, Nobuyuki*; Kawase, Shoichiro*; Takeuchi, Satoshi*; Sumikama, Toshiyuki*; Koyama, Shumpei*; Sakurai, Hiroyoshi*; Watanabe, Yukinobu*; Nakayama, Shinsuke; et al.
Communications Physics (Internet), 2(1), p.78_1 - 78_6, 2019/07
Times Cited Count：4 Percentile：25.43(Physics, Multidisciplinary)Searching for effective pathways for the production of proton- and neutron-rich isotopes through an optimal combination of reaction mechanism and energy is one of the main driving forces behind experimental and theoretical nuclear reaction studies as well as for practical applications in nuclear transmutation of radioactive waste. We report on a study on incomplete fusion induced by deuteron, which contains one proton and one neutron with a weak binding energy and is easily broken up. This reaction study was achieved by measuring directly the cross sections for both proton and deuteron for 
Pd at 50 MeV/u via inverse kinematics technique. The results provide direct experimental evidence for the onset of a cross-section enhancement at high energy, indicating the potential of incomplete fusion induced by loosely-bound nuclei for creating proton-rich isotopes and nuclear transmutation of radioactive waste.
Idomura, Yasuhiro; Ina, Takuya*; Yamashita, Susumu; Onodera, Naoyuki; Yamada, Susumu; Imamura, Toshiyuki*
Proceedings of 9th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2018) (Internet), p.17 - 24, 2018/11
Times Cited Count：0 Percentile：100A communication avoiding (CA) multigrid preconditioned conjugate gradient method (CAMGCG) is applied to the pressure Poisson equation in a multiphase CFD code JUPITER, and its computational performance and convergence property are compared against CA Krylov methods. In the JUPITER code, the CAMGCG solver has robust convergence properties regardless of the problem size, and shows both communication reduction and convergence improvement, leading to higher performance gain than CA Krylov solvers, which achieve only the former. The CAMGCG solver is applied to extreme scale multiphase CFD simulations with 
billion DOFs, and it is shown that compared with a preconditioned CG solver, the number of iterations is reduced to , and speedup is achieved with keeping excellent strong scaling up to 8,000 nodes on the Oakforest-PACS.
Kamikubota, Norihiko*; Yamada, Shuei*; Sato, Kenichiro*; Kikuzawa, Nobuhiro; Yamamoto, Noboru*; Yoshida, Susumu*; Nemoto, Hiroyuki*
Proceedings of 16th International Conference on Accelerator and Large Experimental Physics Control Systems (ICALEPCS 2017) (Internet), p.1470 - 1473, 2018/01
no abstracts in English
Yamada, Susumu; Imamura, Toshiyuki*; Machida, Masahiko
Lecture Notes in Computer Science 10776, p.243 - 256, 2018/00
Times Cited Count：0 Percentile：100no abstracts in English
Idomura, Yasuhiro; Ina, Takuya*; Mayumi, Akie; Yamada, Susumu; Imamura, Toshiyuki*
Lecture Notes in Computer Science 10776, p.257 - 273, 2018/00
A preconditioned Chebyshev basis communication-avoiding conjugate gradient method (P-CBCG) is applied to the pressure Poisson equation in a multiphase thermal-hydraulic CFD code JUPITER, and its computational performance and convergence properties are compared against a preconditioned conjugate gradient (P-CG) method and a preconditioned communication-avoiding conjugate gradient (P-CACG) method on the Oakforest-PACS, which consists of 8,208 KNLs. The P-CBCG method reduces the number of collective communications with keeping the robustness of convergence properties. Compared with the P-CACG method, an order of magnitude larger communication-avoiding steps are enabled by the improved robustness. It is shown that the P-CBCG method is 
and 
faster than the P-CG and P-CACG methods at 2,000 processors, respectively.
Yamada, Susumu; Imamura, Toshiyuki*; Machida, Masahiko
Parallel Computing is Everywhere, p.27 - 36, 2018/00
no abstracts in English
Idomura, Yasuhiro; Ina, Takuya*; Mayumi, Akie; Yamada, Susumu; Matsumoto, Kazuya*; Asahi, Yuichi*; Imamura, Toshiyuki*
Proceedings of 8th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2017), p.7_1 - 7_8, 2017/11
A communication-avoiding generalized minimal residual (CA-GMRES) method is applied to the gyrokinetic toroidal five dimensional Eulerian code GT5D, and its performance is compared against the original code with a generalized conjugate residual (GCR) method on the JAEA ICEX (Haswell), the Plasma Simulator (FX100), and the Oakforest-PACS (KNL). The CA-GMRES method has 
higher arithmetic intensity than the GCR method, and thus, is suitable for future Exa-scale architectures with limited memory and network bandwidths. In the performance evaluation, it is shown that compared with the GCR solver, its computing kernels are accelerated by 
, and the cost of data reduction communication is reduced from 
to 
of the total cost at 1,280 nodes.
Yamada, Susumu; Ina, Takuya*; Sasa, Narimasa; Idomura, Yasuhiro; Machida, Masahiko; Imamura, Toshiyuki*
Proceedings of 2017 IEEE International Parallel & Distributed Processing Symposium Workshops (IPDPSW) (Internet), p.1418 - 1425, 2017/08
Times Cited Count：3 Percentile：21.83no abstracts in English
Kobayashi, Keita*; Okumura, Masahiko; Yamada, Susumu; Machida, Masahiko; Aoki, Hideo*
Physical Review B, 94(21), p.214501_1 - 214501_7, 2016/12
Times Cited Count：24 Percentile：16.9(Materials Science, Multidisciplinary)no abstracts in English
Mayumi, Akie; Idomura, Yasuhiro; Ina, Takuya; Yamada, Susumu; Imamura, Toshiyuki*
Proceedings of 7th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2016) (Internet), p.17 - 24, 2016/11
The left-preconditioned communication avoiding conjugate gradient (LP-CA-CG) method is applied to the pressure Poisson equation in the multiphase CFD code JUPITER. The arithmetic intensity of the LP-CA-CG method is analyzed, and is dramatically improved by loop splitting for inner product operations and for three term recurrence operations. Two LP-CA-CG solvers with block Jacobi preconditioning and with underlap preconditioning are developed. It is shown that on the K computer, the LP-CA-CG solvers with block Jacobi preconditioning is faster, because the performance of local point-to-point communications scales well, and the convergence property becomes worse with underlap preconditioning. The LP-CA-CG solver shows good strong scaling up to 30,000 nodes, where the LP-CA-CG solver achieved higher performance than the original CG solver by reducing the cost of global collective communications by 69%.
