COMBINATION OF REGIONS 201
not give an accurate indication of the behaviour of the system towards
infinity but the effect of the far region on the domain of interest is
introduced (see Bettess
Diffraction and refraction of waves by a
parabolic shoal, surmounted by a cylindrical island
Figure 7.18 shows the geometry of the problem and the element mesh
used. The finite element mesh is enclosed by a number of infinite
elements, the shape function in the radial direction is given by
(for time dependence e
). Here P(r) is a polynomial in r and L is the
so-called decay length, and κ is the wavenumber corresponding to the
frequency ω. This shape function satisfies the Sommerfeld radiation
L was chosen so that near the island the decay e~
matched the decay of the first term of the general series solution,
(/cr), where H
is the Hankel function of the zeroth order of the
second kind. See Bettess and Zienkiewicz
for a fuller explanation of
the above theory.
Figure 7.19, from Bettess and Zienkiewicz
, shows a comparison
of the elevations at the island itself compared with the analytic
solution given in Homas
and Vastano and Reid.
The agreement is
very good especially at moderate wavelengths.
The weakness in this method is that to determine the parameter L
some knowledge of the exact solution is required before it can be used.
However, the method is very efficient if this information is available.
7.5 COMBINATION OF FINITE AND BOUNDARY ELEMENTS
The idea of combining both techniques can be attributed to Wexler
who started to use integral equation solutions to represent the
unbounded field problem early in the 1970s, the advantage being that
this allowed for the use of appropriate conditions to represent the
infinite domain. The integral equation technique is also of interest
when regions of high stress or potential gradients exist, but finite
elements are adequate for other parts of a body and may be simpler to