Self-Gravitational Force Calculation of Second Order Accuracy for Infinitesimally Thin Gaseous Disks in Polar Coordinates

Abstract

Investigating the evolution of disk galaxies and the dynamics of proto-stellar disks can involve the use of both a hydrodynamical and a Poisson solver. These systems are usually approximated as infinitesimally thin disks using two- dimensional Cartesian or polar coordinates. In Cartesian coordinates, the calcu- lations of the hydrodynamics and self-gravitational forces are relatively straight- forward for attaining second order accuracy. However, in polar coordinates, a second order calculation of self-gravitational forces is required for matching the second order accuracy of hydrodynamical schemes. We present a direct algorithm for calculating self-gravitational forces with second order accuracy without artifi- cial boundary conditions. The Poisson integral in polar coordinates is expressed in a convolution form and the corresponding numerical complexity is nearly lin- ear using a fast Fourier transform. Examples with analytic solutions are used to verify that the truncated error of this algorithm is of second order. The kernel integral around the singularity is applied to modify the particle method. The use of a softening length is avoided and the accuracy of the particle method is significantly improved.

35 pages, 7 figures, accepted for publication in ApJS

References in corpus (9)

• Type I planetary migration in a self-gravitating disk
• Treating gravity in thin disk simulations
• Hydrodynamical Simulations of the Barred Spiral Galaxy NGC 1097
• Feathering Instability of Spiral Arms. I: Formulation of the Problem
• Embedded Star Formation in S4G Galaxy Dust Lanes
• Feathering Instability of Spiral Arms. II. Parameter Study
• On The Orbital Evolution of Jupiter Mass Protoplanet Embedded in A Self-Gravity Disk
• Effects of Spiral Arms on Star Formation in Nuclear Rings of Barred-spiral Galaxies
• A Fast Potential and Self-Gravity Solver for Non-Axisymmetric Disks

Cited by in corpus (4)

• A Fast Poisson Solver of Second-Order Accuracy for Isolated Systems in Three-Dimensional Cartesian and Cylindrical Coordinates
• Self-Gravitational Force Calculation of Infinitesimally Thin Gaseous Disks on Nested Grids
• Self-gravitational Force Calculation of Second Order Accuracy Using Multigrid Method on Nested Grids
• Self-gravitational Force Calculation of High-Order Accuracy for Infinitesimally Thin Gaseous Disks