Higher-order spectral element method for the stationary Stokes interface problem in two dimensions

This article presents a higher-order spectral element method for the two-dimensional Stokes interface problem involving a piecewise constant viscosity coefficient. The proposed numerical formulation is based on least-squares formulation. The mesh is aligned with the interface, and the interface is completely resolved using blending element functions. The higher-order spectral element functions are nonconforming, and the same-order approximation is used for both velocity and pressure variables. The interface conditions are added to the minimizing functional in appropriate Sobolev norms. Stability and error estimates are proven. The proposed method is shown to be exponentially accurate in both velocity and pressure variables. The theoretical estimates are validated through various numerical examples.