Numerical Techniques for Threedimensional Smoothed Particle Hydrodynamics Simulations: Applications to Accretion Disks
Abstract
Numerical techniques are described for threedimensional fluid systems in the absence of selfgravity using the Lagrangian method of smoothed particle hydrodynamics (SPH). In particular, we present an efficient method for locating nearest neighbors that uses an ancillary Eulerian grid and conserves memory by partitioning the computational space into manageable layers. Further savings in both memory and computational time are achieved by using interparticle distances that are discretized with respect to small integral increments of the smoothing length. We also present a time integration algorithm using multiple time steps which guarantees that all particles are always synchronous in phase space to a least firstorder accuracy with respect to the individual time steps.
These techniques are used to simulate an accretion disk in a low mass ratio (M_{2}/M_{1} = 0.08) binary system with the ideal gas law, low adiabatic gamma (γ = 1.01), and excluding radiation effects and magnetic fields. The results agree qualitatively with the ShukuraSunyaev αdisk model but overestimate the radial temperature profile by a factor of ∼10, indicating that radiation effects must be included for a complete model.
 Publication:

The Astrophysical Journal
 Pub Date:
 August 1995
 DOI:
 10.1086/176010
 Bibcode:
 1995ApJ...448..822S
 Keywords:

 ACCRETION;
 ACCRETION DISKS;
 HYDRODYNAMICS;
 METHODS: NUMERICAL;
 STARS: BINARIES: CLOSE