Multiresolution-based adaptive central high resolution schemes for modeling of nonlinear propagating fronts

Hassan Yousefi, Timon Rabczuk

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

A class of wavelet-based methods, which combine adaptive grids together with second-order central and central-upwind high-resolution schemes, is introduced for accurate solution of first order hyperbolic systems of conservation laws and related equations. It is shown that the proposed class recovers stability which otherwise may be lost when the underlying central schemes are utilized over non-uniform grids. For simulations on irregular grids, both full-discretized and semi-discretized forms are derived; in particular, the effect of using certain shifts of the center of the computational allow to use standard slope limiters on non-uniform cells. Thereafter, the questions of numerical entropy production, local truncation errors, and the Total Variation Diminishing (TVD) criterion for scalar equations, are investigated. Central schemes are sensitives for irregular grids. To cure this drawback in the present context of multiresolution, the adapted grid is locally modified around high-gradient zones by adding new neighboring points. For an adapted point of resolution scale j, the new points are locally added in both resolution j and all successive coarser resolutions. It is shown that this simple grid modification improves stability of numerical solutions; this, however, comes at the expense that cell-centered central and central-upwind high resolution schemes, may not satisfy the TVD property. Finally, we present numerical simulations for both scalar and systems of non-linear conservation laws which confirm the simplicity and efficiency of the proposed method. These simulations demonstrate the high accuracy, the entropy production of wavelet-based algorithms which can effectively detect high gradient zones of shock waves, rarefaction regions, and contact discontinuities.

Original languageEnglish
Pages (from-to)172-195
Number of pages24
JournalEngineering Analysis with Boundary Elements
Volume103
DOIs
Publication statusPublished - 2019 Jun 1

Fingerprint

High-resolution Schemes
Central Schemes
Multiresolution
Irregular Grids
Conservation
Upwind Scheme
Entropy
Entropy Production
Diminishing
Total Variation
Limiters
Wavelets
Modeling
Scalar
Shock waves
Gradient
Grid
Contact Discontinuity
Non-uniform Grid
Hyperbolic Systems of Conservation Laws

Keywords

  • Adapted grid
  • Central high-resolution schemes
  • Discontinuous fronts
  • Interpolating wavelet
  • Nonlinear waves
  • Second order accuracy

ASJC Scopus subject areas

  • Analysis
  • Engineering(all)
  • Computational Mathematics
  • Applied Mathematics

Cite this

Multiresolution-based adaptive central high resolution schemes for modeling of nonlinear propagating fronts. / Yousefi, Hassan; Rabczuk, Timon.

In: Engineering Analysis with Boundary Elements, Vol. 103, 01.06.2019, p. 172-195.

Research output: Contribution to journalArticle

@article{1b2ab84d94d04bca9dc2fbc5f224ef6c,
title = "Multiresolution-based adaptive central high resolution schemes for modeling of nonlinear propagating fronts",
abstract = "A class of wavelet-based methods, which combine adaptive grids together with second-order central and central-upwind high-resolution schemes, is introduced for accurate solution of first order hyperbolic systems of conservation laws and related equations. It is shown that the proposed class recovers stability which otherwise may be lost when the underlying central schemes are utilized over non-uniform grids. For simulations on irregular grids, both full-discretized and semi-discretized forms are derived; in particular, the effect of using certain shifts of the center of the computational allow to use standard slope limiters on non-uniform cells. Thereafter, the questions of numerical entropy production, local truncation errors, and the Total Variation Diminishing (TVD) criterion for scalar equations, are investigated. Central schemes are sensitives for irregular grids. To cure this drawback in the present context of multiresolution, the adapted grid is locally modified around high-gradient zones by adding new neighboring points. For an adapted point of resolution scale j, the new points are locally added in both resolution j and all successive coarser resolutions. It is shown that this simple grid modification improves stability of numerical solutions; this, however, comes at the expense that cell-centered central and central-upwind high resolution schemes, may not satisfy the TVD property. Finally, we present numerical simulations for both scalar and systems of non-linear conservation laws which confirm the simplicity and efficiency of the proposed method. These simulations demonstrate the high accuracy, the entropy production of wavelet-based algorithms which can effectively detect high gradient zones of shock waves, rarefaction regions, and contact discontinuities.",
keywords = "Adapted grid, Central high-resolution schemes, Discontinuous fronts, Interpolating wavelet, Nonlinear waves, Second order accuracy",
author = "Hassan Yousefi and Timon Rabczuk",
year = "2019",
month = "6",
day = "1",
doi = "10.1016/j.enganabound.2019.03.002",
language = "English",
volume = "103",
pages = "172--195",
journal = "Engineering Analysis with Boundary Elements",
issn = "0955-7997",
publisher = "Elsevier Limited",

}

TY - JOUR

T1 - Multiresolution-based adaptive central high resolution schemes for modeling of nonlinear propagating fronts

AU - Yousefi, Hassan

AU - Rabczuk, Timon

PY - 2019/6/1

Y1 - 2019/6/1

N2 - A class of wavelet-based methods, which combine adaptive grids together with second-order central and central-upwind high-resolution schemes, is introduced for accurate solution of first order hyperbolic systems of conservation laws and related equations. It is shown that the proposed class recovers stability which otherwise may be lost when the underlying central schemes are utilized over non-uniform grids. For simulations on irregular grids, both full-discretized and semi-discretized forms are derived; in particular, the effect of using certain shifts of the center of the computational allow to use standard slope limiters on non-uniform cells. Thereafter, the questions of numerical entropy production, local truncation errors, and the Total Variation Diminishing (TVD) criterion for scalar equations, are investigated. Central schemes are sensitives for irregular grids. To cure this drawback in the present context of multiresolution, the adapted grid is locally modified around high-gradient zones by adding new neighboring points. For an adapted point of resolution scale j, the new points are locally added in both resolution j and all successive coarser resolutions. It is shown that this simple grid modification improves stability of numerical solutions; this, however, comes at the expense that cell-centered central and central-upwind high resolution schemes, may not satisfy the TVD property. Finally, we present numerical simulations for both scalar and systems of non-linear conservation laws which confirm the simplicity and efficiency of the proposed method. These simulations demonstrate the high accuracy, the entropy production of wavelet-based algorithms which can effectively detect high gradient zones of shock waves, rarefaction regions, and contact discontinuities.

AB - A class of wavelet-based methods, which combine adaptive grids together with second-order central and central-upwind high-resolution schemes, is introduced for accurate solution of first order hyperbolic systems of conservation laws and related equations. It is shown that the proposed class recovers stability which otherwise may be lost when the underlying central schemes are utilized over non-uniform grids. For simulations on irregular grids, both full-discretized and semi-discretized forms are derived; in particular, the effect of using certain shifts of the center of the computational allow to use standard slope limiters on non-uniform cells. Thereafter, the questions of numerical entropy production, local truncation errors, and the Total Variation Diminishing (TVD) criterion for scalar equations, are investigated. Central schemes are sensitives for irregular grids. To cure this drawback in the present context of multiresolution, the adapted grid is locally modified around high-gradient zones by adding new neighboring points. For an adapted point of resolution scale j, the new points are locally added in both resolution j and all successive coarser resolutions. It is shown that this simple grid modification improves stability of numerical solutions; this, however, comes at the expense that cell-centered central and central-upwind high resolution schemes, may not satisfy the TVD property. Finally, we present numerical simulations for both scalar and systems of non-linear conservation laws which confirm the simplicity and efficiency of the proposed method. These simulations demonstrate the high accuracy, the entropy production of wavelet-based algorithms which can effectively detect high gradient zones of shock waves, rarefaction regions, and contact discontinuities.

KW - Adapted grid

KW - Central high-resolution schemes

KW - Discontinuous fronts

KW - Interpolating wavelet

KW - Nonlinear waves

KW - Second order accuracy

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

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

U2 - 10.1016/j.enganabound.2019.03.002

DO - 10.1016/j.enganabound.2019.03.002

M3 - Article

VL - 103

SP - 172

EP - 195

JO - Engineering Analysis with Boundary Elements

JF - Engineering Analysis with Boundary Elements

SN - 0955-7997

ER -