Related papers: Forward and adjoint calculations of gravitational …
New exact solutions are derived for the gravitational potential inside and outside a homogeneous torus as rapidly converging series of toroidal harmonics. The approach consists of splitting the inter- nal potential into a known solution to…
Astrophysical accretion discs that carry a significant mass compared with their central object are subject to the effect of self-gravity. In the context of circumstellar discs, this can, for instance, cause fragmentation of the disc gas,…
We present an accurate and efficient method to calculate the gravitational potential of an isolated system in three-dimensional Cartesian and cylindrical coordinates subject to vacuum (open) boundary conditions. Our method consists of two…
Adjoint field methods are both elegant and efficient for calculating sensitivity information required across a wide range of physics-based inverse problems. Here we provide a unified approach to the derivation of such methods for problems…
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…
We demonstrate the high accuracy of the density splitting method to compute the gravitational potential and field in the plane of razor-thin, axially symmetric discs, as preliminarily outlined in Pierens & Hure (2004). Because residual…
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…
Accurate gravity field calculations are necessary for landing on planets, moons, asteroids, minimoons, or other irregularly shaped bodies, but current methods become increasingly inaccurate and slow near the surface. We present high…
We present a fast, direct and adaptive Poisson solver for complex two-dimensional geometries based on potential theory and fast multipole acceleration. More precisely, the solver relies on the standard decomposition of the solution as the…
We introduce an accurate and efficient method for a class of nonlocal potential evaluations with free boundary condition, including the 3D/2D Coulomb, 2D Poisson and 3D dipolar potentials. Our method is based on a Gaussian-sum approximation…
We find a new solution to calculate the orbital periastron advance of a test body subject to a central gravitational force field, for relativistic theories and models beyond Einstein. This analitycal formula has general validity that…
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…
For a fundamental solution of Laplace's equation on the $R$-radius $d$-dimensional hypersphere, we compute the azimuthal Fourier coefficients in closed form in two and three dimensions. We also compute the Gegenbauer polynomial expansion…
Over the last two decades, several fast, robust, and high-order accurate methods have been developed for solving the Poisson equation in complicated geometry using potential theory. In this approach, rather than discretizing the partial…
We have developed an adaptive multigrid code for solving the Poisson equation in gravitational simulations. Finer rectangular subgrids are adaptively created in locations where the density exceeds a local level-dependent threshold. We…
We present an efficient and accurate algorithm for solving the Poisson equation in spherical polar coordinates with a logarithmic radial grid and open boundary conditions. The method employs a divide-and-conquer strategy, decomposing the…
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 introduce a technique to compute exact anelastic sensitivity kernels in the time domain using parsimonious disk storage. The method is based on a reordering of the time loop of time-domain forward/adjoint wave propagation solvers…
We solve Poisson's equation by combining a spectral-element method with a mapped infinite-element method. We focus on problems in geostatics and geodynamics, where Earth's gravitational field is determined by Poisson's equation inside the…
Aims. Our aim is to derive a fast and accurate method for computing the gravitational potential of astrophysical objects with high contrasts in density, for which nested or adaptive meshes are required. Methods. We present an extension of…