×

Multi-symplectic integrators: numerical schemes for Hamiltonian PDEs that conserve symplecticity. (English) Zbl 0984.37104

Summary: The symplectic numerical integration of finite-dimensional Hamiltonian systems is a well established subject and has led to a deeper understanding of existing methods as well as to the development of new very efficient and accurate schemes, e.g., for rigid body, constrained, and molecular dynamics. The numerical integration of infinite-dimensional Hamiltonian systems or Hamiltonian PDEs is much less explored. In this Letter, we suggest a new theoretical framework for generalizing symplectic numerical integrators for ODEs to Hamiltonian PDEs in \(\mathbb{R}^2\): time plus one space dimension. The central idea is that symplecticity for Hamiltonian PDEs is directional: the symplectic structure of the PDE is decomposed into distinct components representing space and time independently. In this setting PDE integrators can be constructed by concatenating uni-directional ODE symplectic integrators. This suggests a natural definition of multi-symplectic integrator as a discretization that conserves a discrete version of the conservation of symplecticity for Hamiltonian PDEs. We show that this approach leads to a general framework for geometric numerical schemes for Hamiltonian PDEs, which have remarkable energy and momentum conservation properties. Generalizations, including development of higher-order methods, application to the Euler equations in fluid mechanics, application to perturbed systems, and extension to more than one space dimension are also discussed.

MSC:

37M15 Discretization methods and integrators (symplectic, variational, geometric, etc.) for dynamical systems
37J05 Relations of dynamical systems with symplectic geometry and topology (MSC2010)
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Abbott, M. B., Computational Hydraulics (1979), Pitman: Pitman London · Zbl 0507.76026
[2] Abbott, M. B.; Basco, D. R., Computational Fluid Dynamics (1989), Longman Scientific
[3] Benettin, G.; Giorgilli, A., J. Stat. Phys., 74, 1117 (1994)
[4] Binz, E.; Śniatycki, J.; Fischer, H., Geometry of Classical Fields (1988), North-Holland · Zbl 0675.53065
[5] Bridges, T. J., Math. Proc. Cambridge Philos. Soc., 121, 147 (1997)
[6] Bridges, T. J., Proc. R. Soc. London A, 453, 1365 (1997)
[7] Bridges, T. J., Eur. J. Mech. B/Fluids, 18, 493 (1999)
[8] Bridges, T. J.; Derks, G., Proc. R. Soc. London A, 455, 2427 (1999)
[9] T.J. Bridges, S. Reich, in preparation (2001); T.J. Bridges, S. Reich, in preparation (2001)
[10] (Deuflhard, P.; Hermans, J.; Leimkuhler, B.; Mark, A. E.; Reich, S.; Skeel, R. D., Lecture Notes in Computational Science and Engineering, 4 (1999), Springer)
[11] García, P. L., Sympos. Math., 14, 219 (1974)
[12] Gotay, M. J., (Francaviglia, M., Mechanics, Analysis and Geometry: 200 Years After Lagrange (1991), North-Holland: North-Holland Amsterdam), 203-235
[13] Hairer, E.; Lubich, Ch., Numer. Math., 76, 441 (1997)
[14] Kijowski, J.; Tulczyjew, W., A Symplectic Framework for Field Theories (1979), Springer · Zbl 0439.58002
[15] Marsden, J. E.; Patrick, G. P.; Shkoller, S., Comm. Math. Phys., 199, 351 (1999)
[16] Marsden, J. E.; Shkoller, S., Math. Proc. Cambridge Philos. Soc., 125, 553 (1999)
[17] McLachlan, R. I., Numer. Math., 66, 465 (1994)
[18] B. Moore, Ph.D. progress report, University of Surrey (2000); B. Moore, Ph.D. progress report, University of Surrey (2000)
[19] Reich, S., SIAM Numer. Anal., 36, 1549 (1999)
[20] Reich, S., J. Comput. Phys., 157, 473 (2000)
[21] Reich, S., BIT, 40, 559 (2000)
[22] Sanz-Serna, J. M.; Calvo, M. P., Numerical Hamiltonian Systems (1994), Chapman and Hall · Zbl 0816.65042
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.