One-dimensional shock-capturing for high-order discontinuous Galerkin methods

SUMMARY Discontinuous Galerkin methods have emerged in recent years as an alternative for nonlinear conservation equations. In particular, their inherent structure (a numerical flux based on a suitable approximate Riemann solver introduces some stabilization) suggests that they are specially adapted...

Full description

Saved in:
Bibliographic Details
Published in:International journal for numerical methods in fluids Vol. 71; no. 6; pp. 737 - 755
Main Authors: Casoni, E., Peraire, J., Huerta, A.
Format: Journal Article Publication
Language:English
Published: Bognor Regis Blackwell Publishing Ltd 28.02.2013
Wiley Subscription Services, Inc
Subjects:
65 Numerical analysis
65E05 Numerical methods in complex analysis (potential theory, etc.)
Algorithms
Anàlisi numèrica
Approximation
artificial viscosity
Classificació AMS
Constraining
Diffusion
discontinuous Galerkin
Elements finits, Mètode dels
Finite volume method
Galerkin methods
high-order approximation
Matemàtiques i estadística
Mathematical analysis
Mathematical models
Mètodes numèrics
Nonlinearity
Numerical methods and algorithms
Reproduction
shock capturing
Studies
Àrees temàtiques de la UPC
ISSN:0271-2091, 1097-0363
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:SUMMARY Discontinuous Galerkin methods have emerged in recent years as an alternative for nonlinear conservation equations. In particular, their inherent structure (a numerical flux based on a suitable approximate Riemann solver introduces some stabilization) suggests that they are specially adapted to capture shocks. However, numerical fluxes are not sufficient to stabilize the solution in the presence of shocks. Thus, slope limiter methods, which are extensions of finite volume methods, have been proposed. These techniques require, in practice, mesh adaption to localize the shock structure. This is is more obvious for large elements typical of high‐order approximations. Here, a new approach based on the introduction of artificial diffusion into the original equations is presented. The order is not systematically decreased to one in the presence of the shock, large high‐order elements can be used, and several linear and nonlinear tests demonstrate the efficiency of the proposed methodology. Copyright © 2012 John Wiley & Sons, Ltd. Figure 1 compares the artificial diffusion technique proposed here with high‐order limiters (described as moments’ in the figures), with the same number of degrees of freedom. The artificial diffusion technique outperforms the high‐order limiting scheme, demonstrating that it seems more beneficial to increase the degree of approximation than to refine the mesh in order to obtain accurate solutions.
Bibliography:Generalitat de Catalunya AGAUR - No. 2009SGR875
ark:/67375/WNG-2KXHRLSV-9
istex:144E5772A03EF09FB90C6A677EBB44FB84551A09
ArticleID:FLD3682
ObjectType-Article-1
SourceType-Scholarly Journals-1
content type line 14
ObjectType-Feature-2
content type line 23
ISSN:0271-2091
1097-0363
DOI:10.1002/fld.3682