Stable high-order absorbing boundary condition based on new continued fraction for scalar wave propagation in unbounded multilayer media

2018 ◽  
Vol 334 ◽  
pp. 111-137 ◽  
Author(s):  
Mi Zhao ◽  
Lihua Wu ◽  
Xiuli Du ◽  
Zilan Zhong ◽  
Chengshun Xu ◽  
...  
2019 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Huifang Li ◽  
Mi Zhao ◽  
Lihua Wu ◽  
Piguang Wang ◽  
Xiuli Du

Purpose The purpose of this paper is to propose a stable high-order absorbing boundary condition (ABC) based on new continued fraction for scalar wave propagation in 2D and 3D unbounded layers. Design/methodology/approach The ABC is obtained based on continued fraction (CF) expansion of the frequency-domain dynamic stiffness coefficient (DtN kernel) on the artificial boundary of a truncated infinite domain. The CF which has been used to the thin layer method in [69] will be applied to the DtN method to develop a time-domain high-order ABC for the transient scalar wave propagation in 2D. Furthermore, a new stable composite-CF is proposed in this study for 3D unbounded layers by nesting the above CF for 2D layer and another CF. Findings The ABS has been transformed from frequency to time domain by using the auxiliary variable technique. The high-order time-domain ABC can couple seamlessly with the finite element method. The instability of the ABC-FEM coupled system is discussed and cured. Originality/value This manuscript establishes a stable high-order time-domain ABC for the scalar wave equation in 2D and 3D unbounded layers, which is based on the new continued fraction. The high-order time-domain ABC can couple seamlessly with the finite element method. The instability of the coupled system is discussed and cured.


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.


Sign in / Sign up

Export Citation Format

Share Document