Perfectly matched layer boundary conditions for frequency‐domain acoustic wave simulation in the mesh‐free discretization. (15th April 2019)

• Record Type: Journal Article
• Authors: Liu, Xin; Liu, Yang; Ren, Zhiming; Li, Bei
• Is Part Of: Geophysical prospecting. Volume 67:Number 7(2019)
• Page Start: 1732
• Page End: 1744
• Publication Date: 2019-04-15
• Subjects: Acoustic -- Wave propagation -- Finite difference -- Modelling -- Numerical study; Prospecting -- Geophysical methods -- Periodicals; 622.15
• DOI: 10.1111/1365-2478.12788
• Languages: English
• ISSNs: 0016-8025

Abstract

ABSTRACT: Mesh‐free discretization, flexibly distributing nodes without computationally expensive meshing process, is able to deal with staircase problem, oversampling and undersampling problems and saves plenty of nodes through distributing nodes suitably with respect to irregular boundaries and model parameters. However, the time‐domain mesh‐free discretization usually exhibits poorer stability than that in regular grid discretization. In order to reach unconditional stability and easy implementation in parallel computing, we develop the frequency‐domain finite‐difference method in a mesh‐free discretization, incorporated with two perfectly matched layer boundary conditions. Furthermore, to maintain the flexibility of mesh‐free discretization, the nodes are still irregularly distributed in the absorbing zone, which complicates the situation of artificial boundary reflections. In this paper, we implement frequency‐domain acoustic wave modelling in a mesh‐free system. First, we present the perfectly matched layer boundary condition to suppress spurious reflections. Moreover, we develop the complex frequency shifted–perfectly matched layer boundary condition to improve the attenuation of grazing waves. In addition, we employ the radial‐basis‐function‐generated finite difference method in the mesh‐free discretization to calculate spatial derivatives. The numerical experiment on a rectangle homogeneous model shows the effectiveness of the perfectly matched layer boundary condition and the complex frequency shifted–perfectly matched layer boundary condition, and the latter one is better than the former one when absorbing large angle incident waves. The experiment on the Marmousi model suggests that the complex frequency shifted–perfectly matched layer boundary condition works well for complicated models.