Symmetric parareal algorithms for Hamiltonian systems - Inria - Institut national de recherche en sciences et technologies du numérique Access content directly
Journal Articles ESAIM: Mathematical Modelling and Numerical Analysis Year : 2013

Symmetric parareal algorithms for Hamiltonian systems

Abstract

The parareal in time algorithm allows to efficiently use parallel computing for the simulation 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 processors. 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 [Lions, Maday and Turinici 2001, Baffico, Bernard, Maday, Turinici and Zerah 2002, Bal and Maday 2002] 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.

Dates and versions

hal-00541166 , version 1 (30-11-2010)

Identifiers

Cite

Xiaoying Dai, Claude Le Bris, Frédéric Legoll, Yvon Maday. Symmetric parareal algorithms for Hamiltonian systems. ESAIM: Mathematical Modelling and Numerical Analysis, 2013, 47 (3), pp.717-742. ⟨10.1051/m2an/2012046⟩. ⟨hal-00541166⟩
289 View
0 Download

Altmetric

Share

Gmail Facebook X LinkedIn More