Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Asahi, Yuichi; Onodera, Naoyuki; Hasegawa, Yuta; Shimokawabe, Takashi*; Shiba, Hayato*; Idomura, Yasuhiro
Boundary-Layer Meteorology, 186(3), p.659 - 692, 2023/03
Times Cited Count：0 Percentile：0.01(Meteorology & Atmospheric Sciences)We develop a Transformer-based deep learning model to predict the plume concentrations in the urban area under uniform flow conditions. Our model has two distinct input layers: Transformer layers for sequential data and convolutional layers in convolutional neural networks (CNNs) for image-like data. Our model can predict the plume concentration from realistically available data such as the time series monitoring data at a few observation stations and the building shapes and the source location. It is shown that the model can give reasonably accurate prediction with orders of magnitude faster than CFD simulations. It is also shown that the exactly same model can be applied to predict the source location, which also gives reasonable prediction accuracy.
Onodera, Naoyuki; Idomura, Yasuhiro; Hasegawa, Yuta; Shimokawabe, Takashi*; Aoki, Takayuki*
Keisan Kogaku Koenkai Rombunshu (CD-ROM), 27, 4 Pages, 2022/06
We have developed a wind simulation code named CityLBM to realize wind digital twins. Mesoscale wind conditions are given as boundary conditions in CityLBM by using a nudging data assimilation method. It is found that conventional approaches with constant nudging coefficients fail to reproduce turbulent intensity in long time simulations, where atmospheric stability conditions change significantly. We propose a dynamic parameter optimization method for the nudging coefficient based on a particle filter. CityLBM was validated against plume dispersion experiments in the complex urban environment of Oklahoma City. The nudging coefficient was updated to reduce the error of the turbulent intensity between the simulation and the observation, and the atmospheric boundary layer was reproduced throughout the day.
Asahi, Yuichi; Onodera, Naoyuki; Hasegawa, Yuta; Shimokawabe, Takashi*; Shiba, Hayato*; Idomura, Yasuhiro
Keisan Kogaku Koenkai Rombunshu (CD-ROM), 27, 5 Pages, 2022/06
We have ported the GPU accelerated Lattice Boltzmann Method code "CityLBM" to AMD MI100 GPU. We present the performance of CityLBM achieved on NVIDIA P100, V100, A100 GPUs and AMDMI100 GPU. Using the host to host MPI communications, the performance on MI100 GPU is around 20% better than on V100 GPU. It has turned out that most of the kernels are successfully accelerated except for interpolation kernels for Adaptive Mesh Refinement (AMR) method.
Asahi, Yuichi; Hatayama, Sora*; Shimokawabe, Takashi*; Onodera, Naoyuki; Hasegawa, Yuta; Idomura, Yasuhiro
Proceedings of 2021 IEEE International Conference on Cluster Computing (IEEE Cluster 2021) (Internet), p.686 - 691, 2021/10
Times Cited Count：2 Percentile：72.38(Computer Science, Hardware & Architecture)We develop a convolutional neural network model to predict the multi-resolution steady flow. Based on the state-of-the-art image-to-image translation model pix2pixHD, our model can predict the high resolution flow field from the set of patched signed distance functions. By patching the high resolution data, the memory requirements in our model is suppressed compared to pix2pixHD.
Onodera, Naoyuki; Idomura, Yasuhiro; Hasegawa, Yuta; Nakayama, Hiromasa; Shimokawabe, Takashi*; Aoki, Takayuki*
Boundary-Layer Meteorology, 179(2), p.187 - 208, 2021/05
Times Cited Count：12 Percentile：75.33(Meteorology & Atmospheric Sciences)A plume dispersion simulation code named CityLBM enables a real time simulation for several km by applying adaptive mesh refinement (AMR) method on GPU supercomputers. We assess plume dispersion problems in the complex urban environment of Oklahoma City (JU2003). Realistic mesoscale wind boundary conditions of JU2003 produced by a Weather Research and Forecasting Model (WRF), building structures, and a plant canopy model are introduced to CityLBM. Ensemble calculations are performed to reduce turbulence uncertainties. The statistics of the plume dispersion field, mean and max concentrations show that ensemble calculations improve the accuracy of the estimation, and the ensemble-averaged concentration values in the simulations over 4 km areas with 2-m resolution satisfied factor 2 agreements for 70% of 24 target measurement points and periods in JU2003.
Onodera, Naoyuki; Idomura, Yasuhiro; Hasegawa, Yuta; Shimokawabe, Takashi*; Aoki, Takayuki*
Keisan Kogaku Koenkai Rombunshu (CD-ROM), 26, 3 Pages, 2021/05
We develop a mixed-precision preconditioner for the pressure Poisson equation in a two-phase flow CFD code JUPITER-AMR. The multi-grid (MG) preconditioner is constructed based on the geometric MG method with a three- stage V-cycle, and a cache-reuse SOR (CR-SOR) method at each stage. The numerical experiments are conducted for two-phase flows in a fuel bundle of a nuclear reactor. The MG-CG solver in single-precision shows the same convergence histories as double-precision, which is about 75% of the computational time in double-precision. In the strong scaling test, the MG-CG solver in single-precision is accelerated by 1.88 times between 32 and 96 GPUs.
Asahi, Yuichi; Hatayama, Sora*; Shimokawabe, Takashi*; Onodera, Naoyuki; Hasegawa, Yuta; Idomura, Yasuhiro
Keisan Kogaku Koenkai Rombunshu (CD-ROM), 26, 4 Pages, 2021/05
We develop a convolutional neural network model to predict the multi-resolution steady flow. Based on the state-of-the-art image-to-image translation model Pix2PixHD, our model can predict the high resolution flow field from the signed distance function. By patching the high resolution data, the memory requirements in our model is suppressed compared to Pix2PixHD.
Onodera, Naoyuki; Idomura, Yasuhiro; Hasegawa, Yuta; Yamashita, Susumu; Shimokawabe, Takashi*; Aoki, Takayuki*
Proceedings of International Conference on High Performance Computing in Asia-Pacific Region (HPC Asia 2021) (Internet), p.120 - 128, 2021/01
Times Cited Count：0 Percentile：0.01(Computer Science, Hardware & Architecture)We develop a multigrid preconditioned conjugate gradient (MG-CG) solver for the pressure Poisson equation in a two-phase flow CFD code JUPITER. The MG preconditioner is constructed based on the geometric MG method with a three-stage V-cycle, and a RB-SOR smoother and its variant with cache-reuse optimization (CR-SOR) are applied at each stage. The numerical experiments are conducted for two-phase flows in a fuel bundle of a nuclear reactor. The MG-CG solvers with the RB-SOR and CR-SOR smoothers reduce the number of iterations to less than 15% and 9% of the original preconditioned CG method, leading to 3.1- and 5.9-times speedups, respectively. The obtained performance indicates that the MG-CG solver designed for the block-structured grid is highly efficient and enables large-scale simulations of two-phase flows on GPU based supercomputers.
Onodera, Naoyuki; Idomura, Yasuhiro; Asahi, Yuichi; Hasegawa, Yuta; Shimokawabe, Takashi*; Aoki, Takayuki*
Dai-34-Kai Suchi Ryutai Rikigaku Shimpojiumu Koen Rombunshu (Internet), 2 Pages, 2020/12
We develop a multigrid preconditioned conjugate gradient (MG-CG) solver for the pressure Poisson equation in a two-phase flow CFD code JUPITER. The code is written in C++ and CUDA to keep the portability on multi-platforms. The main kernels of the CG solver achieve reasonable performance as 0.4 
0.75 of the roofline performances, and the performances of the MG-preconditioner are also reasonable on NVIDIA GPU and Intel CPU. However, the performance degradation of the SpMV kernel on ARM is significant. It is confirmed that the optimization does not work if any functions are included in the loop.
Onodera, Naoyuki; Idomura, Yasuhiro; Ali, Y.*; Yamashita, Susumu; Shimokawabe, Takashi*; Aoki, Takayuki*
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.210 - 215, 2020/10
This paper presents a GPU-based Poisson solver on a block-based adaptive mesh refinement (block-AMR) framework. The block-AMR method is essential for GPU computation and efficient description of the nuclear reactor. In this paper, we successfully implement a conjugate gradient method with a state-of-the-art multi-grid preconditioner (MG-CG) on the block-AMR framework. GPU kernel performance was measured on the GPU-based supercomputer TSUBAME3.0. The bandwidth of a vector-vector sum, a matrix-vector product, and a dot product in the CG kernel gave good performance at about 60% of the peak performance. In the MG kernel, the smoothers in a three-stage V-cycle MG method are implemented using a mixed precision RB-SOR method, which also gave good performance. For a large-scale Poisson problem with 
cells, the developed MG-CG method reduced the number of iterations to less than 30% and achieved 
2.5 speedup compared with the original preconditioned CG method.
Onodera, Naoyuki; Idomura, Yasuhiro; Ali, Y.*; Shimokawabe, Takashi*; Aoki, Takayuki*
Keisan Kogaku Koenkai Rombunshu (CD-ROM), 25, 4 Pages, 2020/06
We have developed the stencil-based CFD code JUPITER for simulating three-dimensional multiphase flows. A GPU-accelerated Poisson solver based on the preconditioned conjugate gradient (P-CG) method with a multigrid preconditioner was developed for the JUPITER with block-structured AMR mesh. All Poisson kernels were implemented using CUDA, and the GPU kernel function is well tuned to achieve high performance on GPU supercomputers. The developed multigrid solver shows good convergence of about 1/7 compared with the original P-CG method, and 3 speed up is achieved with strong scaling test from 8 to 216 GPUs on TSUBAME 3.0.
Tomiyasu, Keisuke*; Ito, Naoko*; Okazaki, Ryuji*; Takahashi, Yuki*; Onodera, Mitsugi*; Iwasa, Kazuaki*; Nojima, Tsutomu*; Aoyama, Takuya*; Ogushi, Kenya*; Ishikawa, Yoshihisa*; et al.
Advanced Quantum Technologies (Internet), 1(3), p.1800057_1 - 1800057_7, 2018/12
Spin-state transition, also known as spin crossover, plays a key role in diverse systems. In theory, the boundary range between the low- and high-spin states is expected to enrich the transition and give rise to unusual physical states. However, no compound that realizes a nearly degenerate critical range as the ground state without requiring special external conditions has yet been experimentally identified. This study reports that the Sc substitution in LaCoO3 destabilizes its nonmagnetic low-spin state and generates an anomalous paramagnetic state accompanied by the enhancement of transport gap and magneto-lattice-expansion as well as the contraction of Co-O distance with the increase of electron site transfer. These phenomena are not well described by the mixture of conventional low- and high-spin states, but by their quantum superposition occurring on the verge of a spin-state transition.
Onodera, Naoyuki; Idomura, Yasuhiro; Ali, Y.*; Shimokawabe, Takashi*
Proceedings of 9th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2018) (Internet), p.9 - 16, 2018/11
Times Cited Count：9 Percentile：94.7(Computer Science, Theory & Methods)We develop a communication reduced multi-time- step (CRMT) algorithm for a Lattice Boltzmann method (LBM) based on a block-structured adaptive mesh refinement (AMR). This algorithm is based on the temporal blocking method, and can improve computational efficiency by replacing a communication bottleneck with additional computation. The proposed method is implemented on an extreme scale airflow simulation code CityLBM, and its impact on the scalability is tested on GPU based supercomputers, TSUBAME and Reedbush. Thanks to the CRMT algorithm, the communication cost is reduced by 
, and weak and strong scalings are improved up to 
GPUs. The obtained performance indicates that real time airflow simulations for about 2km square area with the wind speed of 
is feasible using 1m resolution.
Shimokawabe, Takashi*; Endo, Toshio*; Onodera, Naoyuki; Aoki, Takayuki*
Proceedings of 2017 IEEE International Conference on Cluster Computing (IEEE Cluster 2017) (Internet), p.525 - 529, 2017/09
Stencil-based applications such as CFD have succeeded in obtaining high performance on GPU supercomputers. The problem sizes of these applications are limited by the GPU device memory capacity, which is typically smaller than the host memory. On GPU supercomputers, a locality improvement technique using temporal blocking method with memory swapping between host and device enables large computation beyond the device memory capacity. Our high-productivity stencil framework automatically applies temporal blocking to boundary exchange required for stencil computation and supports automatic memory swapping provided by a MPI/CUDA wrapper library. The framework-based application for the airflow in an urban city maintains 80% performance even with the twice larger than the GPU memory capacity and have demonstrated good weak scalability on the TSUBAME 2.5 supercomputer.
Esaka, Fumitaka; Watanabe, Kazuo; Onodera, Takashi; Lee, C. G.; Magara, Masaaki; Sakurai, Satoshi; Usuda, Shigekazu
Applied Surface Science, 255(4), p.1512 - 1515, 2008/12
Times Cited Count：9 Percentile：40.55(Chemistry, Physical)Isotope ratio analysis of individual uranium particles in environmental samples taken at nuclear facilities is one of the key techniques for nuclear safeguards. For accurate analysis, we developed a combined technique of single particle manipulation, which was carried out in scanning electron microscope (SEM), and SIMS. This technique gives information on particle size as well as isotope ratios for each particle. In the present study, the particle size dependence of the precision of uranium isotope ratio is investigated by using the SEM-SIMS technique. In the analysis of individual particles of a CRM U050 (5% enriched uranium) sample, even if the particle diameter was 0.5 
m, the analysis with the precision of 
U/
U isotope ratio within 5% was possible. In the presentation, the results obtained for the other uranium particles having different isotopic composition will also be presented.
Magara, Masaaki; Usuda, Shigekazu; Sakurai, Satoshi; Shinohara, Nobuo; Esaka, Fumitaka; Kokubu, Yoko; Suzuki, Daisuke; Yasuda, Kenichiro; Lee, C. G.; Inagawa, Jun; et al.
Kaku Busshitsu Kanri Gakkai (INMM) Nihon Shibu Dai-29-Kai Nenji Taikai Rombunshu (CD-ROM), 6 Pages, 2008/00
JAEA has been developing, under the auspices of the Ministry of Education, Culture, Sports, Science and Technology of Japan, analytical techniques for ultra-trace amounts of nuclear materials in the environmental samples in order to contribute to the strengthened safeguards system. In January 2003, JAEA was qualified as a member of the IAEA network analytical laboratories (NWAL) for environmental sample analysis. Since then, JAEA has conducted the analysis of domestic and the IAEA samples. Two techniques, bulk and particle analyses, are available for the environmental samples and the latter method generally provides more detailed information about history of nuclear materials in a facility. However, isotope ratios of uranium are measured in the particle analysis at present and it is wished to develop analytical method for isotope ratios of plutonium in individual particles. We commence the development of the plutonium particle analysis and the consideration of age-dating for plutonium particles through the atomic ratio of Pu-241 and Am-241.
Iguchi, Kazunari; Esaka, Konomi; Lee, C. G.; Inagawa, Jun; Esaka, Fumitaka; Onodera, Takashi; Fukuyama, Hiroyasu; Suzuki, Daisuke; Sakurai, Satoshi; Watanabe, Kazuo; et al.
Radiation Measurements, 40(2-6), p.363 - 366, 2005/11
Times Cited Count：11 Percentile：59.93(Nuclear Science & Technology)In particle analysis for safeguards environmental samples, the fission track technique is very important to detect sub-micrometer particles containing uranium. In the technique the authors developed, the particles were recovered onto the polycarbonate membrane filter. The filter was dissolved in solvent and dried to form a thin film of detector, in which the particles were confined. After thermal neutron irradiation and etching, the particles of interest in the detector were easily identified with fission tracks, and were picked up for isotope ratio analysis. It was found, however, that the particles in the vicinity of the detector surface may fall off during the etching process. Therefore, optimization of the etching condition is required. In this work, the effects of etching time and enrichment of uranium in particles were investigated. Preliminary results suggest that etching time should be shorter with the increase in the enrichment.
Sakurai, Satoshi; Magara, Masaaki; Usuda, Shigekazu; Watanabe, Kazuo; Esaka, Fumitaka; Hirayama, Fumio; Lee, C. G.; Yasuda, Kenichiro; Kono, Nobuaki; Inagawa, Jun; et al.
Proceedings of International Conference on Nuclear Energy System for Future Generation and Global Sustainability (GLOBAL 2005) (CD-ROM), 6 Pages, 2005/10
no abstracts in English
Magara, Masaaki; Usuda, Shigekazu; Sakurai, Satoshi; Watanabe, Kazuo; Esaka, Fumitaka; Hirayama, Fumio; Lee, C. G.; Yasuda, Kenichiro; Kono, Nobuaki; Inagawa, Jun; et al.
Dai-26-Kai Kaku Busshitsu Kanri Gakkai (INMM) Nihon Shibu Nenji Taikai Rombunshu, p.157 - 164, 2005/00
JAERI has conducted the analysis of domestic and the IAEA samples. JAERI is developing the analytical techniques to improve the analytical ability for the safeguards environmental samples. For bulk analysis, study is focused on the improvement of reliability of isotope ratio measurements by ICP-MS. New chemical separation techniques are under development and a desolvation module is introduced to reduce the polyatomic interferences. In particle analysis, the sample preparation procedure for SIMS method is modified to measure the 
U/U and 
U/U ratios for individual particles. We are also developing fission track-TIMS method to measure uranium isotope ratios in particles of sub-micrometer size. A screening instrument of X-ray fluorescent analysis is equipped to measure elemental distribution on a swipe surface.
Magara, Masaaki; Usuda, Shigekazu; Sakurai, Satoshi; Watanabe, Kazuo; Esaka, Fumitaka; Hirayama, Fumio; Lee, C. G.; Yasuda, Kenichiro; Kono, Nobuaki; Inagawa, Jun; et al.
Proceedings of INMM 46th Annual Meeting (CD-ROM), 8 Pages, 2005/00
JAERI has been developing analytical techniques for ultra-trace amounts of nuclear materials in the environmental samples in order to contribute to the strengthened safeguards system. Development of essential techniques for bulk and particle analysis of the environmental swipe sample has been established as an ultra-trace analytical method of uranium and plutonium. In January 2003, JAERI was qualified as a member of the IAEA network analytical laboratories for environmental samples. Since then, JAERI has conducted the analysis of domestic and the IAEA samples. From Japanese fiscal year 2003, the second phase of the project was started for the development of advanced techniques, such as analyzing minor actinides and fission products as well as uranium and plutonium, particle analysis using fission-track technique, more efficient particle analysis using ICP-TOFMS and screening by X-ray fluorescent analysis. This paper deals with the progress in the development of the new techniques, applications and future perspective.
