Poisson-Bracket-SWE.MWR2004.i1520-0469-61-16-2016

Poisson-Bracket-SWE.MWR2004.i1520-0469-61-16-2016 - 2016...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
2016 V OLUME 61 JOURNAL OF THE ATMOSPHERIC SCIENCES q 2004 American Meteorological Society Poisson-Bracket Approach to the Construction of Energy- and Potential-Enstrophy- Conserving Algorithms for the Shallow-Water Equations R ICK SALMON Scripps Institution of Oceanography, University of California, San Diego, La Jolla, California (Manuscript received 11 September 2003, in final form 27 February 2004) ABSTRACT Arakawa and Lamb discovered a finite-difference approximation to the shallow-water equations that exactly conserves finite-difference approximations to the energy and potential enstrophy of the fluid. The Arakawa– Lamb (AL) algorithm is a stunning and important achievement—stunning, because in the shallow-water case, neither energy nor potential enstrophy is a simple quadratic, and important because the simultaneous conservation of energy and potential enstrophy is known to prevent the spurious cascade of energy to high wavenumbers. However, the method followed by AL is somewhat ad hoc, and it is difficult to see how it might be generalized to other systems. In this paper, the AL algorithm is rederived and greatly generalized in a way that should permit still further generalizations. Beginning with the Hamiltonian formulation of shallow-water dynamics, its two essential in- gredients—the Hamiltonian functional and the Poisson-bracket operator—are replaced by finite-difference ap- proximations that maintain the desired conservation laws. Energy conservation is maintained if the discrete Poisson bracket retains the antisymmetry property of the exact bracket, a trivial constraint. Potential enstrophy is conserved if a set of otherwise arbitrary coefficients is chosen in such a way that a very large quadratic form contains only diagonal terms. Using a symbolic manipulation program to satisfy the potential-enstrophy con- straint, it is found that the energy- and potential-enstrophy-conserving schemes corresponding to a stencil of 25 grid points contain 22 free parameters. The AL scheme corresponds to the vanishing of all free parameters. No parameter setting can increase the overall accuracy of the schemes beyond second order, but 19 of the free parameters may be independently adjusted to yield a scheme with fourth-order accuracy in the vorticity equation. 1. Introduction In a remarkable paper, Arakawa and Lamb (1981, hereafter AL) discovered a finite-difference scheme for the shallow-water equations that—apart from errors as- sociated with the finite time step—exactly conserves finite-difference analogs of the energy and potential en- strophy of the flow. The practical importance of retain- ing energy- and enstrophy-conservation laws in nu- merical models had been appreciated since the earlier work of Arakawa (1966) on two-dimensional, nondiv- ergent flow; but whereas the derivation of Arakawa’s Jacobian for nondivergent flow may be transparently understood in a variety of ways (see, e.g., Salmon and Talley 1989), the AL algorithm for shallow-water dy-
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Image of page 2
This is the end of the preview. Sign up to access the rest of the document.

This note was uploaded on 01/08/2012 for the course MPO 662 taught by Professor Iskandarani,m during the Spring '08 term at University of Miami.

Page1 / 21

Poisson-Bracket-SWE.MWR2004.i1520-0469-61-16-2016 - 2016...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online