Abstract:We develop two unfitted finite element methods for the Stokes equations based on -conforming finite elements. Both cut finite element methods exhibit optimal convergence order for the velocity, pointwise divergence-free velocity fields, and well-posed linear systems, independently of the position of the boundary relative to the computational mesh. The first method is a cut finite element discretization of the Stokes equations based on the Brezzi–Douglas–Marini (BDM) elements and involves interior penalty terms to enforce tangential continuity of the velocity at interior edges in the mesh. The second method is a cut finite element discretization of a three-field formulation of the Stokes problem involving the vorticity, velocity, and pressure and uses the Raviart–Thomas (RT) space for the velocity. We present mixed ghost penalty stabilization terms for both methods so that the resulting discrete problems are stable and the divergence-free property of the -conforming elements is preserved also on unfitted meshes. In both methods boundary conditions are imposed weakly. We show that imposing Dirichlet boundary conditions weakly introduces additional challenges; (1) The divergence-free property of the RT and the BDM finite elements may be lost depending on how the normal component of the velocity field at the boundary is imposed. (2) Pressure robustness is affected by how well the boundary condition is satisfied and may not hold even if the incompressibility condition holds pointwise. We study two approaches of weakly imposing the normal component of the velocity at the boundary; we either use a penalty parameter and Nitsche's method or a Lagrange multiplier method. We show that appropriate conditions on the velocity space has to be imposed when Nitsche's method or penalty is used. Pressure robustness can hold with both approaches by reducing the error at the boundary but the price we pay is seen in the condition numbers of the resulting linear systems, independent of if the mesh is fitted or unfitted to the boundary.

Face‐centred finite volume methods for Stokes flows with variable viscosity

A quasi-Lagrangian finite element method for the Navier-Stokes equations in a time-dependent domain

Fully consistent lowest-order finite element methods for generalised Stokes flows with variable viscosity

A Volume of Fluid Method Algorithm for Simulation of Surface Tension Dominant Two-Phase Flows

A perturbed twofold saddle point-based mixed finite element method for the Navier-Stokes equations with variable viscosity

A face-centred finite volume method for laminar and turbulent incompressible flows

Finite difference methods for 3D viscous incompressible flows in the vorticity-vector potential formulation on nonstaggered grids

Finite Difference Schemes for Incompressible Flows in Vorticity Formulations

Finite Element Iterative Methods for the 3D Steady Navier--Stokes Equations

Galerkin/least-square Finite-Element Methods for Steady Viscoelastic Flows

Finite element setting for fluid flow simulations with natural enforcement of capillary effects

Finite element methods for multicomponent convection-diffusion

Non-Conforming Multiscale Finite Element Method for Stokes Flows in Heterogeneous Media. Part I: Methodologies and Numerical Experiments

Divergence-free cut finite element methods for Stokes flow

A numerical stabilization framework for viscoelastic fluid flow using the finite volume method on general unstructured meshes

Stable and locally mass- and momentum-conservative control-volume finite-element schemes for the Stokes problem

Numerical methods for 3-D unsteady viscous flow and their applications

Discretisation of an Oldroyd-B viscoelastic fluid flow using a Lie derivative formulation

Improved Crouzeix-Raviart scheme for the Stokes and Navier-Stokes problem