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


Contact
Martin.Losch [ at ] awi.de

Abstract

Numerical convergence properties of a recently developed Jacobian-free Newton-Krylov (JFNK) solver are compared to the ones of the widely used EVP model when solving the sea ice momentum equation with a Viscous-Plastic (VP) formulation. To do so, very accurate reference solutions are produced with an independent Picard solver with an advective time step of 10 s and a tight nonlinear convergence criterion on 10, 20, 40, and 80 km grids. Approximate solutions with the JFNK and EVP solvers are obtained for advective time steps of 10, 20 and 30 min. Because of an artificial elastic term, the EVP model permits an explicit time-stepping scheme with a relatively large subcycling time step. The elastic waves excited during the subcycling are intended to damp out and almost entirely disappear such that the approximate solution should be close to the VP solution. Results show that residual elastic waves cause the EVP approximate solution to have notable differences with the reference solution and that these differences get more important as the grid is refined. Compared to the reference solution, additional shear lines and zones of strong convergence/divergence are seen in the EVP approximate solution. The approximate solution obtained with the JFNK solver is very close to the reference solution for all spatial resolutions tested.



Item Type
Article
Authors
Divisions
Programs
Publication Status
Published
Eprint ID
30470
DOI 10.1016/j.jcp.2012.05.024

Cite as
Lemieux, J. F. , Knoll, D. , Tremblay, B. , Holland, D. M. and Losch, M. (2012): 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 , Journal of Computational Physics, 231 (17), pp. 5926-5944 . doi: 10.1016/j.jcp.2012.05.024


Download
[thumbnail of JCOMP-D-12-00289R1_textonly.pdf]
Preview
PDF
JCOMP-D-12-00289R1_textonly.pdf - Accepted Version

Download (3MB) | Preview
Cite this document as:

Share
Add to AnyAdd to TwitterAdd to FacebookAdd to LinkedinAdd to PinterestAdd to Email


Citation

Research Platforms
N/A

Campaigns
N/A


Actions
Edit Item Edit Item