1887
Volume 67, Issue 7
  • E-ISSN: 1365-2478

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.

Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.12788
2019-04-15
2024-04-25
Loading full text...

Full text loading...

References

  1. BérengerJ.P.1994. A perfectly matched layer for the absorption of electromagnetic waves. Journal of Computational Physics114, 185–200.
    [Google Scholar]
  2. BérengerJ.P.2002a. Application of the CFS PML to the absorption of evanescent waves in waveguides. IEEE Microwave and Wireless Components Letters12, 218–220.
    [Google Scholar]
  3. BérengerJ.P.2002b. Numerical reflection from FDTD PMLs: A comparison of the split PML with the unsplit and CFS PMLs. IEEE Transactions on Antennas and Propagation50, 258–265.
    [Google Scholar]
  4. ChaljubE., KomatitschD., J. P.Vilotte, Y.Capdeville, B.Valette and G.Festa. 2007. Spectral‐element analysis in seismology. Advances in Geophysics, 48, 365–419.
    [Google Scholar]
  5. ChewW.C. and LiuQ.1996. Perfectly matched layers for elastodynamics. A new absorbing boundary condition. Journal of Computational Acoustics 4, 341–359.
    [Google Scholar]
  6. ChewW.C. and WeedonW.H.1994. A 3‐D perfectly matched medium from modified Maxwell's equations with stretched coordinates. Microwave and Optical Technology Letters7, 599–604.
    [Google Scholar]
  7. ChuC. and StoffaP.L.2010. Frequency‐domain modeling using implicit spatial finite difference operators. 80th SEG meeting, Denver, Colorado, USA, Expanded Abstracts, 3076–3080.
  8. CollinoF. and MonkP.1998. Optimizing the perfectly matched layer. Computer Methods in Applied Mechanics and Engineering164, 157–171.
    [Google Scholar]
  9. CollinoF. and TsogkaC.2001. Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media. Geophysics66, 294–307.
    [Google Scholar]
  10. FestaG. and VilotteJ.P.2005. The Newmark scheme as velocity‐stress time‐staggering: An efficient PML implementation for spectral‐element simulations of elastodynamics. Geophysical Journal International161, 789–812.
    [Google Scholar]
  11. FlyerN., FornbergB., BayonaV. and BarnettG.A.2016. On the role of polynomials in RBF‐FD approximations: I. Interpolation and accuracy. Journal of Computational Physics321, 21–38.
    [Google Scholar]
  12. FornbergB. and FlyerN.2015a. Solving PDEs with radial basis function. Acta Numerica24, 215–258.
    [Google Scholar]
  13. FornbergB. and FlyerN.2015b. Fast generation of 2‐D node distributions for mesh‐free PDE discretizations. Computers & Mathematics with Applications67, 531–544.
    [Google Scholar]
  14. FornbergB. and LehtoE.2011. Stabilization of RBF‐generated finite difference methods for convective PDEs. Journal of Computational Physics. 230, 2270–2285.
    [Google Scholar]
  15. FornbergB., LehtoE. and PowellC.2013. Stable calculation of Gaussian‐based RBF‐FD stencils. Computers & Mathematics with Applications65, 627–637.
    [Google Scholar]
  16. FriesT. and MatthiesH.2003. Classification and Overview of Mesh‐Free Methods. Braunschweig University of TechnologyInformatikbericht Nr. 2003‐3, Institute of Scientific Computing, Technical University Braunschweig, Brunswick, Germany.
    [Google Scholar]
  17. GuoH., WangS., GuoN. and ChenW.2012. Wave equation simulation by finite‐element method with perfectly matched layer. Advanced Materials Research524, 96–100.
    [Google Scholar]
  18. KomatitschD. and MartinR.2007. An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation. Geophysics72, SM155–SM167.
    [Google Scholar]
  19. KomatitschD. and TrompJ.2003. A perfectly matched layer absorbing boundary condition for the second‐order seismic wave equation. Geophysics Journal International154, 146–153.
    [Google Scholar]
  20. KuzuogluM. and MittraR.1996. Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers. IEEE Microwave and Guided Wave Letters6, 447–449.
    [Google Scholar]
  21. LarssonE. and FornbergB.2005. Theoretical and computational aspects of multivariate interpolation with increasingly flat radial basis functions. Computers & Mathematics with Applications49, 103–130.
    [Google Scholar]
  22. LeeG.H., ChungH.J. and ChoiC.K.2003. Adaptive crack propagation analysis with the element‐free Galerkin method. International Journal for Numerical Methods in Engineering56, 331–350.
    [Google Scholar]
  23. LiB., LiuY., SenM.K. and RenZ.M.2017. Time‐space‐domain mesh‐free finite difference based on least squares for 2D acoustic‐wave modeling. Geophysics82, T143–T157.
    [Google Scholar]
  24. LiuQ. and TaoJ.P.1997. The perfectly matched layer for acoustic waves in absorptive media. Journal of the Acoustical Society of America102, 2072–2082.
    [Google Scholar]
  25. LiuY. and SenM.K.2010. A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation. Geophysics75, A1–A6.
    [Google Scholar]
  26. LombardB. and J.Piraux. 2004. Numerical treatment of two‐dimensional interfaces for acoustic and elastic waves. Journal of Computational Physics195, 90–116.
    [Google Scholar]
  27. MaY., YuJ. and WangY.2014. An efficient complex‐frequency shifted‐perfectly matched layer for second‐order acoustic wave equation. International Journal for Numerical Methods in Engineering97, 130–148.
    [Google Scholar]
  28. MartinB., FornbergB. and St‐CyrA.2015. Seismic modeling with radial‐basis‐function‐generated finite differences. Geophysics80, T137–T146.
    [Google Scholar]
  29. MartinR. and KomatitschD.2009. An unsplit convolutional perfectly matched layer technique improved at grazing incidence for the viscoelastic wave equation. Geophysical Journal International179, 333–344.
    [Google Scholar]
  30. MinD.J., ShinC., KwonB.D. and ChungS.2000. Improved frequencydomain elastic wave modeling using weighted‐averaging difference operators. Geophysics65, 884–895.
    [Google Scholar]
  31. MoczoP., J.Kristek, V.Vavryčuk, R. J.Archuleta and L.Halada. 2002. 3D heterogeneous staggered‐grid finite‐difference modeling of seismic motion with volume harmonic and arithmetic averaging of elastic moduli and densities. Bulletin of Seismological Society of America. 92, 3042–3066.
    [Google Scholar]
  32. MoczoP., J.Kristek and M.Gális. 2014. The Finite‐Difference Modelling of Earthquake Motions: Waves and Ruptures, pp. 1–365. Cambridge University Press.
    [Google Scholar]
  33. PrattR.G. and WorthingtonM.H.1990a. Inverse theory applied to multi‐source cross‐hole tomography I: Acoustic wave‐equation method. Geophysical Prospecting38, 287–310.
    [Google Scholar]
  34. PrattR.G. and WorthingtonM.H.1990b. Inverse theory applied to multi‐source cross‐hole tomography II: Elastic wave‐equation method. Geophysical Prospecting38, 311–330.
    [Google Scholar]
  35. RamadanO.2003. Auxiliary differential equation formulation: An efficient implementation of the perfectly matched layer. IEEE Microwave and Wireless Components Letters13, 69–71.
    [Google Scholar]
  36. RenéM.2011. An efficient finite element time‐domain formulation for the elastic second‐order wave equation: A non‐split complex‐frequency shifted convolutional PML. International Journal for Numerical Methods in Engineering88, 951–973.
    [Google Scholar]
  37. RenZ. and LiuY.2013. A hybrid absorbing boundary condition for frequency‐domain finite‐difference modelling. Journal of Geophysics and Engineering10, 054003.
    [Google Scholar]
  38. RodenJ.A. and GedneyS.D.2000. Convolution PML (CPML): An efficient FDTD implementation of the CFS‐PML for arbitrary media. Microwave and Optical Technology Letters27, 334–339.
    [Google Scholar]
  39. ShiR., WangS. and ZhaoJ.2012. An unsplit complex‐frequency‐shifted PML based on matched Z‐transform for FDTD modelling of seismic wave equations. Journal of Geophysical Engineering9, 218–229.
    [Google Scholar]
  40. ShinC. and SohnH.1998. A frequency‐space 2‐D scalar wave extrapolator using extended 25‐point finite‐difference operator. Geophysics63, 289–296.
    [Google Scholar]
  41. TakekawaJ. and MikadaH.2016. An absorbing boundary condition for acoustic‐wave propagation using a mesh‐free method. Geophysics81, T145–T154.
    [Google Scholar]
  42. TakekawaJ., MikadaH. and ImamuraN.2015. A mesh‐free method with arbitrary‐order accuracy for acoustic wave propagation. Computers & Geosciences78, 15–25.
    [Google Scholar]
  43. TeixeiraF.L. and ChewW.C.1999a. Unified analysis of perfectly matched layers using differential forms. Microwave and Optical Technology Letters20, 124–126.
    [Google Scholar]
  44. TeixeiraF.L. and ChewW.C.1999b. Differential forms, metrics, and the reflectionless absorption of electromagnetic waves. Journal of Electromagnetic Waves and Applications13, 665–686.
    [Google Scholar]
  45. WenterodtC. and EstorffO.2009. Dispersion analysis of the mesh‐free radial point interpolation method for the Helmholtz equation. International Journal for Numerical Methods in Engineering77, 1670–1689.
    [Google Scholar]
  46. VirieuxJ., CalandraH. and PlessixR. É.2011. A review of the spectral, pseudo‐spectral, finite‐difference and finite‐element modelling techniques for geophysical imaging. Geophysical Prospecting59, 794–813.
    [Google Scholar]
  47. WenterodtC. and EstorffO.2011. Optimized mesh‐free methods for acoustics. Computer Methods in Applied Mechanics and Engineering200, 2223–2236.
    [Google Scholar]
  48. WittkeJ. and TezkanB.2014. Mesh‐free magnetotelluric modeling. Geophysical Journal International198, 1255–1268.
    [Google Scholar]
  49. XueZ., FomelS. and SunJ.2018. Increasing resolution of reverse‐time migration using time‐shift gathers. Geophysical Prospecting66, 726–735.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.12788
Loading
/content/journals/10.1111/1365-2478.12788
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): Acoustic; Finite difference; Modelling; Numerical study; Wave propagation

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error