Related papers: An FFT-based Solution Method for the Poisson Equat…
We propose a collocation method based on multivariate polynomial splines over triangulation or tetrahedralization for the numerical solution of partial differential equations. We start with a detailed explanation of the method for the…
One of the challenges when simulating astrophysical flows with self-gravity is to compute the gravitational forces. In contrast to the hyperbolic hydrodynamic equations, the gravity field is described by an elliptic Poisson equation. We…
This work presents a multigrid preconditioned high order immersed finite difference solver to accurately and efficiently solve the Poisson equation on complex 2D and 3D domains. The solver employs a low order Shortley-Weller multigrid…
Two long-standing problems with the post-Newtonian approximation for isolated slowly-moving systems in general relativity are: (i) the appearance at high post-Newtonian orders of divergent Poisson integrals, casting a doubt on the soundness…
We compute the structure of a self-gravitating torus with polytropic equation-of-state (EOS) rotating in an imposed centrifugal potential. The Poisson-solver is based on isotropic multigrid with optimal covering factor (fluid…
A new method is presented for solving Poisson's equation inside an open-ended rectangular pipe. The method uses Fast Fourier Transforms (FFTs) to perform mixed convolutions and correlations of the charge density with the Green function.…
When solving the Poisson equation by the finite element method, we use one degree of freedom for interpolation by the given Laplacian - the right hand side function in the partial differential equation. The finite element solution is the…
Investigating the evolution of disk galaxies and the dynamics of proto-stellar disks can involve the use of both a hydrodynamical and a Poisson solver. These systems are usually approximated as infinitesimally thin disks using two-…
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…
High order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution $\Delta x$) with only a moderate increase to computational expense. Significant…
An expansion of a density field or particle distribution in basis functions which solve the Poisson equation both provides an easily parallelized n-body force algorithm and simplifies perturbation theories. The expansion converges quickly…
We present a generalized algorithm based on a spherical harmonics expansion method for efficient computation of the three-dimensional gravitational potential on a multi-patch grid in spherical geometry. Instead of solving for the…
We present an explicit solver of the three-dimensional screened and unscreened Poisson's equation which combines accuracy, computational efficiency and versatility. The solver, based on a mixed plane-wave / interpolating scaling function…
The self-similar asymptotics for solutions to the drift-diffusion equation with fractional dissipation, coupled to the Poisson equation, is analyzed in the whole space. It is shown that in the subcritical and supercritical cases, the…
Many problems in beam physics and plasma physics require the solution of Poisson's equation with free-space boundary conditions. The algorithm proposed by Hockney and Eastwood is a popular scheme to solve this problem numerically, used by…
We consider the numerical solution of Poisson's equation on structured grids using geometric multigrid with nonstandard coarse grids and coarse level operators. We are motivated by the problem of developing high-order accurate numerical…
It is well known that the usual mixed method for solving the biharmonic eigenvalue problem by decomposing the operator into two Laplacians may generate spurious eigenvalues on non-convex domains. To overcome this difficulty, we adopt a…
A new approximation method for inverting the Poisson's equation is presented for a continuously distributed and finite-sized source in an unbound domain. The advantage of this image multipole method arises from its ability to place the…
In multi-phase fluid flow, fluid-structure interaction, and other applications, partial differential equations (PDEs) often arise with discontinuous coefficients and singular sources (e.g., Dirac delta functions). These complexities arise…
Singular charge sources in terms of Dirac delta functions present a well-known numerical challenge for solving Poisson's equation. For a sharp interface between inhomogeneous media, singular charges could be analytically treated by…