High order edge sensors with ℓ 1 regularization for enhanced discontinuous Galerkin methods

Jan Glaubitz, Anne Gelb

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

This paper investigates the use of ℓ 1 regularization for solving hyperbolic conservation laws based on high order discontinuous Galerkin (DG) approximations. We first use the polynomial annihilation method to construct a high order edge sensor which enables us to flag “troubled” elements. The DG approximation is enhanced in these troubled regions by activating ℓ 1 regularization to promote sparsity in the corresponding jump function of the numerical solution. The resulting ℓ 1 optimization problem is efficiently implemented using the alternating direction method of multipliers. By enacting ℓ 1 regularization only in troubled cells, our method remains accurate and efficient, as no additional regularization or expensive iterative procedures are needed in smooth regions. We present results for the inviscid Burgers' equation as well as a nonlinear system of conservation laws using a nodal collocation-type DG method as a solver.

Original languageEnglish (US)
Pages (from-to)A1304-A1330
JournalSIAM Journal on Scientific Computing
Volume41
Issue number2
DOIs
StatePublished - Jan 1 2019
Externally publishedYes

Fingerprint

Discontinuous Galerkin Method
Galerkin methods
Conservation
Regularization
Higher Order
Sensor
Discontinuous Galerkin
Galerkin Approximation
Sensors
Nonlinear systems
Polynomials
Method of multipliers
Alternating Direction Method
Systems of Conservation Laws
Hyperbolic Conservation Laws
Iterative Procedure
Burgers Equation
Annihilation
Collocation
Sparsity

Keywords

  • Discontinuity sensor
  • Discontinuous Galerkin
  • Hyperbolic conservation laws
  • Polynomial annihilation
  • Shock capturing
  • ℓ regularization

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Cite this

High order edge sensors with ℓ 1 regularization for enhanced discontinuous Galerkin methods . / Glaubitz, Jan; Gelb, Anne.

In: SIAM Journal on Scientific Computing, Vol. 41, No. 2, 01.01.2019, p. A1304-A1330.

Research output: Contribution to journalArticle

@article{2d86ce17e7274cf4bec020b8d10f70a7,
title = "High order edge sensors with ℓ 1 regularization for enhanced discontinuous Galerkin methods",
abstract = "This paper investigates the use of ℓ 1 regularization for solving hyperbolic conservation laws based on high order discontinuous Galerkin (DG) approximations. We first use the polynomial annihilation method to construct a high order edge sensor which enables us to flag “troubled” elements. The DG approximation is enhanced in these troubled regions by activating ℓ 1 regularization to promote sparsity in the corresponding jump function of the numerical solution. The resulting ℓ 1 optimization problem is efficiently implemented using the alternating direction method of multipliers. By enacting ℓ 1 regularization only in troubled cells, our method remains accurate and efficient, as no additional regularization or expensive iterative procedures are needed in smooth regions. We present results for the inviscid Burgers' equation as well as a nonlinear system of conservation laws using a nodal collocation-type DG method as a solver.",
keywords = "Discontinuity sensor, Discontinuous Galerkin, Hyperbolic conservation laws, Polynomial annihilation, Shock capturing, ℓ regularization",
author = "Jan Glaubitz and Anne Gelb",
year = "2019",
month = "1",
day = "1",
doi = "10.1137/18M1195280",
language = "English (US)",
volume = "41",
pages = "A1304--A1330",
journal = "SIAM Journal of Scientific Computing",
issn = "1064-8275",
publisher = "Society for Industrial and Applied Mathematics Publications",
number = "2",

}

TY - JOUR

T1 - High order edge sensors with ℓ 1 regularization for enhanced discontinuous Galerkin methods

AU - Glaubitz, Jan

AU - Gelb, Anne

PY - 2019/1/1

Y1 - 2019/1/1

N2 - This paper investigates the use of ℓ 1 regularization for solving hyperbolic conservation laws based on high order discontinuous Galerkin (DG) approximations. We first use the polynomial annihilation method to construct a high order edge sensor which enables us to flag “troubled” elements. The DG approximation is enhanced in these troubled regions by activating ℓ 1 regularization to promote sparsity in the corresponding jump function of the numerical solution. The resulting ℓ 1 optimization problem is efficiently implemented using the alternating direction method of multipliers. By enacting ℓ 1 regularization only in troubled cells, our method remains accurate and efficient, as no additional regularization or expensive iterative procedures are needed in smooth regions. We present results for the inviscid Burgers' equation as well as a nonlinear system of conservation laws using a nodal collocation-type DG method as a solver.

AB - This paper investigates the use of ℓ 1 regularization for solving hyperbolic conservation laws based on high order discontinuous Galerkin (DG) approximations. We first use the polynomial annihilation method to construct a high order edge sensor which enables us to flag “troubled” elements. The DG approximation is enhanced in these troubled regions by activating ℓ 1 regularization to promote sparsity in the corresponding jump function of the numerical solution. The resulting ℓ 1 optimization problem is efficiently implemented using the alternating direction method of multipliers. By enacting ℓ 1 regularization only in troubled cells, our method remains accurate and efficient, as no additional regularization or expensive iterative procedures are needed in smooth regions. We present results for the inviscid Burgers' equation as well as a nonlinear system of conservation laws using a nodal collocation-type DG method as a solver.

KW - Discontinuity sensor

KW - Discontinuous Galerkin

KW - Hyperbolic conservation laws

KW - Polynomial annihilation

KW - Shock capturing

KW - ℓ regularization

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

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

U2 - 10.1137/18M1195280

DO - 10.1137/18M1195280

M3 - Article

AN - SCOPUS:85065583896

VL - 41

SP - A1304-A1330

JO - SIAM Journal of Scientific Computing

JF - SIAM Journal of Scientific Computing

SN - 1064-8275

IS - 2

ER -