Grid-characteristic method on unstructured tetrahedral meshes

The goal of this paper is to develop a grid-characteristic method intended for high-performance computer systems and implemented on unstructured tetrahedral hierarchical meshes with the use of a multiple time step and high-order interpolation, including interpolation with a limiter, piecewise parabolic interpolation, and monotone interpolation. The method is designed for simulating complex three-dimensional dynamical processes in heterogeneous media. It involves accurately stated contact conditions and produces physically correct solutions of problems in seismology and seismic exploration. Hierarchical meshes make it possible to take into account numerous inhomogeneous inclusions (cracks, cavities, etc.) and to solve problems in a real-life formulation. The grid-characteristic method enables the use of a multiple time step. As a result, the computation time is considerably reduced and the efficiency of the method is raised. The method is parallelized on a computer cluster with an optimal use of system resources.