On numerical methods for Hamiltonian PDEs and a collocation method for the Vlasov-Maxwell equations (Q2564465)

The purpose of this paper is to examine how the conserved quantities -- and particularly the Hamiltonian itself -- of Hamiltonian partial differential equations can be preserved in weighted residual and collocation methods by examining the properties of the Poisson bracket. By using weighted residual methods and exploiting the Hamiltonian structure of the equations it is shown at least some of the conservation laws are preserved in a method of lines (continuous in time). In particular, the Hamiltonian can always be exactly preserved as a constant of the motion. Other conservation laws, in particular linear and quadratic Casimirs and momenta, can sometimes be conserved too, depending on the details of the equations under consideration and the form of discretization employed. Collocation methods also offer automatic conservation of linear and quadratic Casimirs. Some standard discretization methods, when applied to Hamiltonian problems are shown to be derived from a numerical approximation to the exact Poisson bracket of the system. A method for the Vlasov-Maxwell equations based on Legendre-Gauss-Lobatto collocation is presented as an example of these ideas.