STRESS–VELOCITY COMPLETE RADIATION BOUNDARY CONDITIONS
Abstract
A new high-order local Absorbing Boundary Condition (ABC) has been recently proposed for use on an artificial boundary for time-dependent elastic waves in unbounded domains, in two dimensions. It is based on the stress–velocity formulation of the elastodynamics problem, and on the general Complete Radiation Boundary Condition (CRBC) approach, originally devised by Hagstrom and Warburton in 2009. The work presented here is a sequel to previous work that concentrated on the stability of the scheme; this is the first known high-order ABC for elastodynamics which is long-time stable. Stability was established both theoretically and numerically. The present paper focuses on the accuracy of the scheme. In particular, two accuracy-related issues are investigated. First, the reflection coefficients associated with the new CRBC for different types of incident and reflected elastic waves are analyzed. Second, various choices of computational parameters for the CRBC, and their effect on the accuracy, are discussed. These choices include the optimal coefficients proposed by Hagstrom and Warburton for the acoustic case, and a simplified formula for these coefficients. A finite difference discretization is employed in space and time. Numerical examples are used to experiment with the scheme and demonstrate the above-mentioned accuracy issues.
