×

Local oscillations in finite difference solutions of hyperbolic conservation laws. (English) Zbl 1198.65169

Summary: 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.

MSC:

65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
65T50 Numerical methods for discrete and fast Fourier transforms
35L40 First-order hyperbolic systems
35L65 Hyperbolic conservation laws
76M20 Finite difference methods applied to problems in fluid mechanics
Full Text: DOI

References:

[1] Michael Breuss, The correct use of the Lax-Friedrichs method, M2AN Math. Model. Numer. Anal. 38 (2004), no. 3, 519 – 540. · Zbl 1077.65089 · doi:10.1051/m2an:2004027
[2] Michael 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), no. 5, 965 – 993. · Zbl 1077.35089 · doi:10.1051/m2an:2005042
[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), no. 3, 297 – 322. With an appendix by B. Keyfitz. · Zbl 0351.76070 · doi:10.1002/cpa.3160290305
[4] Peter Lax and Burton Wendroff, Systems of conservation laws, Comm. Pure Appl. Math. 13 (1960), 217 – 237. · Zbl 0152.44802 · doi:10.1002/cpa.3160130205
[5] Philippe G. LeFloch and Jian-Guo Liu, Generalized monotone schemes, discrete paths of extrema, and discrete entropy conditions, Math. Comp. 68 (1999), no. 227, 1025 – 1055. · Zbl 0915.35069
[6] Randall J. LeVeque, Numerical methods for conservation laws, 2nd ed., Lectures in Mathematics ETH Zürich, Birkhäuser Verlag, Basel, 1992. · Zbl 0847.65053
[7] K. W. Morton and D. F. Mayers, Numerical solution of partial differential equations, 2nd ed., Cambridge University Press, Cambridge, 2005. An introduction. · Zbl 1126.65077
[8] Eitan Tadmor, Numerical viscosity and the entropy condition for conservative difference schemes, Math. Comp. 43 (1984), no. 168, 369 – 381. · Zbl 0587.65058
[9] Huazhong Tang and Gerald Warnecke, A note on (2\?+1)-point conservative monotone schemes, M2AN Math. Model. Numer. Anal. 38 (2004), no. 2, 345 – 357. · Zbl 1075.65113 · doi:10.1051/m2an:2004016
[10] J. W. Thomas, Numerical partial differential equations: finite difference methods, Texts in Applied Mathematics, vol. 22, Springer-Verlag, New York, 1995. · Zbl 0831.65087
[11] R. F. Warming and B. J. Hyett, The modified equation approach to the stability and accuracy analysis of finite-difference methods, J. Computational Phys. 14 (1974), 159 – 179. · Zbl 0291.65023
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.