An absorbing boundary condition for acoustic-wave propagation using a mesh-free method

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. T145-T154 ◽  
Author(s):  
Junichi Takekawa ◽  
Hitoshi Mikada

We have developed an absorbing boundary condition for acoustic-wave propagation using a mesh-free method without sacrificing the flexibility of the mesh-free framework. When we simulate acoustic-wave propagation using a numerical method, artificial reflections from model edges induced by a truncated computational domain should be avoided. Although many absorbing boundary conditions have been developed, most of them have been based on a regular latticed alignment of grids or nodes, and the efficiency of such absorbing boundary conditions for irregular arrangement of grids or nodes has not been examined yet. We have studied the artificial reflections generated at the boundaries of a model for a mesh-free method, and we have proposed a novel approach for suppressing the artifacts. The method uses a hybrid approach with a transition zone, in which the wavefield is estimated by a weighted average of solutions from the one- and two-way wave equations. Numerical experiments indicate that the proposed method can provide good performance in suppression of the artificial edge reflections even for irregular distributions of calculation points in the vicinity of model edges.

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 ◽  
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.


2015 ◽  
Vol 78 ◽  
pp. 15-25 ◽  
Author(s):  
Junichi Takekawa ◽  
Hitoshi Mikada ◽  
Naoto Imamura

Author(s):  
Kadir Akbudak ◽  
Hatem Ltaief ◽  
Vincent Etienne ◽  
Rached Abdelkhalak ◽  
Thierry Tonellot ◽  
...  

The aim of this study is to design and implement an asynchronous computational scheme for solving the acoustic wave propagation equation with absorbing boundary conditions (ABCs) in the context of seismic imaging applications. While the convolutional perfectly matched layer (CPML) is typically used for ABCs in the oil and gas industry, its formulation further stresses memory accesses and decreases the arithmetic intensity at the physical domain boundaries. The challenges with CPML are twofold: (1) the strong, inherent data dependencies imposed on the explicit time-stepping scheme render asynchronous time integration cumbersome and (2) the idle time is further exacerbated by the load imbalance introduced among processing units. In fact, the CPML formulation of the ABCs requires expensive synchronization points, which may hinder the parallel performance of the overall asynchronous time integration. In particular, when deployed in conjunction with the multicore-optimized wavefront diamond temporal blocking (MWD-TB) approach for the inner domain points, it results in a major performance slow down. To relax CPML’s synchrony and mitigate the resulting load imbalance, we embed CPML’s calculation into MWD-TB’s inner loop and carry on the time integration with fine-grained computations in an asynchronous, holistic way. This comes at the price of storing transient results to alleviate dependencies from critical data hazards while maintaining the numerical accuracy of the original scheme. Performance and scalability results on various x86 architectures demonstrate the superiority of MWD-TB with CPML support against the standard spatial blocking on various grid sizes. To our knowledge, this is the first practical study that highlights the consolidation of CPML ABCs with asynchronous temporal blocking stencil computations.


Geophysics ◽  
1998 ◽  
Vol 63 (1) ◽  
pp. 231-243 ◽  
Author(s):  
Shunhua Cao ◽  
Stewart Greenhalgh

Four types of boundary conditions: Dirichlet, Neumann, transmitting, and modified transmitting, are derived by combining the damped wave equation with corresponding boundary conditions. The Dirichlet attenuating boundary condition is the easiest to implement. For an appropriate choice of attenuation parameter, it can achieve a boundary reflection coefficient of a few percent in a one‐wavelength wide zone. The Neumann‐attenuating boundary condition has characteristics similar to the Dirichlet attenuating boundary condition, but it is numerically more difficult to implement. Both the transmitting boundary condition and the modified transmitting boundary condition need an absorbing boundary condition at the termination of the attenuating region. The modified transmitting boundary condition is the most effective in the suppression of boundary reflections. For multidimensional modeling, there is no perfect absorbing boundary condition, and an approximate absorbing boundary condition is used.


1995 ◽  
Vol 85 (3) ◽  
pp. 954-956
Author(s):  
Ningya Cheng ◽  
Chuen Hon Cheng

Abstract In this short note, we derive the differential form of Liao's multi-transmitting formula. The reflection coefficient of the multi-transmitting formula in the acoustic case is obtained in closed form. These formulas are compared with Clayton-Engquist absorbing boundary condition and are shown to have a very close relationship.


Sign in / Sign up

Export Citation Format

Share Document