Open Source Finite Element Software for Ice Sheet, Glaciers and Ice Flow Modelling

Elmer/Ice is a full-Stokes, finite element, ice sheet / ice flow model. The aim of this website is to present the capabilities of Elmer/Ice and to distribute course materials and tutorials.
Elmer/Ice is an add-on package to Elmer, which is a multi-physics FEM suite mainly developed by CSC-IT Center for Science Ltd., Espoo, Finland. Initially started by CSC,  IGE and ILTS, currently multiple institutions and individuals contribute to the development of Elmer/Ice.

Best of both worlds: SIA and full-Stokes

In a recent article in the Journal of Computational Physics, Ahlkrona et al. present a method (ISCAL - Ice  Sheet  Coupled Approximation Levels) implemented in Elmer/Ice that uses an error estimate to identify areas where the shallow ice approximation (SIA) can and cannot be deployed. According to this estimate, those regions, where the shallowness assumption prevails can be switched to a - compared to the full-Stokes solution - extremely cheap SIA (simply by replacing the matrix entries in the system matrix of the linearized and discretized full-Stokes problem. In their article the authors demonstrate a reduction of CPU time by ISCAL on an in size reduced real world case of up to 75% compared to a full-Stokes simulation run of the same problem.

Ahlkrona, J., P. Lötstedt, N. Kirchner, and T. Zwinger, 2016. Dynamically coupling the non-linear Stokes equations with the shallow ice approximation in glaciology: Description and first applications of the ISCAL method. J. Comp. Phys., 308, 1-19, doi:10.1016/j.jcp.2015.12.025.

Elmer/Ice project © 2018 -- Conception : iGrafic