Elastodynamics of a non-ideal interface: Application to crack and fracture scattering

J.M. Carcione,
Osservatorio Geofisico Sperimentale, P. O. Box 2011, 34016 Trieste, Italy

Abstract:

This work introduces a numerical algorithm for solving wave propagation in the presence of an imperfect contact between two isotropic, elastic and heterogeneous media. Non-ideal interfaces of general type can be modeled as boundary discontinuities of the displacement u and its first time derivative (the particle velocity v). The stress field is continuous, and the quantity tex2html_wrap_inline29, where the brackets denote discontinuities across the interface, is equal to the corresponding stress component. The specific stiffness tex2html_wrap_inline31 introduces frequency-dependence and phase changes in the interface response. On the other hand, the specific viscosity tex2html_wrap_inline33 is related to the energy loss. It is shown here that, in the velocity stress formulation of the wave equation, such a model is described by Maxwell relaxation-like functions. I compute the reflection and transmission coefficients in terms of the corresponding incident propagation angle and complex moduli, together with the energy dissipated at the interface. This analysis characterizes the properties of the non-ideal interface. The numerical method is based on a domain decomposition technique that assigns a different mesh to each side of the interface. As stated above, the effects of the interface on wave propagation are modeled through the boundary conditions, that require a special boundary treatment based on characteristic variables. The algorithm solves the velocity-stress w ave equations and two additional first-order differential e quations (in two-dimensional space) in the displacement discontinuity. For each mesh, the spatial derivatives normal to the interface are solved by the Chebyshev method, and the spatial derivatives parallel to the interface are computed with the Fourier method. The algorithm allows general material variability. The modeling is applied to the problems of crack and fracture scattering.

AGU Index Terms: 0900 Exploration geophysics; Mineralogy, Petrology, and Rock Chemistry; 0935 Seismic methods; 5102 Acoustic properties; 5144 Wave attenuation
Keywords/Free Terms: Crack, fracture, scattering, grid method.

JGR-Solid Earth 96JB02658
Vol. 101 , No. B12 , p. 28,177


© 1996 AGU