Frequency-domain meshless solver for acoustic wave equation using a stable radial basis-finite difference (RBF-FD) algorithm with hybrid kernels

2017 ◽  
Author(s):  
Pankaj Mishra ◽  
Sankar Nath ◽  
Gregory Fasshauer ◽  
Mrinal Sen
Keyword(s):  
Wave Equation ◽  
Finite Difference ◽  
Acoustic Wave ◽  
Frequency Domain ◽  
Acoustic Wave Equation ◽  
Radial Basis

Adaptive 27-point finite difference frequency domain method for wave simulation of 3D acoustic wave equation

Geophysics ◽  
2021 ◽  
pp. 1-39
Author(s):  
Wenhao Xu ◽  
Bangyu Wu ◽  
Yang Zhong ◽  
Jinghuai Gao ◽  
Qing Huo Liu
Keyword(s):  
Wave Equation ◽  
Finite Difference ◽  
Acoustic Wave ◽  
Frequency Domain ◽  
Difference Frequency ◽  
Lookup Table ◽  
Acoustic Wave Equation ◽  
Wave Simulation ◽  
Finite Difference Frequency Domain ◽  
Point Stencil

The finite-difference frequency-domain (FDFD) method has important applications in the wave simulation of various wave equations. To promote the accuracy and efficiency for wave simulation with the FDFD method, we have developed a new 27-point FDFD stencil for 3D acoustic wave equation. In the developed stencil, the FDFD coefficients not only depend on the ratios of cell sizes in the x-, y-, and z-directions, but we also depend on the spatial sampling density (SD) in terms of the number of wavelengths per grid. The corresponding FDFD coefficients can be determined efficiently by making use of the plane-wave expression and the lookup table technique. We also develop a new way for designing an adaptive FDFD stencil by directly adding some correction terms to the conventional second-order FDFD stencil, which is simpler to use and easier to generalize. Corresponding dispersion analysis indicates that, compared to the optimal 27-point stencil derived from the average-derivative method (ADM), the developed adaptive 27-point stencil can reduce the required SD from approximately 4 to 2.2 points per wavelength (PPW) for a cubic mesh and to 2.7 PPW for a general cuboid mesh. Numerical examples of a 3D homogeneous model and SEG/EAGE salt-dome model indicate that the developed stencil is more accurate than the ADM 27-point stencil for cubic and general cuboid meshes, while requiring similar CPU time and computational memory as the ADM 27-point stencil for direct and iterative solvers.

Sign in / Sign up

Export Citation Format

Share Document