Finite difference schemes for modelling the propagation of axisymmetric elastic longitudinal waves

An efficient finite difference shock-capturing scheme for the solution of direct seismic problems is constructed. Problem formulation is based on equations of the dynamics of elastic medium with axial symmetry. When implementating the scheme on multiprocessor computing systems, the two-cyclic splitting method with respect to spatial variables is used. One-dimensional systems of equations that arise in the context of splitting procedure are represented as subsystems for longitudinal, transverse and torsional waves. The case of longitudinal waves is considered in this paper. The results of simulations with the use of explicit grid-characteristic schemes and implicit schemes of the "predictor–corrector" type with controllable dissipation of energy are compared with exact solutions that describe propagation of monochromatic waves.