Remote Access Mathematics of Computation
Green Open Access

Mathematics of Computation

ISSN 1088-6842(online) ISSN 0025-5718(print)

 
 

 

Local oscillations in finite difference solutions of hyperbolic conservation laws


Authors: Jiequan Li, Huazhong Tang, Gerald Warnecke and Lumei Zhang
Journal: Math. Comp. 78 (2009), 1997-2018
MSC (2000): Primary 65M06, 65T50; Secondary 35L40, 35L65, 76M20
DOI: https://doi.org/10.1090/S0025-5718-09-02219-4
Published electronically: January 28, 2009
MathSciNet review: 2521276
Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: It was generally expected that monotone schemes are oscillation-free for hyperbolic conservation laws. However, recently local oscillations were observed and usually understood to be caused by relative phase errors. In order to further explain this, we first investigate the discretization of initial data that trigger the chequerboard mode, the highest frequency mode. Then we proceed to use the discrete Fourier analysis and the modified equation analysis to distinguish the dissipative and dispersive effects of numerical schemes for low frequency and high frequency modes, respectively. It is shown that the relative phase error is of order $ {\mathcal O}(1)$ for the high frequency modes $ u_j^n=\lambda^n_k e^{i\xi j}$, $ \xi\approx \pi$, but of order $ {\mathcal O}(\xi^2)$ for low frequency modes ( $ \xi \approx 0$). In order to avoid numerical oscillations, the relative phase errors should be offset by numerical dissipation of at least the same order. Numerical damping, i.e. the zero order term in the corresponding modified equation, is important to dissipate the oscillations caused by the relative phase errors of high frequency modes. This is in contrast to the role of numerical viscosity, the second order term, which is the lowest order term usually present to suppress the relative phase errors of low frequency modes.


References [Enhancements On Off] (What's this?)

  • 1. M. Breuss, The correct use of the Lax-Friedrichs method, M2AN Math. Model. Numer. Anal. 38(2004), 519-540. MR 2075758 (2006e:65137)
  • 2. M. Breuss, An analysis of the influence of data extrema on some first and second order central approximations of hyperbolic conservation laws, M2AN Math. Model. Numer. Anal., 39(2005), 965-993. MR 2178569 (2007f:35183)
  • 3. A. Harten, J. M. Hyman and P. D. Lax, On finite-difference approximations and entropy conditions for shocks, Comm. Pure Appl. Math., 29(1976), 297-322. MR 0413526 (54:1640)
  • 4. P. Lax and B. Wendroff, Systems of conservation laws, Comm. Pure Appl. Math., 13(1960), 217-237. MR 0120774 (22:11523)
  • 5. P.D. LeFloch and J.-G. Liu, Generalized monotone schemes, discrete and paths of extrema, and discrete entropy conditions, Math. Comp., 68(1999), 1025-1055. MR 1627801 (99i:65092)
  • 6. R. LeVeque, Numerical Methods for Conservation Laws, Lectures in Mathematics ETH Zurich, 2nd ed. Basel, Birkhäuser Verlag, 1992. MR 1153252 (92m:65106)
  • 7. K.W. Morton and D.F. Mayers, Numerical Solution of Partial Differential Equations, 2nd edition, Cambridge University Press, 2005. MR 2153063 (2006a:65003)
  • 8. E. Tadmor, Numerical viscosity and the entropy condition for conservative difference schemes. Math. Comp., 43(1984), 369-381. MR 758189 (86g:65163)
  • 9. H.-Z. Tang and G. Warnecke, A note on $ (2K+1)$-point conservative monotone schemes, M2AN Math. Model. Numer. Anal., 38(2004), 345-357. MR 2069150 (2005f:65111)
  • 10. J.W. Thomas, Numerical Partial Differential Equations: Finite Difference Methods, Springer-Verlag, 1995. MR 1367964 (97a:65001)
  • 11. R. F. Warming and B. J. Hyett, The modified equation approach to the stability and accuracy of finite difference methods, J. Comput. Phys., 14(1974), 159-179. MR 0339526 (49:4284)

Similar Articles

Retrieve articles in Mathematics of Computation with MSC (2000): 65M06, 65T50, 35L40, 35L65, 76M20

Retrieve articles in all journals with MSC (2000): 65M06, 65T50, 35L40, 35L65, 76M20


Additional Information

Jiequan Li
Affiliation: School of Mathematics, Capital Normal University, Beijing 100037, People’s Republic of China
Email: jiequan@mail.cnu.edu.cn

Huazhong Tang
Affiliation: LMAM, School of Mathematical Sciences, Peking University, Beijing 100871, People’s Republic of China
Email: hztang@math.pku.edu.cn

Gerald Warnecke
Affiliation: Institut für Analysis und Numerik, Otto-von-Guericke-Universität, PSF 4120, 39016 Magdeburg, F.R. Germany
Email: warnecke@ovgu.de

Lumei Zhang
Affiliation: The high school attached to the Central University for Nationalities, Beijing 100081, People’s Republic of China
Email: zhanglumei99@163.com

DOI: https://doi.org/10.1090/S0025-5718-09-02219-4
Keywords: Finite difference schemes, high and low frequency modes, oscillations, chequerboard modes, numerical damping, numerical viscosity, relative phase error, modified equation analysis, discrete Fourier analysis.
Received by editor(s): March 14, 2008
Received by editor(s) in revised form: September 13, 2008
Published electronically: January 28, 2009
Article copyright: © Copyright 2009 American Mathematical Society

American Mathematical Society