The Discrete Helmholtz Operators

Linear Helmholtz or Helmholtz problems obviously give rise to a range of integral equation formulations, for example depending on whether the Helmholtz field lies in an interior or exterior domain. However, the integral equations that arise in all problems in the same spatial dimensions contain similar integral operators. For example a computational method for evaluating the discrete from of the integral operators in an three-dimensional exterior Helmholtz problem can also be used in any other three-dimensional Helmholtz problem.

For each dimensional space it is possible to develop a module for computing the discrete form of the integral operators that is common to the interior, exterior and modal analysis subroutines. The purpose of this section is to show how the discrete forms of the integral operators in the three spatial dimensions are computed and to introduce the Fortran subroutines H2LC, H3LC and H3ALC have been developed in order to carry this out. The naming of the subroutines is such that the H represents Helmholtz, the 2, 3, 3A identifies the dimensionality, the LC represents the linear boundary approximation with a constant function approximation on each panel.

A set of test runs for the subroutines are also available: H2LC_T, H3LC_T and H3ALC_T.

The subroutines were originally published in the paper " Fortran Codes for Computing the Discrete Helmholtz Integral Operators". The subroutines were applied to acoustics problems in the book " The Boundary Element Method in Acoustics" and the particular associated codes can be downloaded from this site www.boundary-element-method.com .