It is well known that isotropic, nondispersive continuous hyperbolic problems become dispersive and anisotropic upon discretization. The purpose of this paper is to conduct a dispersion analysis of the nondissipative numerical approximations to plane wave motions in isotropic elastic solids. The discrete formulations considered are: an explicit, second-order accurate ®nite difference scheme, a consistent mass matrix formulation with linear quadrilateral elements and the corresponding lumped mass matrix formulation. Dispersion relation is derived for each of these formulations. In the context of the ®nite difference scheme, expressions for group velocity for both the shear and longitudinal waves are derived and the effect of using meshes of unequal size in x and y directions is studied. Results from numerical experiments con®rming the predictions of analysis are also presented.
IntroductionDispersion relation is the relation between the frequency and wavenumber. The ratio of the frequency and wavenumber is the phase speed whereas the gradient of the frequency with respect to the wavenumber is the group speed. For a nondispersive problem, the group speed does not depend on the frequency and is identical to the phase speed. For a dispersive problem on the other hand, the group speed depends on the frequency and in general is different from the phase speed.It is well known that when a nondispersive continuous problem is discretized, the discrete model becomes dispersive. In the case of one-dimensional problems, dispersion leads to waves of different lengths propagating at different speeds. In the case of two and three dimensions, not only do different wavelengths travel at different speeds but they also travel in wrong directions. That is, the discrete models become anisotropic even if the continuous problem is isotropic.A survey of the literature shows that a large body of work has been done on analyzing the dispersion errors induced due to numerical approximations. Much of the work related to one-dimensional and two-dimensional scalar wave equations can be found in the monograph by Vichnevetsky [10]. Dispersion analysis of a three-dimensional scalar wave equation was conducted by Abboud and Pinsky [1]. The excellent review article by Trefethen [9] surveys the role of group velocity in ®nite difference approximations to one-dimensional and two-dimensional scalar wave equations. Recently, Monk and Parrott [8] studied the dispersion errors introduced in various ®nite element approximations to Maxwell's equations.However, surprisingly, there appears to be very little work done on the dispersion analysis of numerical approximations to the wave propagation problems in isotropic, elastic solids in two and three dimensions. Only the one-dimensional elastic wave propagation problem appears to have received the attention of several researchers previously (see [3,7] and the references cited therein). In addition, most of this work appears to consider only one type of wave propagating through the body.In the case of two and th...