Wavelet and finite element solutions for the Neumann problem using fictitious domains

Roland GLOWINSKI, T. W. Pan, Raymond O. Wells, Xiaodong Zhou

Research output: Contribution to journalJournal articlepeer-review

46 Citations (Scopus)


This paper presents a new fictitious domain formulation for the solution of a strongly elliptic boundary value problem with Neumann boundary conditions for a bounded domain in a finite-dimensional Euclidean space with a smooth (possibly only Lipschitz) boundary. This extends the domain to a larger rectangular domain with periodic boundary conditions for which fast solvers are available. The extended solution converges on the original domain in the appropriate function spaces as the penalty parameter approaches zero. Both wavelet-Galerkin and finite elements numerical approximation schemes are developed using this methodology. The convergence rates of both schemes are comparable, and the use of finite elements requires a parameterization of the boundary, while the wavelet-Galerkin method admits an implicit description of the boundary in terms of a wavelet representation of the boundary measure defined as the distributional gradient of the characteristic function of the interior. The accuracy of both methods is investigated and compared, both theoretically and for numerical test cases. The conclusion is that the methods are comparable, and that the wavelet method allows the use of more general boundaries which are not explicitly parametrized, which would be of greater advantage in higher dimensions (the numerical tests are carried out in two dimensions).

Original languageEnglish
Pages (from-to)40-51
Number of pages12
JournalJournal of Computational Physics
Issue number1
Publication statusPublished - Jun 1996

Scopus Subject Areas

  • Numerical Analysis
  • Modelling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'Wavelet and finite element solutions for the Neumann problem using fictitious domains'. Together they form a unique fingerprint.

Cite this