Self-gravitational force calculation of infinitesimally thin gaseous disks

Abstract

A thin gaseous disk has often been investigated in the context of various phenomena in galaxies, which point to the existence of starburst rings and dense circumnuclear molecular disks. The effect of self-gravity of the gas in the 2D disk can be important in confronting observations and numerical simulations in detail. For use in such applications, a new method for the calculation of the gravitational force of a 2D disk is presented. Instead of solving the complete potential function problem, we calculate the force in infinite planes in Cartesian and polar coordinates by a reproducing kernel method. Under the limitation of a 2D disk, we specifically represent the force as a double summation of a convolution of the surface density and a fundamental kernel and employ a fast Fourier transform technique. In this method, the entire computational complexity can be reduced from $O(N^2\times N^2)$ to $O((N\times {\log }_2(N{)}^2)$, where $N$ is the number of zones in one dimension. This approach does not require softening. The proposed method is similar to a spectral method, but without the necessity of imposing a periodic boundary condition. We further show this approach is of near second order accuracy for a smooth surface density in a Cartesian coordinate system.

16 pages

References in corpus (1)

• On The Orbital Evolution of Jupiter Mass Protoplanet Embedded in A Self-Gravity Disk

Cited by in corpus (4)

• Spiral Density Waves in M81. II. Hydrodynamic Simulations for Gas Response to Stellar Spiral Density Waves
• 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
• Direct Calculation of Self-Gravitational Force for Infinitesimally Thin Gaseous Disks Using Adaptive Mesh Refinement