Jouvet, G. and Gräser, C. (2013) An adaptive Newton multigrid method for a model of marine ice sheets. Journal of Computational Physics, 252 . pp. 419437. ISSN 00219991

PDF
1MB 
Official URL: http://dx.doi.org/10.1016/j.jcp.2013.06.032
Abstract
In this paper, we consider a model for the time evolution of threedimensional marine ice sheets. This model combines the Shallow Ice Approximation (SIA) for the ice deformation, the Shallow Shelf Approximation (SSA) for the basal sliding, and the mass conservation principle. At each time step, we solve a scalar pLaplace minimizationtype problem with obstacle (SIA), a vectorial pLaplace minimizationtype problem (SSA) and a transport equation (mass conservation). The two minimization problems are solved using a truncated nonsmooth Newton multigrid method while the transport equation is solved using a vertexcentred finite volume method. Our approach is combined to an heuristic mesh adaptive refinement procedure to face the large gradients of the solution that are expected between the ice sheet and the ice shelf. As applications, we present some simulations of the Marine Ice Sheet Model Intercomparison Project MISMIP (2D and 3D) and validate our results against an analytic solution (2D) and other participant model results (3D). Further numerical results show that the convergence of our Newton multigrid method is insensitive to local refinements making our overall adaptive strategy fully efficient.
Item Type:  Article 

Uncontrolled Keywords:  Marine ice sheets, Grounding line, pLaplace, Newton multigrid method, Adaptive method 
Subjects:  Mathematical and Computer Sciences > Mathematics > Numerical Analysis 
Divisions:  Department of Mathematics and Computer Science > Institute of Mathematics 
ID Code:  1794 
Deposited By:  Ekaterina Engel 
Deposited On:  17 Feb 2016 10:01 
Last Modified:  03 Mar 2017 14:41 
Repository Staff Only: item control page