A Direct Method for the Solution of Poisson's Equation with Neumann Boundary Conditions on a Staggered Grid of Arbitrary Size
Abstract
A method based on cyclic reduction is described for the solution of the discrete Poisson equation on a rectangular twodimensional staggered grid with an arbitrary number of grid points in each direction. Neumann boundary conditions are assumed in one direction and any boundary condition may be used in the other direction. The coefficients of the equation can be functions of the latter direction so that, e.g., nonequidistant grid spacings or nonCartesian coordinates can be used. Poisson's equation with these boundary conditions describes, e.g., the pressure field of an incompressible fluid flow within rigid boundaries. Numerical results are reported for a FORTRAN subroutine using the method. For an M × N grid the operation count is proportional to MN log _{2}N, and about MN storage locations are required.
 Publication:

Journal of Computational Physics
 Pub Date:
 February 1976
 DOI:
 10.1016/00219991(76)900620
 Bibcode:
 1976JCoPh..20..171S