A numerical method for solving variable coefficient elliptic equation with interfaces

  • Authors:

  • Songming Hou;Xu-Dong Liu

  • Affiliations:

  • Department of Mathematics, University of California Santa Barbara, Santa Barbara, CA 93106, USA;Department of Mathematics, University of California Santa Barbara, Santa Barbara, CA 93106, USA

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2005

Quantified Score

Hi-index 31.54

Visualization

Abstract

A new 2nd order accurate numerical method on non-body-fitting grids is proposed for solving the variable coefficient elliptic equation in disjoint subdomains @W^+/- separated by interfaces @C. 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 C^2 and solutions are C^2 on the closures of the subdomains. The method can handle the problems when the solutions and/or the interfaces are weaker than C^2. For example, u@?H^2(@W^+/-), @C 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.