{"title":"Dispersion-dissipation analysis of quadrilateral- and triangular-based discontinuous Galerkin methods for the visco-acoustic wave equation","authors":"Rubing Han , Jiandong Huang , Xijun He","doi":"10.1016/j.jappgeo.2025.105931","DOIUrl":null,"url":null,"abstract":"<div><div>The discontinuous Galerkin method (DGM) has been extensively applied to numerically discretize acoustic and elastic wave equations. However, few studies are focused on viscous media. In this study, we conduct a comprehensive numerical dispersion-dissipation analysis of DGM for the visco-acoustic wave equation and simulate wave propagation in viscous media. The plane-wave analysis is based on the standard-linear-solid (SLS) -based model. A weighted Runge- Kutta (WRK) time scheme, Legendre polynomials, and a local Lax-Friedrichs flux are employed. The fully discrete analyses are implemented in quadrilateral and triangular elements. We consider two types of triangular elements, the quality factors and the number of mechanisms. Our results show that the quality factor leads to the numerical dispersion, the value of which is constant within a range of small sampling rate, but does not cause numerical dissipation. The visco-acoustic wave equation for different mechanisms for viscosity also introduces an inherent numerical-dispersion value that is not affected by the sampling rate. Meanwhile, we find that the stability of the quadrilateral mesh is stronger than that of the triangular mesh under the same order of basis functions. Several numerical experiments are provided to validate some theoretical findings. Seismic waves go through attenuation and phase distortion during propagation, and the numerical results indicate that the DGM is suitable for simulating seismic wave propagation in the SLS-based model.</div></div>","PeriodicalId":54882,"journal":{"name":"Journal of Applied Geophysics","volume":"243 ","pages":"Article 105931"},"PeriodicalIF":2.1000,"publicationDate":"2025-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Applied Geophysics","FirstCategoryId":"89","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S092698512500312X","RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"GEOSCIENCES, MULTIDISCIPLINARY","Score":null,"Total":0}
引用次数: 0
Abstract
The discontinuous Galerkin method (DGM) has been extensively applied to numerically discretize acoustic and elastic wave equations. However, few studies are focused on viscous media. In this study, we conduct a comprehensive numerical dispersion-dissipation analysis of DGM for the visco-acoustic wave equation and simulate wave propagation in viscous media. The plane-wave analysis is based on the standard-linear-solid (SLS) -based model. A weighted Runge- Kutta (WRK) time scheme, Legendre polynomials, and a local Lax-Friedrichs flux are employed. The fully discrete analyses are implemented in quadrilateral and triangular elements. We consider two types of triangular elements, the quality factors and the number of mechanisms. Our results show that the quality factor leads to the numerical dispersion, the value of which is constant within a range of small sampling rate, but does not cause numerical dissipation. The visco-acoustic wave equation for different mechanisms for viscosity also introduces an inherent numerical-dispersion value that is not affected by the sampling rate. Meanwhile, we find that the stability of the quadrilateral mesh is stronger than that of the triangular mesh under the same order of basis functions. Several numerical experiments are provided to validate some theoretical findings. Seismic waves go through attenuation and phase distortion during propagation, and the numerical results indicate that the DGM is suitable for simulating seismic wave propagation in the SLS-based model.
期刊介绍:
The Journal of Applied Geophysics with its key objective of responding to pertinent and timely needs, places particular emphasis on methodological developments and innovative applications of geophysical techniques for addressing environmental, engineering, and hydrological problems. Related topical research in exploration geophysics and in soil and rock physics is also covered by the Journal of Applied Geophysics.