2017 Annual Meeting

Presentation information

Oral presentation

III. Fission Energy Engineering » 307-1 Computational Science and Engineering

[1G01-07] Numerical simulation and Visualization of Thermal Fluid Dynamics and Nuclear Fuel

Mon. Mar 27, 2017 10:00 AM - 11:55 AM Room G (16-207 Building No.16)

Chair: Nakajima Norihiro (JAEA)

11:00 AM - 11:15 AM

[1G05] Development of preconditioned communication avoiding CG solver for multiphase CFD code JUPITER

*Akie Mayumi1, Yasuhiro Idomura1, Takuya Ina1, Susumu Yamada1, Toshiyuki Imamura2 (1. JAEA, 2. AICS, Riken)

Keywords:Multiphase CFD code, Krylov subspace solver, Communication avoiding algorithm

The left-preconditioned communication avoiding conjugate gradient (LP-CA-CG) method is applied to the pressure Poisson equation in the multiphase CFD code JUPITER. Two LP- CA-CG solvers with block Jacobi preconditioning and with underlap preconditioning, in which point Jacobi preconditioning approximation is applied to boundary regions, are developed. The former is developed based on a hybrid CA approach, in which CA is applied only to global collective communications for inner product operations. The latter is a full CA approach, in which CA is applied also to local point-to-point communications in sparse matrix-vector (SpMV) operations and preconditioning. It is shown that on the K computer, the former hybrid approach 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 becomes 2x faster than the original CG solver by reducing the cost of global collective communications by 69 percent.