On the convergence of the modified elastic-viscous-plastic method for solving the sea ice momentum equation

Most dynamic sea ice models for climate type simulations are based on the viscous-plastic (VP) rheology. The resulting stiff system of partial differential equations for ice velocity is either solved implicitly at great computational cost, or explicitly with added pseudo-elasticity (elastic-viscous-plastic, EVP). A recent modification of the EVP approach seeks to improve the convergence of the EVP method by re-interpreting it as a pseudotime VP solver. The question of convergence of this modified EVP method is revisited here and it is shown that convergence is reached provided the stability requirements are satisfied and the number of pseudotime iterations is sufficiently high. Only in this limit, the VP and the modified EVP solvers converge to the same solution. Related questions of the impact of mesh resolution and incomplete convergence are also addressed.

[1]  Jean-François Lemieux,et al.  Numerical convergence of viscous‐plastic sea ice models , 2009 .

[2]  Sylvain Bouillon,et al.  The elastic–viscous–plastic method revisited , 2013 .

[3]  D. Menemenlis,et al.  On the formulation of sea-ice models. Part 1: Effects of different solver implementations and parameterizations , 2010 .

[4]  Paul F. Tupper,et al.  Improving the numerical convergence of viscous-plastic sea ice models with the Jacobian-free Newton-Krylov method , 2010, J. Comput. Phys..

[5]  William D. Hibler,et al.  On an efficient numerical method for modeling sea ice dynamics , 1997 .

[6]  W. Hibler.  A Dynamic Thermodynamic Sea Ice Model , 1979 .

[7]  Jean-François Lemieux,et al.  A parallel Jacobian-free Newton-Krylov solver for a coupled sea ice-ocean model , 2014, J. Comput. Phys..

[8]  J. Schröter,et al.  Ocean circulation and sea ice distribution in a finite element global sea ice–ocean model , 2009 .

[9]  J. Peraire,et al.  Finite Element Flux-Corrected Transport (FEM-FCT) for the Euler and Navier-Stokes equations , 1987 .

[10]  E. Hunke,et al.  An Elastic–Viscous–Plastic Model for Sea Ice Dynamics , 1996 .

[11]  Jens Schröter,et al.  Finite-Element Sea Ice Model (FESIM), version 2 , 2015 .

[12]  Sergey Danilov,et al.  On Solving the Momentum Equations of Dynamic Sea , 2011 .

[13]  Elizabeth C. Hunke,et al.  Viscous–Plastic Sea Ice Dynamics with the EVP Model: Linearization Issues , 2001 .

[14]  Jens Schröter,et al.  The Finite Element Sea Ice-Ocean Model (FESOM) v.1.4: formulation of an ocean general circulation model , 2014 .

[15]  David M. Holland,et al.  A comparison of the Jacobian-free Newton-Krylov method and the EVP model for solving the sea ice momentum equation with a viscous-plastic formulation: A serial algorithm study , 2012, J. Comput. Phys..