# zbMATH — the first resource for mathematics

##### Examples
 Geometry Search for the term Geometry in any field. Queries are case-independent. Funct* Wildcard queries are specified by * (e.g. functions, functorial, etc.). Otherwise the search is exact. "Topological group" Phrases (multi-words) should be set in "straight quotation marks". au: Bourbaki & ti: Algebra Search for author and title. The and-operator & is default and can be omitted. Chebyshev | Tschebyscheff The or-operator | allows to search for Chebyshev or Tschebyscheff. "Quasi* map*" py: 1989 The resulting documents have publication year 1989. so: Eur* J* Mat* Soc* cc: 14 Search for publications in a particular source with a Mathematics Subject Classification code (cc) in 14. "Partial diff* eq*" ! elliptic The not-operator ! eliminates all results containing the word elliptic. dt: b & au: Hilbert The document type is set to books; alternatively: j for journal articles, a for book articles. py: 2000-2015 cc: (94A | 11T) Number ranges are accepted. Terms can be grouped within (parentheses). la: chinese Find documents in a given language. ISO 639-1 language codes can also be used.

##### Operators
 a & b logic and a | b logic or !ab logic not abc* right wildcard "ab c" phrase (ab c) parentheses
##### Fields
 any anywhere an internal document identifier au author, editor ai internal author identifier ti title la language so source ab review, abstract py publication year rv reviewer cc MSC code ut uncontrolled term dt document type (j: journal article; b: book; a: book article)
A numerical method for solving variable coefficient elliptic equation with interfaces. (English) Zbl 1061.65123
Summary: A new second-order accurate numerical method on non-body-fitting grids is proposed for solving the variable coefficient elliptic equation in disjoint subdomains $\Omega^{\pm}$ separated by interfaces $\Gamma$. 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 second-order accurate in solution and first-order accurate in its gradient in $L^\infty$ 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\in H^2(\Omega^{\pm})$, $\Gamma$ 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.

##### MSC:
 65N30 Finite elements, Rayleigh-Ritz and Galerkin methods, finite methods (BVP of PDE) 35J25 Second order elliptic equations, boundary value problems
##### Keywords:
finite element methods; numerical experiments
Full Text:
##### References:
 [1] Zhiming, Chen; Jun, Zou: Finite element methods and their convergence for elliptic and parabolic interface problems. Numerische Mathematik 79, 175-202 (1998) · Zbl 0909.65085 [2] Fedkiw, R.; Aslam, T.; Merriman, B.; Osher, S.: A non-oscillatory Eulerian approach to interfaces in multimaterial flows (The ghost fluid method). J. comput. Phys. 152, No. 2, 457-492 (1999) · Zbl 0957.76052 [3] Gibou, F.; Fedkiw, R.; Cheng, L. -T.; Kang, M. A.: Second order accurate symmetric discretization of the Poisson equation on irregular domains. J. comput. Phys. 176, 1-23 (2002) · Zbl 0996.65108 [4] P. Grisvard, Elliptic problems in nonsmooth domains - Monographs and Studies in Mathematics, Pitman Advanced Publishing Program. ISSN 0743-0329, 1985 · Zbl 0695.35060 [5] Leveque, R. J.; Li, Z.: The immersed interface method for elliptic equations with discontinuous coefficients and singular sources. SIAM J. Numer. anal. 31, 1019 (1994) · Zbl 0811.65083 [6] Li, Z.: A fast iterative algorithm for elliptic interface problems. SIAM J. Numer. anal. 35, No. 1, 230-254 (1998) · Zbl 0915.65121 [7] Z. Li, T. Lin, X. Wu, New Cartesian grid methods for interface problems using the finite element formulation, Numerische Mathematik, 9661-98, (2003). Preprint. NCSU CRSC-TR99-12 · Zbl 1055.65130 [8] Liu, Xu-Dong; Fedkiw, Ronald P.; Kang, Myungjoo: A boundary condition capturing method for Poisson’s equation on irregular domains. J. comput. Phys. 160, No. 1, 151-178 (2000) · Zbl 0958.65105 [9] Liu, X. -D.; Sideris, T.: Convergence of the ghost fluid method for elliptic equations with interfaces. Math. comp. 72 (2003) · Zbl 1027.65140 [10] Mayo, A.: The fast solution of Poisson’s and the biharmonic equations in irregular domains. SIAM J. Numer. anal. 21, No. 2, 285-299 (1984) · Zbl 1131.65303 [11] Mayo, A.: Fast high order accurate solutions of Laplace’s equation on irregular domains. SIAM J. Sci. stat. Comput. 6, No. 1, 144-157 (1985) · Zbl 0559.65082 [12] J. Necas, Introduction to the theory of nonlinear elliptic equations, Teubner-Texte zur Mathematik. Band 52, ISSN 0138-502X, 1983 [13] Peskin, C.: Numerical analysis of blood flow in the heart. J. comput. Phys. 25, 220-252 (1977) · Zbl 0403.76100 [14] Peskin; C.; ; Printz, B.: Improved volume conservation in the computation of flows with immersed elastic boundaries. J. comput. Phys. 105, 33-46 (1993) · Zbl 0762.92011 [15] Sussman, M.; Smereka, P.; Osher, S.: A level set approach for computing solutions to incompressible two-phase flow. J. comput. Phys. 114, 146-154 (1994) · Zbl 0808.76077 [16] Wan, Justin W. L.; Liu, Xu-Dong: A boundary condition capturing multigrid approach to irregular boundary problems. SIAM J. Sci. comput. 25, No. 6, 1982-2003 (2004) · Zbl 1107.65032