Statistical mechanics of Arakawa's discretizations

  • Authors:
  • Svetlana Dubinkina;Jason Frank

  • Affiliations:
  • CWI, P.O. Box 94079, 1090 GB Amsterdam, The Netherlands;CWI, P.O. Box 94079, 1090 GB Amsterdam, The Netherlands

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2007

Quantified Score

Hi-index 31.45

Visualization

Abstract

The results of statistical analysis of simulation data obtained from long time integrations of geophysical fluid models greatly depend on the conservation properties of the numerical discretization used. This is illustrated for quasi-geostrophic flow with topographic forcing, for which a well established statistical mechanics exists. Statistical mechanical theories are constructed for the discrete dynamical systems arising from three discretizations due to Arakawa [Arakawa, Computational design for long-term numerical integration of the equations of fluid motion: two-dimensional incompressible flow. Part I. J. Comput. Phys. 1 (1966) 119-143] which conserve energy, enstrophy or both. Numerical experiments with conservative and projected time integrators show that the statistical theories accurately explain the differences observed in statistics derived from the discretizations.