A New Approach to Handle Curved Meshes in the Hybrid High-Order Method

Abstract

We present here a novel approach to handling curved meshes in polytopal methods within the framework of hybrid high-order methods. The hybrid high-order method is a modern numerical scheme for the approximation of elliptic PDEs. An extension to curved meshes allows for the strong enforcement of boundary conditions on curved domains and for the capture of curved geometries that appear internally in the domain e.g. discontinuities in a diffusion coefficient. The method makes use of non-polynomial functions on the curved faces and does not require any mappings between reference elements/faces. Such an approach does not require the faces to be polynomial and has a strict upper bound on the number of degrees of freedom on a curved face for a given polynomial degree. Moreover, this approach of enriching the space of unknowns on the curved faces with non-polynomial functions should extend naturally to other polytopal methods. We show the method to be stable and consistent on curved meshes and derive optimal error estimates in $$L^2$$ L 2 and energy norms. We present numerical examples of the method on a domain with curved boundary and for a diffusion problem such that the diffusion tensor is discontinuous along a curved arc.

Design and analysis of the Extended Hybrid High-Order method for the Poisson problem

“…Given enrichment spaces satisfying Assumption 2, an enriched method is designed and shown to be well-posed and achieve optimal error estimates in energy norm. Moreover, the method of Yemm [20] of defining polynomials on curved faces is followed, leading to the scheme proposed in this paper being valid for highly generic curved meshes.…”