Abstract
Region transformation methods have been used by many authors to integrate singular multivariate integrands, the most notable being the Duffy transformation that maps squares onto triangles and cubes onto pyramids. In this paper, a new method that generalizes and improves the existing ones is provided. We show that behind these variable transformation methods there is an underlying linear structure involving the shape functions of the integration domains. This structure can be modified with further composition of the original transformation with regular mappings on the unit square in order to improve the performance of the numerical integration. The proposed method is neutral in the sense that the transformation composition is suitable for different types of singularities, without making any particular assumption on the regular part of the integrand. Moreover, this is achieved without a significant increase in the computational cost and making use of standard Gaussian quadrature rules only. We illustrate the efficiency of the proposed methods with comparative numerical examples. This technique proves to be both efficient and practical and has a better convergence than those based on classical transformations such as power or trigonometric maps.
Similar content being viewed by others
References
Belytschko, T., Gracie, R., Ventura, G.: A review of extended/generalized finite element methods for material modelling. Model. Simul. Mater. Sci. Eng. 17(4), 043001–043001 (2009)
Davis, P.J., Rabinowitz, P.: Methods of Numerical Integration, 2nd edn. Academic Press, London (1984)
Duffy, M.G.: Quadrature over a pyramid or cube of integrands with a singularity at a vertex. SIAM J. Numer. Anal. 19(6), 1260–1262 (1982)
Fairweather, G., Rizzo, F.J., Shippy, D.J.: Computation of double integrals in the boundary integral equation method. In: Vichnevetsky, R., Stepleman, R.S. (eds.) Advances in Computer Methods for Partial Differential Equations - III, vol. 36, pp. 331–334. IMACS Publ. Brussels, Belgium (1979)
Laborde, P., Pommier, J., Renard, Y., Salaün, M.: High-order extended finite element method for cracked domains. Int. J. Numer. Methods Eng. 64, 354–381 (2005)
Eisenberg, M.A., Malvern, L.E.: On Finite element integration in natural co-ordinates. Int. J. Numer. Methods Eng., Short Communications, no. 7.4: 574–575 (1973)
Gautschi, W.: Numerical integration over the square in the presence of algebraic/logarithmic singularities with an application to aerodynamics.Numer. Algoritm. 61, 275–290 (2012)
Mousavi, S.E., Sukumar, N.: Generalized Duffy transformation for integrating vertex singularities. Comput. Mech. 45(2–3), 127–140 (2010)
Mousavi, S.E., Sukumar, N.: Generalized Gaussian quadrature rules for discontinuities and crack singularities in the extended finite element method. Comput. Methods Appl. Mech. Eng. 199, 3237–3249 (2010)
Mustard, D., Lyness, J.N., Blatt, J.M.: Numerical quadrature in n dimensions. Comput. J. 6, 75–85 (1963)
Nagarajan, A., Mukherjee, S.: A mapping method for numerical evaluation of two-dimensional integrals with 1/r singularity. Comput. Mech. 12, 19–26 (1993)
Park, K., Pereira, J.P., Duarte, C.A., Paulino, G.H.: Integration of singular enrichment functions in the generalized/ extended finite element method for three-dimensional problems. Int. J. Numer. Methods Eng. 78(10), 1220–1257 (2009)
Rathod, H.T., Venkatesh, B.: Gauss Legendre - Gauss Jacobi Quadrature Rules over a Tetrahedral Region. Int. J. Math. Anal. 5(4), 189–198 (2011)
Sag, T.W., Szekeres, G.: Numerical evaluation of high-dimensional integrals. Math. Comput. 18, 245–253 (1964)
Stroud, A. H.: Approximate Calculation of Multiple Integrals. Prentice-Hall, Englewood Cliffs, New Jersey (1971)
Stroud, A.H., Secrest, D.: Gaussian Quadrature Formulas. Prentice-Hall Inc. (1966)
Tracey, D.M.: Finite elements for determination of crack tip elastic stress intensity factors. Eng. Fract. Mech. 3, 255–265 (1971)
Ungar, A.A.: Barycentric Calculus in Euclidean and Hyperbolic Geometry: A Comparative Introduction. World Scientific, Singapore (2010)
Author information
Authors and Affiliations
Corresponding author
Additional information
This work has been partially supported by Plan Nacional I+D+i,(MTM2011-23976), Spain
Rights and permissions
About this article
Cite this article
Cano, A., Moreno, C. A new method for numerical integration of singular functions on the plane. Numer Algor 68, 547–568 (2015). https://doi.org/10.1007/s11075-014-9860-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-014-9860-x