High-order full discretization for anisotropic wave equations

Two-dimensional linear wave equation in anisotropic media, on a rectangular domain with initial conditions and periodic boundary conditions, is considered. The energy of the problem is contemplated. The space discretization is reached by means of finite differences on a uniform grid, paying attention to the mixed derivative of the equation. The discrete energy of the semi-discrete problem is introduced. For the time integration of the system of ordinary differential equations obtained, a fourth order exponential splitting method, which is a geometric integrator, is proposed. This time integrator is efficient and easy to implement. The stability condition for time step and space step ratio is deduced. Numerical experiments displaying the good behavior in the long time integration and the efficiency of the numerical solution are provided.