Mechanical Quadrature Method and Splitting Extrapolation for Solving Dirichlet Boundary Integral Equation of Helmholtz Equation on Polygons

摘要

We study the numerical solution of Helmholtz equation with Dirichlet boundary condition. Based on the potential theory, the problem can be converted into a boundary integral equation. We propose the mechanical quadrature method (MQM) using specific quadrature rule to deal with weakly singular integrals. Denote byhmthe mesh width of a curved edgeΓm  (m=1,…,d)of polygons. Then, the multivariate asymptotic error expansion of MQM accompanied withO(hm3)for all mesh widthshmis obtained. Hence, once discrete equations with coarse meshes are solved in parallel, the higher accuracy order of numerical approximations can be at leastO(hmax⁡5)by splitting extrapolation algorithm (SEA). A numerical example is provided to support our theoretical analysis.