Related papers: A Fast Potential and Self-Gravity Solver for Non-A…
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…
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 new fast solver to calculate fixed-boundary plasma equilibria in toroidally axisymmetric geometries. By combining conformal mapping with Fourier and integral equation methods on the unit disk, we show that high-order accuracy…
A fast $N$-body code has been developed for simulating a stellar disk embedded in a live dark matter halo. In generating its Poisson solver, a self-consistent field (SCF) code which inherently possesses perfect scalability is incorporated…
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 calculate ab initio the gravitational potential energy per unit area for a gravitationally coupled multi-component galactic disk of stars and gas, which is given as the integration over vertical density distribution, vertical…
We introduce a fast algorithm for computing volume potentials - that is, the convolution of a translation invariant, free-space Green's function with a compactly supported source distribution defined on a uniform grid. The algorithm relies…
We present a novel numerical method for solving the elliptic partial differential equation problem for the electrostatic potential with piecewise constant conductivity. We employ an integral equation approach for which we derive a system of…
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…
In this article we propose a new efficient strategy to construct exact solutions of Einstein gravities with a minimally coupled self-interacting scalar field. The strategy is to use the symmetry of the equations of motion (EOMs) to give a…
Large classes of materials systems in physics and engineering are governed by magnetic and electrostatic interactions. Continuum or mesoscale descriptions of such systems can be cast in terms of integral equations, whose direct…
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…
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 describe a parallel version of our tree-code for the simulation of self-gravitating systems in Astrophysics. It is based on a dynamic and adaptive method for the domain decomposition, which exploits the hierarchical data arrangement used…
This paper describes the implementation of the direct solution method (DSM) using radial spectral elements for the calculation of synthetic seismograms in self-gravitating, spherically symmetric, non-rotating, anelastic, and transversely…
We describe a short, reproducible workflow for applying finite differences on nonuniform grids determined by a positive weight function g. The grid is obtained by equidistribution, mapping uniform computational coordinates $\xi\in[0,1]$ to…
The approximate computation of all gravitational forces between $N$ interacting particles via the fast multipole method (FMM) can be made as accurate as direct summation, but requires less than $\mathcal{O}(N)$ operations. FMM groups…
Near-future surveys promise a dramatic improvement in the number and precision of astrometric, photometric and spectroscopic measurements of stars in the Milky Way's disk. We examine the impact of such surveys on our understanding of the…
We present a simple and effective multigrid-based Poisson solver of second-order accuracy in both gravitational potential and forces in terms of the one, two and infinity norms. The method is especially suitable for numerical simulations…
In 2D simulations of thin gaseous disks with embedded planets or self-gravity the gravitational potential needs to be smoothed to avoid singularities in the numerical evaluation of the gravitational potential or force. In order to correctly…