The problem of computing the properties of the acoustic field exterior to a vibrating surface for the complete wavenumber range by the boundary element method is considered. A particular computational method based on the Kussmaul formulation is described. The method is derived through approximating the surface by a set of planar triangles and approximating the surface functions by a constant on each element. The method is successfully applied to test problems and to the Ricardo crankcase simulation rig.