Capacitance matrix methods for the Helmholtz equation on general three-dimensional regions

Journal Article · Sun Jul 01 04:00:00 EDT 1979 · Math. Comput.; (United States)

DOI:https://doi.org/10.2307/2006065· OSTI ID:6091381

O'Leary, D P; Widlund, O

Capacitance matrix methods provide techniques for extending the use of fast Poisson solvers to arbitrary bounded regions. These techniques are further studied and developed with a focus on the three-dimensional case. A discrete analogue of classical potential theory is used as a guide in the design of rapidly convergent iterative methods. Algorithms and programmings aspects of the methods are also explored in detail. Several conjugate gradient methods are discussed for the solution of the capacitance matrix equation. A fast Poisson solver is developed which is numerically very stable even for indefinite Helmholtz equations. Variants thereof allow substantial savings in primary storage for problems on very fine meshes. Numerical results show that accurate solutions can be obtained at a cost which is proportional to that of the fast Helmholtz solver in use.

• Research Organization: National Sience Foundation Grant, University of Michigan
• OSTI ID: 6091381
• Journal Information: Math. Comput.; (United States), Journal Name: Math. Comput.; (United States) Vol. 33:147; ISSN MCMPA
• Country of Publication: United States
• Language: English

Related Subjects

99 GENERAL AND MISCELLANEOUS
990200* -- Mathematics & Computers
ALGORITHMS
BOUNDARY CONDITIONS
DIFFERENTIAL EQUATIONS
EQUATIONS
ITERATIVE METHODS
MATHEMATICAL LOGIC
MATRICES
THREE-DIMENSIONAL CALCULATIONS