BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jean-Marc Huré

Jean-Marc Huré contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.IMastro-ph.EPphysics.class-phastro-ph.SRastro-phastro-ph.COmath-phmath.MP

Trust snapshot

Quick read

Trust 21 - Emerging
12works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Nested spheroidal figures of equilibrium I. Approximate solutions for rigid rotations

We discuss the equilibrium conditions for a body made of two homogeneous components separated by oblate spheroidal surfaces and in relative motion. While exact solutions are not permitted for rigid rotation (unless a specific ambient pressure), approximations can be obtained for configurations involving a small confocal parameter. The problem then admits two families of solutions, depending on the pressure along the common interface (constant or quadratic with the cylindrical radius). We give in both cases the pressure and the rotation rates as a function of the fractional radius, ellipticities and mass-density jump. Various degrees of flattening are allowed but there are severe limitations for global rotation, as already known from classical theory (e.g. impossibility of confocal and coelliptical solutions, gradient of ellipticity outward). States of relative rotation are much less constrained, but these require a mass-density jump. This analytical approach compares successfully with the numerical solutions obtained from the Self-Consistent-Field method. Practical formula are derived in the limit of small ellipticities appropriate for slowly-rotating star/planet interiors.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nested spheroidal figures of equilibrium II. Generalization to L layers

We present a vectorial formalism to determine the approximate solutions to the problem of a composite body made of $L$ homogeneous, rigidly rotating layers bounded by spheroidal surfaces. The method is based on the 1st-order expansion of the gravitational potential over confocal parameters, thereby generalizing the method described in Paper I for $L=2$. For a given relative geometry of the ellipses and a given set of mass-density jumps at the interfaces, the sequence of rotation rates and interface pressures is obtained analytically by recursion. A wide range of equilibria result when layers rotate in an asynchronous manner, although configurations with a negative oblateness gradient are more favorable. In contrast, states of global rotation (all layers move at the same rate), found by solving a linear system of $L-1$ equations, are much more constrained. In this case, we mathematically demonstrate that confocal and coelliptical configurations are not permitted. Approximate formula for small ellipticities are derived. These results reinforce and prolongate known results and classical theorems restricted to small elliptiticities. Comparisons with the numerical solutions computed from the Self-Consistent-Field method are successful.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A truly Newtonian softening length for disc simulations

The softened point mass model is commonly used in simulations of gaseous discs including self-gravity while the value of associated length λremains, to some degree, controversial. This ``parameter'' is however fully constrained when, in a discretized disc, all fluid cells are demanded to obey Newton's law. We examine the topology of solutions in this context, focusing on cylindrical cells more or less vertically elongated. We find that not only the nominal length depends critically on the cell's shape (curvature, radial extension, height), but it is either a real or an imaginary number. Setting λas a fraction of the local disc thickness -- as usually done -- is indeed not the optimal choice. We then propose a novel prescription valid irrespective of the disc properties and grid spacings. The benefit, which amounts to 2-3 more digits typically, is illustrated in a few concrete cases. A detailed mathematical analysis is in progress.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-gravity in curved mesh elements

The local character of self-gravity along with the number of spatial dimensions are critical issues when computing the potential and forces inside massive systems like stars and disks. This appears from the discretisation scale where each cell of the numerical grid is a self-interacting body in itself. There is apparently no closed-form expression yet giving the potential of a three-dimensional homogeneous cylindrical or spherical cell, in contrast with the Cartesian case. By using Green's theorem, we show that the potential integral for such polar-type 3D sectors -- initially, a volume integral with singular kernel -- can be converted into a regular line-integral running over the lateral contour, thereby generalising a formula already known under axial symmetry. It therefore is a step towards the obtention of another potential/density pair. The new kernel is a finite function of the cell's shape (with the simplest form in cylindrical geometry), and mixes incomplete elliptic integrals, inverse trigonometric and hyperbolic functions. The contour integral is easy to compute; it is valid in the whole physical space, exterior and interior to the sector itself and works in fact for a wide variety of shapes of astrophysical interest (e.g. sectors of tori or flared discs). This result is suited to easily providing reference solutions, and to reconstructing potential and forces in inhomogeneous systems by superposition. The contour integrals for the 3 components of the acceleration vector are explicitely given.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-gravity in thin discs and edge effects: an extension of Paczynski's approximation

As hydrostatic equilibrium of gaseous discs is partly governed by the gravity field, we have estimated the component caused by a vertically homogeneous disc, with a special attention for the outer regions where self-gravity classically appears. The accuracy of the integral formula is better than 1%, whatever the disc thickness, radial extension and radial density profile. At order zero, the field is even algebraic for thin discs and writes $- 4 πG Σ(R) f_{edge} (R)$ at disc surface, thereby correcting Paczynski's formula by a multiplying factor $f_{edge} \gtrsim 1/2$, which depends on the relative distance to the edges and the local disc thickness. For very centrally condensed discs however, this local contribution can be surpassed by action of mass stored in the inner regions, possibly resulting in $f_{edge} \gg 1$. A criterion setting the limit between these two regimes is derived. These result are robust in the sense that the details of vertical stratification are not critical. We briefly discuss how hydrostatic equilibrium is impacted. In particular, the disc flaring should not reverse in the self-gravitating region, which contradicts what is usually obtained from Paczynski's formula. This suggests that i) these outer regions are probably not fully shadowed by the inner ones (important when illuminated by a central star), and ii) the flared shape of discs does not firmly prove the absence or weakness of self-gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exact, singularity-free recasting of the Newtonian potential in continuous media

The gravitational potential is a key function involved in many astrophysical problems. Its evaluation inside continuous media from Newton's law is known to be challenging because of the diverging kernel 1/|r-r'|. This difficulty is generally treated with avoidance techniques (e.g. multipole expansions, softening length) themselves not without drawbacks. In this article, we present a new path that basically fixes the point-mass singularity problem in systems with, at least, two dimensions. It consists of recasting the gravitational potential in an equivalent integro-differential form, namely the cross-derivative of a "hyperpotential" (i.e., an auxiliary scalar function). In contrast with the potential, the hyperpotential is the convolution of the mass density with a finite amplitude kernel. We show that closed-form expressions for this new kernel can be directly deduced from the potential of homogeneous sheets. We then give a few formulae appropriate to the Cartesian, cylindrical and spherical coordinate systems, including axial symmetry. The method is essentially not limited, either on the geometry of the source or on the distribution, and its implementation is straightforward. Several tests based upon simple quadrature/differentiation schemes are presented (the homogeneous rectangular sheet, cuboid and disk, the Maclaurin disk and a truncated Lane-Emden solution). Compared with a direct summation, the extra computational cost is low and the gain is real: no truncated series, no free parameter, and a relative accuracy better than 1% for typically 16 nodes per spatial direction using the most basic numerical schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A key-formula to compute the gravitational potential of inhomogeneous discs in cylindrical coordinates

We have established the exact expression for the gravitational potential of a homogeneous polar cell - an elementary pattern used in hydrodynamical simulations of gravitating discs. This formula, which is a closed-form, works for any opening angle and radial extension of the cell. It is valid at any point in space, i.e. in the plane of the distribution (inside and outside) as well as off-plane, thereby generalizing the results reported by Durand (1953) for the circular disc. The three components of the gravitational acceleration are given. The mathematical demonstration proceeds from the "incomplete version of Durand's formula" for the potential (based on complete elliptic integrals). We determine first the potential due to the circular sector (i.e. a pie-slice sheet), and then deduce that of the polar cell (from convenient radial scaling and subtraction). As a by-product, we generate an integral theorem stating that "the angular average of the potential of any circular sector along its tangent circle is 2/PI times the value at the corner". A few examples are presented. For numerical resolutions and cell shapes commonly used in disc simulations, we quantify the importance of curvature effects by performing a direct comparison between the potential of the polar cell and that of the Cartesian (i.e. rectangular) cell having the same mass. Edge values are found to deviate roughly like 2E-3 x N/256 in relative (N is the number of grid points in the radial direction), while the agreement is typically four orders of magnitude better for values at the cell's center. We also produce a reliable approximation for the potential, valid in the cell's plane, inside and close to the cell. Its remarkable accuracy, about 5E-4 x N/256 in relative, is sufficient to estimate the cell's self-acceleration.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A substitute for the singular Green kernel in the Newtonian potential of celestial bodies

The "point mass singularity" inherent in Newton's law for gravitation represents a major difficulty in accurately determining the potential and forces inside continuous bodies. Here we report a simple and efficient analytical method to bypass the singular Green kernel 1/|r-r'| inside the source without altering the nature of the interaction. We build an equivalent kernel made up of a "cool kernel", which is fully regular (and contains the long-range -GM/r asymptotic behavior), and the gradient of a "hyperkernel", which is also regular. Compared to the initial kernel, these two components are easily integrated over the source volume using standard numerical techniques. The demonstration is presented for three-dimensional distributions in cylindrical coordinates, which are well-suited to describing rotating bodies (stars, discs, asteroids, etc.) as commonly found in the Universe. An example of implementation is given. The case of axial symmetry is treated in detail, and the accuracy is checked by considering an exact potential/surface density pair corresponding to a flat circular disc. This framework provides new tools to keep or even improve the physical realism of models and simulations of self-gravitating systems, and represents, for some of them, a conclusive alternative to softened gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

AGN disks and black holes on the weighting scales

We exploit our formula for the gravitational potential of finite size, power-law disks to derive a general expression linking the mass of the black hole in active galactic nuclei (AGN), the mass of the surrounding disk, its surface density profile (through the power index s), and the differential rotation law. We find that the global rotation curve v(R) of the disk in centrifugal balance does not obey a power law of the cylindrical radius R (except in the confusing case s = -2 that mimics a Keplerian motion), and discuss the local velocity index. This formula can help to understand how, from position-velocity diagrams, mass is shared between the disk and the black hole. To this purpose, we have checked the idea by generating a sample of synthetic data with different levels of Gaussian noise, added in radius. It turns out that, when observations are spread over a large radial domain and exhibit low dispersion (standard deviation less than 10% typically), the disk properties (mass and s-parameter) and black hole mass can be deduced from a non linear fit of kinematic data plotted on a (R, Rv 2)-diagram. For a deviation higher than 10%, masses are estimated fairly well from a linear regression (corresponding to the zeroth-order treatment of the formula), but the power index s is no longer accessible. We have applied the model to 7 AGN disks whose rotation has already been probed through water maser emission. For NGC3393 and UGC3789, the masses seem well constrained through the linear approach. For IC1481, the power-law exponent s can even be deduced. Because the model is scale-free, it applies to any kind of star/disk system. Extension to disks around young stars showing deviation from Keplerian motion is thus straightforward.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Newtonian potential of thin disks

The one-dimensional, ordinary differential equation (ODE) by Huré & Hersant (2007) that satisfies the midplane gravitational potential of truncated, flat power-law disks is extended to the whole physical space. It is shown that thickness effects (i.e. non-flatness) can be easily accounted for by implementing an appropriate "softening length" $λ$. The solution of this "softened ODE" has the following properties: i) it is regular at the edges (finite radial accelerations), ii) it possesses the correct long-range properties, iii) it matches the Newtonian potential of a geometrically thin disk very well, and iv) it tends continuously to the flat disk solution in the limit $λ\rightarrow 0$. As illustrated by many examples, the ODE, subject to exact Dirichlet conditions, can be solved numerically with efficiency for any given colatitude at second-order from center to infinity using radial mapping. This approach is therefore particularly well-suited to generating grids of gravitational forces in order to study particles moving under the field of a gravitating disk as found in various contexts (active nuclei, stellar systems, young stellar objects). Extension to non-power-law surface density profiles is straightforward through superposition. Grids can be produced upon request.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A local prescription for the softening length in self-gravitating gaseous discs

In 2D-simulations of self-gravitating gaseous discs, the potential is often computed in the framework of &#34;softened gravity&#34; initially designed for N-body codes. In this special context, the role of the softening length LAMBDA is twofold: i) to avoid numerical singularities in the integral representation of the potential (i.e., arising when the relative separation vanishes), and ii) to acount for stratification of matter in the direction perpendicular to the disc mid-plane. So far, most studies have considered LAMBDA 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 LAMBDA such that the gravitational potential of a flat disc coincides at order zero with that of a geometically thin disc of the same surface density. Our prescription for LAMBDA, 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 LAMBDA : i) is independant of the numerical mesh, ii) is always a fraction of the local thickness H, iii) goes through a minimum at the singularity (i.e., at null separation), and iv) is such that 0.13 < LAMBDA/H < 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).

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

The global structure of thin, stratified &#34;alpha&#34;-discs and the reliability of the one layer approximation

We report the results of a systematic comparison between the vertically averaged model and the vertically explicit model of steady state, Keplerian, optically thick &#34;alpha&#34;-discs. The simulations have concerned discs currently found in three different systems: dwarf novae, young stellar objects and active galactic nuclei. In each case, we have explored four decades of accretion rates and almost the whole disc area (except the narrow region where the vertically averaged model has degenerate solutions). We find that the one layer approach gives a remarkably good estimate of the main physical quantities in the disc, and specially the temperature at the equatorial plane which is accurate to within 30% for cases considered. The major deviations (by a factor < 4) are observed on the disc half-thickness. The sensitivity of the results to the &#34;alpha&#34;-parameter value has been tested for 0.001 < alpha < 0.1 and appears to be weak. This study suggests that the ``precision&#39;&#39; of the vertically averaged model which is easy to implement should be sufficient in practice for many astrophysical applications.

0 citations0 reviewsNo code linkedOpen access