Abstract:A numerical scheme is presented for solving the Helmholtz equation with Dirichlet or Neumann boundary conditions on piecewise smooth open curves, where the curves may have corners and multiple junctions. Existing integral equation methods for smooth open curves rely on analyzing the exact singularities of the density at endpoints for associated integral operators, explicitly extracting these singularities from the densities in the formulation, and using global quadrature to discretize the boundary integral equation. Extending these methods to handle curves with corners and multiple junctions is challenging because the singularity analysis becomes much more complex, and constructing high-order quadrature for discretizing layer potentials with singular and hypersingular kernels and singular densities is nontrivial. The proposed scheme is built upon the following two observations. First, the single-layer potential operator and the normal derivative of the double-layer potential operator serve as effective preconditioners for each other locally. Second, the recursively compressed inverse preconditioning (RCIP) method can be extended to address "implicit" second-kind integral equations. The scheme is high-order, adaptive, and capable of handling corners and multiple junctions without prior knowledge of the density singularity. It is also compatible with fast algorithms, such as the fast multipole method. The performance of the scheme is illustrated with several numerical examples.

A Fully Discrete High-Order Fast Multiscale Galerkin Method for Solving Boundary Integral Equations in a Domain with Corners

A Fast Multiscale Galerkin Method for Solving a Boundary Integral Equation in a Domain with Corners

A Fast Fourier-Galerkin Method for Solving Boundary Integral Equations on Torus-Shaped Surfaces

An Improved MOD Discontinuous Galerkin Method for Solving Time-Domain Electric Field Integral Equation

Discontinuous Galerkin Method Using Laguerre Polynomials for Solving a Time-Domain Electric Field Integral Equation

A Fast Fourier-Galerkin Method Solving Boundary Integral Equations for the Helmholtz Equation with Exponential Convergence

A Fully Discrete Fast Fourier–Galerkin Method Solving a Boundary Integral Equation for the Biharmonic Equation

The Helmholtz Dirichlet and Neumann problems on piecewise smooth open curves

Boundary treatment for high-order IMEX Runge-Kutta local discontinuous Galerkin schemes for multidimensional nonlinear parabolic PDEs

Discontinuous Galerkin Methods for the Laplace-Beltrami Operator on Point Cloud

Fast direct solvers for integral equations in complex three-dimensional domains

A discontinuous Galerkin discretization of elliptic problems with improved convergence properties using summation by parts operators

Discontinuous Galerkin method based on the reduced space for the nonlinear convection-diffusion-reaction equation

Discontinuous Galerkin Approximations to Elliptic and Parabolic Problems with a Dirac Line Source

A selective immersed discontinuous Galerkin method for elliptic interface problems

Multi-Level Solution Strategies for Implicit Time Discontinuous Galerkin Discretizations

A hybrid interpolation ACA accelerated method for parabolic boundary integral operators

A Hybrid Discontinuous Galerkin Scheme for Multi-scale Kinetic Equations

Parallel Skeletonization for Integral Equations in Evolving Multiply-Connected Domains

Efficient Exponential Integrator Finite Element Method for Semilinear Parabolic Equations

Discontinuous Galerkin Method Based on the Reduced Space for the Nonlinear Convection–Diffusion–Reaction Equation