Theory for the Exterior Field Calculation: The Helmholtz-Kirchhoff Integral
The Acoustics Module has functionality to evaluate the acoustic pressure field in the exterior field region of the model (outside the computational domain). This section gives the relevant definitions and mathematical background as well as some general advice for analyzing the exterior field. Details about how to use the exterior field functionality is described in Exterior Field Calculation and in the modeling section Evaluating the Acoustic Field in the Exterior: Near- and Far-Field.
The Helmholtz-Kirchhoff Integral Representation
In many cases, solving the acoustic Helmholtz equation everywhere in the domain where results are requested is neither practical nor necessary. For homogeneous media, the solution anywhere outside a closed surface containing all sources and scatterers can be written as a boundary integral in terms of quantities evaluated on the surface. To evaluate this Helmholtz-Kirchhoff   integral, it is necessary to know both Dirichlet and Neumann values on the surface. Applied to acoustics, this means that if the pressure and its normal derivative (which is related to the normal velocity) is known on a closed surface, the acoustic field can be calculated at any point outside.
In general, the solution p to Helmholtz’ equation
in the homogeneous domain exterior to a closed surface, S, can be explicitly expressed in terms of the values of  p and its normal derivative on S:
Here the coordinate vector r parameterizes S. The unit vector n is the outward normal to the exterior infinite domain; thus, n points into  the domain that S encloses. The function G (Rr) is a Green’s function satisfying
This essentially means that the Green’s function, seen as a function of  r, is an outgoing traveling wave excited by a simple source at R. In 3D, the Green’s function is therefore:
In 2D, the Green’s function contains a Hankel function instead of the exponential:
Inserting the 3D Green’s function in the general representation formula gives:
(2-23)
while in 2D, the Hankel function leads to a slightly different expression:
(2-24)
For axially symmetric geometries, the full 3D integral must be evaluated. The Acoustics Module uses an adaptive numerical quadrature in the azimuthal direction on a fictitious revolved geometry in addition to the standard mesh-based quadrature in the rz-plane.
The default in the Exterior Field Calculation feature is to evaluate the full Helmholtz-Kirchhoff integral given in Equation 2-23 and Equation 2-24.
The Far-Field Limit
The full Helmholtz-Kirchhoff integral gives the pressure at any point at a finite distance from the source surface, but the numerical integration tends to lose accuracy at very large distances. At the same time, in some applications, the quantity of interest is the far-field radiation pattern, which can be defined as the limit of r | p | when r goes to infinity in a given direction.
Taking the limit of Equation 2-23 when |R| goes to infinity and ignoring the rapidly oscillating phase factor, the far field, pfar is defined as
The relevant quantity is | pfar| rather than pfar because the phase of the latter is undefined. For the same reason, only the direction of  R is important, not its magnitude.
Because Hankel functions asymptotically approach an exponential, the limiting 2D integral is remarkably similar to that in the 3D case:
For axially symmetric geometries, the azimuthal integral of the limiting 3D case can be handled analytically, which leads to a rather complicated expression but avoids the numerical quadrature required in the general case. For zero azimuthal mode number m = 0, the expression is:
(2-25)
In this integral, r and z are the radial and axial components of r, while R and Z are the radial and axial components of R.
To evaluate the pressure in the far-field limit according to the equations in this section, set the Type of integral option to the Far-field integral approximation for r  option in the Exterior Field Calculation section in the Settings window for the feature. See Exterior Field Calculation.
The Elkernel Element
These integrals can be implemented as integration coupling variables in COMSOL Multiphysics. However, such an approach is very inefficient because then the simple structure of the integration kernels cannot be exploited. In the Acoustics Module, convolution integrals of this type are therefore evaluated in optimized codes that hides all details from the user.