- Journal Home
- Volume 36 - 2024
- Volume 35 - 2024
- Volume 34 - 2023
- Volume 33 - 2023
- Volume 32 - 2022
- Volume 31 - 2022
- Volume 30 - 2021
- Volume 29 - 2021
- Volume 28 - 2020
- Volume 27 - 2020
- Volume 26 - 2019
- Volume 25 - 2019
- Volume 24 - 2018
- Volume 23 - 2018
- Volume 22 - 2017
- Volume 21 - 2017
- Volume 20 - 2016
- Volume 19 - 2016
- Volume 18 - 2015
- Volume 17 - 2015
- Volume 16 - 2014
- Volume 15 - 2014
- Volume 14 - 2013
- Volume 13 - 2013
- Volume 12 - 2012
- Volume 11 - 2012
- Volume 10 - 2011
- Volume 9 - 2011
- Volume 8 - 2010
- Volume 7 - 2010
- Volume 6 - 2009
- Volume 5 - 2009
- Volume 4 - 2008
- Volume 3 - 2008
- Volume 2 - 2007
- Volume 1 - 2006
Commun. Comput. Phys., 16 (2014), pp. 1201-1238.
Published online: 2014-11
Cited by
- BibTex
- RIS
- TXT
The solution of systems of hyperbolic conservation laws remains an interesting and challenging task due to the diversity of physical origins and complexity of the physical situations. The present work introduces the use of the partial differential equation (PDE) transform, paired with the Fourier pseudospectral method (FPM), as a new approach for hyperbolic conservation law problems. The PDE transform, based on the scheme of adaptive high order evolution PDEs, has recently been applied to decompose signals, images, surfaces and data to various target functional mode functions such as trend, edge, texture, feature, trait, noise, etc. Like wavelet transform, the PDE transform has controllable time-frequency localization and perfect reconstruction. A fast PDE transform implemented by the fast Fourier Transform (FFT) is introduced to avoid stability constraint of integrating high order PDEs. The parameters of the PDE transform are adaptively computed to optimize the weighted total variation during the time integration of conservation law equations. A variety of standard benchmark problems of hyperbolic conservation laws is employed to systematically validate the performance of the present PDE transform based FPM. The impact of two PDE transform parameters, i.e., the highest order and the propagation time, is carefully studied to deliver the best effect of suppressing Gibbs' oscillations. The PDE orders of 2-6 are used for hyperbolic conservation laws of low oscillatory solutions, while the PDE orders of 8-12 are often required for problems involving highly oscillatory solutions, such as shock-entropy wave interactions. The present results are compared with those in the literature. It is found that the present approach not only works well for problems that favor low order shock capturing schemes, but also exhibits superb behavior for problems that require the use of high order shock capturing methods.
}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.031113.160514a}, url = {http://global-sci.org/intro/article_detail/cicp/7078.html} }The solution of systems of hyperbolic conservation laws remains an interesting and challenging task due to the diversity of physical origins and complexity of the physical situations. The present work introduces the use of the partial differential equation (PDE) transform, paired with the Fourier pseudospectral method (FPM), as a new approach for hyperbolic conservation law problems. The PDE transform, based on the scheme of adaptive high order evolution PDEs, has recently been applied to decompose signals, images, surfaces and data to various target functional mode functions such as trend, edge, texture, feature, trait, noise, etc. Like wavelet transform, the PDE transform has controllable time-frequency localization and perfect reconstruction. A fast PDE transform implemented by the fast Fourier Transform (FFT) is introduced to avoid stability constraint of integrating high order PDEs. The parameters of the PDE transform are adaptively computed to optimize the weighted total variation during the time integration of conservation law equations. A variety of standard benchmark problems of hyperbolic conservation laws is employed to systematically validate the performance of the present PDE transform based FPM. The impact of two PDE transform parameters, i.e., the highest order and the propagation time, is carefully studied to deliver the best effect of suppressing Gibbs' oscillations. The PDE orders of 2-6 are used for hyperbolic conservation laws of low oscillatory solutions, while the PDE orders of 8-12 are often required for problems involving highly oscillatory solutions, such as shock-entropy wave interactions. The present results are compared with those in the literature. It is found that the present approach not only works well for problems that favor low order shock capturing schemes, but also exhibits superb behavior for problems that require the use of high order shock capturing methods.