Accuracy and conservation properties in numerical integration: the case of the Korteweg-de Vries equation

When numerically integrating time-dependent differential equations, it is often recommended to employ methods that preserve some of the invariant quantities (mass, energy, etc.) of the problem being considered. This recommendation is usually justified on the grounds that conservation of invariant quantities may ensure that the numerical solution possesses some important qualitative features. However there are cases where schemes that preserve invariants are also advantageous in that they possess favourable error propagation mechanisms that render them superior from a quantitative point of view. In the present paper we consider the Korteweg-de Vries equation as a case study. We show rigorously that, for soliton problems and at leading order, the error of conservative schemes consists of a phase error that grows linearly with time plus a complementary term that is bounded in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document} $H^1$\end{document} norm uniformly in time. For ‘general’, nonconservative schemes the error involves a linearly growing amplitude error, a quadratically growing phase error and a complementary term that grows linearly in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document} $H^1$\end{document} norm. Numerical experiments are presented.