We conduct long simulations with a Hamiltonian particle-mesh method for ideal fluid flow, to determine the statistical mean vorticity field. Lagrangian and Eulerian statistical models are proposed for the discrete dynamics, and these are compared against numerical experiments. The observed results are in excellent agreement with theoretical models, as well as the continuum statistical mechanical theory for ideal flow developed by Eillis, Haven & Turkington (Nonlinearity, 2002). In particular the results verify that the apparently trivial conservation of potential vorticity along particle paths using the HPM method significantly influences the mean state. As a side note, the numerical experiments show that a nonzero fourth moment of potential vorticity can influence the statistical mean.
conservative discretizations, statistical mechanics, geometric numerical integration, quasigeostrophic flow, geophysical fluid dynamics
Meteorology and atmospheric physics (msc 86A10)
Modelling, Analysis and Simulation [MAS]
The investigations were in part supported by the Research Council for Earth and Life Sciences (ALW) with financial aid from the Netherlands Organization for Scientific Research (NWO).
Computational Dynamics

Dubinkina, S, & Frank, J.E. (2009). Statistical relevance of vorticity conservation with the Hamiltonian particle-mesh method. Modelling, Analysis and Simulation [MAS]. CWI.