zbMATH — the first resource for mathematics

Statistical mechanics of Arakawa’s discretizations. (English) Zbl 1130.82005
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 exist. Statistical mechanical theories are constructed for the discrete dynamical systems arising from three discretizations due to A. Arakawa [J. Comput. Phys. 1, 119–143 (1966; Zbl 0147.44202)], which conserve energy, enstrophy or both. Numerical experiments with conservative and projected time integrators show that the statistical theories accurately explain the difference observed in statistics derived from the discretizations.
The authors of the present paper first briefly recall the quasi-geostrophic potential vorticity equation and its conservation properties. Then, they recall Arakawa’s three discretizations, their conservation properties, and they prove that all of these define divergence-free vector fields. Generalizations to Arakawa’s discretizations with the Nambu bracket approach are considered. Equilibrium statistical mechanical theories are developed for the three discretizations for an ideal fluid in vorticity-stream function form. After that the results are extended to the cases in which only one of the quantities (energy or enstrophy) is conserved. Time integration aspects (symmetric implicit midpoint method) are discussed in detail. Finally, numerical experiments confirming the statistical predictions are presented. The results of the paper make a strong argument for the use of conservative discretizations in weather and climate simulations.

82-08 Computational methods (statistical mechanics) (MSC2010)
82-01 Introductory exposition (textbooks, tutorial papers, etc.) pertaining to statistical mechanics
Full Text: DOI
[1] Abramov, R.V.; Majda, A.J., Statistically relevant conserved quantities for truncated quasigeostrophic flow, Proc. natl. acad. sci. USA, 100, 7, 3841-3846, (2003), (electronic) · Zbl 1063.86001
[2] Arakawa, A., Computational design for long-term numerical integration of the equations of fluid motion: two-dimensional incompressible flow. part I., J. comput. phys., 1, 119-143, (1966) · Zbl 0147.44202
[3] Carnevale, G.F.; Frederiksen, J.S., Nonlinear stability and statistical mechanics of flow over topography, J. fluid mech., 175, 157-181, (1987) · Zbl 0651.76020
[4] Frank, J.; Gottwald, G.; Reich, S., A Hamiltonian particle-mesh method for the rotating shallow-water equations, (), 131-142 · Zbl 1062.76045
[5] Frank, J.; Reich, S., Conservation properties of smoothed particle hydrodynamics applied to the shallow water equation, Bit, 43, 1, 41-55, (2003) · Zbl 1023.76037
[6] Frank, J.; Reich, S., The Hamiltonian particle-mesh method for the spherical shallow water equations, Atmos. sci. lett., 5, 89-95, (2004)
[7] Hairer, E.; Lubich, C.; Wanner, G., Geometric numerical integration, Springer series in computational mathematics, vol. 31, (2006), Springer-Verlag Berlin, structure-preserving algorithms for ordinary differential equations
[8] Hairer, E.; Nørsett, S.P.; Wanner, G., Solving ordinary differential equations. I, Springer series in computational mathematics, vol. 8, (1993), Springer-Verlag Berlin, nonstiff problems · Zbl 0789.65048
[9] Kraichnan, R.H., Statistical dynamics of two-dimensional flow, J. fluid mech., 67, 155-175, (1975) · Zbl 0297.76042
[10] Leimkuhler, B.; Reich, S., Simulating Hamiltonian dynamics, Cambridge monographs on applied and computational mathematics, vol. 14, (2004), Cambridge University Press Cambridge · Zbl 1069.65139
[11] Lorenz, E.N., Energy and numerical weather prediction, Tellus, 12, 364-373, (1960)
[12] Majda, A.J.; Wang, X., Non-linear dynamics and statistical theories for basic geophysical flows, (2006), Cambridge University Press Cambridge · Zbl 1141.86001
[13] McLachlan, R., Symplectic integration of Hamiltonian wave equations, Numer. math., 66, 4, 465-492, (1994) · Zbl 0831.65099
[14] McLachlan, R.I., Explicit Lie-Poisson integration and the Euler equations, Phys. rev. lett., 71, 19, 3043-3046, (1993) · Zbl 0972.65509
[15] Morrison, P.J., Hamiltonian description of the ideal fluid, Rev. modern phys., 70, 2, 467-521, (1998) · Zbl 1205.37093
[16] Nambu, Y., Generalized Hamiltonian dynamics, Phys. rev. D, 7, 3, 2405-2412, (1973) · Zbl 1027.70503
[17] Névir, P.; Blender, R., A Nambu representation of incompressible hydrodynamics using helicity and enstrophy, J. phys. A, 26, 22, L1189-L1193, (1993) · Zbl 0812.76005
[18] Pedlosky, J., Geophysical fluid dynamics, (2005), Springer
[19] Salmon, R., Lectures on geophysical fluid dynamics, (1998), Oxford University Press New York
[20] Salmon, R., Poisson-bracket approach to the construction of energy- and potential-enstrophy-conserving algorithms for the shallow-water equations, J. atmos. sci., 61, 16, 2016-2036, (2004)
[21] Salmon, R., A general method for conserving quantities related to potential vorticity in numerical models, Nonlinearity, 18, 5, R1-R16, (2005) · Zbl 1213.76143
[22] Salmon, R.; Holloway, G.; Hendershott, M.C., The equilibrium statistical mechanics of simple quasi-geostrophic models, J. fluid mech., 75, 691-703, (1976) · Zbl 0342.76015
[23] Salmon, R.; Talley, L.D., Generalizations of arakawa’s Jacobian, J. comput. phys., 83, 2, 247-259, (1989) · Zbl 0672.76002
[24] Sanz-Serna, J.M.; Calvo, M.P., Numerical Hamiltonian problems, Applied mathematics and mathematical computation, vol. 7, (1994), Chapman & Hall London · Zbl 0816.65042
[25] Stuart, A.M.; Humphries, A.R., Dynamical systems and numerical analysis, Cambridge monographs on applied and computational mathematics, vol. 2, (1996), Cambridge University Press Cambridge · Zbl 0869.65043
[26] Zeitlin, V., Finite-mode analogs of 2D ideal hydrodynamics: coadjoint orbits and local canonical structure, Physica. D, 49, 3, 353-362, (1991) · Zbl 0724.76004
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.