Overview Statistic: PDF-Downloads (blue) and Frontdoor-Views (gray)
  • search hit 7 of 15
Back to Result List

The midpoint scheme and variants for Hamiltonian systems: advantages and pitfalls

Please always quote using this URN: urn:nbn:de:0297-zib-3201
  • The (implicit) midpoint scheme, like higher order Gauss-collocation schemes, is algebraically stable and symplectic, and it preserves quadratic integral invariants. It may appear particularly suitable for the numerical solution of highly oscillatory Hamiltonian systems, such as those arising in molecular dynamics or structural mechanics, because there is no stability restriction when it is applied to a simple harmonic oscillator. Although it is well-known that the midpoint scheme may also exhibit instabilities in various stiff situations, one might still hope for good results when resonance-type instabilities are avoided. In this paper we investigate the suitability of the midpoint scheme for highly oscillatory, frictionless mechanical systems, where the step-size $k$ is much larger than the system's small parameter $\varepsilon$, in case that the solution remains bounded as $\varepsilon \rightarrow 0$. We show that in general one must require that $k^2/\varepsilon$ be small enough, or else, even the errors in slowly varying quantities like the energy may grow undesirably (especially when fast and slow modes are tightly coupled) or, worse, the computation may yield misleading information. In some cases this may already happen when $k = O(\varepsilon )$. The same holds for higher order collocation at Gaussian points. The encountered restrictions on $k$ are still better than the corresponding ones for explicit schemes.

Download full text files

Export metadata

Additional Services

Share in Twitter Search Google Scholar Statistics - number of accesses to the document
Metadaten
Author:Uri Ascher, Sebastian Reich
Document Type:ZIB-Report
Date of first Publication:1997/11/05
Series (Serial Number):ZIB-Report (SC-97-51)
ZIB-Reportnumber:SC-97-51
Published in:Appeared in: SIAM J. Scient. Comp. 21 (1999) 1045-1065
Accept ✔
Diese Webseite verwendet technisch erforderliche Session-Cookies. Durch die weitere Nutzung der Webseite stimmen Sie diesem zu. Unsere Datenschutzerklärung finden Sie hier.