Seismic waves in heterogeneous material: subcell resolution of the discontinuous Galerkin method

Details Seismic waves in heterogeneous material: subcell resolution of the discontinuous Galerkin method Seismic waves in heterogeneous material: subcell resolution of the discontinuous Galerkin method

Jul 2010

adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2010geoji.182..250c&link_type=abstract

Geophysical Journal International, Volume 182, Issue 1, pp. 250-264.

Statistics

Computation

1

Numerical Solutions, Numerical Approximations And Analysis, Computational Seismology, Wave Propagation

Scientific paper

We present an important extension of the arbitrary high-order discontinuous Galerkin (DG) finite-element method to model 2-D elastic wave propagation in highly heterogeneous material. In this new approach we include space-variable coefficients to describe smooth or discontinuous material variations inside each element using the same numerical approximation strategy as for the velocity-stress variables in the formulation of the elastic wave equation. The combination of the DG method with a time integration scheme based on the solution of arbitrary accuracy derivatives Riemann problems still provides an explicit, one-step scheme which achieves arbitrary high-order accuracy in space and time. Compared to previous formulations the new scheme contains two additional terms in the form of volume integrals. We show that the increasing computational cost per element can be overcompensated due to the improved material representation inside each element as coarser meshes can be used which reduces the total number of elements and therefore computational time to reach a desired error level. We confirm the accuracy of the proposed scheme performing convergence tests and several numerical experiments considering smooth and highly heterogeneous material. As the approximation of the velocity and stress variables in the wave equation and of the material properties in the model can be chosen independently, we investigate the influence of the polynomial material representation on the accuracy of the synthetic seismograms with respect to computational cost. Moreover, we study the behaviour of the new method on strong material discontinuities, in the case where the mesh is not aligned with such a material interface. In this case second-order linear material approximation seems to be the best choice, with higher-order intra-cell approximation leading to potential instable behaviour. For all test cases we validate our solution against the well-established standard fourth-order finite difference and spectral element method.

Affiliated with

Also associated with

No associations

Rating

Seismic waves in heterogeneous material: subcell resolution of the discontinuous Galerkin method does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Seismic waves in heterogeneous material: subcell resolution of the discontinuous Galerkin method, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Seismic waves in heterogeneous material: subcell resolution of the discontinuous Galerkin method will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-957120