Japan Geoscience Union Meeting 2019

Presentation information

[J] Oral

M (Multidisciplinary and Interdisciplinary) » M-GI General Geosciences, Information Geosciences & Simulations

[M-GI35] Development of computational sciences on planetary formation, evolution and surface environment

Tue. May 28, 2019 1:45 PM - 3:15 PM 301B (3F)

convener:Yoshi-Yuki Hayashi(Department of Planetology/CPS, Graduate School of Science, Kobe University), Masaki Ogawa(Division of General Systems Studies, Graduate School of Arts and Sciences, University of Tokyo), Shigeru Ida(Earth-Life Science Institute, Tokyo Institute of Technology), Kanya Kusano(Institute for Space-Earth Environmental Research, Nagoya University), Chairperson:Hiroki Kashimura

2:30 PM - 2:45 PM

[MGI35-04] Development of N-body Simulation Code for Planetary System Formation With Particle-Particle Particle-Tree Scheme

*Yota Ishigaki1,4, Junko Kominami2, Junichiro Makino3, Masaki Fujimoto4 (1.University of Tokyo, 2.Tokyo Institute of Technology, 3.RIKEN AICS, 4.Institite of Space and Astronautical Science, Japan Aerospace Exploration Agency)

Keywords:N-body simulation, planetary system formation

In general, the planetary system is thought to be formed from protoplanetary disks surrounding the central star. In particular, it is accepted that terrestrial planets and the cores of gas giant planets are formed by accretion of km size objects (planetesimals) in the protoplanetary disk. The accretion process of the planets has been mainly discussed using the gravitational multi-body simulations (N-body simulations) of planetesimal systems. Various unresolved problems have been pointed out in the conventional standard planetary system formation theory.In recent years, new planetary system formation models have been proposed to solve these problems.However, there is no model with no problem even in these new models. Since N-body simulations with sufficient number of particles has not been performed, the generalized planetary system formation theory has not been constructed yet.

We have developed a new N-body simulation code with particle-particle particle-tree (P3T) scheme for planetary system formation, GPLUM. The code GPLUM uses a fourth-order Hermite scheme to calculate gravitational interactions between particles within a cut-off radius and a Barnes-Hut tree scheme for gravitational interactions from particles beyond. The conventional simulation codes with P3T scheme has the bottleneck that the calculation speed decreases when the mass ratio among the particles becomes large. We have solved it by implementing an algorithm which determines the cut-off radius based on mass and orbital semi-major axis of particles individually for each particle. The performance of GPLUM is significantly improved for the simulations of particle systems with mass distribution.

By improving the performance of the N-body simulation code, we made it possible to carry out numerous global simulation with various parameters. GPLUM allows us to perform N-body simulations with ~106-107 particles. By using GPLUM, we will perform N-body simulations with wide range and high resolution and investigate various parameters by perform parameter studies with N-body simulations.