Absorbing boundary condition for the elastic wave equation

Geophysics ◽  
1988 ◽  
Vol 53 (5) ◽  
pp. 611-624 ◽  
Author(s):  
C. J. Randall

Extant absorbing boundary conditions for the elastic wave equation are generally effective only for waves nearly normally incident upon the boundary. High reflectivity is exhibited for waves traveling obliquely to the boundary. In this paper, a new and efficient absorbing boundary condition for two‐dimensional and three‐dimensional finite‐difference calculations of elastic wave propagation is presented. Compressional and shear components of the incident vector displacement fields are separated by calculating intermediary scalar potentials, allowing the use of Lindman’s boundary condition for scalar fields, which is highly absorbing for waves incident at any angle. The elastic medium is assumed to be homogeneous in the region immediately adjacent to the boundary. The reflectivity matrix of the resulting absorbing boundary for elastic waves is calculated, including the effects of finite‐difference truncation error. For effectively all angles of incidence, reflectivities are much smaller than those of the commonly employed paraxial absorbing boundaries, and the boundary condition is stable for any physical Poisson’s ratio. The nearly complete absorption predicted by the reflectivity matrix calculations, even at near grazing incidence, is demonstrated in a finite‐difference application.

Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 296-301 ◽  
Author(s):  
Chengbin Peng ◽  
M. Nafi Toksöz

Absorbing boundary conditions are widely used in numerical modeling of wave propagation in unbounded media to reduce reflections from artificial boundaries (Lindman, 1975; Clayton and Engquist, 1977; Reynolds, 1978; Liao et al., 1984; Cerjan et al., 1985; Randall, 1988; Higdon, 1991). We are interested in a particular absorbing boundary condition that has maximum absorbing ability with a minimum amount of computation and storage. This is practical for 3-D simulation of elastic wave propagation by a finite‐difference method. Peng and Toksöz (1994) developed a method to design a class of optimal absorbing boundary conditions for a given operator length. In this short note, we give a brief introduction to this technique, and we compare the optimal absorbing boundary conditions against those by Reynolds (1978) and Higdon (1991) using examples of 3-D elastic finite‐difference modeling on an nCUBE-2 parallel computer. In the Appendix, we also give explicit formulas for computing coefficients of the optimal absorbing boundary conditions.


Geophysics ◽  
1989 ◽  
Vol 54 (9) ◽  
pp. 1141-1152 ◽  
Author(s):  
C. J. Randall

This paper describes an absorbing boundary condition for finite‐difference modeling of elastic wave propagation in two and three dimensions. The boundary condition is particularly effective for obliquely incident waves, typically quite troublesome for absorbing boundaries. Analytical predictions of the boundary reflection coefficients of a few percent or less for angles of incidence up to 89° are verified in example finite‐difference applications. The algorithm is appropriate for use in a velocity‐stress finite‐difference (vs‐fd) formulation. It is computationally simpler than a similar absorbing boundary given previously for the standard displacement formulation. A second algorithm is presented which may be advantageous when the boundary of interest is exposed to strong evanescent waves. Both algorithms require that the adjacent elastic medium be locally homogeneous.


2014 ◽  
Vol 2014 ◽  
pp. 1-14 ◽  
Author(s):  
Ke-Yang Chen

Elastic wave equation simulation offers a way to study the wave propagation when creating seismic data. We implement an equivalent dual elastic wave separation equation to simulate the velocity, pressure, divergence, and curl fields in pure P- and S-modes, and apply it in full elastic wave numerical simulation. We give the complete derivations of explicit high-order staggered-grid finite-difference operators, stability condition, dispersion relation, and perfectly matched layer (PML) absorbing boundary condition, and present the resulting discretized formulas for the proposed elastic wave equation. The final numerical results of pure P- and S-modes are completely separated. Storage and computing time requirements are strongly reduced compared to the previous works. Numerical testing is used further to demonstrate the performance of the presented method.


Geophysics ◽  
1998 ◽  
Vol 63 (3) ◽  
pp. 1017-1023 ◽  
Author(s):  
Antonio Simone ◽  
Stig Hestholm

The problem of artificial reflections from grid boundaries in the numerical discretization of elastic and acoustic wave equations has long plagued geophysicists. Even if modern computers have made it possible to extend the synthetics over more wavelengths (equivalent to larger propagation distances), efficient absorption methods are still needed to minimize interference from unwanted reflections from the numerical grid boundaries. In this study, we examine applicabilities and stabilities of the optimal absorbing boundary condition (OABC) of Peng and Toksöz (1994, 1995) for 2-D and 3-D acoustic and elastic wave modeling. As a basis for comparison, we use exponential damping (ED) (Cerjan et. al., 1985), in which velocities and stresses are multiplied by progressively decreasing terms when approaching the boundaries of the numerical grid.


Geophysics ◽  
1986 ◽  
Vol 51 (7) ◽  
pp. 1499-1501 ◽  
Author(s):  
Kenneth D. Mahrer

One of the persistent problems with numerical solutions to the elastic wave equation is the finite size of the numerical grid. As with a physical body, the grid boundaries will reflect incident energy. If not eliminated or reduced substantially, these reflections will invade the grid interior and interfere with the desired solution. One method for eliminating reflections is creating a large and/or expanding grid. This method may be impractical since it can be quite costly in both computer time and memory. Another method is making the grid boundary “transparent” to outgoing energy. This method is ideally done by designing absorbing or nonreflecting boundaries which are mathematically equivalent to a one‐way, or outgoing, elastic wave equation only. In practice, an outgoing elastic wave equation is an approximation since the wave equation is not generally separable into outgoing and incoming parts. Two absorbing boundary condition approximations commonly used are those from Reynolds (Reynolds, 1978) and Clayton and Engquist, (Clayton and Engquist, 1977).


2009 ◽  
Vol 6 (2) ◽  
pp. 113-121 ◽  
Author(s):  
Zhen Qin ◽  
Minghui Lu ◽  
Xiaodong Zheng ◽  
Yao Yao ◽  
Cai Zhang ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document