Trace Finite Element Methods for PDEs on Surfaces

In this paper we consider a class of unfitted finite element methods for discretization of partial differential equations on surfaces. In this class of methods known as the Trace Finite Element Method (TraceFEM), restrictions or traces of background surface-independent finite element functions are used to approximate the solution of a PDE on a surface. We treat equations on steady and time-dependent (evolving) surfaces. Higher order TraceFEM is explained in detail. We review the error analysis and algebraic properties of the method. The paper navigates through the known variants of the TraceFEM and the literature on the subject

A trace finite element method for a class of coupled bulk-interface transport problems

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

A cell-centred finite volume approximation for second order partial derivative operators with full matrix on unstructured meshes in any space dimension

Finite volume methods for problems involving second order operators with full diffusion matrix can be used thanks to the definition of a discrete gradient for piecewise constant functions on unstructured meshes satisfying an orthogonality condition. This discrete gradient is shown to satisfy a strong convergence property on the interpolation of regular functions, and a weak one on functions bounded for a discrete $H^1$ norm. To highlight the importance of both properties, the convergence of the finite volume scheme on a homogeneous Dirichlet problem with full diffusion matrix is proven, and an error estimate is provided. Numerical tests show the actual accuracy of the method