Open AccessA hybrid scheme for seismic modelling based on Galerkin method
Author(s) -
Liu Tao,
Hu Tianyue,
Sen Mrinal K.,
Yang Jinhua,
Wang Runqiu,
Wei Jianxin,
Wang Shangxu
Publication year - 2011
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.2011.05094.x
Subject(s) - discretization , finite element method , mass matrix , matrix (chemical analysis) , diagonal , mathematics , dispersion (optics) , spectral element method , spectral method , discontinuous galerkin method , finite difference method , finite difference , computer science , extended finite element method , mathematical analysis , geometry , structural engineering , physics , materials science , nuclear physics , neutrino , composite material , engineering , optics
SUMMARY Purely numerical methods for modelling seismic wave propagation are now fairly popular. In this paper, we report on the development of a hybrid finite element‐finite difference method. Our method first employs semi‐discretization of the finite element method in a part of the spatial domain (in the z ‐direction for 2‐D situation) to obtain a wave equation in weak form, and then uses the finite difference method to solve it. This offers significant advantages in applying non‐uniform grids with improved accuracy. To improve the computational efficiency, we introduce the spectral element method to replace the traditional finite element method, which leads to a diagonal mass matrix with sufficient accuracy. After that, we carry out detailed analyses of the dispersion behaviours for both uniform and non‐uniform cases; the dispersion curves demonstrate high precision of our method in numerical modelling, especially in dealing with non‐uniform grids. Some examples are presented to demonstrate performance of this method and confirm the analytic results.