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 journalArticleAcademicpeer-review

10 Citations (Scopus)
227 Downloads (Pure)

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(10(5))) degrees of freedom such as multi-layer quasi-geostrophic and shallow-water models and relatively low-resolution (e.g., 4 degrees 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. (C) 2009 Elsevier Ltd. All rights reserved.

Original languageEnglish
Pages (from-to)287-297
Number of pages11
JournalOcean modelling
Volume30
Issue number4
DOIs
Publication statusPublished - 2009

Keywords

  • Bifurcation
  • Newton-Krylov methods
  • Implicit ocean model
  • Parallel implementation
  • Dynamical system
  • Trilinos
  • LOW-FREQUENCY VARIABILITY
  • NEWTON-KRYLOV METHODS
  • SHALLOW-WATER MODELS
  • WIND-DRIVEN
  • MULTIPLE EQUILIBRIA
  • DOUBLE-GYRE
  • THERMOSOLUTAL CONVECTION
  • THERMOHALINE CONVECTION
  • LINEAR-SYSTEMS
  • SPIN-UP

Cite this