Controlling Artificial Viscosity in SPH simulations of accretion disks

Abstract

We test the operation of two methods for selective application of Artificial Viscosity (AV) in SPH simulations of Keplerian Accretion Disks, using a ring spreading test to quantify effective viscosity, and a correlation coefficient technique to measure the formation of unwanted prograde alignments of particles. Neither the Balsara Switch nor Time Dependent Viscosity work effectively, as they leave AV active in areas of smooth shearing flow, and do not eliminate the accumulation of alignments of particles in the prograde direction. The effect of both switches is periodic, the periodicity dependent on radius and unaffected by the density of particles. We demonstrate that a very simple algorithm activates AV only when truly convergent flow is detected and reduces the unwanted formation of prograde alignments. The new switch works by testing whether all the neighbours of a particle are in Keplerian orbit around the same point, rather than calculating the divergence of the velocity field, which is very strongly affected by Poisson noise in the positions of the SPH particles.

8 pages, 5 figures

References in corpus (3)

• The properties of brown dwarfs and low-mass hydrogen-burning stars formed by disc fragmentation
• Pseudo-viscous modeling of self-gravitating discs and the formation of low mass ratio binaries
• Radiative transfer and the energy equation in SPH simulations of star formation

Cited by in corpus (7)

• Asteroseismic determination of obliquities of the exoplanet systems Kepler-50 and Kepler-65
• SEREN - A new SPH code for star and planet formation simulations
• Episodic accretion, protostellar radiative feedback, and their role in low-mass star formation
• Convergence of simulations of self-gravitating accretion discs II: Sensitivity to the implementation of radiative cooling and artificial viscosity
• SPH Simulations of Black Hole Accretion: A Step to Model Black Hole Feedback in Galaxies
• Tendency of spherically imploding plasma liners formed by merging plasma jets to evolve toward spherical symmetry
• Examining the Accuracy of Astrophysical Disk Simulations With a Generalized Hydrodynamical Test Problem