Self-gravitational force calculation of infinitesimally thin gaseous disks
Autor: | Ronald E. Taam, Chien-Chang Yen, Kang C. Jea, Ken Huai-Che Yeh |
---|---|
Rok vydání: | 2012 |
Předmět: |
Numerical Analysis
Physics and Astronomy (miscellaneous) Applied Mathematics Fast Fourier transform Mathematical analysis FOS: Physical sciences Mathematical Physics (math-ph) Computer Science Applications law.invention Convolution Computational Mathematics Classical mechanics Kernel (image processing) law Modeling and Simulation Periodic boundary conditions Cartesian coordinate system Poisson's equation Polar coordinate system Spectral method Mathematical Physics Astrophysics::Galaxy Astrophysics Mathematics |
Zdroj: | Journal of Computational Physics. 231:8246-8261 |
ISSN: | 0021-9991 |
DOI: | 10.1016/j.jcp.2012.08.003 |
Popis: | 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. Comment: 16 pages |
Databáze: | OpenAIRE |
Externí odkaz: |