Fixing EVP — Solving the Momentum Equations of Dynamic Sea Ice Models

Most dynamic sea ice models for climate type simulations are based on the viscous-plastic (VP) rheology. The most common methods of solving the momentum equations with VP-rheologies and elliptic yields curves include implicit methods with different flavors of Picard (or fixed point iterative) solver...

Full description

Bibliographic Details
Main Author: Losch, Martin
Format: Conference Object
Language:unknown
Published: 2018
Subjects:
Online Access:https://epic.awi.de/id/eprint/48415/
https://hdl.handle.net/10013/epic.13f66618-ef3d-4727-85fa-7009fc377630
Description
Summary:Most dynamic sea ice models for climate type simulations are based on the viscous-plastic (VP) rheology. The most common methods of solving the momentum equations with VP-rheologies and elliptic yields curves include implicit methods with different flavors of Picard (or fixed point iterative) solvers, Newton methods, and different variants of the explicit Elastic-Viscous-Plastic (EVP) solver. The EVP solver, in particular, is now used in many Earth system models and forced ice-ocean models, because it is easier to implement than implicit methods and generally cheaper on massively parallel computer. In its original form, however, it leads to noisy solutions that do not converge to the solution of the momentum equations. Simple and easy to implement modifications recover the full momentum equations. The improved EVP-solutions correspond to those obtained with implicit methods, still at comparatively low cost. At very high resolution, the effect of the correction is particularly striking.