A General Approach for Producing Hamiltonian Numerical Schemes for Fluid Equations

arXivmath/0501468MaRDI QIDQ6475007FDOQ6475007

Authors: C. J. Cotter 

Publication date: 26 January 2005

Abstract: Given a fluid equation with reduced Lagrangian $l$ which is a functional of velocity $MMu$ and advected density $D$ given in Eulerian coordinates, we give a general method for semidiscretising the equations to give a canonical Hamiltonian system; this system may then be integrated in time using a symplectic integrator. The method is Lagrangian, with the variables being a set of Lagrangian particle positions and their associated momenta. The canonical equations obtained yield a discrete form of Euler-Poincar'e equations for $l$ when projected onto the grid, with a new form of discrete calculus to represent the gradient and divergence operators. Practical symplectic time integrators are suggested for a large family of equations which include the shallow-water equations, the EP-Diff equations and the 3D compressible Euler equations, and we illustrate the technique by showing results from a numerical experiment for the EP-Diff equations.

Mathematics Subject Classification ID

Particle methods and lattice-gas methods (76M28)

This page was built for publication: A General Approach for Producing Hamiltonian Numerical Schemes for Fluid Equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6475007)