Computational design for long-term numerical integration of the equations of fluid motion: Two-dimensional incompressible flow. Part I

Abstract

The integral constraints on quadratic quantities of physical importance, such as conservation of mean kinetic energy and mean square vorticity, will not be maintained in finite difference analogues of the equation of motion for two-dimensional incompressible flow, unless the finite difference Jacobian expression for the advection term is restricted to a form which properly represents the interaction between grid points, as derived in this paper. It is shown that the derived form of the finite difference Jacobian prevents nonlinear computational instability and thereby permits long-term numerical integrations.