A novel fast solver for Poisson equation with the Neumann boundary condition

In this paper we present a novel fast method to solve Poisson equation in an arbitrary two dimensional region with Neumann boundary condition. The basic idea is to solve the original Poisson problem by a two-step procedure: the first one finds the electric displacement field $\mathbf{D}$ and the second one involves the solution of potential $\phi$. The first step exploits loop-tree decomposition technique that has been applied widely in integral equations within the computational electromagnetics (CEM) community. We expand the electric displacement field in terms of a tree basis. Then, coefficients of the tree basis can be found by the fast tree solver in O(N) operations. Such obtained solution, however, fails to expand the exact field because the tree basis is not completely curl free. Despite of this, the accurate field could be retrieved by carrying out a procedure of divergence free field removal. Subsequently, potential distribution $\phi$ can be found rapidly at the second stage with another fast approach of O(N) complexity. As a result, the method dramatically reduces solution time comparing with traditional FEM with iterative method. Numerical examples including electrostatic simulations are presented to demonstrate the efficiency of the proposed method.

[1]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[2]  D. A. Dunnett.  Classical Electrodynamics , 2020, Nature.

[3]  William L. Briggs,et al.  A multigrid tutorial , 1987 .

[4]  J. Jackson.  Classical Electrodynamics, 3rd Edition , 1998 .

[5]  W. Chew,et al.  Integral equation solution of Maxwell's equations from zero frequency to microwave frequencies , 2000 .

[6]  J. Bladel,et al.  Electromagnetic Fields , 1985 .

[7]  Peter Arbenz,et al.  A fast parallel Poisson solver on irregular domains applied to beam dynamics simulations , 2009, J. Comput. Phys..

[8]  Leslie Greengard,et al.  A free-space adaptive fmm-based pde solver in three dimensions , 2011 .

[9]  An FFT Based Fast Poisson Solver on Spherical Shells , 2011 .

[10]  A. Brigo,et al.  The Poisson–Boltzmann equation for biomolecular electrostatics: a tool for structural biology , 2002, Journal of molecular recognition : JMR.

[11]  Ming-Chih Lai,et al.  Fast direct solvers for Poisson equation on 2D polar and spherical geometries , 2002 .

[12]  Leslie Greengard,et al.  A New Fast-Multipole Accelerated Poisson Solver in Two Dimensions , 2001, SIAM J. Sci. Comput..

[13]  Roger F. Harrington,et al.  An E-Field solution for a conducting surface small or comparable to the wavelength , 1984 .

[14]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[15]  Weng Cho Chew,et al.  Integral Equation Methods for Electromagnetic and Elastic Waves , 2007, Synthesis Lectures on Computational Electromagnetics.

[16]  Wayne H. Schubert,et al.  Multigrid methods for elliptic problems: a review , 1986 .

[17]  L. Greengard,et al.  A Fast Poisson Solver for Complex Geometries , 1995 .

[18]  G. Vecchi,et al.  A multiresolution method of moments for triangular meshes , 2005, IEEE Transactions on Antennas and Propagation.

[19]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[20]  L. Greengard,et al.  A Direct Adaptive Poisson Solver of Arbitrary Order Accuracy , 1996 .

[21]  Eftychios Sifakis,et al.  A parallel multigrid Poisson solver for fluids simulation on large grids , 2010, SCA '10.