Picture of a black hole

Strathclyde Open Access research that creates ripples...

The Strathprints institutional repository is a digital archive of University of Strathclyde's Open Access research outputs. Strathprints provides access to thousands of research papers by University of Strathclyde researchers, including by Strathclyde physicists involved in observing gravitational waves and black hole mergers as part of the Laser Interferometer Gravitational-Wave Observatory (LIGO) - but also other internationally significant research from the Department of Physics. Discover why Strathclyde's physics research is making ripples...

Strathprints also exposes world leading research from the Faculties of Science, Engineering, Humanities & Social Sciences, and from the Strathclyde Business School.

Discover more...

Stability of numerical schemes on staggered grids

Oishi, C.M. and Cuminato, J.A. and Yuan, J.Y. and McKee, S. (2008) Stability of numerical schemes on staggered grids. Numerical Linear Algebra with Applications, 15 (10). pp. 945-967. ISSN 1070-5325

Full text not available in this repository. (Request a copy from the Strathclyde author)

Abstract

This paper considers the stability of explicit, implicit and Crank-Nicolson schemes for the one-dimensional heat equation on a staggered grid. Furthermore, we consider the cases when both explicit and implicit approximations of the boundary conditions are employed. Why we choose to do this is clearly motivated and arises from solving fluid flow equations with free surfaces when the Reynolds number can be very small, in at least parts of the spatial domain. A comprehensive stability analysis is supplied: a novel result is the precise stability restriction on the Crank-Nicolson method when the boundary conditions are approximated explicitly, that is, at t=nt rather than t=(n+1)t. The two-dimensional Navier-Stokes equations were then solved by a marker and cell approach for two simple problems that had analytic solutions. It was found that the stability results provided in this paper were qualitatively very similar, thereby providing insight as to why a Crank-Nicolson approximation of the momentum equations is only conditionally stable.