A tailored solver for bifurcation analysis of ocean-climate models

Arie de Niet, Fred Wubs*, Arjen Terwisscha van Scheltinga, Henk A. Dijkstra

*Corresponding author for this work

Research output: Contribution to journalArticleAcademicpeer-review

24 Citations (Scopus)
388 Downloads (Pure)

Abstract

In this paper, we present a new linear system solver for use in a fully-implicit ocean model. The new solver allows to perform bifurcation analysis of relatively high-resolution primitive-equation ocean-climate models. It is based on a block-ILU approach and takes special advantage of the mathematical structure of the governing equations. In implicit models Jacobian matrices have to be constructed. Analytical construction is hard for complicated but more realistic representations of mixing. This is overcome by evaluating the Jacobian in part numerically. The performance of the new implicit ocean model is demonstrated using (i) a high-resolution model of the wind-forced double-gyre flow problem in a (relatively small) midlatitude spherical basin, and (ii) a medium-resolution model of thermohaline and wind-driven flows in an Atlantic size single-hemispheric basin. (c) 2007 Elsevier Inc. All rights reserved.

Original languageEnglish
Pages (from-to)654-679
Number of pages26
JournalJournal of computational physics
Volume227
Issue number1
DOIs
Publication statusPublished - 10-Nov-2007

Keywords

  • LOW-FREQUENCY VARIABILITY
  • SHALLOW-WATER MODELS
  • FULLY-IMPLICIT MODEL
  • WIND-DRIVEN
  • DOUBLE-GYRE
  • THERMOHALINE CIRCULATION
  • TRACER TRANSPORTS
  • ALGORITHM
  • MATRICES
  • GMRES

Cite this