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.
- Curved domains
- Discontinuous Galerkin methods
- Elliptic equations
- Immersed boundary methods