A High Order Method for the Approximation of Integrals Over Implicitly Defined Hypersurfaces

International audienceWe introduce a novel method to compute approximations of integrals over implicitly defined hyper-surfaces. The new method is based on a weak formulation in L 2 (0, 1), that uses the coarea formula to circumvent an explicit integration over the hypersurfaces. As such it is possible to use standard quadrature rules in the spirit of hp/spectral finite element methods, and the expensive computation of explicit hypersurface parametrizations is avoided. We derive error estimates showing that high order convergence can be achieved provided the integrand and the hypersurface defining function are sufficiently smooth. The theoretical results are supplemented by numerical experiments including an application for plasma modeling in nuclear fusion