Ce travail de thèse qui a pour objet la génération et l'étude des couches absorbantes dans les problèmes impliquant la dynamique des sols, est divisé en trois parties essentielles. La première consiste à proposer une méthode de dimensionnement des couches absorbantes par l'amortissement de Rayleigh afin de simuler des problèmes de propagation d'ondes dans les milieux infinis. Cette méthode repose sur une analyse mathématique du problème de propagation d'ondes dans un milieu caractérisé par la matrice de Rayleigh, qui nous permet, d'une part, d'établir des conditions de minimisation des réflexions parasites aux interfaces, et d'autre part, de proposer une simple relation de dimensionnement du domaine absorbant basée sur la notion de décrément logarithmique. On se propose dans la deuxième partie d'appliquer une stratégie de couplage des schémas temporels pour des problèmes de propagation d'ondes dans les milieux infinis 1D et 2D. L'approche proposée est d'intégrer le domaine d'étude par un schéma explicite et le domaine absorbant par un schéma implicite, et d'évaluer le potentiel de cette méthode en faisant varier les rapports de pas de temps entre les sous domaines. Une attention particulière est accordée au cas 1D pour lequel l'effet de la finesse du maillage définie par le nombre d'éléments finis par longueur d'onde est également analysé. Par ailleurs, l'évolution du temps de calcul en fonction du rapport entre les pas de temps est étudiée afin d'estimer les gains réalisés par rapport à un calcul de référence où le problème global est intégré uniquement avec un schéma explicite. La dernière partie est dédiée à l'étude des couches amortissantes de type PML ("Perfectly Matched Layer") dans le cadre des couplages hybrides multi-pas de temps. Cette partie est introduite par une étude de stabilité des schémas temporels dans le cas d'une PML en 1D. La couche absorbante PML est intégrée selon un schéma implicite en adoptant des pas de temps plus importants que le domaine d'intérêt intégré selon un schéma explicite. Bien que cette méthodologie de couplage s'avère très efficace pour la reproduction des milieux infinis, les études paramétriques montrent une sensibilité à la taille du pas de temps plus forte que celle exhibée par les couches amortissantes de Rayleigh. / This thesis which deals with the study of absorbing layers for soil dynamics problems, is divided into three essential parts. The first part aims to propose a design method of absorbing layers by the Rayleigh damping to simulate wave propagation problems in infinite media. This method is based on a mathematical analysis of the wave propagation problem in a media characterized by a Rayleigh damping matrix, which allows us, firstly, to establish conditions for minimizing spurious waves at the interfaces, and another hand, to provide a simple design relationship for the absorbing domain based on the notion of the logarithmic decrement. The second part aims to apply the multi-time step strategy for wave propagation problems in 1D and 2D infinite media. The proposed approach is to integrate the physical domain by an explicit scheme and the absorbing domain by an implicit scheme and to evaluate the potential of this method by varying the time step ratio between subdomains. Special attention is given to the 1D case for which the effect of the mesh fineness, defined by the number of finite elements per wavelength, is also analyzed. Furthermore, the evolution of computing time depending on the time ratio is studied in order to estimate the gains made with respect to a reference computation achieved by a full explicit integration. The last part is dedicated to the study of the Perfectly Matched Layer (PML) as part of hybrid couplings multi-time step. This section is introduced by a stability study of temporal scheme for 1D cases. The absorbing layer PML is integrated by an implicit scheme with a time step larger than that of the domain of interest. Although this coupling methodology is very effective for the reproduction of infinite media, parametric studies show a sensitivity to the time ratio greater than that exhibited by the Rayleigh damping layers.

Cette thèse étudie le raffinement local de maillage à la fois en espace et en temps pour l’équation de l’elastodynamique du second ordre pour le calcul haute performance. L’objectif est de mettre en place des méthodes numériques pour traiter des hétérogénéités de petite taille ayant un impact important sur la propagation des ondes. Nous utilisons une approche par éléments finis de Galerkin discontinus avec pénalisation pour leur flexibilité et facilité de parallélisation. La formulation éléments finis que nous proposons a pour particularité d’être élasto-acoustique, pour pouvoir prendre en compte des hétérogénéités acoustiques de petite taille. Par ailleurs, nous proposons un terme de pénalisation optimisé qui est mieux adapté à l’équation de l’élastodynamique, conduisant en particulier à une meilleure condition CFL. Nous avons aussi amélioré une formulation PML du second ordre pour laquelle nous avons proposé une nouvelle discrétisation temporelle qui rend la formulation plus stable. En tirant parti de la p-adaptivité et des maillages non-conformes des méthodes de Galerkin discontinues combiné à une méthode de pas de temps local, nous avons grandement réduit le coût du raffinement local. Ces méthodes ont été implémentées en C++, en utilisant des techniques de template metaprogramming, au sein d’un code parallèle à mémoire distribuée (MPI) et partagée (OpenMP). Enfin, nous montrons le potentiel de notre approche sur des cas tests de validation et sur des cas plus réalistes avec des milieux présentant des hydrofractures. / This thesis studies local mesh refinement both in time and space for the second order elastodynamic equation in a high performance computing context. The objective is to develop numerical methods to treat small heterogeneities that have global impact on wave propagation. We use an internal penalty discontinuous Galerkin finite element approach for its flexibity and parallelization capabilities. The elasto-acoustic finite element formulation we discuss is elasto-acoustic in order to handle local acoustic heterogeneities. We also propose an optimized penalty term more suited to the elastodynamic equation that results in better CFL condition. We improve a second order PML formulation with an original time discretization that results in a more stable formulation. Using the p-adaptivity and nonconforming mesh capabilities of discontinuous Galerkin methods combined with a local time stepping method, we greatly reduce the high computational cost of local refinements. These methods have been implemented in C++, using template metaprogramming, in a distributed memory (MPI) and shared memory (OpenMP) parallel code. Finally, we show the potential of our methods on validation test cases and on more realistic test cases with medium including hydrofractures.