Solving dirichlet boundary-value problems on curved domains by extensions from subdomains

Bernardo Cockburn, Manuel Solano

Research output: Contribution to journalArticlepeer-review

21 Scopus citations

Abstract

We present a technique for numerically solving Dirichlet boundary-value problems for second-order elliptic equations on domains Ω with curved boundaries. This is achieved by using suitably defined extensions from polyhedral subdomains Dh; the problem of dealing with curved boundaries is thus reduced to the evaluations of simple line integrals. The technique is independent of the representation of the boundary, of the space dimension, and allows for the use of only polyhedral elements. We apply this technique to the hybridizable discontinuous Galerkin method and provide extensive numerical experiments showing that the convergence properties of the resulting method are the same as those for the case in which Ω = Dh whenever the distance of Dh to ∂Ω is of order h. In particular, we find that when using polynomial approximations of degree k on the polyhedral subdomains Dh, both the scalar and vector unknowns converge with the optimal order of k + 1 in the whole domain Ω for k ≥ 0. Moreover, we also find that for k ≥ 1, both a postprocessing of the scalar unknown in the polyhedral subdomain Dh as well as its extension from it converge with order k + 2. The extension also converges with order 2 for k = 0.

Original languageEnglish (US)
JournalSIAM Journal on Scientific Computing
Volume34
Issue number1
DOIs
StatePublished - May 28 2012

Keywords

  • Curved domains
  • Discontinuous Galerkin methods
  • Elliptic equations
  • Immersed boundary methods

Fingerprint Dive into the research topics of 'Solving dirichlet boundary-value problems on curved domains by extensions from subdomains'. Together they form a unique fingerprint.

Cite this