Vladimir Yushutin

Maxim A. Olshanskii, Annalisa Quaini, Arnold Reusken et al.

We consider a Stokes problem posed on a 2D surface embedded in a 3D domain. The equations describe an equilibrium, area-preserving tangential flow of a viscous surface fluid and serve as a model problem in the dynamics of material interfaces. In this paper, we develop and analyze a Trace finite element method (TraceFEM) for such a surface Stokes problem. TraceFEM relies on finite element spaces defined on a fixed, surface-independent background mesh which consists of shape-regular tetrahedra. Thus, there is no need for surface parametrization or surface fitting with the mesh. The TraceFEM treated here is based on $P_1$ bulk finite elements for both the velocity and the pressure. In order to enforce the velocity vector field to be tangential to the surface we introduce a penalty term. The method is straightforward to implement and has an $O(h^2)$ geometric consistency error, which is of the same order as the approximation error due to the $P_1$--$P_1$ pair for velocity and pressure. We prove stability and optimal order discretization error bounds in the surface $H^1$ and $L^2$ norms. A series of numerical experiments is presented to illustrate certain features of the proposed TraceFEM.

Maxim A. Olshanskii, Vladimir Yushutin

The paper introduces a finite element method for the incompressible Navier--Stokes equations posed on a closed surface $Γ\subset\R^3$. The method needs a shape regular tetrahedra mesh in $\mathbb{R}^3$ to discretize equations on the surface, which can cut through this mesh in a fairly arbitrary way. Stability and error analysis of the fully discrete (in space and in time) scheme is given. The tangentiality condition for the velocity field on $Γ$ is enforced weakly by a penalty term. The paper studies both theoretically and numerically the dependence of the error on the penalty parameter. Several numerical examples demonstrate convergence and conservation properties of the finite element method.

Vladimir Yushutin, Annalisa Quaini, Sheereen Majd et al.

Conservative and non-conservative phase-field models are considered for the numerical simulation of lateral phase separation and coarsening in biological membranes. An unfitted finite element method is devised for these models to allow for a flexible treatment of complex shapes in the absence of an explicit surface parametrization. For a set of biologically relevant shapes and parameter values, the paper compares the dynamic coarsening produced by conservative and non-conservative numerical models, its dependence on certain geometric characteristics and convergence to the final equilibrium