A high-order discontinuous Galerkin Method using a mixture of Gauss-Legendre and Gauss-Lobatto quadratures for improved efficiency - Ecole Centrale de Marseille Accéder directement au contenu
Pré-Publication, Document De Travail Année : 2022

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

Résumé

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.
Fichier principal
Vignette du fichier
DG_NodalSEM.pdf (3.99 Mo) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)

Dates et versions

hal-03695573 , version 1 (15-06-2022)

Licence

Paternité

Identifiants

  • HAL Id : hal-03695573 , version 1

Citer

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⟩
128 Consultations
104 Téléchargements

Partager

Gmail Facebook X LinkedIn More