Method of virtual sources using on-surface radiation conditions for the Helmholtz equation

by   Sebastian Acosta, et al.

We develop a novel method of virtual sources to formulate boundary integral equations for exterior wave propagation problems. However, by contrast to classical boundary integral formulations, we displace the singularity of the Green's function by a small distance h>0. As a result, the discretization can be performed on the actual physical boundary with continuous kernels so that any naive quadrature scheme can be used to approximate integral operators. Using on-surface radiation conditions, we combine single- and double-layer potential representations of the solution to arrive at a well-conditioned system upon discretization. The virtual displacement parameter h controls the conditioning of the discrete system. We provide mathematical guidance to choose h, in terms of the wavelength and mesh refinements, in order to strike a balance between accuracy and stability. Proof-of-concept implementations are presented, including piecewise linear and isogeometric element formulations in two- and three-dimensional settings. We observe exceptionally well-behaved spectra, and solve the corresponding systems using matrix-free GMRES iterations. The results are compared to analytical solutions for canonical problems. We conclude that the proposed method leads to accurate solutions and good stability for a wide range of wavelengths and mesh refinements.


page 5

page 7

page 9

page 10


Boundary integral formulations for transient linear thermoelasticity with combined-type boundary conditions

We study boundary integral formulations for an interior/exterior initial...

Time-Domain Boundary Integral Methods in Linear Thermoelasticity

This paper is concerned with the application of time domain boundary int...

Towards optimal boundary integral formulations of the Poisson-Boltzmann equation for molecular electrostatics

The Poisson-Boltzmann equation offers an efficient way to study electros...

A Single-Layer Dual-Mesh Boundary Element Method for Multiscale Electromagnetic Modeling of Penetrable Objects in Layered Media

A surface integral representation of Maxwell's equations allows the effi...

Accuracy and Conditioning of Surface-Source Based Near-Field to Far-Field Transformations

The conditioning and accuracy of various inverse surface-source formulat...

Combined field-only boundary integral equations for PEC electromagnetic scattering problem in spherical geometries

We analyze the well posedness of certain field-only boundary integral eq...

SLIM: A Well-Conditioned Single-Source Boundary Element Method for Modeling Lossy Conductors in Layered Media

The boundary element method (BEM) enables the efficient electromagnetic ...

Please sign up or login with your details

Forgot password? Click here to reset