Chapter 4
The Interior Helmholtz Problem

The underlying problem addressed in this Chapter is that of computing the Helmholtz field within an enclosed region and subject to a specified boundary condition. Let the closed region be denoted D with boundary S, as illustrated in Figure 4.1.

Fig 4.1. The domain of the interior Helmholtz problem.

Following the analysis of Subsection 1.3.1, the problem is equivalent to the solution of the Helmholtz equation

Ñ2 f(p) + k2 f(p) = 0    (p Î D) .
The boundary condition is assumed to take a general form
a(p) f(p) + b(p) v(p) = f(p)
(4.1)
where a, b and f are complex-valued functions defined on the boundary.

In this application the finite element method is now an established computational technique (ref. [68], [69]) and it may well also be often more efficient than the BEM . However, the BEM allows much more flexibility when the geometry of the domain is complicated and is a more natural method to apply if the domain is to be coupled with neighbouring domains in a wider computational method.

The volume of published research into the problem of determining the Helmholtz field within an enclosure or cavity by the boundary element method is minute in comparison to that of the corresponding exterior problem. There are two important reasons for this. The first is that the interior problem can be solved much more straightforwardly by the finite element method than the exterior problem, hence there is no pressing need for an alternative method. The other reason is that the development of the BEM for the exterior problem has been beset by difficulties; extensive research has been required to achieve reliable methods. The BEM for the interior problem is relatively straightforward.

The boundary element method for the solution of the interior Helmholtz boundary-value problem have been developed in Bell et al [9], Bernard et al [11], and Kipp and Bernard [39]. Further analysis or applications of the method are described in Seybert and Cheng [78], Cheng et al [19] and Kopuz and Lalor [53]. In this Chapter the application of the BEM to the interior Helmholtz problem is developed further so that the solution with a general boundary, boundary condition and incident field can be obtained. The subroutines HIBEM2, HIBEM3 and HIBEMA [33] for solving the two-, three- and axisymmetric three-dimensional problems are described and demonstrated.