Directly Simulation of Second Order Hyperbolic Systems in Second Order Form via the Regularization Concept

Hassan Yousefi, Seyed Shahram Ghorashi, Timon Rabczuk

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

We present an efficient and robust method for stress wave propagation problems (second order hyperbolic systems) having discontinuities directly in their second order form. Due to the numerical dispersion around discontinuities and lack of the inherent dissipation in hyperbolic systems, proper simulation of such problems are challenging. The proposed idea is to denoise spurious oscillations by a post-processing stage from solutions obtained from higher-order grid-based methods (e.g., high-order collocation or finite-difference schemes). The denoising is done so that the solutions remain higher-order (here, second order) around discontinuities and are still free from spurious oscillations. For this purpose, improved Tikhonov regularization approach is advised. This means to let data themselves select proper denoised solutions (since there is no pre-Assumptions about regularized results). The improved approach can directly be done on uniform or non-uniform sampled data in a way that the regularized results maintenance continuous derivatives up to some desired order. It is shown how to improve the smoothing method so that it remains conservative and has local estimating feature. To confirm effectiveness of the proposed approach, finally, some one and two dimensional examples will be provided. It will be shown how both the numerical (artificial) dispersion and dissipation can be controlled around discontinuous solutions and stochastic-like results.

Original languageEnglish
Pages (from-to)86-135
Number of pages50
JournalCommunications in Computational Physics
Volume20
Issue number1
DOIs
Publication statusPublished - 2016 Jun 22
Externally publishedYes

Fingerprint

hyperbolic systems
discontinuity
dissipation
simulation
stress waves
oscillations
collocation
smoothing
maintenance
wave propagation
estimating
grids

Keywords

  • Numerical (artificial) dispersion
  • Second order hyperbolic systems
  • Tikhonov regularization

ASJC Scopus subject areas

  • Physics and Astronomy (miscellaneous)

Cite this

Directly Simulation of Second Order Hyperbolic Systems in Second Order Form via the Regularization Concept. / Yousefi, Hassan; Ghorashi, Seyed Shahram; Rabczuk, Timon.

In: Communications in Computational Physics, Vol. 20, No. 1, 22.06.2016, p. 86-135.

Research output: Contribution to journalArticle

Yousefi, Hassan ; Ghorashi, Seyed Shahram ; Rabczuk, Timon. / Directly Simulation of Second Order Hyperbolic Systems in Second Order Form via the Regularization Concept. In: Communications in Computational Physics. 2016 ; Vol. 20, No. 1. pp. 86-135.
@article{96b9e97cd53f40be8c3a9ed805db174b,
title = "Directly Simulation of Second Order Hyperbolic Systems in Second Order Form via the Regularization Concept",
abstract = "We present an efficient and robust method for stress wave propagation problems (second order hyperbolic systems) having discontinuities directly in their second order form. Due to the numerical dispersion around discontinuities and lack of the inherent dissipation in hyperbolic systems, proper simulation of such problems are challenging. The proposed idea is to denoise spurious oscillations by a post-processing stage from solutions obtained from higher-order grid-based methods (e.g., high-order collocation or finite-difference schemes). The denoising is done so that the solutions remain higher-order (here, second order) around discontinuities and are still free from spurious oscillations. For this purpose, improved Tikhonov regularization approach is advised. This means to let data themselves select proper denoised solutions (since there is no pre-Assumptions about regularized results). The improved approach can directly be done on uniform or non-uniform sampled data in a way that the regularized results maintenance continuous derivatives up to some desired order. It is shown how to improve the smoothing method so that it remains conservative and has local estimating feature. To confirm effectiveness of the proposed approach, finally, some one and two dimensional examples will be provided. It will be shown how both the numerical (artificial) dispersion and dissipation can be controlled around discontinuous solutions and stochastic-like results.",
keywords = "Numerical (artificial) dispersion, Second order hyperbolic systems, Tikhonov regularization",
author = "Hassan Yousefi and Ghorashi, {Seyed Shahram} and Timon Rabczuk",
year = "2016",
month = "6",
day = "22",
doi = "10.4208/cicp.101214.011015a",
language = "English",
volume = "20",
pages = "86--135",
journal = "Communications in Computational Physics",
issn = "1815-2406",
publisher = "Global Science Press",
number = "1",

}

TY - JOUR

T1 - Directly Simulation of Second Order Hyperbolic Systems in Second Order Form via the Regularization Concept

AU - Yousefi, Hassan

AU - Ghorashi, Seyed Shahram

AU - Rabczuk, Timon

PY - 2016/6/22

Y1 - 2016/6/22

N2 - We present an efficient and robust method for stress wave propagation problems (second order hyperbolic systems) having discontinuities directly in their second order form. Due to the numerical dispersion around discontinuities and lack of the inherent dissipation in hyperbolic systems, proper simulation of such problems are challenging. The proposed idea is to denoise spurious oscillations by a post-processing stage from solutions obtained from higher-order grid-based methods (e.g., high-order collocation or finite-difference schemes). The denoising is done so that the solutions remain higher-order (here, second order) around discontinuities and are still free from spurious oscillations. For this purpose, improved Tikhonov regularization approach is advised. This means to let data themselves select proper denoised solutions (since there is no pre-Assumptions about regularized results). The improved approach can directly be done on uniform or non-uniform sampled data in a way that the regularized results maintenance continuous derivatives up to some desired order. It is shown how to improve the smoothing method so that it remains conservative and has local estimating feature. To confirm effectiveness of the proposed approach, finally, some one and two dimensional examples will be provided. It will be shown how both the numerical (artificial) dispersion and dissipation can be controlled around discontinuous solutions and stochastic-like results.

AB - We present an efficient and robust method for stress wave propagation problems (second order hyperbolic systems) having discontinuities directly in their second order form. Due to the numerical dispersion around discontinuities and lack of the inherent dissipation in hyperbolic systems, proper simulation of such problems are challenging. The proposed idea is to denoise spurious oscillations by a post-processing stage from solutions obtained from higher-order grid-based methods (e.g., high-order collocation or finite-difference schemes). The denoising is done so that the solutions remain higher-order (here, second order) around discontinuities and are still free from spurious oscillations. For this purpose, improved Tikhonov regularization approach is advised. This means to let data themselves select proper denoised solutions (since there is no pre-Assumptions about regularized results). The improved approach can directly be done on uniform or non-uniform sampled data in a way that the regularized results maintenance continuous derivatives up to some desired order. It is shown how to improve the smoothing method so that it remains conservative and has local estimating feature. To confirm effectiveness of the proposed approach, finally, some one and two dimensional examples will be provided. It will be shown how both the numerical (artificial) dispersion and dissipation can be controlled around discontinuous solutions and stochastic-like results.

KW - Numerical (artificial) dispersion

KW - Second order hyperbolic systems

KW - Tikhonov regularization

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

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

U2 - 10.4208/cicp.101214.011015a

DO - 10.4208/cicp.101214.011015a

M3 - Article

AN - SCOPUS:84975867888

VL - 20

SP - 86

EP - 135

JO - Communications in Computational Physics

JF - Communications in Computational Physics

SN - 1815-2406

IS - 1

ER -