A second-order-accurate symmetric discretization of the Poisson equation on irregular domains

In this paper, we consider the variable coefficient Poisson equation with Dirichlet boundary conditions on an irregular domain and show that one can obtain second-order accuracy with a rather simple discretization. Moreover, since our discretization matrix is symmetric, it can be inverted rather quickly as opposed to the more complicated nonsymmetric discretization matrices found in other second-order-accurate discretizations of this problem. Multidimensional computational results are presented to demonstrate the second-order accuracy of this numerical method. In addition, we use our approach to formulate a second-order-accurate symmetric implicit time discretization of the heat equation on irregular domains. Then we briefly consider Stefan problems.

[1]  M. Muir.  Physical Chemistry , 1888, Nature.

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

[3]  Gene H. Golub,et al.  Matrix computations , 1983 .

[4]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[5]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[6]  J. Sethian,et al.  Crystal growth and dendritic solidification , 1992 .

[7]  R. LeVeque,et al.  Analysis of a one-dimensional model for the immersed boundary method , 1992 .

[8]  R. Almgren.  Variational algorithms and pattern formation in dendritic solidification , 1993 .

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

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

[11]  D. Juric,et al.  A Front-Tracking Method for Dendritic Solidification , 1996 .

[12]  A. Schmidt.  Computation of Three Dimensional Dendrites with Finite Elements , 1996 .

[13]  A. Karma,et al.  Quantitative phase-field modeling of dendritic growth in two and three dimensions , 1996 .

[14]  S. Osher,et al.  A Simple Level Set Method for Solving Stefan Problems , 1997, Journal of Computational Physics.

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

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

[17]  Ronald Fedkiw,et al.  Regular Article: The Ghost Fluid Method for Deflagration and Detonation Discontinuities , 1999 .

[18]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

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

[20]  Ronald Fedkiw,et al.  An Isobaric Fix for the Overheating Problem in Multimaterial Compressible Flows , 1999 .

[21]  J. A. Sethian,et al.  Fast Marching Methods , 1999, SIAM Rev..

[22]  W. Shyy,et al.  Computation of Solid-Liquid Phase Fronts in the Sharp Interface Limit on Fixed Grids , 1999 .

[23]  W. P. Jones,et al.  Analysis of the cell-centred finite volume method for the diffusion equation , 2000 .

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

[25]  Dantzig,et al.  Computation of dendritic microstructures using a level set method , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[26]  Alain Karma,et al.  Multiscale finite-difference-diffusion-Monte-Carlo method for simulating dendritic solidification , 2000 .

[27]  Ronald Fedkiw,et al.  A Boundary Condition Capturing Method for Multiphase Incompressible Flow , 2000, J. Sci. Comput..

[28]  R. Fedkiw,et al.  A boundary condition capturing method for incompressible flame discontinuities , 2001 .