“…Typically, XFEM and related methods add, locally, special functions to the FE approximation, in order to better approximate particular features of a problem. For problems with discontinuities in the primal dependent variable (cracks [29,28,13], fluid-structure interaction [26]), discontinuous functions are added. To handle discontinuities in the dual unknown (multi-materials [41], solidification [16], biofilms [35], multi-phase flow [15], Stokes flow [47]), functions with discontinuous derivatives enhance the standard basis.…”