Space-time FLAVORS: finite difference, multisymlectic, and pseudospectral integrators for multiscale PDEs

We present a new class of integrators for stiff PDEs. These integrators are generalizations of FLow AVeraging integratORS (FLAVORS) for stiff ODEs and SDEs introduced in [Tao, Owhadi and Marsden 2010] with the following properties: (i) Multiscale: they are based on flow averaging and have a computational cost determined by mesoscopic steps in space and time instead of microscopic steps in space and time; (ii) Versatile: the method is based on averaging the flows of the given PDEs (which may have hidden slow and fast processes). This bypasses the need for identifying explicitly (or numerically) the slow variables or reduced effective PDEs; (iii) Nonintrusive: A pre-existing numerical scheme resolving the microscopic time scale can be used as a black box and easily turned into one of the integrators in this paper by turning the large coefficients on over a microscopic timescale and off during a mesoscopic timescale; (iv) Convergent over two scales: strongly over slow processes and in the sense of measures over fast ones; (v) Structure-preserving: for stiff Hamiltonian PDEs (possibly on manifolds), they can be made to be multi-symplectic, symmetry-preserving (symmetries are group actions that leave the system invariant) in all variables and variational

Multi-symplectic discretisation of wave map equations

We present a new multi-symplectic formulation of constrained Hamiltonian partial differential equations, and we study the associated local conservation laws. A multi-symplectic discretisation based on this new formulation is exemplified by means of the Euler box scheme. When applied to the wave map equation, this numerical scheme is explicit, preserves the constraint and can be seen as a generalisation of the Shake algorithm for constrained mechanical systems. Furthermore, numerical experiments show excellent conservation properties of the numerical solutions

On the multi-symplectic structure of the Serre-Green-Naghdi equations

In this short note, we present a multi-symplectic structure of the Serre-Green-Naghdi (SGN) equations modelling nonlinear long surface waves in shallow water. This multi-symplectic structure allow the use of efficient finite difference or pseudo-spectral numerical schemes preserving exactly the multi-symplectic form at the discrete level.Comment: 10 pages, 1 figure, 30 references. Other author's papers can be downloaded at http://www.denys-dutykh.com

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