×

A decomposed immersed interface method for variable coefficient elliptic equations with non-smooth and discontinuous solutions. (English) Zbl 1052.65100

Summary: A second order accurate finite difference method is presented for solving two-dimensional variable coefficient elliptic equations on Cartesian grids, in which the coefficients, the source term, the solution and its derivatives may be non-smooth or discontinuous across an interface. A correction term is introduced to the standard central difference stencil so that the numerical discretization is well-defined across the interface.
We also propose a new method to approximate the correction term as part of the iterative procedure. The method is easy to implement since the correction term only needs to be added to the right-hand-side of the system. Therefore, the coefficient matrix remains symmetric and diagonally dominant, allowing for most standard solvers to be used. Numerical examples show good agreements with exact solutions, and the order of accuracy is comparable with other immersed interface methods.

MSC:

65N06 Finite difference methods for boundary value problems involving PDEs
35J25 Boundary value problems for second-order elliptic equations
35R05 PDEs with low regular coefficients and/or low regular data
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Adams, L.; Li, Z., The immersed interface/multigrid methods for interface problems, SIAM J. Sci. Comput., 24, 2, 463-479 (2002) · Zbl 1014.65099
[2] Beyer, R. P.; Leveque, R. J., Analysis of a one-dimensional model for the immersed boundary method, SIAM J. Numer. Anal., 29, 2, 332-364 (1992) · Zbl 0762.65052
[3] Brackbill, J. U.; Kothe, D. B.; Zemach, C., A continuum method for modeling surface tension, J. Comput. Phys., 100, 335-354 (1992) · Zbl 0775.76110
[4] Chang, Y. C.; Hou, T. Y.; Merriman, B.; Osher, S., A level set formulation of Eulerian interface capturing methods for incompressible fluid flows, J. Comput. Phys., 124, 449-464 (1996) · Zbl 0847.76048
[5] Fedkiw, R. P.; Aslam, T.; Merriman, B.; Osher, S., A non-oscillatory Eulerian approach to interfaces in multimaterial flows (The Ghost Fluid Method), J. Comput. Phys., 152, 457-492 (1999) · Zbl 0957.76052
[6] Huang, H.; Li, Z., Convergence analysis of the immersed interface method, IMA J. Numer. Anal., 19, 583-608 (1999) · Zbl 0940.65114
[7] Kang, M.; Fedkiw, R. P.; Liu, X.-D., A boundary condition capturing method for multiphase incompressible flow, J. Sci. Comput., 15, 3, 323-360 (2000) · Zbl 1049.76046
[8] Leveque, R. J.; Li, Z., The immersed interface method for elliptic equations with discontinuous coefficients and singular sources, SIAM J. Numer. Anal., 31, 4, 1019-1044 (1994) · Zbl 0811.65083
[9] Li, Z., A note on immersed interface method for three-dimensional elliptic equations, Comput. Math. Applic, 31, 3, 9-17 (1996) · Zbl 0876.65074
[10] Li, Z., A fast iterative algorithm for elliptic interface problems, SIAM J. Numer. Anal., 35, 1, 230-254 (1998) · Zbl 0915.65121
[11] Li, Z.; Lai, M.-C., The immersed interface method for the Navier-Stokes equations with singular forces, J. Comput. Phys, 171, 822-842 (2001) · Zbl 1065.76568
[12] Li, Z.; Ito, K., Maximum principle preserving schemes for interface problems with discontinuous coefficients, SIAM J. Sci. Comput., 23, 1, 330-361 (2001) · Zbl 1001.65115
[13] Liu, X.-D.; Fedkiw, R. P.; Kang, M., A boundary condition capturing method for Poisson’s equation on irregular domains, J. Comput. Phys., 160, 151-178 (2000) · Zbl 0958.65105
[14] Osher, S.; Sethian, J. A., Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations, J. Comput. Phys., 79, 12-49 (1988) · Zbl 0659.65132
[15] Peskin, C. S., Numerical analysis of blood flow in heart, J. Comput. Phys., 25, 220-252 (1977) · Zbl 0403.76100
[16] Sussman, M.; Smereka, P.; Osher, S., A level set approach for computing solutions to incompressible two-phase flow, J. Comput. Phys., 114, 146-159 (1994) · Zbl 0808.76077
[17] Unverdi, S. O.; Tryggvason, G., A front-tracking method for viscous, incompressible multi-fluid flows, J. Comput. Phys., 100, 25-37 (1992) · Zbl 0758.76047
[18] Wiegmann, A.; Bube, K. P., The explicit-jump immersed interface method: finite difference methods for PDEs with piecewise smooth solutions, SIAM J. Numer. Anal., 37, 3, 827-862 (2000) · Zbl 0948.65107
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.