A numerical method for solving variable coefficient elliptic equation with interfaces

A new 2nd order accurate numerical method on non-body-fitting grids is proposed for solving the variable coefficient elliptic equation in disjoint subdomains Ω± separated by interfaces Γ. The variable coefficients, the source term, and hence the solution itself and its derivatives may be discontinuous across the interfaces. Jump conditions in solution and its co-normal derivative at interface are prescribed. Instead of smooth, the interfaces are only required to be Lipschitz continuous as submanifold. A weak formulation is developed, the existence, uniqueness and regularity of the solutions are studied. The numerical method is derived by discretizing the weak formulation. The method is different from traditional finite element methods. Extensive numerical experiments are presented and show that the method is 2nd order accurate in solution and 1st order accurate in its gradient in L∞ norm if the interface is C2 and solutions are C2 on the closures of the subdomains. The method can handle the problems when the solutions and/or the interfaces are weaker than C2. For example, u ∈ H2(Ω±), Γ is Lipschitz continuous and their singularities coincide, see Example 18 in Section 4. The accuracies of the method under various circumstances are listed in Table 19.