Dynamically coupling the non-linear Stokes equations with the shallow ice approximation in glaciology: description and first applications of the ISCAL method

DOI10.1016/J.JCP.2015.12.025MaRDI QIDQ729439zbMATH OpenOpenAlexWikidataFDO

Authors Josefin Ahlkrona, Nina Kirchner, Thomas Zwinger, Per Lötstedt

Publication date 20 December 2016

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1503.07743

zbMATH Keywords

error estimationfinite element methodshallow ice approximationfull Stokes equationsice sheet modeling

Mathematics Subject Classification ID

Non-Newtonian fluids (76A05) Stokes and related (Oseen, etc.) flows (76D07) Finite element methods applied to problems in fluid mechanics (76M10) Computational methods for problems pertaining to geophysics (86-08) Glaciology (86A40)

Abstract: We propose and implement a new method, called the Ice Sheet Coupled Approximation Levels (ISCAL) method, for simulation of ice sheet flow in large domains under long time-intervals. The method couples the exact, full Stokes (FS) equations with the Shallow Ice Approximation (SIA). The part of the domain where SIA is applied is determined automatically and dynamically based on estimates of the modeling error. For a three dimensional model problem where the number of degrees of freedom is comparable to a real world application, ISCAL performs almost an order of magnitude faster with a low reduction in accuracy compared to a monolithic FS. Furthermore, ISCAL is shown to be able to detect rapid dynamic changes in the flow. Three different error estimations are applied and compared. Finally, ISCAL is applied to the Greenland Ice Sheet, proving ISCAL to be a potential valuable tool for the ice sheet modeling community.

Recommendations

• A new fully nonisothermal coupled model for the simulation of ice sheet flow
• Finite elements numerical solution of a coupled profile--velocity--temperature shallow ice sheet approximation model.
• A fluidity-based first-order system least-squares method for ice sheets
• Coupled models and parallel simulations for three-dimensional full-Stokes ice sheet modeling
• A multilayer ice-flow model generalising the shallow shelf approximation

Cites work

• scientific article; zbMATH DE number 780776 (Why is no real title available?)
• A computational framework for infinite-dimensional Bayesian inverse problems. II: stochastic Newton MCMC with application to ice sheet flow inverse problems
• A numerical study of scaling relations for non-Newtonian thin-film flows with applications in ice sheet modelling
• Achieving textbook multigrid efficiency for hydrostatic ice sheet flow
• Parallel iterative methods for sparse linear systems
• Thin-film flows with wall slip: an asymptotic analysis of higher order glacier flow models

Cited in

(7)

• A meshfree approach to non-Newtonian free surface ice flow: application to the Haut Glacier d'Arolla
• Mechanical error estimators for shallow ice flow models
• A multilayer ice-flow model generalising the shallow shelf approximation
• Accurate and stable time stepping in ice sheet modeling
• Selection of a method for solving nonlinear equations in shallow-water icing model implementation
• Anisotropic radial basis function methods for continental size ice sheet simulations
• A fluidity-based first-order system least-squares method for ice sheets

Describes a project that uses

Uses Software

• Elmer/Ice
• ISSM

This page was built for publication: Dynamically coupling the non-linear Stokes equations with the shallow ice approximation in glaciology: description and first applications of the ISCAL method

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q729439)