Numerical approximations of singular source terms in differential equations. (English) Zbl 1115.76392

Summary: Singular terms in differential equations pose severe challenges for numerical approximations on regular grids. Regularization of the singularities is a very useful technique for their representation on the grid. We analyze such techniques for the practically preferred case of narrow support of the regularizations, extending our earlier results for wider support. The analysis also generalizes existing theory for one dimensional problems to multi-dimensions. New high order multidimensional techniques for differential equations and numerical quadrature are introduced based on the analysis and numerical results are presented. We also show that the common use of distance functions in level-set methods to extend one dimensional regularization to higher dimensions may produce \(O(1)\) errors.


76M25 Other numerical methods (fluid mechanics) (MSC2010)
74S30 Other numerical methods in solid mechanics (MSC2010)
65N99 Numerical methods for partial differential equations, boundary value problems
Full Text: DOI


[1] Beale, J.T.; Majda, A., Vortex methods II: higher order accuracy in two and three dimensions, Math. comput, 39, 29-52, (1982) · Zbl 0488.76025
[2] Beyer, R.P.; LeVeque, R.J., Analysis of a one-dimensional model for the immersed boundary method, SIAM J. num. anal, 29, 332-364, (1992) · Zbl 0762.65052
[3] Chorin, A.J., Numerical study of slightly viscous flow, J. fluid mech, 57, 785-796, (1973)
[4] Juric, D.; Tryggvason, G., A front-tracking method for dendritic solidification, J. comput. phys, 123, 127-148, (1996) · Zbl 0843.65093
[5] Kreiss, H.O.; Thomée, V.; Widlund, O., Smoothing of initial data and rates of convergence for parabolic difference equations, Commun. pur. appl. math, 23, 241-259, (1970) · Zbl 0188.41001
[6] G. Ledfelt, A thin wire sub cell model for arbitrary oriented wires for the fd – td method, in: G. Kristensson (Ed.), Proceedings of the EMB 98 - Electromagnetic Computations for Analysis and Design of Complex Systems, 1998, pp. 148-155
[7] LeVeque, R.J.; Li, Z.L., Immersed interface methods for Stokes flow with elastic boundaries or surface tension, SIAM J. sci. comput, 18, 709-735, (1997) · Zbl 0879.76061
[8] Monaghan, J.J., Extrapolating B splines for interpolation, J. comput. phys, 60, 253-262, (1985) · Zbl 0588.41005
[9] Osher, S.J.; Fedkiw, R.P., Level set methods and dynamic implicit surfaces, (2002), Springer Verlag Berlin
[10] Peskin, C.S., Numerical analysis of blood flow in the heart, J. comput. phys, 25, 220-252, (1977) · Zbl 0403.76100
[11] Peskin, C.S., The immersed boundary method, Acta numer, 11, 479-517, (2002) · Zbl 1123.74309
[12] Raviart, P.A., An analysis of particle methods, (), 253-262
[13] Sethian, J.A., Level set methods and fast marching methods. evolving interfaces in computational geometry, fluid mechanics, computer vision and materials science, (1999), Cambridge University Press Cambridge · Zbl 0973.76003
[14] Tornberg, A.K., Multi-dimensional quadrature of singular and discontinuous functions, Bit, 42, 644-669, (2002) · Zbl 1021.65010
[15] Tornberg, A.K.; Engquist, B., Regularization techniques for numerical approximation of PDEs with singularities, J. sci. comput, 19, 527-552, (2003) · Zbl 1035.65085
[16] Tornberg, A.K.; Engquist, B., The segment projection method for interface tracking, Commun. pur. appl. math, 56, 47-79, (2003) · Zbl 1205.76205
[17] Tryggvason, G.; Bunner, B.; Esmaeeli, A.; Juric, D.; Al-Rawahi, N.; Tauber, W.; Han, J.; Nas, S.; Jan, Y.J., A front-tracking method for the computations of multiphase flow, J. comput. phys, 169, 708-759, (2001) · Zbl 1047.76574
[18] Waldén, J., On the approximation of singular source terms in differential equations, Numer. meth. part D E, 15, 503-520, (1999) · Zbl 0938.65112
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. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.