Japan Geoscience Union Meeting 2016

Presentation information


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

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

Tue. May 24, 2016 5:15 PM - 6:30 PM Poster Hall (International Exhibition Hall HALL6)

Convener:*Junichiro Makino(RIKEN AICS), Yoshi-Yuki Hayashi(Department of Planetology/CPS, Graduate School of Science, Kobe University), Shigeru Ida(Department of Earth and Planetary Science, Graduate School of Science and Technology, Tokyo Institute of Technology), Yuri Aikawa(Center for Computational Sciences, University of Tsukuba), Masaki Ogawa(Division of General Systems Studies, Graduate School of Arts and Sciences, University of Tokyo), Masayuki Umemura(Center for Computational Sciences, University of Tsukuba)

5:15 PM - 6:30 PM

[MGI22-P03] Comprehensive tests of artificial viscosities, their switches and derivative operators used in Smoothed Particle Hydrodynamics

*Natsuki Hosono1, Takayuki R Saitoh2, Junichiro Makino1,2 (1.RIKEN, Advanced Institute for Computational Science, 2.EARTH-LIFE SCIENCE INSTITUTE)

Keywords:numerical hydrodynamics

In the field of astrophysical and planetary science, hydrodynamical numerical simulations for rotating disk play important role.
So far, Smoothed Particle Hydrodynamics (SPH) has been widely applied for such simulations.
It, however, has been known that with SPH, a cold and thin Kepler disk breaks up due to the unphysical angular momentum transfer.
There are two possible reasons for the breaking up of the disk; the artificial viscosity (AV) and the numerical error in the evaluation of pressure gradient.
However, which one is dominant has been still unclear.
Thus, we performed a systematic survey of how the lifetime of a cold disk varies depending on known implementations of AV and various switchs.
As a result, we found that the angular momentum transfer due to AV at the inner edge triggers the breaking up of the disk in the case of Monaghan (1997)'s AV.
We also found that with the classical von-Neumann-Richtmyer-Landshoff type AV with a high order derivative estimate the disk survives for more than $100$ orbits.