Accéder directement au contenu Accéder directement à la navigation
Pré-publication, Document de travail

A high-order discontinuous Galerkin Method using a mixture of Gauss-Legendre and Gauss-Lobatto quadratures for improved efficiency

Stéphanie Chaillat 1 Régis Cottereau 2 Ruben Sevilla 3 
1 POEMS - Propagation des Ondes : Étude Mathématique et Simulation
Inria Saclay - Ile de France, UMA - Unité de Mathématiques Appliquées, CNRS - Centre National de la Recherche Scientifique : UMR7231
Abstract : In discontinuous Galerkin spectral element methods (DGSEM), the two most common approaches to numerically integrate the terms of the weak form are either using Gauss-Legendre or Gauss-Lobatto quadratures. The former yields more accurate results but at a higher computational cost, so that a priori it is not clear whether one approach is more efficient that the other. In this paper, it is shown (theoretically for a particular case and numerically for the general case) that using Gauss-Lobatto quadrature for the convection matrix actually introduces a negligible error. In contrast, using Gauss-Lobatto quadratures for the evaluation of the jump term in the element faces introduces a sizeable error. This leads to the proposal of a new DG approach, where the convection matrix is evaluated using Gauss-Lobatto quadratures, whereas the face mass matrices are integrated using Gauss-Legendre quadratures. For elements with constant Jacobian and constant coefficients, a formal proof shows that no numerical integration error is actually introduced in the evaluation of the residual, even though both the mass and the convection matrices are not computed exactly with Gauss-Lobatto quadratures. For elements with non-constant Jacobian and/or non-constant coefficients, the impact of numerical integration error on the overall error is evaluated through a series of numerical tests, showing that this is also negligible. In addition, the computational cost associated to the matrix-vector products required to evaluate the residual is evaluated precisely for the different cases considered. The proposed approach is particularly attractive in the most general case, since the use of Gauss-Lobatto quadratures significantly speeds-up the evaluation of the residual.
Type de document :
Pré-publication, Document de travail
Liste complète des métadonnées

https://hal.archives-ouvertes.fr/hal-03695573
Contributeur : Régis Cottereau Connectez-vous pour contacter le contributeur
Soumis le : mercredi 15 juin 2022 - 06:43:19
Dernière modification le : samedi 25 juin 2022 - 03:26:47

Fichier

DG_NodalSEM.pdf
Fichiers produits par l'(les) auteur(s)

Licence


Distributed under a Creative Commons Paternité 4.0 International License

Identifiants

  • HAL Id : hal-03695573, version 1

Citation

Stéphanie Chaillat, Régis Cottereau, Ruben Sevilla. A high-order discontinuous Galerkin Method using a mixture of Gauss-Legendre and Gauss-Lobatto quadratures for improved efficiency. 2022. ⟨hal-03695573⟩

Partager

Métriques

Consultations de la notice

50

Téléchargements de fichiers

17