Accurate Solution and Gradient Computation for Elliptic Interface Problems with Variable Coefficients

A new augmented method is proposed for elliptic interface problems with a piecewise variable coefficient that has a finite jump across a smooth interface. The main motivation is not only to get a second order accurate solution but also a second order accurate gradient from each side of the interface. The key of the new method is to introduce the jump in the normal derivative of the solution as an augmented variable and re-write the interface problem as a new PDE that consists of a leading Laplacian operator plus lower order derivative terms near the interface. In this way, the leading second order derivatives jump relations are independent of the jump in the coefficient that appears only in the lower order terms after the scaling. An upwind type discretization is used for the finite difference discretization at the irregular grid points near or on the interface so that the resulting coefficient matrix is an M-matrix. A multi-grid solver is used to solve the linear system of equations and the GMRES iterative method is used to solve the augmented variable. Second order convergence for the solution and the gradient from each side of the interface has also been proved in this paper. Numerical examples for general elliptic interface problems have confirmed the theoretical analysis and efficiency of the new method.

[1]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[2]  Semyon Tsynkov,et al.  The Method of Difference Potentials for the Helmholtz Equation Using Compact High Order Schemes , 2012, Journal of Scientific Computing.

[3]  J. R. Cannon,et al.  A two phase Stefan problem with flux boundary conditions , 1971 .

[4]  Zhilin Li.  A Fast Iterative Algorithm for Elliptic Interface Problems , 1998 .

[5]  Weiwei Sun,et al.  Error estimation of a class of quadratic immersed finite element methods for elliptic interface problems , 2007 .

[6]  Liwei Shi,et al.  A Numerical Method for Solving 3D Elasticity Equations with Sharp-Edged Interfaces , 2013 .

[7]  Zhilin Li.  The immersed interface method using a finite element formulation , 1998 .

[8]  Ralf Massjung,et al.  An Unfitted Discontinuous Galerkin Method Applied to Elliptic Interface Problems , 2012, SIAM J. Numer. Anal..

[9]  Danny C. Sorensen,et al.  The Approximation and Computation of a Basis of the Trace Space H1/2 , 2007, J. Sci. Comput..

[10]  Sarah L. Mitchell,et al.  Finite-difference methods with increased accuracy and correct initialization for one-dimensional Stefan problems , 2009, Appl. Math. Comput..

[11]  Yekaterina Epshteyn,et al.  On the Solution of the Elliptic Interface Problems by Difference Potentials Method , 2015 .

[12]  Tao Lin,et al.  New Cartesian grid methods for interface problems using the finite element formulation , 2003, Numerische Mathematik.

[13]  P. Bassanini,et al.  Elliptic Partial Differential Equations of Second Order , 1997 .

[14]  Tao Lin,et al.  Partially Penalized Immersed Finite Element Methods For Elliptic Interface Problems , 2015, SIAM J. Numer. Anal..

[15]  Xu-dong Liu,et al.  A numerical method for solving variable coefficient elliptic equation with interfaces , 2005 .

[16]  Zhilin Li,et al.  Fast and accurate numerical approaches for Stefan problems and crystal growth , 1999 .

[17]  Yekaterina Epshteyn,et al.  High-order difference potentials methods for 1D elliptic type models , 2015 .

[18]  Ronald Fedkiw,et al.  The immersed interface method. Numerical solutions of PDEs involving interfaces and irregular domains , 2007, Math. Comput..

[19]  Maximum Norm Regularity of Periodic Elliptic Difference Operators With Variable Coefficients , 2015 .

[20]  Petr Kloucek,et al.  Generalized Harmonic Functions and the Dewetting of Thin Films , 2007 .

[21]  P. Hansbo,et al.  An unfitted finite element method, based on Nitsche's method, for elliptic interface problems , 2002 .

[22]  J. Zou,et al.  Finite element methods and their convergence for elliptic and parabolic interface problems , 1998 .

[23]  Wenjun Ying,et al.  A kernel-free boundary integral method for elliptic boundary value problems , 2007, J. Comput. Phys..

[24]  Xiaoming He,et al.  The convergence of the bilinear and linear immersed finite element solutions to interface problems , 2012 .

[25]  Randall J. LeVeque,et al.  Finite difference methods for ordinary and partial differential equations - steady-state and time-dependent problems , 2007 .

[26]  Tao Lin,et al.  Linear and bilinear immersed finite elements for planar elasticity interface problems , 2012, J. Comput. Appl. Math..

[27]  Kazufumi Ito,et al.  Maximum Principle Preserving Schemes for Interface Problems with Discontinuous Coefficients , 2001, SIAM J. Sci. Comput..

[28]  Eftychios Sifakis,et al.  A second order virtual node method for elliptic problems with interfaces and irregular domains in three dimensions , 2012, J. Comput. Phys..

[29]  Vincent Martin,et al.  Modeling Fractures and Barriers as Interfaces for Flow in Porous Media , 2005, SIAM J. Sci. Comput..

[30]  Ivan Yotov,et al.  Coupling Fluid Flow with Porous Media Flow , 2002, SIAM J. Numer. Anal..

[31]  Yekaterina Epshteyn,et al.  High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces , 2017 .

[32]  Haijun Wu,et al.  An unfitted $hp$-interface penalty finite element method for elliptic interface problems , 2010, 1007.2893.

[33]  Ivo Babuska,et al.  The finite element method for elliptic equations with discontinuous coefficients , 1970, Computing.

[34]  Do Y. Kwak,et al.  Optimal convergence analysis of an immersed interface finite element method , 2010, Adv. Comput. Math..

[35]  L. Zhilin,et al.  ANALYSIS AND NUMERICAL METHODS FOR SOME CRACK PROBLEMS , 2011 .

[36]  James H. Bramble,et al.  A finite element method for interface problems in domains with smooth boundaries and interfaces , 1996, Adv. Comput. Math..

[37]  Xiaoming He,et al.  Immersed finite element methods for elliptic interface problems with non-homogeneous jump conditions , 2011 .

[38]  T. Hou,et al.  Removing the Cell Resonance Error in the Multiscale Finite Element Method via a Petrov-Galerkin Formulation , 2004 .

[39]  Shan Zhao,et al.  High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources , 2006, J. Comput. Phys..

[40]  R. M. Furzeland,et al.  A Comparative Study of Numerical Methods for Moving Boundary Problems , 1980 .

[41]  Kye T. Wee,et al.  An Analysis of a Broken P1-Nonconforming Finite Element Method for Interface Problems , 2009, SIAM J. Numer. Anal..

[42]  Zhilin Li,et al.  On convergence of the immersed boundary method for elliptic interface problems , 2014, Math. Comput..

[43]  L. Wahlbin.  Superconvergence in Galerkin Finite Element Methods , 1995 .

[44]  Junping Wang,et al.  A weak Galerkin finite element method for second-order elliptic problems , 2011, J. Comput. Appl. Math..

[45]  Wolfgang Hackbusch,et al.  Elliptic Differential Equations: Theory and Numerical Treatment , 2017 .

[46]  Zhimin Zhang,et al.  A New Finite Element Gradient Recovery Method: Superconvergence Property , 2005, SIAM J. Sci. Comput..

[47]  T. Hou,et al.  Removing the stiffness from interfacial flows with surface tension , 1994 .

[48]  Michael Pruitt,et al.  Large time step maximum norm regularity of L-stable difference methods for parabolic equations , 2014, Numerische Mathematik.

[49]  Weiwei Sun,et al.  Quadratic immersed finite element spaces and their approximation capabilities , 2006, Adv. Comput. Math..

[50]  Crack jump conditions for elliptic problems , 1999 .

[51]  Xingzhou Yang,et al.  The immersed interface method for elasticity problems with interfaces , 2002 .

[52]  Hongkai Zhao,et al.  A semi-implicit augmented IIM for Navier-Stokes equations with open and traction boundary conditions , 2013 .

[53]  Analytic solutions of a multi-interface transmission problem and crack approximation , 2000 .

[54]  Houde Han.  The numerical solutions of interface problems by infinite element method , 1982 .

[55]  Fernando A. Morales,et al.  The narrow fracture approximation by channeled flow , 2010 .

[56]  Jie Liu,et al.  Open and traction boundary conditions for the incompressible Navier-Stokes equations , 2009, J. Comput. Phys..

[57]  P. Kloucek,et al.  Spectral Solutions of Self-adjoint Elliptic Problems with Immersed Interfaces , 2011 .

[58]  N. K. Kulman,et al.  Method of difference potentials and its applications , 2001 .

[59]  Yekaterina Epshteyn,et al.  High-order accurate difference potentials methods for parabolic problems , 2015 .

[60]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[61]  R. Fedkiw,et al.  A Boundary Condition Capturing Method for Poisson's Equation on Irregular Domains , 2000 .

[62]  J. Thomas Beale,et al.  ON THE ACCURACY OF FINITE DIFFERENCE METHODS FOR ELLIPTIC PROBLEMS WITH INTERFACES , 2006 .

[63]  Stanley Osher,et al.  A Hybrid Method for Moving Interface Problems with Application to the Hele-Shaw Flow , 1997 .

[64]  Giles Auchmuty,et al.  Spectral Characterization of the Trace Spaces Hs(∂Ω) , 2006, SIAM J. Math. Anal..

[65]  Giles Auchmuty,et al.  Steklov Eigenproblems and the Representation of Solutions of Elliptic Boundary Value Problems , 2005 .

[66]  Xiaoming He,et al.  Approximation capability of a bilinear immersed finite element space , 2008 .

[67]  Bo Li,et al.  Immersed-Interface Finite-Element Methods for Elliptic Interface Problems with Nonhomogeneous Jump Conditions , 2007, SIAM J. Numer. Anal..