9:15 AM - 9:30 AM

# [MGI28-02] Development of a particle-based parallel code for mantle convection with the variable inertia method

Keywords:Mantle, Numerical simulation, Parallel computing

In order to handle the extreme conditions of the mantle, most of the previous studies ignore the inertia term based on the fact that $Pr\sim10^{24}$ and solve the approximated governing equation implicitly with (1) global matrix inversion or (2) iterative method. The first approach requires the frequent global communication which is not suitable for the current massively-parallel computers. The second approach is faster compared to the first one. However, in some cases, it requires many iterations because of the strongly temperature dependent viscosity, resulting in the slowdown of computation. We thus adopt a different approach; we adopt the variable inertia method (VIM; Takeyama et al. 2017). VIM is an explicit method and thus it is fit for the current cutting-edge supercomputers. In a naive implementation, explicit method is difficult to solve long-term evolution of mantle because of its too tight CFL condition. To overcome this, VIM increases the inertia term so that we can obtain the numerical solution within reasonable computational time. In addition to this, VIM introduces two factors to alter the viscosity and thermal diffusion terms in order to keep the original Rayleigh number of the mantle. According to Takeyama et al. 2017, as long as the modified non-dimensional numbers are being reasonable ranges (Pr > 10, M < 1, and Re < 1), we obtain the almost identical results with the conventional formulation.

We implement VIM to a parallel code, VIM++. The parallel implementation of VIM++ is assisted by a framework for developing particle simulator (FDPS; Iwasawa et al. 2016) developed by the particle simulator team in AICS/RIKEN. With this, our code can scale from one to million CPU cores without changing codes. In the current version of VIM++, we adopted density-independent smoothed particle hydrodynamics (DISPH; Saitoh & Makino 2013, Hosono et al. 2013) as the hydro scheme, but we will implement a novel hydro scheme, consistent particle hydrodynamics in strong form (CPHSF; Yamamoto & Makino 2017), which has space high-order. CPHSF can solve problems of SPH induced by its low accuracy and will improve the numerical modeling of mantle convection with particle-based methods.

In our presentation, we explain the overview of VIM++, details of schemes, and preliminary results of two- and three-dimensional convection simulations obtained by VIM++.