A high-order Nyström discretization scheme for boundary integral equations defined on rotationally symmetric surfaces

DOI10.1016/J.JCP.2012.02.008zbMATH Open1250.65146arXiv1201.0077OpenAlexW2035372130MaRDI QIDQ447587FDOQ447587

Authors: P. Young, S. Hao, P. G. Martinsson 

Publication date: 4 September 2012

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: A scheme for rapidly and accurately computing solutions to boundary integral equations (BIEs) on rotationally symmetric surfaces in R^3 is presented. The scheme uses the Fourier transform to reduce the original BIE defined on a surface to a sequence of BIEs defined on a generating curve for the surface. It can handle loads that are not necessarily rotationally symmetric. Nystrom discretization is used to discretize the BIEs on the generating curve. The quadrature is a high-order Gaussian rule that is modified near the diagonal to retain high-order accuracy for singular kernels. The reduction in dimensionality, along with the use of high-order accurate quadratures, leads to small linear systems that can be inverted directly via, e.g., Gaussian elimination. This makes the scheme particularly fast in environments involving multiple right hand sides. It is demonstrated that for BIEs associated with the Laplace and Helmholtz equations, the kernel in the reduced equations can be evaluated very rapidly by exploiting recursion relations for Legendre functions. Numerical examples illustrate the performance of the scheme; in particular, it is demonstrated that for a BIE associated with Laplace's equation on a surface discretized using 320,800 points, the set-up phase of the algorithm takes 1 minute on a standard laptop, and then solves can be executed in 0.5 seconds.

Full work available at URL: https://arxiv.org/abs/1201.0077

Recommendations

• Relative \(\pi\)-blocks of \(\pi\)-separable groups. II
• An explicit kernel-split panel-based Nyström scheme for integral equations on axially symmetric surfaces
• Fast solution of boundary integral equations with the generalized Neumann kernel
• Numerical solution of the Helmholtz equation in 2D and 3D using a high-order Nyström discretization
• High-order accurate methods for Nyström discretization of integral equations on smooth curves in the plane

zbMATH Keywords

Gaussian eliminationnumerical examplesscatteringboundary integral equationdirect solverLegendre functionquadraturebody of revolutionhigh order discretizationGaussian ruleLaplace and Helmholtz equationsrotationally symmetric surfacesNyström discretization

Mathematics Subject Classification ID

Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Boundary element methods for boundary value problems involving PDEs (65N38)

Cites Work

• Title not available (Why is that?)
• Numerical quadratures for singular and hypersingular integrals
• Title not available (Why is that?)
• Numerical Methods for Special Functions
• Hybrid Gauss-Trapezoidal Quadrature Rules
• The Numerical Solution of Integral Equations of the Second Kind
• Title not available (Why is that?)
• High-Order Corrected Trapezoidal Quadrature Rules for Singular Functions
• A fast, high-order algorithm for the solution of surface scattering problems: Basic implementation, tests, and applications
• A fast direct solver for boundary integral equations in two dimensions
• Treatment of singular integular integrals in the boundary integral equation method for axisymmetric elastostatics problems
• On the condition number of boundary integral operators for the exterior Dirichlet problem for the Helmholtz equation
• A boundary integral approach to potential and elasticity problems for axisymmetric bodies with arbitrary boundary conditions
• A boundary element method for axisymmetric potential problems with non‐axisymmetric boundary conditions using fast Fourier transform
• Locally corrected Nyström method for EM scattering by bodies of revolution

Cited In (24)

• An efficient full-wave solver for eddy currents
• Comparison of integral equations for the Maxwell transmission problem with general permittivities
• The spectra of harmonic layer potential operators on domains with rotationally symmetric conical points
• Inverse scattering reconstruction of a three dimensional sound-soft axis-symmetric impenetrable object ^*
• A high-order wideband direct solver for electromagnetic scattering from bodies of revolution
• An FFT-accelerated direct solver for electromagnetic scattering from penetrable axisymmetric objects
• Robust integral formulations for electromagnetic scattering from three-dimensional cavities
• A fast solver for the narrow capture and narrow escape problems in the sphere
• A high-accuracy mode solver for acoustic scattering by a periodic array of axially symmetric obstacles
• Debye sources, Beltrami fields, and a complex structure on Maxwell fields
• A fast Fourier-Galerkin method solving boundary integral equations for the Helmholtz equation with exponential convergence
• A fast solver for elastic scattering from axisymmetric objects by boundary integral equations
• High-order accurate methods for Nyström discretization of integral equations on smooth curves in the plane
• On the numerical evaluation of the singular integrals of scattering theory
• Determination of normalized electric eigenfields in microwave cavities with sharp edges
• On the efficient evaluation of the azimuthal Fourier components of the Green's function for Helmholtz's equation in cylindrical coordinates
• A direct solver with \(O(N)\) complexity for integral equations on one-dimensional domains
• Convergence of the boundary integral method for interfacial Stokes flow
• An Extended Charge-Current Formulation of the Electromagnetic Transmission Problem
• An integral equation-based numerical solver for Taylor states in toroidal geometries
• An explicit kernel-split panel-based Nyström scheme for integral equations on axially symmetric surfaces
• Relative \(\pi\)-blocks of \(\pi\)-separable groups. II
• An efficient and highly accurate solver for multi-body acoustic scattering problems involving rotationally symmetric scatterers
• Efficient numerical solution of acoustic scattering from doubly-periodic arrays of axisymmetric objects

This page was built for publication: A high-order Nyström discretization scheme for boundary integral equations defined on rotationally symmetric surfaces

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q447587)