The numerical solution of the biharmonic equation by conformal mapping

Raymond H. CHAN*, Thomas K. DELILLO, Mark A. HORN

*Corresponding author for this work

Research output: Journal PublicationsJournal Article (refereed)peer-review

16 Citations (Scopus)

Abstract

The solution to the biharmonic equation in a simply connected region Ω in the plane is computed in terms of the Goursat functions. The boundary conditions are conformally transplanted to the disk with a numerical conformal map. A linear system is obtained for the Taylor coefficients of the Goursat functions. The coefficient matrix of the linear system can be put in the form I + K, where K is the discretization of a compact operator. K can be thought of as the composition of a block Hankel matrix with a diagonal matrix. The compactness leads to clustering of eigenvalues, and the Hankel structure yields a matrix-vector multiplication cost of O(N log N). Thus, if the conjugate gradient method is applied to the system, then superlinear convergence will be obtained. Numerical results are given to illustrate the spectrum clustering and superlinear convergence.

Original languageEnglish
Pages (from-to)1571-1582
Number of pages12
JournalSIAM Journal on Scientific Computing
Volume18
Issue number6
DOIs
Publication statusPublished - Nov 1997
Externally publishedYes

Keywords

  • Biharmonic equation
  • Hankel matrices
  • Numerical conformal mapping

Fingerprint

Dive into the research topics of 'The numerical solution of the biharmonic equation by conformal mapping'. Together they form a unique fingerprint.

Cite this