A transparent boundary technique for numerical modeling of elastic waves

Geophysics ◽  
1999 ◽  
Vol 64 (3) ◽  
pp. 963-966 ◽  
Author(s):  
Jianlin Zhu

In numerical modeling of wave motions, strong reflections from artificial model boundaries may contaminate or mask true reflections from the interior model interfaces. Hence, developing a kind of exterior model boundary transparent to the outgoing waves is of critical importance. Among proposed solutions, e.g., Smith (1974), Kausel and Tassoulas (1981), and Higdon (1991), the most widely used may be the Clayton and Engquist (1977) method of absorbing boundary conditions, based on paraxial approximations for acoustic and elastic‐wave equations. However, absorbing boundary conditions make the reflection coefficients zero only for normal incidence, and suppression of reflected S-waves (Clayton and Engquist, 1977) becomes poorer as the ratio of P- to S-wave velocity ([Formula: see text]) becomes larger.

1977 ◽  
Vol 67 (6) ◽  
pp. 1529-1540 ◽  
Author(s):  
Robert Clayton ◽  
Björn Engquist

abstract Boundary conditions are derived for numerical wave simulation that minimize artificial reflections from the edges of the domain of computation. In this way acoustic and elastic wave propagation in a limited area can be efficiently used to describe physical behavior in an unbounded domain. The boundary conditions are based on paraxial approximations of the scalar and elastic wave equations. They are computationally inexpensive and simple to apply, and they reduce reflections over a wide range of incident angles.


2013 ◽  
Vol 21 (02) ◽  
pp. 1250028 ◽  
Author(s):  
IGOR SHEVCHENKO ◽  
MANFRED KALTENBACHER ◽  
BARBARA WOHLMUTH

In this work, new absorbing boundary conditions (ABCs) for a wave equation with a temperature-dependent speed of sound are proposed. Based on the theory of pseudo-differential calculus, first- and second-order ABCs for the one- and two-dimensional wave equations are derived. Both boundary conditions are local in space and time. The well-posedness of the wave equation with the developed ABCs is shown through the reduction of the original problem to an equivalent one for which the uniqueness and existence of the solution has already been established. Although the second-order ABC is more accurate, the numerical realization is more challenging. Here we use a Lagrange multiplier approach which fits into the abstract framework of saddle point formulations and yields stable results. Numerical examples illustrating stability, accuracy and flexibility of the ABCs are given. As a test setting, we perform computations for a high-intensity focused ultrasound (HIFU) application, which is a typical thermo-acoustic multi-physics problem.


Sign in / Sign up

Export Citation Format

Share Document