2009/09/03 by J.-M. Huré, Jean-Marc Huré, A. Pierens +1 · 1 citation
Mathematics · Physics and Astronomy · #Astrophysical Phenomena and Observations #Astrophysics and Star Formation Studies #Classical mechanics #Context (archaeology) #Gravitation #Gravitational potential #Gravitational singularity #Lambda #Mathematical analysis #Mathematical physics #Mathematics #Physics #Quantum mechanics #Singularity #Stellar, planetary, and galactic studies #Stratification (seeds) #astro-ph.IM
paper · pdf · doi:10.1051/0004-6361/200912348
accepted in A&A, 7 pages, 7 figures, web link for the F90 code and on-line calculations : http://www.obs.u-bordeaux1.fr/radio/JMHure/intro2single.php
openalex publication_date 2009/09/03 · arxiv created 2009/09/04 · arxiv updated 2015/05/14 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/05
In 2D-simulations of self-gravitating gaseous discs, the potential is often computed in the framework of “softened gravity” initially designed for <i>N<i/>-body codes. In this special context, the role of the softening length <i>λ<i/> is twofold: i) to avoid numerical singularities in the integral representation of the potential (i.e., arising when the separation ); and ii) to account for stratification of matter in the direction perpendicular to the disc mid-plane. So far, most studies have considered <i>λ<i/> as a free parameter and various values or formulae have been proposed without much mathematical justification. In this paper, we demonstrate by means of a rigorous calculus that it is possible to define <i>λ<i/> such that the gravitational potential of a flat disc coincides at order zero with that of a geometrically thin disc of the same surface density. Our prescription for <i>λ<i/>, valid in the local, axisymmetric limit, has the required properties i) and ii). It is mainly an analytical function of the radius and disc thickness, and is sensitive to the vertical stratification. For mass density profiles considered (namely, profiles expandable over even powers of the altitude), we find that <i>λ<i/>: i) is independant of the numerical mesh, ii) is always a fraction of the local thickness <i>H<i/>; iii) goes through a minimum at the singularity (i.e., at null separation); and iv) is such that 0.13 <i>λ<i/>/<i>H<i/> 0.29 typically (depending on the separation and on density profile). These results should help us to improve the quality of 2D- and 3D-simulations of gaseous discs in several respects (physical realism, accuracy, and computing time).