On a class of ocean model instabilities that may occur when applying small time steps, implicit methods, and low viscosities

The horizontal grid size in ocean models may be gradually reduced with increasing computer power. This allows simulations with less viscosity. The hyperbolic, free wave nature of the problem will become more dominant. However, as viscosities are reduced, instabilities may more easily occur. To avoid...

Full description

Saved in:
Bibliographic Details
Published in:Ocean modelling (Oxford) 2004, Vol.7 (1), p.135-144
Main Authors: Heimsund, Bjørn-Ove, Berntsen, Jarle
Format: Article
Language:English
Subjects:
C-grid
Ocean modelling
Predictor–corrector
ROMS
Stability analysis
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The horizontal grid size in ocean models may be gradually reduced with increasing computer power. This allows simulations with less viscosity. The hyperbolic, free wave nature of the problem will become more dominant. However, as viscosities are reduced, instabilities may more easily occur. To avoid such instabilities, it is well known from the literature on numerical methods for ordinary differential equations that implicit methods have attractive stability properties. Therefore, it may be tempting to add implicitness into the time stepping also in low viscosity ocean modelling to improve the robustness. However, when using methods with implicit features and low viscosity, it may happen that models are stable for longer time steps, but become unstable as the time step is reduced. In the present paper, an explanation for this phenomenon is given. It is shown that the explanation may be valid when solving the shallow water equations on a C-grid.
ISSN:1463-5003
1463-5011
DOI:10.1016/S1463-5003(03)00041-6