Structure-preserving parametric finite element methods for two-phase Stokes flow based on Lagrange multiplier approaches

We present a novel formulation for parametric finite element methods to approximate two-phase Stokes flow. The new formulation is based on the classical Stokes equation in the bulk and a novel choice of interface conditions with additional Lagrange multipliers. This new Lagrange multiplier approach ensures that the numerical methods exactly preserve two physical structures of two-phase Stokes flow at the fully discrete level: (i) the energy-decaying and (ii) the volume-preserving properties. Moreover, different types of higher-order time discretization methods are employed, including the Crank--Nicolson method and the second-order backward differentiation formula approach. The resulting schemes are nonlinear and can be efficiently solved by using the Newton method with a decoupling technique. Extensive numerical experiments demonstrate that our methods achieve the desired temporal accuracy while preserving the two physical structures of the two-phase Stokes system.