A high-order integral algorithm for highly singular PDE solutions in Lipschitz domains

We present a new algorithm, based on integral equation formulations, for the solution of constant-coefficient elliptic partial differential equations (PDE) in closed two-dimensional domains with non-smooth boundaries; we focus on cases in which the integral-equation solutions as well as physically meaningful quantities (such as, stresses, electric/magnetic fields, etc.) tend to infinity at singular boundary points (corners). While, for simplicity, we restrict our discussion to integral equations associated with the Neumann problem for the Laplace equation, the proposed methodology applies to integral equations arising from other types of PDEs, including the Helmholtz, Maxwell, and linear elasticity equations. Our numerical results demonstrate excellent convergence as discretizations are refined, even around singular points at which solutions tend to infinity. We demonstrate the efficacy of this algorithm through applications to solution of Neumann problems for the Laplace operator over a variety of domains—including domains containing extremely sharp concave and convex corners, with angles as small as π/100 and as large as 199π/100.

[1]  M. Bleszynski,et al.  AIM: Adaptive integral method for solving large‐scale electromagnetic scattering and radiation problems , 1996 .

[2]  G. Fix,et al.  On the use of singular functions with finite element approximations , 1973 .

[3]  V. Maz'ya,et al.  The asymptotic form of the solutions of the integral equations of potential theory in the neighbourhood of the corner points of a contour , 1984 .

[4]  Martin Costabel,et al.  Quadrature and Collocation Methods for the Double Layer Potential on Polygons , 1993 .

[5]  Rainer Kress,et al.  A Nyström method for boundary integral equations in domains with corners , 1990 .

[6]  Mikhail Borsuk,et al.  Chapter 10 - Sharp estimates of solutions to the Robin boundary value problem for elliptic non divergence second order equations in a neighborhood of the conical point , 2006 .

[7]  Ernst P. Stephan,et al.  Boundary Integral Operators in Countably Normed Spaces , 1998 .

[8]  W. D. Evans,et al.  PARTIAL DIFFERENTIAL EQUATIONS , 1941 .

[9]  Neil M. Wigley,et al.  On a Method to Subtract off a Singularity at a Corner for the Dirichlet or Neumann Problem , 1969 .

[10]  Leszek Demkowicz,et al.  On an h-type mesh-refinement strategy based on minimization of interpolation errors☆ , 1985 .

[11]  John Strain,et al.  Locally Corrected Multidimensional Quadrature Rules for Singular Functions , 1995, SIAM J. Sci. Comput..

[12]  Giovanni Monegato,et al.  A polynomial collocation method for the numerical solution of weakly singular and singular integral equations on non‐smooth boundaries , 2003 .

[13]  V. Rokhlin.  Diagonal Forms of Translation Operators for the Helmholtz Equation in Three Dimensions , 1993 .

[14]  Youngmok Jeon,et al.  The collocation method for mixed boundary value problems on domains with curved polygonal boundaries , 1997 .

[15]  I. Babuska,et al.  Numerical treatment of vertex singularities and intensity factors for mixed boundary value problems for the Laplace equation in R 3 , 1994 .

[16]  Dan Givoli,et al.  A finite element method for domains with corners , 1992 .

[17]  R. Kress.  Linear Integral Equations , 1989 .

[18]  Norbert Heuer,et al.  A p‐adaptive algorithm for the BEM with the hypersingular operator on the plane screen , 2002 .

[19]  Xiaonan Wu,et al.  A Finite-Element Method for Laplace- and Helmholtz-Type Boundary Value Problems with Singularities , 1997 .

[20]  V. Kondratiev,et al.  Elliptic Boundary Value Problems of Second Order in Piecewise Smooth Domains , 2006 .

[21]  Ivan G. Graham,et al.  A high-order algorithm for obstacle scattering in three dimensions , 2004 .

[22]  Åke Björck,et al.  Numerical methods for least square problems , 1996 .

[23]  Państwowe Wydawnictwo Naukowe.,et al.  Mathematical models and methods in mechanics , 1985 .

[24]  V. A. Kondrat'ev,et al.  Boundary problems for elliptic equations in domains with conical or angular points , 1967 .

[25]  M. Costabel,et al.  The normal dervative of the double layer potential on polygons and galerkin approximation , 1983 .

[26]  P. Tong,et al.  Singular finite elements for the fracture analysis of V‐notched plate , 1980 .

[27]  David Jerison,et al.  The Neumann problem on Lipschitz domains , 1981 .

[28]  K. Atkinson.  The Numerical Solution of Integral Equations of the Second Kind , 1997 .

[29]  D. B. Huntera,et al.  A quadrature formula of Clenshaw–Curtis type for the Gegenbauer weight-function , 2005 .

[30]  P. Grisvard.  Singularities in Boundary Value Problems , 1992 .

[31]  Lloyd N. Trefethen,et al.  Is Gauss Quadrature Better than Clenshaw-Curtis? , 2008, SIAM Rev..

[32]  Lloyd N. Trefethen,et al.  The kink phenomenon in Fejér and Clenshaw–Curtis quadrature , 2007, Numerische Mathematik.

[33]  Ivan G. Graham,et al.  High-order methods for linear functionals of solutions of second kind integral equations , 1988 .

[34]  F. John.  Partial differential equations , 1967 .

[35]  Johannes Elschner,et al.  A discrete collocation method for Symm's integral equation on curves with corners , 1996 .

[36]  P. Grisvard.  Elliptic Problems in Nonsmooth Domains , 1985 .

[37]  A. Maue,et al.  Zur Formulierung eines allgemeinen Beugungs-problems durch eine Integralgleichung , 1949 .

[38]  G. Verchota.  Layer potentials and regularity for the Dirichlet problem for Laplace's equation in Lipschitz domains , 1984 .

[39]  Ivo Babuška,et al.  The post‐processing approach in the finite element method—Part 2: The calculation of stress intensity factors , 1984 .

[40]  T. Hughes,et al.  Techniques for developing ‘special’ finite element shape functions with particular reference to singularities , 1980 .

[41]  Xiaonan Wu,et al.  Discrete boundary conditions for elasticity problems with singularities , 2003 .

[42]  William E. Smith,et al.  Product integration with the Clenshaw-Curtis points: Implementation and error estimates , 1980 .

[43]  Michael Taylor,et al.  The ̄∂-Neumann Problem , 1996 .

[44]  M. Costabel,et al.  Singularities of Electromagnetic Fields¶in Polyhedral Domains , 2000 .

[45]  O. Bruno,et al.  A fast, high-order algorithm for the solution of surface scattering problems: basic implementation, tests, and applications , 2001 .

[46]  I. Babuška,et al.  Direct and inverse error estimates for finite elements with mesh refinements , 1979 .

[47]  E. Stephan,et al.  The hp-Version of the Boundary Element Method on Polygons , 1996 .

[48]  Neil M. Wigley,et al.  Stress intensity factors and improved convergence estimates at a corner , 1987 .

[49]  M. Stern.  Families of consistent conforming elements with singular derivative fields , 1979 .

[50]  G. Fix,et al.  On the accuracy of least squares methods in the presence of corner singularities , 1984 .

[51]  Matthias Maischak,et al.  The Hp-version of the Boundary Element Method in Ir the Basic Approximation Results , 1996 .

[52]  Jiming Song,et al.  Fast Illinois solver code (FISC) , 1998 .

[53]  Martin Costabel,et al.  Boundary integral equations for mixed boundary value problems in polygonal domains and Galerkin approximation , 1985 .

[54]  J. Roßmann,et al.  Elliptic Boundary Value Problems in Domains with Point Singularities , 2002 .

[55]  J. Whiteman.  The Mathematics of Finite Elements and Applications. , 1983 .