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

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–Kryl...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics 2014-01, Vol.257, p.901-911
Main Authors: Losch, Martin, Fuchs, Annika, Lemieux, Jean-François, Vanselow, Anna
Format: Article
Language:English
Subjects:
Climate models
Jacobian-free Newton–Krylov solver
Mathematical analysis
Mathematical models
Nonlinear dynamics
Numerical sea ice modeling
Parallel implementation
Preconditioning
Reduction
Rheology
Sea ice
Sea ice dynamics
Solvers
Vector implementation
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary: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.
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2013.09.026