Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • 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
Stabilized equal low-order finite elements in ice sheet modeling - accuracy and robustness
Stockholm University, Faculty of Science, Department of Physical Geography.ORCID iD: 0000-0003-4310-4873
Number of Authors: 22018 (English)In: Computational Geosciences, ISSN 1420-0597, E-ISSN 1573-1499, Vol. 22, no 4, p. 951-974Article in journal (Refereed) Published
Abstract [en]

We investigate the accuracy and robustness of one of the most common methods used in glaciology for finite element discretization of the oe-Stokes equations: linear equal order finite elements with Galerkin least-squares (GLS) stabilization on anisotropic meshes. 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 reasons, other stabilization techniques, in particular the interior penalty method, result in better accuracy and are less sensitive to the choice of 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.

Place, publisher, year, edition, pages
2018. Vol. 22, no 4, p. 951-974
Keywords [en]
Finite element method, Galerkin least-squares, p-Stokes, Ice-sheet modeling, Anisotropic mesh
National Category
Earth and Related Environmental Sciences
Identifiers
URN: urn:nbn:se:su:diva-159029DOI: 10.1007/s10596-017-9713-5ISI: 000438821700001OAI: oai:DiVA.org:su-159029DiVA, id: diva2:1245151
Available from: 2018-09-04 Created: 2018-09-04 Last updated: 2020-01-23Bibliographically approved

Open Access in DiVA

No full text in DiVA

Other links

Publisher's full text

Search in DiVA

By author/editor
Helanow, Christian
By organisation
Department of Physical Geography
In the same journal
Computational Geosciences
Earth and Related Environmental Sciences

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 13 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • 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