Sven Gross

Thomas Ludescher, Sven Gross, Arnold Reusken

We consider discrete Poisson interface problems resulting from linear unfitted finite elements, also called cut finite elements (CutFEM). Three of these unfitted finite element methods known from the literature are studied. All three methods rely on Nitsche s method to incorporate the interface conditions. The main topic of the paper is the development of a multigrid method, based on a novel prolongation operator for the unfitted finite element space and an interface smoother that is designed to yield robustness for large jumps in the diffusion coefficients. Numerical results are presented which illustrate efficiency of this multigrid method and demonstrate its robustness properties with respect to variation of the mesh size, location of the interface and contrast in the diffusion coefficients.

Sven Gross, Maxim A. Olshanskii, Arnold Reusken

In this paper we study a system of advection-diffusion equations in a bulk domain coupled to an advection-diffusion equation on an embedded surface. Such systems of coupled partial differential equations arise in, for example, the modeling of transport and diffusion of surfactants in two-phase flows. The model considered here accounts for adsorption-desorption of the surfactants at a sharp interface between two fluids and their transport and diffusion in both fluid phases and along the interface. The paper gives a well-posedness analysis for the system of bulk-surface equations and introduces a finite element method for its numerical solution. The finite element method is unfitted, i.e., the mesh is not aligned to the interface. The method is based on taking traces of a standard finite element space both on the bulk domains and the embedded surface. The numerical approach allows an implicit definition of the surface as the zero level of a level-set function. Optimal order error estimates are proved for the finite element method both in the bulk-surface energy norm and the $L^2$-norm. The analysis is not restricted to linear finite elements and a piecewise planar reconstruction of the surface, but also covers the discretization with higher order elements and a higher order surface reconstruction.