Finite volume formulation of the MIB method for elliptic interface problems

The matched interface and boundary (MIB) method has a proven ability for delivering the second order accuracy in handling elliptic interface problems with arbitrarily complex interface geometries. However, its collocation formulation requires relatively high solution regularity. Finite volume method (FVM) has its merit in dealing with conservation law problems and its integral formulation works well with relatively low solution regularity. We propose an MIB-FVM to take the advantages of both MIB and FVM for solving elliptic interface problems. We construct the proposed method on Cartesian meshes with vertex-centered control volumes. A large number of numerical experiments are designed to validate the present method in both two dimensional (2D) and three dimensional (3D) domains. It is found that the proposed MIB-FVM achieves the second order convergence for elliptic interface problems with complex interface geometries in both $L_{\infty}$ and $L_2$ norms.

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

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

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

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

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

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

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

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

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

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

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

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

[13]  Yang Xiang,et al.  Discrete singular convolution and its application to the analysis of plates with internal supports. Part 1: Theory and algorithm , 2002 .

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

[33]  C. Engwer,et al.  An unfitted finite element method using discontinuous Galerkin , 2009 .

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

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

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

[37]  Guo-Wei Wei,et al.  Second order method for solving 3D elasticity equations with complex interfaces , 2014, J. Comput. Phys..

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

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

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

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

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

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

[44]  B. Fornberg.  CALCULATION OF WEIGHTS IN FINITE DIFFERENCE FORMULAS∗ , 1998 .

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

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

[47]  Guo-Wei Wei,et al.  Matched interface and boundary method for elasticity interface problems , 2014, J. Comput. Appl. Math..

[48]  Bengt Fornberg,et al.  Classroom Note: Calculation of Weights in Finite Difference Formulas , 1998, SIAM Rev..

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

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

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

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

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

[54]  Guo-Wei Wei,et al.  A Galerkin formulation of the MIB method for three dimensional elliptic interface problems , 2014, Comput. Math. Appl..

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

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

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

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

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

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

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

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

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

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

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

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

[67]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[68]  Chang-Ock Lee,et al.  A discontinuous Galerkin method for elliptic interface problems with application to electroporation , 2009 .

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

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

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