Hybrid absorbing boundary condition for piecewise smooth curved boundary in 2D acoustic finite difference modelling

Flexible computational domains and their corresponding irregular absorbing boundary conditions (ABCs) have previously been shown to improve the efficiency of finite difference (FD) modelling. However, these proposed ABCs for irregular boundaries are based on irregular grid methods. Although they can improve geometric flexibility of FD modelling, irregular grid methods are still complex to implement due to computationally expensive meshing process. To avoid complex mesh generation, FDs in mesh-free discretisation have been developed for nontrivial geometric settings, where scattered nodes can be placed suitably with respect to irregular boundaries and arbitrarily shaped anomalies without coordinate mapping or mesh-element forming. Radial-basis-function-generated FD (RBF-FD) has been proven successful in modelling seismic wave propagation based on mesh-free discretisation. Using RBF-FD, we develop a hybrid ABC for piecewise smooth curved boundary in 2D acoustic wave modelling based on straightforward expanding strategies for sampled boundary and corner nodes. The whole irregular computational domain is combined by an objective zone and a transition zone separated by piecewise smooth curved boundary. Nodal distribution in the objective zone is adaptive to the model structure through special alignment and varying density of nodes, while the transition zone are designed for solving one way wave equation more easily and accurately in the presence of piecewise smooth curved boundary. Modelling examples for homogenous and heterogeneous models with differently shaped computational domains demonstrate the effectiveness of our proposed method. Flexible computational domains with irregular absorbing boundaries can improve the efficiency of seismic modelling. We propose a hybrid absorbing boundary condition for piecewise smooth curved boundary in mesh-free discretisation based on radial-basis-function-generated finite difference modelling. Modelling examples in both homogeneous and heterogeneous media with flexibly shaped computational domain demonstrate the effectiveness of our proposed method.

[1]  Hongwei Gao,et al.  Irregular perfectly matched layers for 3D elastic wave modeling , 2011 .

[2]  Ying Rao,et al.  Seismic waveform simulation with pseudo-orthogonal grids for irregular topographic models , 2013 .

[3]  Bengt Fornberg,et al.  Seismic modeling with radial-basis-function-generated finite differences , 2015 .

[4]  Natasha Flyer,et al.  A 3-D RBF-FD solver for modeling the atmospheric global electric circuit with topography (GEC-RBFFD v1.0) , 2015 .

[5]  Jon Louis Bentley,et al.  Multidimensional binary search trees used for associative searching , 1975, CACM.

[6]  R. Higdon Absorbing boundary conditions for elastic waves , 1991 .

[7]  B. Engquist,et al.  Absorbing boundary conditions for acoustic and elastic wave equations , 1977, Bulletin of the Seismological Society of America.

[8]  Bengt Fornberg,et al.  A Stable Algorithm for Flat Radial Basis Functions on a Sphere , 2007, SIAM J. Sci. Comput..

[9]  D. Komatitsch,et al.  An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation , 2007 .

[10]  Zhang Jianfeng,et al.  Elastic wave modelling in 3D heterogeneous media: 3D grid method , 2002 .

[11]  Heung-Jin Chung,et al.  Adaptive crack propagation analysis with the element‐free Galerkin method , 2003 .

[12]  R. Kosloff,et al.  Absorbing boundaries for wave propagation problems , 1986 .

[13]  Wei Zhang,et al.  Unsplit complex frequency-shifted PML implementation using auxiliary differential equations for seismic wave modeling , 2010 .

[14]  B. Fornberg,et al.  Radial Basis Function-Generated Finite Differences: A Mesh-Free Method for Computational Geosciences , 2015 .

[15]  Bengt Fornberg,et al.  Solving PDEs with radial basis functions * , 2015, Acta Numerica.

[16]  Mrinal K. Sen,et al.  A hybrid absorbing boundary condition for elastic staggered-grid modelling , 2012 .

[17]  Stephen D. Gedney,et al.  Convolution PML (CPML): An efficient FDTD implementation of the CFS–PML for arbitrary media , 2000 .

[18]  X. Zhang,et al.  The AWWE-based hybrid absorbing boundary condition for finite-difference modeling and its application in reverse-time migration , 2015 .

[19]  Junichi Takekawa,et al.  A mesh-free method with arbitrary-order accuracy for acoustic wave propagation , 2015, Comput. Geosci..

[20]  Bengt Fornberg,et al.  Fast generation of 2-D node distributions for mesh-free PDE discretizations , 2015, Comput. Math. Appl..

[21]  Bengt Fornberg,et al.  On the role of polynomials in RBF-FD approximations: I. Interpolation and accuracy , 2016, J. Comput. Phys..

[22]  FornbergBengt,et al.  On the role of polynomials in RBF-FD approximations , 2016 .

[23]  Hitoshi Mikada,et al.  An absorbing boundary condition for acoustic-wave propagation using a mesh-free method , 2016 .

[24]  Mrinal K. Sen,et al.  A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation , 2010 .

[25]  Hongwei Gao,et al.  Implementation of perfectly matched layers in an arbitrary geometrical boundary for elastic wave modelling , 2008 .

[26]  Kim B. Olsen,et al.  On the implementation of perfectly matched layers in a three‐dimensional fourth‐order velocity‐stress finite difference scheme , 2003 .

[27]  Qing Huo Liu,et al.  The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media , 2001 .

[28]  Jinhai Zhang,et al.  Comparison of artificial absorbing boundaries for acoustic wave equation modelling , 2017 .

[29]  Erik Lehto,et al.  A guide to RBF-generated finite differences for nonlinear transport: Shallow water simulations on a sphere , 2012, J. Comput. Phys..

[30]  H. Matthies,et al.  Classification and Overview of Meshfree Methods , 2004 .

[31]  Henri Calandra,et al.  A review of the spectral, pseudo‐spectral, finite‐difference and finite‐element modelling techniques for geophysical imaging , 2011 .

[32]  Roland Martin,et al.  A High-Order Time and Space Formulation of the Unsplit Perfectly Matched Layer for the Seismic Wave Equation Using Auxiliary Differential Equations (ADE-PML) , 2010 .

[33]  Manuel Kindelan,et al.  Radial basis function interpolation in the limit of increasingly flat basis functions , 2016, J. Comput. Phys..

[34]  Zhang Jianfeng,et al.  P–SV-wave propagation in heterogeneous media: grid method , 1999 .

[35]  James S. Sochacki,et al.  Absorbing boundary conditions and surface waves , 1987 .

[36]  A. U.S.,et al.  Stable Computation of Multiquadric Interpolants for All Values of the Shape Parameter , 2003 .

[37]  Bradley Martin,et al.  Development of mesh-less computational algorithms for seismic exploration , 2013 .

[38]  J. Carcione,et al.  Seismic modeling , 1942 .

[39]  Aria Abubakar,et al.  Application of the Nearly Perfectly Matched Layer In Acoustic Wave Modeling , 2006 .

[40]  Yang Liu,et al.  A hybrid absorbing boundary condition for frequency-domain finite-difference modelling , 2013 .

[41]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[42]  M. Gunzburger,et al.  Boundary conditions for the numerical solution of elliptic equations in exterior regions , 1982 .

[43]  S. Gedney,et al.  An Auxiliary Differential Equation Formulation for the Complex-Frequency Shifted PML , 2010, IEEE Transactions on Antennas and Propagation.

[44]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[45]  Bengt Fornberg,et al.  On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs , 2017, J. Comput. Phys..

[46]  Bengt Fornberg,et al.  Stable Computation of Multiquadric Interpolants for All Values of the Shape Parameter , 2004 .

[47]  M. Guddati,et al.  Highly accurate absorbing boundary conditions for wide-angle wave equations , 2006 .

[48]  J. Sochacki Absorbing boundary conditions for the elastic wave equations , 1988 .

[49]  Laurence Eymard,et al.  Use of ship mean data for validating model and satellite flux fields during the FETCH experiment , 2003 .

[50]  Otto von Estorff,et al.  Dispersion analysis of the meshfree radial point interpolation method for the Helmholtz equation , 2009 .