Performance comparison of HDG and classical DG method for the simulation of seismic wave propagation in harmonic domain

International audienceIn the most widely used methods for Seismic Imaging, we have to solve 2N wave equations at each iteration of the selected process, where N is the number of sources and is usually large (about 1000). The efficiency of the inverse solver is thus directly related to the efficiency of the numerical method used to solve the wave equation. Seismic imaging can be performed in the time domain or in the frequency domain regime. We focus here on the second setting. The drawback of time domain is that it requires to store the solution at each time step of the forward simulation. The difficulties related to frequency domain inversion lie in the solution of huge linear systems, which cannot be achieved today when considering realistic 3D elastic media, even with the progress of high-performance computing. In this context, the goal is to develop new forward solvers that reduce the number of degrees of freedom without hampering the accuracy of the numerical solution.We consider here discontinuous Galerkin (DG) methods which are more convenient than finite difference methods to handle the topography of the subsurface. Moreover, they are more adapted than continuous Galerkin (CG) methods to deal with hp-adaptivity. This last characteristics is crucial to adapt the mesh to the different regions of the subsurface which is generally highly heterogeneous. Nevertheless, the main drawback of classical DG methods is that they are expensive because they require a large number of degrees of freedom as compared to CG methods on a given mesh. In this work we consider a new class of DG method, the hybridizable DG (HDG) method. Instead of solving a linear system involving the degrees of freedom of all volumic cells of the mesh, the principle of HDG consists in introducing a Lagrange multiplier representing the trace of the numerical solution on each face of the mesh. Hence, it reduces the number of unknowns of the global linear systems and the volumic solution is recovered thanks to a local computation on each element. We compare the performances of the HDG method with those of classical nodal DG methods and we present our first results using a HDG method for the first-order form of the elastic wave propagation equations for 2D realistic test-cases

A non-overlapping Schwarz algorithm for the HDG method

In this paper, we present two non-overlapping Schwarz algorithms for the hybridizable discontinuous Galerkin (HDG) method. The first algorithm is based on the Neumann-Neumann method. The second one is an iterative algorithm uses both trace and flux interface unknowns on interfaces between subdomains. Numerical results are provided to verify the validity of our algorithms

Performances analysis of a Hybridizable discontinuous Galerkin solver for the 3D Helmholtz equations in geophysical context

International audienc

Unified analysis of HDG methods using scalar and vector hybrid variables

In this paper, hybridizable discontinuous Galerkin (HDG) methods using scalar and vector hybrid variables for steady-state diffusion problems are considered. We propose a unified framework to analyze the methods, where both the hybrid variables are treated as double-valued functions. If either of them is single valued, the well-posedness is ensured under some assumptions on approximation spaces. Moreover, we prove that all methods are superconvergent, based on the so-called $M$-decomposition theory. Numerical results are presented to validate our theoretical results.Comment: 16 page