Self-gravitational Force Calculation of High-Order Accuracy for Infinitesimally Thin Gaseous Disks
arXiv:1904.07400 · doi:10.3847/1538-4365/ab1a33
Abstract
Self-gravitational force calculation for infinitesimally thin disks is important for studies on the evolution of galactic and protoplanetary disks. Although high-order methods have been developed for hydrodynamic and magneto-hydrodynamic equations, high-order improvement is desirable for solving self-gravitational forces for thin disks. In this work, we present a new numerical algorithm that is of linear complexity and of high-order accuracy. This approach is fast since the force calculation is associated with a convolution form, and the fast calculation can be achieved using Fast Fourier Transform. The nice properties, such as the finite supports and smoothness, of B-splines are exploited to stably interpolate a surface density and achieve a high-order accuracy in forces. Moreover, if the mass distribution of interest is exclusively confined within a calculation domain, the method does not require artificial boundary values to be specified before the force calculation. To validate the proposed algorithm, a series of numerical tests, ranging from 1st- to 3rd-order implementations, are performed and the results are compared with analytic expressions derived for 3rd- and 4th-order generalized Maclaurin disks. We conclude that the improvement on the numerical accuracy is significant with the order of the method, with only little increase of the complexity of the method.
arXiv comment18 pages, 9 tables, 2 figures, accepted for publication in ApJS
References in corpus (17)
- A Dynamical Model for the Formation of Gas Rings and Episodic Starbursts Near Galactic Centres
- Emergence of Protoplanetary Disks and Successive Formation of Gaseous Planets by Gravitational Instability
- A Dynamical Model for Gas Flows, Star Formation, and Nuclear Winds in Galactic Centres
- Gaseous Structures in Barred Galaxies: Effects of the Bar Strength
- Effects of Gas on Formation and Evolution of Stellar Bars and Nuclear Rings in Disk Galaxies
- Physical properties of the circumnuclear starburst ring in the barred Galaxy NGC 1097
- 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
- Potential-density pairs for a family of finite disks
- Effects of Spiral Arms on Star Formation in Nuclear Rings of Barred-spiral Galaxies
- Self-gravitational force calculation of infinitesimally thin gaseous disks
- Self-Gravitational Force Calculation of Infinitesimally Thin Gaseous Disks on Nested Grids
- Self-Gravitational Force Calculation of Second Order Accuracy for Infinitesimally Thin Gaseous Disks in Polar Coordinates
- Gap formation in a self-gravitating disk and the associated migration of the embedded giant planet