New variational and multisymplectic formulations of the Euler-Poincar\'e equation on the Virasoro-Bott group using the inverse map

We derive a new variational principle, leading to a new momentum map and a new multisymplectic formulation for a family of Euler--Poincar\'e equations defined on the Virasoro-Bott group, by using the inverse map (also called `back-to-labels' map). This family contains as special cases the well-known Korteweg-de Vries, Camassa-Holm, and Hunter-Saxton soliton equations. In the conclusion section, we sketch opportunities for future work that would apply the new Clebsch momentum map with $2$-cocycles derived here to investigate a new type of interplay among nonlinearity, dispersion and noise.Comment: 19 page

The design of conservative finite element discretisations for the vectorial modified KdV equation

We design a consistent Galerkin scheme for the approximation of the vectorial modified Korteweg-de Vries equation. We demonstrate that the scheme conserves energy up to machine precision. In this sense the method is consistent with the energy balance of the continuous system. This energy balance ensures there is no numerical dissipation allowing for extremely accurate long time simulations free from numerical artifacts. Various numerical experiments are shown demonstrating the asymptotic convergence of the method with respect to the discretisation parameters. Some simulations are also presented that correctly capture the unusual interactions between solitons in the vectorial setting

The design of conservative finite element discretisations for the vectorial modified KdV equation

We design a consistent Galerkin scheme for the approximation of the vectorial modified Korteweg–de Vries equation with periodic boundary conditions. We demonstrate that the scheme conserves energy up to solver tolerance. In this sense the method is consistent with the energy balance of the continuous system. This energy balance ensures there is no numerical dissipation allowing for extremely accurate long time simulations free from numerical artifacts. Various numerical experiments are shown demonstrating the asymptotic convergence of the method with respect to the discretisation parameters. Some simulations are also presented that correctly capture the unusual interactions between solitons in the vectorial setting.</p

An explicit finite difference scheme for the Camassa-Holm equation

We put forward and analyze an explicit finite difference scheme for the Camassa-Holm shallow water equation that can handle general $H^1$ initial data and thus peakon-antipeakon interactions. Assuming a specified condition restricting the time step in terms of the spatial discretization parameter, we prove that the difference scheme converges strongly in $H^1$ towards a dissipative weak solution of Camassa-Holm equation.Comment: 45 pages, 6 figure

A numerical study of variational discretizations of the Camassa-Holm equation

We present two semidiscretizations of the Camassa-Holm equation in periodic domains based on variational formulations and energy conservation. The first is a periodic version of an existing conservative multipeakon method on the real line, for which we propose efficient computation algorithms inspired by works of Camassa and collaborators. The second method, and of primary interest, is the periodic counterpart of a novel discretization of a two-component Camassa-Holm system based on variational principles in Lagrangian variables. Applying explicit ODE solvers to integrate in time, we compare the variational discretizations to existing methods over several numerical examples.Comment: 45 pages, 14 figure

Multi-symplectic integration of the Camassa–Holm equation

The Camassa–Holm equation is rich in geometric structures, it is completely integrable, bi-Hamiltonian, and it represents geodesics for a certain metric in the group of diffeomorphism. Here two new multi-symplectic formulations for the Camassa–Holm equation are presented, and the associated local conservation laws are shown to correspond to certain well-known Hamiltonian functionals. The multi-symplectic discretisation of each formulation is exemplified by means of the Euler box scheme. Numerical experiments show that the schemes have good conservative properties, and one of them is designed to handle the conservative continuation of peakon-antipeakon collisions.