A discontinuous Galerkin method for elliptic interface problems with application to electroporation

We solve elliptic interface problems using a discontinuous Galerkin (DG) method, for which discontinuities in the solution and in its normal derivatives are prescribed on an interface inside the domain. Standard ways to solve interface problems with finite element methods consist in enforcing the prescribed discontinuity of the solution in the finite element space. Here, we show that the DG method provides a natural framework to enforce both discontinuities weakly in the DG formulation, provided the triangulation of the domain is fitted to the interface. The resulting discretization leads to a symmetric system that can be efficiently solved with standard algorithms. The method is shown to be optimally convergent in the L2-norm. We apply our method to the numerical study of electroporation, a widely used medical technique with applications to gene therapy and cancer treatment. Mathematical models of electroporation involve elliptic problems with dynamic interface conditions. We discretize such problems into a sequence of elliptic interface problems that can be solved by our method. We obtain numerical results that agree with known exact solutions. Copyright © 2008 John Wiley & Sons, Ltd.

[1]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[2]  Douglas N. Arnold,et al.  Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems , 2001, SIAM J. Numer. Anal..

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

[4]  Igor Mozolevski,et al.  Discontinuous Galerkin finite element approximation of the two-dimensional Navier-Stokes equations in stream-function formulation , 2006 .

[5]  Peter Hansbo,et al.  Piecewise divergence‐free discontinuous Galerkin methods for Stokes flow , 2006 .

[6]  Xu-Dong Liu,et al.  Convergence of the ghost fluid method for elliptic equations with interfaces , 2003, Math. Comput..

[7]  PAUL CASTILLO,et al.  Performance of Discontinuous Galerkin Methods for Elliptic PDEs , 2002, SIAM J. Sci. Comput..

[8]  Airton Ramos,et al.  A new computational approach for electrical analysis of biological tissues. , 2003, Bioelectrochemistry.

[9]  D. Arnold.  An Interior Penalty Finite Element Method with Discontinuous Elements , 1982 .

[10]  Ilaria Perugia,et al.  Local discontinuous Galerkin methods for elliptic problems , 2001 .

[11]  Zhangxin Chen,et al.  Stability and convergence of mixed discontinuous finite element methods for second-order differential problems , 2003, J. Num. Math..

[12]  Boris Rubinsky,et al.  ELECTROPORATION: BIO-ELECTROCHEMICAL MASS TRANSFER AT THE NANO SCALE , 2000 .

[13]  Damijan Miklavčič,et al.  Time course of transmembrane voltage induced by time-varying electric fields—a method for theoretical analysis and its application , 1998 .

[14]  R. Plonsey,et al.  The transient subthreshold response of spherical and cylindrical cell models to extracellular stimulation , 1992, IEEE Transactions on Biomedical Engineering.

[15]  Ilaria Perugia,et al.  An A Priori Error Analysis of the Local Discontinuous Galerkin Method for Elliptic Problems , 2000, SIAM J. Numer. Anal..

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

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

[18]  Ja. A. Roitberg,et al.  a Theorem on Homeomorphisms for Elliptic Systems and its Applications , 1969 .

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

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