Local oscillations in finite difference solutions of hyperbolic conservation laws

Li, Jiequan; Tang, Huazhong; Warnecke, Gerald; Zhang, Lumei

doi:10.1090/S0025-5718-09-02219-4

Local oscillations in finite difference solutions of hyperbolic conservation laws
HTML articles powered by AMS MathViewer

by Jiequan Li, Huazhong Tang, Gerald Warnecke and Lumei Zhang;

Math. Comp. 78 (2009), 1997-2018

DOI: https://doi.org/10.1090/S0025-5718-09-02219-4

Published electronically: January 28, 2009

PDF | Request permission

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

Michael Breuss, The correct use of the Lax-Friedrichs method, M2AN Math. Model. Numer. Anal. 38 (2004), no. 3, 519–540. MR 2075758, DOI 10.1051/m2an:2004027
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. MR 2178569, DOI 10.1051/m2an:2005042
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. MR 413526, DOI 10.1002/cpa.3160290305
Peter Lax and Burton Wendroff, Systems of conservation laws, Comm. Pure Appl. Math. 13 (1960), 217–237. MR 120774, DOI 10.1002/cpa.3160130205
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. MR 1627801, DOI 10.1090/S0025-5718-99-01062-5
Randall J. LeVeque, Numerical methods for conservation laws, 2nd ed., Lectures in Mathematics ETH Zürich, Birkhäuser Verlag, Basel, 1992. MR 1153252, DOI 10.1007/978-3-0348-8629-1
K. W. Morton and D. F. Mayers, Numerical solution of partial differential equations, 2nd ed., Cambridge University Press, Cambridge, 2005. An introduction. MR 2153063, DOI 10.1017/CBO9780511812248
Eitan Tadmor, Numerical viscosity and the entropy condition for conservative difference schemes, Math. Comp. 43 (1984), no. 168, 369–381. MR 758189, DOI 10.1090/S0025-5718-1984-0758189-X
Huazhong Tang and Gerald Warnecke, A note on $(2K+1)$-point conservative monotone schemes, M2AN Math. Model. Numer. Anal. 38 (2004), no. 2, 345–357. MR 2069150, DOI 10.1051/m2an:2004016
J. W. Thomas, Numerical partial differential equations: finite difference methods, Texts in Applied Mathematics, vol. 22, Springer-Verlag, New York, 1995. MR 1367964, DOI 10.1007/978-1-4899-7278-1
R. F. Warming and B. J. Hyett, The modified equation approach to the stability and accuracy analysis of finite-difference methods, J. Comput. Phys. 14 (1974), 159–179. MR 339526, DOI 10.1016/0021-9991(74)90011-4