A Fully Implicit Model of the Three-Dimensional Thermohaline Ocean Circulation

Henk A. Dijkstra, Hakan Oksuzoglu, Fred. W. Wubs, Eugen F.F. Botta

Research output: Contribution to journalArticleAcademicpeer-review

34 Citations (Scopus)
193 Downloads (Pure)

Abstract

In this paper, a fully implicit numerical model of the three-dimensional thermohaline ocean circulation is presented. With this numerical model it is possible to follow branches of steady states in parameter space and monitor their linear stability. Also, transient flows can be computed allowing much larger time steps than those possible with explicit schemes. By using recently developed solvers for linear systems of equations and for generalized eigenvalue problems, results for reasonable spatial resolution can be obtained. Bifurcation diagrams and transient flows are computed for typical flows in a single hemispheric basin situation, with focus on (i) the performance of the methodology and (ii) the new type of information which can be obtained on these flows.
Original languageEnglish
Pages (from-to)685-715
Number of pages31
JournalJournal of computational physics
Volume173
Issue number2
DOIs
Publication statusPublished - 1-Nov-2001

Keywords

  • stability of geophysical flows
  • numerical bifurcation theory
  • iterative linear systems solvers
  • MULTILEVEL ILU PRECONDITIONER
  • MULTIPLE EQUILIBRIA
  • ATMOSPHERIC CO2
  • VARIABILITY
  • MATRICES
  • SENSITIVITY
  • CONVECTION
  • STABILITY
  • FLOWS
  • WIND

Cite this