Fortran codes for computing the discrete Helmholtz integral operators

Abstract

In this paper Fortran subroutines for the evaluation of the discrete form of the Helmholtz integral operators L k, M k, M k t and N k for two-dimensional, three-dimensional and three-dimensional axisymmetric problems are described. The subroutines are useful in the solution of Helmholtz problems via boundary element and related methods. The subroutines have been designed to be easy to use, reliable and efficient. The subroutines are also flexible in that the quadrature rule is defined as a parameter and the library functions (such as the Hankel, exponential and square root functions) are called from external routines. The subroutines are demonstrated on test problems arising from the solution of the Neumann problem exterior to a closed boundary via the Burton and Miller equation.

Available Versions of this Item

• Fortran codes for computing the discrete Helmholtz integral operators. (deposited 08 Jul 2013 14:11) [Currently Displayed]
  • Fortran codes for computing the discrete Helmholtz integral operators. (deposited 01 Jun 2016 08:48)