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

Jonas Thies*, Fred Wubs, Henk A. Dijkstra

*Corresponding author for this work

Research output: Contribution to journalArticleScientificpeer-review

10 Citations (Scopus)


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.

Original languageEnglish
Pages (from-to)287-297
Number of pages11
JournalOcean Modelling
Issue number4
Publication statusPublished - 2009
Externally publishedYes


  • Bifurcation
  • Dynamical system
  • Implicit ocean model
  • Newton-Krylov methods
  • Parallel implementation
  • Trilinos


Dive into the research topics of 'Bifurcation analysis of 3D ocean flows using a parallel fully-implicit ocean model'. Together they form a unique fingerprint.

Cite this