A nonconforming saddle point least squares approach for elliptic interface problems
For researchers solving elliptic interface problems, this work extends a previous conforming method to non-conforming spaces, offering an incremental improvement in flux approximation.
The paper presents a non-conforming saddle point least squares method for elliptic interface problems with discontinuous coefficients, achieving quasi-optimal flux approximations using local projections and gradient recovery. Numerical results in 2D and 3D support the method.
We present a non-conforming least squares method for approximating solutions of second order elliptic problems with discontinuous coefficients. The method is based on a general Saddle Point Least Squares (SPLS) method introduced in previous work based on conforming discrete spaces. The SPLS method has the advantage that a discrete $\inf-\sup$ condition is automatically satisfied for standard choices of test and trial spaces. We explore the SPLS method for non-conforming finite element trial spaces which allow higher order approximation of the fluxes. For the proposed iterative solvers, inversion at each step requires bases only for the test spaces. We focus on using projection trial spaces with local projections that are easy to compute. The choice of the local projections for the trial space can be combined with classical gradient recovery techniques to lead to quasi-optimal approximations of the global flux. Numerical results for 2D and 3D domains are included to support the proposed method.