Diese Seite ist aus Gründen der Barrierefreiheit optimiert für aktuelle Browser. Sollten Sie einen älteren Browser verwenden, kann es zu Einschränkungen der Darstellung und Benutzbarkeit der Website kommen!
Geophysics Homepage
Search:
Log in
print

Arbitrary High Order Finite Volume Schemes for Seismic Wave Propagation on Unstructured Meshes in 2D and 3D

Dumbser, M., M. Käser, and J. de la Puente (2007), Arbitrary High Order Finite Volume Schemes for Seismic Wave Propagation on Unstructured Meshes in 2D and 3D, Geophysical Journal International, 171(2), 665-694, doi:10.1111/j.1365-246X.2007.03421.x.

Abstract
We present a new numerical method to solve the heterogeneous anelastic seismic wave equations with arbitrary high order of accuracy in space and time on unstructured triangular and tetrahedral meshes in two and three space dimensions, respectively. Using the velocity-stress formulation provides a linear hyperbolic system of equations with source terms that is completed by additional equations for the anelastic functions including the strain history of the material. These additional equations result from the rheological model of the generalized Maxwell body and permit the incorporation of realistic attenuation properties of viscoelastic material accounting for the behaviour of elastic solids and viscous fluids. The proposed method relies on the Finite Volume (FV) approach where cell-averaged quantities are evolved in time by computing numerical fluxes at the element interfaces. The basic ingredient of the numerical flux function is the solution of Generalized Riemann Problems at the element interfaces according to the ADER approach of Toro et al., where the initial data is piecewise polynomial instead of piecewise constant as it was in the original first order FV scheme developed by Godunov. The ADER approach automatically produces a scheme of uniform high order of accuracy in space and time. The high order polynomials in space, needed as input for the numerical flux function, are obtained using a reconstruction operator acting on the cell averages. This reconstruction operator uses some techniques originally developed in the Discontinuous Galerkin (DG) Finite Element framework, namely hierarchical orthogonal basis functions in a reference element. In particular, in this article we pay special attention to underline the differences as well as the points in common with the ADER-DG schemes previously developed by the authors. The numerical convergence analysis demonstrates that the proposed Finite Volume schemes
provide very high order of accuracy even on unstructured tetrahedral meshes while computational cost for a desired accuracy can be reduced when applying higher order reconstructions. Applications to a series of well-acknowledged elastic and anelastic test cases and comparisons with analytic and numerical reference solutions, obtained by different well-established numerical methods, confirm the performance of the proposed method. Therefore, the development of the highly accurate ADER-FV approach for tetrahedral meshes including viscoelastic material provides a novel, flexible and efficient numerical technique to approach three-dimensional wave propagation problems including realistic attenuation and complex geometry.
BibTeX
@article{id603,
  author = {M. Dumbser and M. K{\"a}ser and J. de la Puente},
  journal = {Geophysical Journal International},
  number = {2},
  pages = {665-694},
  title = {{Arbitrary High Order Finite Volume Schemes for Seismic Wave Propagation on Unstructured Meshes in 2D and 3D}},
  volume = {171},
  year = {2007},
  doi = {10.1111/j.1365-246X.2007.03421.x},
}
EndNote
%0 Journal Article
%A Dumbser, M.
%A Käser, M.
%A de la Puente, J.
%D 2007
%N 2
%V 171
%J Geophysical Journal International
%P 665-694
%T Arbitrary High Order Finite Volume Schemes for Seismic Wave Propagation on Unstructured Meshes in 2D and 3D
ImprintPrivacy PolicyContact
Printed 06. Dec 2019 01:53