Solutions of the axi-symmetric Poisson equation from elliptic integrals

Free access article

Issue		A&A Volume 434, Number 1, April IV 2005


Page(s)		1 - 15
Section		Astrophysical processes
DOI		http://dx.doi.org/10.1051/0004-6361:20034194

A&A 434, 1-15 (2005)
DOI: 10.1051/0004-6361:20034194

I. Numerical splitting methods

J.-M. Huré^{1, 2}

¹ LUTh (CNRS UMR 8102), Observatoire de Paris-Meudon, Place Jules Janssen, 92195 Meudon Cedex, France
e-mail: Jean-Marc.Hure@obspm.fr
² Université Paris 7 Denis Diderot, 2 place Jussieu, 75251 Paris Cedex 05, France

(Received 14 August 2003 / Accepted 10 December 2004)

Abstract
In a series of two papers, we present numerical, integral-based methods to compute accurately the self-gravitating field and potential induced by tri-dimensional, axially symmetric fluids, with a special regard for tori, discs and rings. This first article is concerned with a fully numerical approach. Complex shapes, small/large aspect ratios, important density gradients and compact/extended systems can be accounted for. Loop singularities in the Poisson integrals are carefully treated from kernel splitting and/or density splitting. Field components are obtained from density splitting: the local density field is separated into a vertically homogeneous contribution for which the integrable singularity is known in a closed form, plus a "residual" contribution which yields a regular integrand. This technique is exact in the vertically homogeneous limit. The potential is computed from double splitting: kernel splitting (to isolate explicitly the singular function), followed by density splitting. In each two directions, numerical quadratures are performed using a dynamical mesh (very efficient for flat and extended systems), combined with a high-order, irregular-spacing scheme. Global performances are demonstrated through a few test-configurations. The accuracy is potentially very high, and can reach the computer precision (for simple geometries and smooth density profiles). In terms of computing time to precision ratio, present methods are asymptotically more efficient than classical finite-difference methods. These can be employed either as a "Poisson-solver" or only to compute a sub-set of Dirichlet/Neumann-type boundary conditions in order to initialize spectral/finite-difference/finite element methods. As an example of an astrophysical application, we determine the structure of an uniformly rotating polytrope belonging to the compressible, Dyson-Wong sequence. In a second paper (Paper II), we show that singularities are integrable analytically for geometrically thin systems where the density is of the form $\rho \propto z^n$ . Improvements and extensions (such as the non-axi-symmetric case) are possible and discussed.

Key words: hydrodynamics -- methods: numerical

What is OpenURL?

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access. An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.

If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
You can define your own OpenURL resolver with your EDPS Account. In this case your choice will be given priority over that of your library.
You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see http://www.openly.com/openurlref/). You should disable this module if you wish to use the OpenURL server that you or your library have defined.