A parallel Jacobian-free Newton-Krylov solver for a coupled sea ice-ocean model

  • Authors:

  • Martin Losch;Annika Fuchs;Jean-François Lemieux;Anna Vanselow

  • Affiliations:

  • Alfred-Wegener-Institut, Helmholtz Zentrum für Polar- und Meeresforschung, Postfach 120161, 27515 Bremerhaven, Germany;Alfred-Wegener-Institut, Helmholtz Zentrum für Polar- und Meeresforschung, Postfach 120161, 27515 Bremerhaven, Germany;Recherche en Prévision Numérique environnementale/Environnement Canada, 2121 route Transcanadienne, Dorval, QC, H9P 1J3, Canada;Universität Oldenburg, Ammerländer Heerstr. 114-118, 26129 Oldenburg, Germany

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2014

Quantified Score

Hi-index 31.45

Visualization

Abstract

The most common representation of sea ice dynamics in climate models assumes that sea ice is a quasi-continuous non-normal fluid with a viscous-plastic rheology. This rheology leads to non-linear sea ice momentum equations that are notoriously difficult to solve. Recently a Jacobian-free Newton-Krylov (JFNK) solver was shown to solve the equations accurately at moderate costs. This solver is extended for massive parallel architectures and vector computers and implemented in a coupled sea ice-ocean general circulation model for climate studies. Numerical performance is discussed along with numerical difficulties in realistic applications with up to 1920 CPUs. The parallel JFNK-solver@?s scalability competes with traditional solvers although the collective communication overhead starts to show a little earlier. When accuracy of the solution is required (i.e. reduction of the residual norm of the momentum equations of more that one or two orders of magnitude) the JFNK-solver is unrivalled in efficiency. The new implementation opens up the opportunity to explore physical mechanisms in the context of large scale sea ice models and climate models and to clearly differentiate these physical effects from numerical artifacts.