A Galerkin formulation of the MIB method for three dimensional elliptic interface problems

We develop a three dimensional (3D) Galerkin formulation of the matched interface and boundary (MIB) method for solving elliptic partial differential equations (PDEs) with discontinuous coefficients, i.e., the elliptic interface problem. The present approach builds up two sets of elements respectively on two extended subdomains which both include the interface. As a result, two sets of elements overlap each other near the interface. Fictitious solutions are defined on the overlapping part of the elements, so that the differentiation operations of the original PDEs can be discretized as if there was no interface. The extra coefficients of polynomial basis functions, which furnish the overlapping elements and solve the fictitious solutions, are determined by interface jump conditions. Consequently, the interface jump conditions are rigorously enforced on the interface. The present method utilizes Cartesian meshes to avoid the mesh generation in conventional finite element methods (FEMs). We implement the proposed MIB Galerkin method with three different elements, namely, rectangular prism element, five-tetrahedron element and six-tetrahedron element, which tile the Cartesian mesh without introducing any new node. The accuracy, stability and robustness of the proposed 3D MIB Galerkin are extensively validated over three types of elliptic interface problems. In the first type, interfaces are analytically defined by level set functions. These interfaces are relatively simple but admit geometric singularities. In the second type, interfaces are defined by protein surfaces, which are truly arbitrarily complex. The last type of interfaces originates from multiprotein complexes, such as molecular motors. Near second order accuracy has been confirmed for all of these problems. To our knowledge, it is the first time for an FEM to show a near second order convergence in solving the Poisson equation with realistic protein surfaces. Additionally, the present work offers the first known near second order accurate method for C1 continuous or H2 continuous solutions associated with a Lipschitz continuous interface in a 3D setting.

[1]  Guido Morgenthal,et al.  An immersed interface method for the vortex-in-cell algorithm , 2007 .

[2]  Boyce E. Griffith,et al.  On the order of accuracy of the immersed boundary method: Higher order convergence rates for sufficiently smooth problems , 2005 .

[3]  Shi Jin,et al.  Robust numerical simulation of porosity evolution in chemical vapor infiltration III: three space dimension , 2003 .

[4]  Lucy T. Zhang,et al.  On computational issues of immersed finite element methods , 2009, J. Comput. Phys..

[5]  Ralf Hiptmair,et al.  Convergence analysis of finite element methods for H(curl; Ω)-elliptic interface problems , 2012, Numerische Mathematik.

[6]  Qiong Zheng,et al.  Second-order Poisson-Nernst-Planck solver for ion transport , 2011, J. Comput. Phys..

[7]  R. Fedkiw,et al.  A fourth order accurate discretization for the Laplace and heat equations on arbitrary domains, with applications to the Stefan problem , 2005 .

[8]  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..

[9]  Weihua Geng,et al.  Treatment of charge singularities in implicit solvent models. , 2007, The Journal of chemical physics.

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

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

[12]  Guo-Wei Wei,et al.  Matched interface and boundary (MIB) method for elliptic problems with sharp-edged interfaces , 2007, J. Comput. Phys..

[13]  Shan Zhao.  Full-Vectorial Matched Interface and Boundary (MIB) Method for the Modal Analysis of Dielectric Waveguides , 2008, Journal of Lightwave Technology.

[14]  Wei Shyy,et al.  COMPUTATIONS OF DROP DYNAMICS WITH THE IMMERSED BOUNDARY METHOD, PART 2: DROP IMPACT AND HEAT TRANSFER , 2003 .

[15]  Guo-Wei Wei,et al.  MIB Galerkin method for elliptic interface problems , 2014, J. Comput. Appl. Math..

[16]  Peng Song,et al.  A weak formulation for solving elliptic interface problems without body fitted grid , 2013, J. Comput. Phys..

[17]  Duan Chen,et al.  Modeling and simulation of electronic structure, material interface and random doping in nano-electronic devices , 2010, J. Comput. Phys..

[18]  C. Peskin.  Numerical analysis of blood flow in the heart , 1977 .

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

[20]  Randall J. LeVeque,et al.  An Immersed Interface Method for Incompressible Navier-Stokes Equations , 2003, SIAM J. Sci. Comput..

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

[22]  Jan Vierendeels,et al.  Flow simulations in rotary volumetric pumps and compressors with the fictitious domain method , 2002 .

[23]  Peter Hansbo,et al.  Interior-penalty-stabilized Lagrange multiplier methods for the finite-element solution of elliptic interface problems , 2010 .

[24]  Wei Shyy,et al.  COMPUTATIONS OF DROP DYNAMICS WITH THE IMMERSED BOUNDARY METHOD, PART 1: NUMERICAL ALGORITHM AND BUOYANCY-INDUCED EFFECT , 2003 .

[25]  Sining Yu,et al.  Three-dimensional matched interface and boundary (MIB) method for treating geometric singularities , 2007, J. Comput. Phys..

[26]  Tianbing Chen,et al.  Piecewise-polynomial discretization and Krylov-accelerated multigrid for elliptic interface problems , 2008, J. Comput. Phys..

[27]  H. Fasel,et al.  A high-order immersed interface method for simulating unsteady incompressible flows on irregular domains , 2005 .

[28]  Ruo Li,et al.  An Edge-based Anisotropic Mesh Refinement Algorithm and Its Application to Interface Problems , 2010 .

[29]  W. Kath,et al.  Modal analysis of circular Bragg fibers with arbitrary index profiles. , 2006, Optics letters.

[30]  Jiangguo Liu,et al.  A matched interface and boundary method for solving multi-flow Navier-Stokes equations with applications to geodynamics , 2012, J. Comput. Phys..

[31]  Guo-Wei Wei,et al.  On the fictitious-domain and interpolation formulations of the matched interface and boundary (MIB) method , 2006, J. Comput. Phys..

[32]  Guo-Wei Wei,et al.  Highly accurate biomolecular electrostatics in continuum dielectric environments , 2008, J. Comput. Chem..

[33]  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 .

[34]  Shun Zhang,et al.  Discontinuous Galerkin Finite Element Methods for Interface Problems: A Priori and A Posteriori Error Estimations , 2011, SIAM J. Numer. Anal..

[35]  Jan S. Hesthaven,et al.  High-order accurate methods in time-domain computational electromagnetics: A review , 2003 .

[36]  Zhilin Li,et al.  The immersed finite volume element methods for the elliptic interface problems , 1999 .

[37]  L. Greengard,et al.  A Fast Poisson Solver for Complex Geometries , 1995 .

[38]  Petter Andreas Berthelsen,et al.  A decomposed immersed interface method for variable coefficient elliptic equations with non-smooth and discontinuous solutions , 2004 .

[39]  M. Burger,et al.  Level set methods for geometric inverse problems in linear elasticity , 2004 .

[40]  R. Kafafy,et al.  Three‐dimensional immersed finite element methods for electric field simulation in composite materials , 2005 .

[41]  Joël Piraux,et al.  How to Incorporate the Spring-Mass Conditions in Finite-Difference Schemes , 2002, SIAM J. Sci. Comput..

[42]  B. Engquist,et al.  Numerical approximations of singular source terms in differential equations , 2004 .

[43]  Shan Zhao,et al.  High-order FDTD methods via derivative matching for Maxwell's equations with material interfaces , 2004 .

[44]  Wei Cai,et al.  An upwinding embedded boundary method for Maxwell's equations in media with material interfaces: 2D case , 2003 .

[45]  Shan Zhao,et al.  WEAK GALERKIN METHODS FOR SECOND ORDER ELLIPTIC INTERFACE PROBLEMS. , 2013, Journal of computational physics.

[46]  Yang Xiang,et al.  Matched interface and boundary (MIB) method for the vibration analysis of plates , 2009 .

[47]  Guo-Wei Wei,et al.  Discrete singular convolution for the solution of the Fokker–Planck equation , 1999 .

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

[49]  D. Zorin,et al.  A fast solver for the Stokes equations with distributed forces in complex geometries , 2004 .

[50]  P. Colella,et al.  A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains , 1998 .

[51]  Anita T. Layton.  Using integral equations and the immersed interface method to solve immersed boundary problems with stiff forces , 2009 .

[52]  Shan Zhao,et al.  Author ' s personal copy High order matched interface and boundary methods for the Helmholtz equation in media with arbitrarily curved interfaces , 2010 .

[53]  Andreas Wiegmann,et al.  The Explicit-Jump Immersed Interface Method: Finite Difference Methods for PDEs with Piecewise Smooth Solutions , 2000, SIAM J. Numer. Anal..

[54]  R. Verzicco,et al.  Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations , 2000 .

[55]  Weihua Geng,et al.  Treatment of geometric singularities in implicit solvent models. , 2007, The Journal of chemical physics.

[56]  James P. Keener,et al.  An Immersed Interface Method for Solving Anisotropic Elliptic Boundary Value Problems in Three Dimensions , 2003, SIAM J. Sci. Comput..

[57]  Zhilin Li,et al.  The Immersed Interface/Multigrid Methods for Interface Problems , 2002, SIAM J. Sci. Comput..

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

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

[60]  Isaac Harari,et al.  Analysis of an efficient finite element method for embedded interface problems , 2010 .

[61]  Duan Chen,et al.  MIBPB: A software package for electrostatic analysis , 2011, J. Comput. Chem..

[62]  Shi Jin,et al.  Robust Numerical Simulation of Porosity Evolution in Chemical Vapor Infiltration , 2002 .

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

[64]  Guo-Wei Wei,et al.  Synchronization-based image edge detection , 2002 .

[65]  David Farrell,et al.  Immersed finite element method and its applications to biological systems. , 2006, Computer methods in applied mechanics and engineering.

[66]  Guo-Wei Wei,et al.  Adaptively deformed mesh based interface method for elliptic equations with discontinuous coefficients , 2012, J. Comput. Phys..

[67]  G. Iaccarino,et al.  Immersed boundary technique for turbulent flow simulations , 2003 .

[68]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[69]  Neelesh A. Patankar,et al.  Erratum to "Immersed finite element method and its applications to biological systems" [Comput. Methods Appl. Mech. Engrg. 195 (2006) 1722-1749] (DOI:10.1016/j.cma.2005.05.049) , 2006 .

[70]  Xiaoming He,et al.  Interior penalty bilinear IFE discontinuous Galerkin methods for elliptic equations with discontinuous coefficient , 2010, J. Syst. Sci. Complex..

[71]  James P. Keener,et al.  Immersed Interface Methods for Neumann and Related Problems in Two and Three Dimensions , 2000, SIAM J. Sci. Comput..

[72]  Martin Schulz,et al.  Two-dimensional modelling of the river Rhine , 2002 .

[73]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[74]  Liqun Wang,et al.  Numerical method for solving matrix coefficient elliptic equation with sharp-edged interfaces , 2010, J. Comput. Phys..

[75]  Wei Shyy,et al.  Multigrid computations with the immersed boundary technique for multiphase flows , 2004 .

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

[77]  M. Lai,et al.  An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced Numerical Viscosity 1 , 2000 .

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

[79]  J. Sethian.  Evolution, implementation, and application of level set and fast marching methods for advancing fronts , 2001 .

[80]  Michael Oevermann,et al.  A Cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces , 2006, J. Comput. Phys..

[81]  Guo-Wei Wei,et al.  Multiscale molecular dynamics using the matched interface and boundary method , 2011, J. Comput. Phys..

[82]  Juan Galvis,et al.  BDDC methods for discontinuous Galerkin discretization of elliptic problems , 2007, J. Complex..

[83]  Juri D. Kandilarov.  Immersed Interface Method for a Reaction-Diffusion Equation with a Moving Own Concentrated Source , 2002, Numerical Methods and Application.

[84]  Isabelle Ramière,et al.  Convergence analysis of the Q1‐finite element method for elliptic problems with non‐boundary‐fitted meshes , 2008 .

[85]  Eftychios Sifakis,et al.  ' s personal copy A second order virtual node method for elliptic problems with interfaces and irregular domains , 2010 .

[86]  Shan Zhao,et al.  Minimal molecular surfaces and their applications , 2008, J. Comput. Chem..

[87]  G. R. Hadley.  High-accuracy finite-difference equations for dielectric waveguide analysis I: uniform regions and dielectric interfaces , 2002 .

[88]  A. Mayo.  The Fast Solution of Poisson’s and the Biharmonic Equations on Irregular Regions , 1984 .