A Fast 3D Poisson Solver of Arbitrary Order Accuracy

E. Braverman, M. Israeli, A. Averbuch, L. Vozovoi

Research output: Contribution to journalArticlepeer-review

42 Scopus citations

Abstract

We present a direct solver for the Poisson and Laplace equations in a 3D rectangular box. The method is based on the application of the discrete Fourier transform accompanied by a subtraction technique which allows reducing the errors associated with the Gibbs phenomenon and achieving any prescribed rate of convergence. The algorithm requiresO(N3logN) operations, whereNis the number of grid points in each direction. We show that our approach allows accurate treatment of singular cases which arise when the boundary function is discontinuous or incompatible with the differential equation.

Original languageEnglish
Pages (from-to)109-136
Number of pages28
JournalJournal of Computational Physics
Volume144
Issue number1
DOIs
StatePublished - 20 Jul 1998
Externally publishedYes

Keywords

  • 3D Poisson solver for Dirichlet problem
  • Corner and edge singularities
  • Fourier method

Fingerprint

Dive into the research topics of 'A Fast 3D Poisson Solver of Arbitrary Order Accuracy'. Together they form a unique fingerprint.

Cite this