Adaptive quadrature by expansion for layer potential evaluation in two dimensions. (English) Zbl 1446.65204
Summary: When solving partial differential equations using boundary integral equation methods, accurate evaluation of singular and nearly singular integrals in layer potentials is crucial. A recent scheme for this is quadrature by expansion (QBX), which solves the problem by locally approximating the potential using a local expansion centered at some distance from the source boundary. In this paper we introduce an extension of the QBX scheme in two dimensions (2D) denoted AQBX – adaptive quadrature by expansion – which combines QBX with an algorithm for automated selection of parameters, based on a target error tolerance. A key component in this algorithm is the ability to accurately estimate the numerical errors in the coefficients of the expansion. Combining previous results for flat panels with a procedure for taking the panel shape into account, we derive such error estimates for arbitrarily shaped boundaries in 2D that are discretized using panel-based Gauss-Legendre quadrature. Applying our scheme to numerical solutions of Dirichlet problems for the Laplace and Helmholtz equations, and also for solving these equations, we find that the scheme is able to satisfy a given target tolerance to within an order of magnitude, making it useful for practical applications. This represents a significant simplification over the original QBX algorithm, in which choosing a good set of parameters can be hard.
MSC:
| 65R20 | Numerical methods for integral equations |
| 65D30 | Numerical integration |
| 65D32 | Numerical quadrature and cubature formulas |
| 35J05 | Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation |
References:
| [1] | L. af Klinteberg and A.-K. Tornberg, {\it A fast integral equation method for solid particles in viscous flow using quadrature by expansion}, J. Comput. Phys., 326 (2016), pp. 420-445, . · Zbl 1373.76037 |
| [2] | L. af Klinteberg and A.-K. Tornberg, {\it Error estimation for quadrature by expansion in layer potential evaluation}, Adv. Comput. Math., 43 (2017), pp. 195-234, . · Zbl 1367.65154 |
| [3] | T. Askham and A. Cerfon, {\it An adaptive fast multipole accelerated Poisson solver for complex geometries}, J. Comput. Phys., 344 (2017), pp. 1-22, . · Zbl 1380.65413 |
| [4] | A. Barnett, B. Wu, and S. Veerapaneni, {\it Spectrally accurate quadratures for evaluation of layer potentials close to the boundary for the 2D Stokes and Laplace equations}, SIAM J. Sci. Comput., 37 (2015), pp. B519-B542, . · Zbl 1433.65323 |
| [5] | A. H. Barnett, {\it Evaluation of layer potentials close to the boundary for Laplace and Helmholtz problems on analytic planar domains}, SIAM J. Sci. Comput., 36 (2014), pp. A427-A451, . · Zbl 1298.65184 |
| [6] | W. Barrett, {\it Convergence properties of Gaussian quadrature formulae}, Comput. J., 3 (1961), pp. 272-277, . · Zbl 0098.31703 |
| [7] | T. Betcke, S. Chandler-Wilde, I. Graham, S. Langdon, and M. Lindner, {\it Condition number estimates for combined potential integral operators in acoustics and their boundary element discretisation}, Numer. Methods Partial Differential Equations, 27 (2011), pp. 31-69, . · Zbl 1272.76217 |
| [8] | J. D. Donaldson and D. Elliott, {\it A unified approach to quadrature rules with asymptotic estimates of their remainders}, SIAM J. Numer. Anal., 9 (1972), pp. 573-602, . · Zbl 0264.65020 |
| [9] | C. L. Epstein, L. Greengard, and A. Klöckner, {\it On the convergence of local expansions of layer potentials}, SIAM J. Numer. Anal., 51 (2013), pp. 2660-2679, . · Zbl 1281.65047 |
| [10] | Z. Gimbutas and L. Greengard, {\it FMMLIB2D}, (2012). |
| [11] | L. Greengard, D. Gueyffier, P.-G. Martinsson, and V. Rokhlin, {\it Fast direct solvers for integral equations in complex three-dimensional domains}, Acta Numer., 18 (2009), pp. 243-275, . · Zbl 1176.65141 |
| [12] | L. Greengard and V. Rokhlin, {\it A new version of the fast multipole method for the Laplace equation in three dimensions}, Acta Numer., 6 (1997), pp. 229-269, . · Zbl 0889.65115 |
| [13] | S. Hao, A. H. Barnett, P. G. Martinsson, and P. Young, {\it High-order accurate methods for Nyström discretization of integral equations on smooth curves in the plane}, Adv. Comput. Math., 40 (2014), pp. 245-272, . · Zbl 1300.65093 |
| [14] | J. Helsing and A. Holst, {\it Variants of an explicit kernel-split panel-based Nyström discretization scheme for Helmholtz boundary value problems}, Adv. Comput. Math., 41 (2015), pp. 691-708, . · Zbl 1319.65117 |
| [15] | J. Helsing and A. Karlsson, {\it An explicit kernel-split panel-based Nyström scheme for integral equations on axially symmetric surfaces}, J. Comput. Phys., 272 (2014), pp. 686-703, . · Zbl 1349.65709 |
| [16] | J. Helsing and R. Ojala, {\it On the evaluation of layer potentials close to their sources}, J. Comput. Phys., 227 (2008), pp. 2899-2921, . · Zbl 1135.65404 |
| [17] | K. L. Ho and L. Greengard, {\it A fast direct solver for structured linear systems by recursive skeletonization}, SIAM J. Sci. Comput., 34 (2012), pp. A2507-A2532, . · Zbl 1259.65062 |
| [18] | A. Klöckner, A. Barnett, L. Greengard, and M. O’Neil, {\it Quadrature by expansion: A new method for the evaluation of layer potentials}, J. Comput. Phys., 252 (2013), pp. 332-349, . · Zbl 1349.65094 |
| [19] | P. Martinsson and V. Rokhlin, {\it A fast direct solver for boundary integral equations in two dimensions}, J. Comput. Phys., 205 (2005), pp. 1-23, . · Zbl 1078.65112 |
| [20] | NIST, {\it Digital Library of Mathematical Functions}, Release 1.0.16, (2017). |
| [21] | M. Rachh, {\it Integral Equation Methods for Problems in Electrostatics, Elastostatics and Viscous Flow}, Ph.D. thesis, New York University, New York, 2015. |
| [22] | M. Rachh, A. Klöckner, and M. O’Neil, {\it Fast algorithms for quadrature by expansion I: Globally valid expansions}, J. Comput. Phys., 345 (2017), pp. 706-731, . · Zbl 1378.65074 |
| [23] | A. Rahimian, A. Barnett, and D. Zorin, {\it Ubiquitous evaluation of layer potentials using quadrature by kernel-independent expansion}, BIT, (2017), . · Zbl 1395.65151 |
| [24] | L. Ricketson, A. Cerfon, M. Rachh, and J. Freidberg, {\it Accurate derivative evaluation for any Grad-Shafranov solver}, J. Comput. Phys., 305 (2016), pp. 744-757, . · Zbl 1349.76258 |
| [25] | Y. Saad and M. H. Schultz, {\it GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems}, SIAM J. Sci. Stat. Comput., 7 (1986), pp. 856-869, . · Zbl 0599.65018 |
| [26] | V. Simoncini and D. B. Szyld, {\it Theory of inexact Krylov subspace methods and applications to scientific computing}, SIAM J. Sci. Comput., 25 (2003), pp. 454-477, . · Zbl 1048.65032 |
| [27] | L. N. Trefethen, {\it Approximation Theory and Approximation Practice}, Other Titles Appl. Math. 128. SIAM, Philadelphia, 2012. · Zbl 1264.41001 |
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.
