Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Galerkin Least-Squares Stabilization in Ice Sheet Modeling - Accuracy, Robustness, and Comparison to other Techniques
Stockholm University, Faculty of Science, Department of Physical Geography. (Geomorphology and glaciology)ORCID iD: 0000-0003-4310-4873
(English)Manuscript (preprint) (Other academic)
Abstract [en]

We investigate the accuracy and robustness of one of the most common methods used in glaciology for the discretization of the p-Stokes equations: equal order finite elements with Galerkin Least-Squares (GLS) stabilization. Furthermore we compare the results to other stabilized methods. We find that the vertical velocity component is more sensitive to the choice of GLS stabilization parameter than horizontal velocity. Additionally, the accuracy of the vertical velocity component is especially important since errors in this component can cause ice surface instabilities and propagate into future ice volume predictions. If the element cell size is set to the minimum edge length and the stabilization parameter is allowed to vary non-linearly with viscosity, the GLS stabilization parameter found in literature is a good choice on simple domains. However, near ice margins the standard parameter choice may result in significant oscillations in the vertical component of the surface velocity. For these cases, other stabilization techniques, such as the interior penalty method, result in better accuracy and are less sensitive to the choice of the stabilization parameter. During this work we also discovered that the manufactured solutions often used to evaluate errors in glaciology are not reliable due to high artificial surface forces at singularities. We perform our numerical experiments in both FEniCS and Elmer/Ice.

Keyword [en]
finite element method, Galerkin Least-Squares, p-Stokes, ice sheet modeling, anisotropic mesh
National Category
Geosciences, Multidisciplinary
Research subject
Computing Science; Geography, Physical Geography
Identifiers
URN: urn:nbn:se:su:diva-141635OAI: oai:DiVA.org:su-141635DiVA: diva2:1087896
Note

This manuscript is also available as a preprint on arXiv.org at

http://arxiv.org/abs/1702.08369

with the arXiv identifier

arXiv:1702.08369

Available from: 2017-04-10 Created: 2017-04-10 Last updated: 2017-04-11
In thesis
1.
The record could not be found. The reason may be that the record is no longer available or you may have typed in a wrong id in the address field.

Open Access in DiVA

No full text

Search in DiVA

By author/editor
Helanow, Christian
By organisation
Department of Physical Geography
Geosciences, Multidisciplinary

Search outside of DiVA

GoogleGoogle Scholar

CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf