Implementation of Hybridizable Discontinuous Galerkin method for time-harmonic anisotropic poroelasticity in two dimensions.

Abstract

International audienceWe apply a Hybridizable Discontinuous Galerkin (HDG) method to numerically solve two-dimensional anisotropic poroelastic wave equations in the frequency domain given by Biot theory. The motivation for choosing HDG method comes from the complexity of the considered equations and the high number of unknowns. The HDG method possesses all the advantages of Discontinuous Galerkin method (hpadaptivity, accuracy, ability to model complex tectonics,...) without a drastic increase in the number of degrees of freedom. After a description of its implementation, we illustrate the accuracy of the proposed method by comparisons with analytical solutions. We then perform a sensitivity analysis of the method as a function of stabilization parameters and frequency, and show in particular that there exists optimal values of these parameters in order to obtain a very good level of accuracy. We also show the ability of the method to reproduce the different types of poroelastic waves including the slow Biot wave