Section: New Results
Numerical methods combining local time stepping and mixed hybrid elements for the terrestrial migration
We work on the development of a software for the Reverse Time Migration of acoustic waves which combines mixed hybrid elements in space and a local time stepping scheme. The discretization in space allows us to take the topography into account, which usually outperforms finite difference schemes. In order to illustrate this point, we have compared the solution obtained with a Dicontinuous Galerkin Finite Element Method (DGFEM) to the one obtained with finite difference methods or with the GSP method and we have shown that DGFEM gives much more accurate results. We have also compared the solution obtained with DGFEM to the one obtained with spectral element methods, which requires the use of quadrilateral (in 2D) of hexaedral (in 3D) cells and is therefore less convenient to deal with a topography. Both solution present the same order of accuracy, which encourage us to use DGFEM to deal with a topography.
Near the topography (or in very thin layers), the cells of the mesh have to be very small compared to the cells far from the topography. In such configuration, it is not useful to user high order elements in the whole domain and we propose to use second order elements in the fine part of the mesh and high order elements in the coarse part. The numerical experiments we have performed show actually that there is no difference between the solutions obtained with meshes composed only of high-order cells and with meshes composed of second order cells near the topography. Moreover the use of second order cells reduces dramatically the computational burden. All these experiments have been performed with a second order time discretization and the next step is to adapt the order of the time discretization to the one of the space discrtization, this will be done during the next year.