Volume 47, Number 3, May-June 2013
|Page(s)||717 - 742|
|Published online||04 March 2013|
Symmetric parareal algorithms for Hamiltonian systems
1 LSEC, Institute of Computational
Mathematics and Scientific/Engineering Computing, Academy of Mathematics and Systems
Science, Chinese Academy of Sciences, Beijing
2 École Nationale des Ponts et Chaussées, 6 et 8 avenue Blaise Pascal, 77455 Marne-La-Vallée Cedex 2, France
3 INRIA Rocquencourt, MICMAC team-project, Domaine de Voluceau, B.P. 105, 78153 Le Chesnay Cedex, France
4 UPMC Univ. Paris 06, UMR 7598, Laboratoire J.-L. Lions, Boîte courrier 187, 75252 Paris Cedex 05, France
5 Division of Applied Mathematics, Brown University, Providence, RI, USA
Revised: 14 February 2012
The parareal in time algorithm allows for efficient parallel numerical simulations of time-dependent problems. It is based on a decomposition of the time interval into subintervals, and on a predictor-corrector strategy, where the propagations over each subinterval for the corrector stage are concurrently performed on the different processors that are available. In this article, we are concerned with the long time integration of Hamiltonian systems. Geometric, structure-preserving integrators are preferably employed for such systems because they show interesting numerical properties, in particular excellent preservation of the total energy of the system. Using a symmetrization procedure and/or a (possibly also symmetric) projection step, we introduce here several variants of the original plain parareal in time algorithm [L. Baffico, et al. Phys. Rev. E 66 (2002) 057701; G. Bal and Y. Maday, A parareal time discretization for nonlinear PDE’s with application to the pricing of an American put, in Recent developments in domain decomposition methods, Lect. Notes Comput. Sci. Eng. 23 (2002) 189–202; J.-L. Lions, Y. Maday and G. Turinici, C. R. Acad. Sci. Paris, Série I 332 (2001) 661–668.] that are better adapted to the Hamiltonian context. These variants are compatible with the geometric structure of the exact dynamics, and are easy to implement. Numerical tests on several model systems illustrate the remarkable properties of the proposed parareal integrators over long integration times. Some formal elements of understanding are also provided.
Mathematics Subject Classification: 65L05 / 65P10 / 65Y05
Key words: Parallel integrators / Hamiltonian dynamics / long-time integration / symmetric algorithms / symmetric projection / geometric integration
© EDP Sciences, SMAI, 2013
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.