High order VEM on curved domains (Q2319592)

Summary: We deal with the virtual element method (VEM) for solving the Poisson equation on a domain \(\Omega\) with curved boundary. Given a polygonal approximation \(\Omega_h\) of the domain \(\Omega\), the standard order \(m\) VEM [\textit{L. Beirão da Veiga} et al., Math. Models Methods Appl. Sci. 24, No. 8, 1541--1573 (2014; Zbl 1291.65336)], for \(m\) increasing, leads to a suboptimal convergence rate. We adapt the approach of \textit{M. Fernando Benedetto} et al. [``A globally conforming method for solving flow in discrete fracture networks using the virtual element method'', Finite Elem. Anal. Des. 109, 23--36 (2016; \url{doi.org/10.1016/j.finel.2015.10.003})] to VEM and we prove that an optimal convergence rate can be achieved by using a suitable correction depending on high order normal derivatives of the discrete solution at the boundary edges of \(\Omega_h\), which, to retain computability, is evaluated after applying the projector \(\Pi^{\nabla}\) onto the space of polynomials. Numerical experiments confirm the theory.