Bifurcation analysis of 3D ocean flows using a parallel fully-implicit ocean model

Journal article (2009)

Authors

J. Thies University Medical Center Groningen
Fred W. Wubs University Medical Center Groningen
Henk A. Dijkstra Universiteit Utrecht
Affiliation
External organisation
Newton-Krylov methods Bifurcation Trilinos Dynamical system Implicit ocean model Parallel implementation
More Info
expand_more

Published Date

2009

Language

English

Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Affiliation

External organisation

Abstract

To understand the physics and dynamics of the ocean circulation, techniques of numerical bifurcation theory such as continuation methods have proved to be useful. Up to now these techniques have been applied to models with relatively few (O (105)) degrees of freedom such as multi-layer quasi-geostrophic and shallow-water models and relatively low-resolution (e.g., 4° horizontal resolution) primitive equation models. In this paper, we present a new approach in which continuation methods are combined with parallel numerical linear system solvers. With this implementation, we show that it is possible to compute steady states versus parameters (and perform fully implicit time integration) of primitive equation ocean models with up to a few million degrees of freedom.