A comparison of SPH artificial viscosities and their impact on the Keplerian disk

Abstract

Hydrodynamical simulations of rotating disk play important roles in the field of astrophysical and planetary science. Smoothed Particle Hydrodynamics (SPH) has been widely used for such simulations. It, however, has been known that with SPH, a cold and thin Kepler disk breaks up due to the unwanted angular momentum transfer. Two possible reasons have been suggested for this breaking up of the disk; the artificial viscosity (AV) and the numerical error in the evaluation of pressure gradient in SPH. Which one is dominant has been still unclear. In this paper, we investigate the reason for this rapid breaking up of the disk. We implemented most of popular formulations of AV and switches and measured the angular momentum transfer due to both AV and the error of SPH estimate of pressure gradient. We found that the angular momentum transfer due to AV at the inner edge triggers the breaking up of the disk. We also found that the classical von-Neumann-Richtmyer-Landshoff type AV with a high order estimate for $\nabla \cdot \vec{v}$ can maintain the disk for $\sim 100$ orbits even when used with the standard formulation of SPH.

arXiv commentAccepted for publication in The Astrophysical Journal. Animations are available at http://v1.jmlab.jp/~hosono/1601.05903/

References in corpus (23)

• The cosmological simulation code GADGET-2
• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• Cosmological SPH simulations: The entropy equation
• Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
• Gasoline: An adaptable implementation of TreeSPH
• A General Class of Lagrangian Smoothed Particle Hydrodynamics Methods and Implications for Fluid Mixing Problems
• Smoothed Particle Hydrodynamics in Astrophysics
• Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
• Inviscid SPH
• A Density Independent Formulation of Smoothed Particle Hydrodynamics
• Resolving mixing in Smoothed Particle Hydrodynamics
• SPHGal: Smoothed Particle Hydrodynamics with improved accuracy for Galaxy simulations
• SPHS: Smoothed Particle Hydrodynamics with a higher order dissipation switch
• MAGMA: a 3D, Lagrangian magnetohydrodynamics code for merger applications
• Astrophysical Weighted Particle Magnetohydrodynamics
• SEREN - A new SPH code for star and planet formation simulations
• Boosting the accuracy of SPH techniques: Newtonian and special-relativistic tests
• VINE -- A numerical code for simulating astrophysical systems using particles I: Description of the physics and the numerical methods
• Momentum transfer across shear flows in Smoothed Particle Hydrodynamic simulations of galaxy formation
• Improving smoothed particle hydrodynamics with an integral approach to calculating gradients
• VINE -- A numerical code for simulating astrophysical systems using particles II: Implementation and performance characteristics
• EvoL: The new Padova T-SPH parallel code for cosmological simulations - I. Basic code: gravity and hydrodynamics
• SPH Simulations of Accretion Disks and Narrow Rings

Cited by in corpus (13)

• Simulating cosmic structure formation with the GADGET-4 code
• Consequences of Giant Impacts on Early Uranus for Rotation, Internal Structure, Debris, and Atmospheric Erosion
• Bifurcation in the history of Uranus and Neptune: the role of giant impacts
• Fortran interface layer of the framework for developing particle simulator FDPS
• The AGORA High-resolution Galaxy Simulations Comparison Project. III: Cosmological zoom-in simulation of a Milky Way-mass halo
• Dynamic Localised Turbulent Diffusion and its Impact on the Galactic Ecosystem
• Improved Performances in Subsonic Flows of an SPH Scheme with Gradients Estimated using an Integral Approach
• Examining the Accuracy of Astrophysical Disk Simulations With a Generalized Hydrodynamical Test Problem
• A Formulation of Consistent Particle Hydrodynamics in Strong Form
• Novel Hydrodynamic Schemes Capturing Shocks and Contact Discontinuities and Comparison Study with Existing Methods
• Dealing with density discontinuities in planetary SPH simulations
• Merging Criteria for Planetesimal Collisions
• Selfgravitating disks in binary systems: an SPH approach -- I. Implementation of the code and reliability tests