Numerical integration of the Ostrovsky equation based on its geometric structures

  • Authors:
  • Yuto Miyatake;Takaharu Yaguchi;Takayasu Matsuo

  • Affiliations:
  • Department of Mathematical Informatics, Graduate School of Information Science and Technology, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan;Department of Computational Science, Graduate School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan;Department of Mathematical Informatics, Graduate School of Information Science and Technology, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2012

Quantified Score

Hi-index 31.45

Visualization

Abstract

We consider structure preserving numerical schemes for the Ostrovsky equation, which describes gravity waves under the influence of Coriolis force. This equation has two associated invariants: an energy function and the L^2 norm. It is widely accepted that structure preserving methods such as invariants-preserving and multi-symplectic integrators generally yield qualitatively better numerical results. In this paper we propose five geometric integrators for this equation: energy-preserving and norm-preserving finite difference and Galerkin schemes, and a multi-symplectic integrator based on a newly found multi-symplectic formulation. A numerical comparison of these schemes is provided, which indicates that the energy-preserving finite difference schemes are more advantageous than the other schemes.