A penalty finite element method for a fluid system posed on embedded surface

TL;DR

This paper develops a finite element method for simulating incompressible fluid flow on embedded surfaces, providing stability, error analysis, and numerical validation for the approach.

Contribution

It introduces a penalty finite element scheme for Navier-Stokes equations on surfaces, including stability, error estimates, and numerical demonstrations.

Findings

Method achieves convergence in numerical tests

Error depends on penalty parameter as analyzed

Conserves physical quantities effectively

Abstract

The paper introduces a finite element method for the incompressible Navier--Stokes equations posed on a closed surface $\Gamma\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 $\Gamma$ 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.