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

DOI10.1016/J.JCP.2012.02.008MaRDI QIDQ447587zbMATH OpenOpenAlexFDO

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

Publication date 4 September 2012

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

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

zbMATH Keywords

body of revolutionboundary integral equationdirect solverGaussian eliminationGaussian rulehigh order discretizationLaplace and Helmholtz equationsLegendre functionnumerical examplesNyström discretizationquadraturerotationally symmetric surfacesscattering

Mathematics Subject Classification ID

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

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.

Recommendations

• Relative -blocks of -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

Cites work

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

Cited in

(30)

• Robust integral formulations for electromagnetic scattering from three-dimensional cavities
• An integral equation-based numerical solver for Taylor states in toroidal geometries
• The spectra of harmonic layer potential operators on domains with rotationally symmetric conical points
• On the efficient evaluation of the azimuthal Fourier components of the Green's function for Helmholtz's equation in cylindrical coordinates
• An efficient full-wave solver for eddy currents
• A fast solver for elastic scattering from axisymmetric objects by boundary integral equations
• 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
• A fast Fourier-Galerkin method solving boundary integral equations for the Helmholtz equation with exponential convergence
• Determination of normalized electric eigenfields in microwave cavities with sharp edges
• Debye sources, Beltrami fields, and a complex structure on Maxwell fields
• On the numerical evaluation of the singular integrals of scattering theory
• An explicit kernel-split panel-based Nyström scheme for integral equations on axially symmetric surfaces
• High-order accurate methods for Nyström discretization of integral equations on smooth curves in the plane
• An extended charge-current formulation of the electromagnetic transmission problem
• Inverse scattering reconstruction of a three dimensional sound-soft axis-symmetric impenetrable object
• Convergence of the boundary integral method for interfacial Stokes flow
• 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
• Relative -blocks of -separable groups. II
• 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
• A fast Fourier-Galerkin method for solving boundary integral equations on non-axisymmetric toroidal surfaces
• A direct solver with O(N) complexity for integral equations on one-dimensional domains
• High-order Nyström/convolution-quadrature solution of time-domain scattering from closed and open Lipschitz boundaries with Dirichlet and Neumann boundary conditions
• A Spectral Solver for Acoustic Scattering by Multiple Quasi-Axisymmetric Structures
• A high-accuracy PML-BIE solver for three-dimensional acoustic wave scattering in axisymmetric layered media
• Boundary integral equation analysis for spheroidal suspensions
• Tolerance-driven close evaluation of the Stokes double layer potential on axisymmetric surfaces
• Comparison of integral equations for the Maxwell transmission problem with general permittivities

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)