Partially Penalized Immersed Finite Element Methods for Elliptic Interface Problems

This article presents new immersed finite element (IFE) methods for solving the popular second order elliptic interface problems on structured Cartesian meshes even if the involved interfaces have nontrivial geometries. These IFE methods contain extra stabilization terms introduced only at interface edges for penalizing the discontinuity in IFE functions. With the enhanced stability due to the added penalty, not only these IFE methods can be proven to have the optimal convergence rate in the H1-norm provided that the exact solution has sufficient regularity, but also numerical results indicate that their convergence rates in both the H1-norm and the L2-norm do not deteriorate when the mesh becomes finer which is a shortcoming of the classic IFE methods in some situations. Trace inequalities are established for both linear and bilinear IFE functions that are not only critical for the error analysis of these new IFE methods, but also are of a great potential to be useful in error analysis for other IFE methods

A CutFEM method for two-phase flow problems

In this article, we present a cut finite element method for two-phase Navier-Stokes flows. The main feature of the method is the formulation of a unified continuous interior penalty stabilisation approach for, on the one hand, stabilising advection and the pressure-velocity coupling and, on the other hand, stabilising the cut region. The accuracy of the algorithm is enhanced by the development of extended fictitious domains to guarantee a well defined velocity from previous time steps in the current geometry. Finally, the robustness of the moving-interface algorithm is further improved by the introduction of a curvature smoothing technique that reduces spurious velocities. The algorithm is shown to perform remarkably well for low capillary number flows, and is a first step towards flexible and robust CutFEM algorithms for the simulation of microfluidic devices