Related papers: Toward an efficient second-order method for comput…
We describe three approaches for computing a gravity signal from a density anomaly. The first approach consists of the classical "summation" technique, whilst the remaining two methods solve the Poisson problem for the gravitational…
We present the exact calculus of the gravitational potential and acceleration along the symmetry axis of a plane, homogeneous, polar cell as a function of mean radius a, radial extension e, and opening angle f. Accurate approximations are…
This study aims to establish an analytical model that reproduces the gravitational field around non-spherical bodies with constant density. Due to the non-spherical geometry of such bodies, their gravitational potential is disturbed…
A new one-dimensional, dynamical model is proposed for geometrically thin, self-gravitating viscous accretion discs. The vertically integrated equations are simplified using the slow accretion limit and the monopole approximation with a…
Globular clusters (GCs) are known to harbor multiple stellar populations. To explain these observations Bastian et al. suggested a scenario in which a second population is formed by the accretion of enriched material onto the low-mass stars…
In the core-accretion model the nominal runaway gas-accretion phase brings most planets to multiple Jupiter masses. However, known giant planets are predominantly Jupiter-mass bodies. Obtaining longer timescales for gas accretion may…
With hundreds of exoplanets detected, it is necessary to revisit giant planets accretion models to explain their mass distribution. In particular, formation of sub-jovian planets remains unclear, given the short timescale for the runaway…
Recently, gas giant planets in nearly circular orbits with large semimajor axes ($a \sim$ 30--1000AU) have been detected by direct imaging. We have investigated orbital evolution in a formation scenario for such planets, based on core…
We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas…
We construct a pseudo-Newtonian potential (PNP) corresponding to a rotating black hole solution in a modified gravity (MGR) framework using a metric-based prescription. The motivation is to enable realistic accretion disk studies in MGR,…
The Poisson equation governing a planet's gravitational field is posed on the unbounded domain, $\mathbb{R}^3$, whereas finite-element computations require bounded meshes. We implement and compare three strategies for handling the infinite…
Many important problems in astrophysics, space physics, and geophysics involve flows of (possibly ionized) gases in the vicinity of a spherical object, such as a star or planet. The geometry of such a system naturally favors numerical…
We propose an efficient algorithm for the evaluation of the potential and its gradient of gravitational/electrostatic $N$-body systems, which we call particle mesh multipole method (PMMM or PM$^3$). PMMM can be understood both as an…
The Fast Multipole Method (FMM) for the Poisson equation is extended to the case of non-axisymmetric problems in an axisymmetric domain, described by cylindrical coordinates. The method is based on a Fourier decomposition of the source into…
We describe a new, adaptive solver for the two-dimensional Poisson equation in complicated geometries. Using classical potential theory, we represent the solution as the sum of a volume potential and a double layer potential. Rather than…
In this work we introduced a new proposal to study the gravitational lensing theory by spherical lenses, starting from its surface mass density $\Sigma(x)$ written in terms of a decreasing function $f$ of a dimensionless coordinate $x$ on…
In many astrophysical problems involving discs (gaseous or particulate) orbiting a dominant central mass, gravitational potential of the disc plays an important dynamical role. Its impact on the motion of external objects, as well as on the…
Accretion discs are composed of ionized gas in motion around a central object. Sometimes, the disc is the source of powerful bipolar jets along its rotation axis. Theoretical models invoke the existence of a bipolar magnetic field crossing…
We present a numerical method for solving the Poisson equation on a nested grid. The nested grid consists of uniform grids having different grid spacing and is designed to cover the space closer to the center with a finer grid. Thus our…
We use the adaptive mesh refinement code RAMSES to model the formation of protoplanetary disks in realistic star formation environments. The resolution scales over up to 29 powers of two ($\sim$ 9 orders of magnitude) covering a range from…