I. Harari et Tjr. Hughes, STUDIES OF DOMAIN-BASED FORMULATIONS FOR COMPUTING EXTERIOR PROBLEMS OF ACOUSTICS, International journal for numerical methods in engineering, 37(17), 1994, pp. 2935-2950
The propagation and decay of acoustic waves in exterior domains in an
essential ingredient in the study of fluid-structure interaction. A st
rategy must be devised to compute solutions over domains which are unb
ounded. Exact impedance conditions at an artificial external boundary
are specified by the DtN method, yielding an equivalent problem that i
s suitable for domain-based computation. The DtN boundary condition is
non-reflective, giving rise to exact (and thereby unique) solutions.
The truncated DtN operator, which is employed in practice, fails to in
hibit the reflection of higher modes, so that non-unique solutions may
occur at their harmonics. Simple expressions determine a sufficient n
umber of terms in the truncated operator for unique solutions at any g
iven wave number. There are three characteristic length scales in the
computational problem: the radius of the artificial boundary, the geom
etry of the body (represented by the internal boundary) and the mesh s
ize. Numerical studies examine the dependence of the conditioning of f
inite element coefficient matrices on the number of terms in the trunc
ated DtN operator vs. the wave number non-dimensionalized by each of t
he length scales. Analytic results regarding the number of terms suffi
cient for unique solutions are confirmed. As long as this criterion is
respected, no upper limit on the allowable wave number is detected. A
local approximation of the boundary conditions restores uniqueness fo
r all wave numbers.