Year: 2021
Author: Jiuyang Liang, Pei Liu, Zhenli Xu
Communications in Computational Physics, Vol. 30 (2021), Iss. 1 : pp. 210–226
Abstract
Poisson's equations in a cuboid are frequently solved in many scientific and engineering applications such as electric structure calculations, molecular dynamics simulations and computational astrophysics. In this paper, a fast and highly accurate algorithm is presented for the solution of the Poisson's equation in a cuboidal domain with boundary conditions of mixed type. This so-called harmonic surface mapping algorithm is a meshless algorithm which can achieve a desired order of accuracy by evaluating a body convolution of the source and the free-space Green's function within a sphere containing the cuboid, and another surface integration over the spherical surface. Numerical quadratures are introduced to approximate the integrals, resulting in the solution represented by a summation of point sources in free space, which can be accelerated by means of the fast multipole algorithm. The complexity of the algorithm is linear to the number of quadrature points, and the convergence rate can be arbitrarily high even when the source term is a piecewise continuous function.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/10.4208/cicp.OA-2020-0052
Communications in Computational Physics, Vol. 30 (2021), Iss. 1 : pp. 210–226
Published online: 2021-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 17
Keywords: Fast algorithm Poisson's equation boundary integral method image charge mixed boundary condition fast multipole method.
Author Details
-
HSMA: An O(N) electrostatics package implemented in LAMMPS
Liang, Jiuyang
Yuan, Jiaxing
Xu, Zhenli
Computer Physics Communications, Vol. 276 (2022), Iss. P.108332
https://doi.org/10.1016/j.cpc.2022.108332 [Citations: 4]