A local projection stabilization method with shock capturing and diagonal mass matrix for solving non-stationary transport dominated problems

Friedhelm Schieweck, Piotr Skrzypacz

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

We consider a time-dependent convection diffusion equation in the transport dominated case. As a stabilization method in space we propose a new variant of Local Projection Stabilization (LPS) which uses special enriched bubble functions such that L 2-orthogonal local basis functions can be constructed. L 2-orthogonal basis functions lead to a diagonal mass matrix which is advantageous for time discretization. We use the discontinuous Galerkin method of polynomial order one for the discretization in time which is superconvergent of order three at the endpoints of the time intervals. In order to avoid the remaining oscillations in the LPS-solution we add for each time step in the space discretization an extra shock capturing term which acts only locally on those mesh cells where an error-indicator is relatively large. The novelty in the shock capturing term is that the scaling factor in front of the additive diffusion term is computed from a low order post-processing error. As a result we obtain both, an oscillation-free discrete solution and the information about the local regions where this solution is still inaccurate due to some smearing. The latter information can be used to create in each time step an adaptively refined space mesh. Whereas the numerical experiments are restricted to one space dimension the proposed ideas work also in the multi-dimensional spatial case. The numerical tests show that the discrete solution with shock capturing is oscillation-free and of optimal accuracy in the regions outside of the shock.

Original languageEnglish
Pages (from-to)221-240
Number of pages20
JournalComputational Methods in Applied Mathematics
Volume12
Issue number2
DOIs
Publication statusPublished - Apr 2012
Externally publishedYes

Fingerprint

Shock Capturing
Stabilization
Projection
Oscillation
Basis Functions
Term
Discretization
Mesh
Bubble Function
Error Indicator
Orthogonal Functions
Orthogonal Basis
Scaling Factor
Discontinuous Galerkin Method
Convection-diffusion Equation
Time Discretization
Galerkin methods
Inaccurate
Post-processing
Shock

Keywords

  • Discontinuous Galerkin time discretization
  • Error indicator
  • Local projection stabilization
  • Post-processing
  • Shock capturing

ASJC Scopus subject areas

  • Numerical Analysis
  • Applied Mathematics
  • Computational Mathematics

Cite this

@article{322cbcf4f66c427589127b22906a42a3,
title = "A local projection stabilization method with shock capturing and diagonal mass matrix for solving non-stationary transport dominated problems",
abstract = "We consider a time-dependent convection diffusion equation in the transport dominated case. As a stabilization method in space we propose a new variant of Local Projection Stabilization (LPS) which uses special enriched bubble functions such that L 2-orthogonal local basis functions can be constructed. L 2-orthogonal basis functions lead to a diagonal mass matrix which is advantageous for time discretization. We use the discontinuous Galerkin method of polynomial order one for the discretization in time which is superconvergent of order three at the endpoints of the time intervals. In order to avoid the remaining oscillations in the LPS-solution we add for each time step in the space discretization an extra shock capturing term which acts only locally on those mesh cells where an error-indicator is relatively large. The novelty in the shock capturing term is that the scaling factor in front of the additive diffusion term is computed from a low order post-processing error. As a result we obtain both, an oscillation-free discrete solution and the information about the local regions where this solution is still inaccurate due to some smearing. The latter information can be used to create in each time step an adaptively refined space mesh. Whereas the numerical experiments are restricted to one space dimension the proposed ideas work also in the multi-dimensional spatial case. The numerical tests show that the discrete solution with shock capturing is oscillation-free and of optimal accuracy in the regions outside of the shock.",
keywords = "Discontinuous Galerkin time discretization, Error indicator, Local projection stabilization, Post-processing, Shock capturing",
author = "Friedhelm Schieweck and Piotr Skrzypacz",
year = "2012",
month = "4",
doi = "10.2478/cmam-2012-0019",
language = "English",
volume = "12",
pages = "221--240",
journal = "Computational Methods in Applied Mathematics",
issn = "1609-4840",
publisher = "Walter de Gruyter GmbH",
number = "2",

}

TY - JOUR

T1 - A local projection stabilization method with shock capturing and diagonal mass matrix for solving non-stationary transport dominated problems

AU - Schieweck, Friedhelm

AU - Skrzypacz, Piotr

PY - 2012/4

Y1 - 2012/4

N2 - We consider a time-dependent convection diffusion equation in the transport dominated case. As a stabilization method in space we propose a new variant of Local Projection Stabilization (LPS) which uses special enriched bubble functions such that L 2-orthogonal local basis functions can be constructed. L 2-orthogonal basis functions lead to a diagonal mass matrix which is advantageous for time discretization. We use the discontinuous Galerkin method of polynomial order one for the discretization in time which is superconvergent of order three at the endpoints of the time intervals. In order to avoid the remaining oscillations in the LPS-solution we add for each time step in the space discretization an extra shock capturing term which acts only locally on those mesh cells where an error-indicator is relatively large. The novelty in the shock capturing term is that the scaling factor in front of the additive diffusion term is computed from a low order post-processing error. As a result we obtain both, an oscillation-free discrete solution and the information about the local regions where this solution is still inaccurate due to some smearing. The latter information can be used to create in each time step an adaptively refined space mesh. Whereas the numerical experiments are restricted to one space dimension the proposed ideas work also in the multi-dimensional spatial case. The numerical tests show that the discrete solution with shock capturing is oscillation-free and of optimal accuracy in the regions outside of the shock.

AB - We consider a time-dependent convection diffusion equation in the transport dominated case. As a stabilization method in space we propose a new variant of Local Projection Stabilization (LPS) which uses special enriched bubble functions such that L 2-orthogonal local basis functions can be constructed. L 2-orthogonal basis functions lead to a diagonal mass matrix which is advantageous for time discretization. We use the discontinuous Galerkin method of polynomial order one for the discretization in time which is superconvergent of order three at the endpoints of the time intervals. In order to avoid the remaining oscillations in the LPS-solution we add for each time step in the space discretization an extra shock capturing term which acts only locally on those mesh cells where an error-indicator is relatively large. The novelty in the shock capturing term is that the scaling factor in front of the additive diffusion term is computed from a low order post-processing error. As a result we obtain both, an oscillation-free discrete solution and the information about the local regions where this solution is still inaccurate due to some smearing. The latter information can be used to create in each time step an adaptively refined space mesh. Whereas the numerical experiments are restricted to one space dimension the proposed ideas work also in the multi-dimensional spatial case. The numerical tests show that the discrete solution with shock capturing is oscillation-free and of optimal accuracy in the regions outside of the shock.

KW - Discontinuous Galerkin time discretization

KW - Error indicator

KW - Local projection stabilization

KW - Post-processing

KW - Shock capturing

UR - http://www.scopus.com/inward/record.url?scp=84868626271&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84868626271&partnerID=8YFLogxK

U2 - 10.2478/cmam-2012-0019

DO - 10.2478/cmam-2012-0019

M3 - Article

VL - 12

SP - 221

EP - 240

JO - Computational Methods in Applied Mathematics

JF - Computational Methods in Applied Mathematics

SN - 1609-4840

IS - 2

ER -