Exploring basal sliding with a fluidity‐based, ice‐sheet model using FOSLS
Summary This paper develops two first‐order system—in this context, first‐order refers to the order of the PDE not to the model—least‐squares, fluidity‐based formulations of a nonlinear Stokes flow model for ice sheets that attempt to overcome the difficulties introduced by unbounded viscosity. One...
Published in: | Numerical Linear Algebra with Applications |
---|---|
Main Authors: | , , |
Other Authors: | , |
Format: | Article in Journal/Newspaper |
Language: | English |
Published: |
Wiley
2018
|
Subjects: | |
Online Access: | http://dx.doi.org/10.1002/nla.2161 https://api.wiley.com/onlinelibrary/tdm/v1/articles/10.1002%2Fnla.2161 https://onlinelibrary.wiley.com/doi/pdf/10.1002/nla.2161 https://onlinelibrary.wiley.com/doi/full-xml/10.1002/nla.2161 https://onlinelibrary.wiley.com/doi/am-pdf/10.1002/nla.2161 |
id |
crwiley:10.1002/nla.2161 |
---|---|
record_format |
openpolar |
spelling |
crwiley:10.1002/nla.2161 2024-06-02T08:08:20+00:00 Exploring basal sliding with a fluidity‐based, ice‐sheet model using FOSLS Allen, Jeffery Manteuffel, Tom Rajaram, Harihar U.S. Department of Energy Lawrence Livermore National Laboratory 2018 http://dx.doi.org/10.1002/nla.2161 https://api.wiley.com/onlinelibrary/tdm/v1/articles/10.1002%2Fnla.2161 https://onlinelibrary.wiley.com/doi/pdf/10.1002/nla.2161 https://onlinelibrary.wiley.com/doi/full-xml/10.1002/nla.2161 https://onlinelibrary.wiley.com/doi/am-pdf/10.1002/nla.2161 en eng Wiley http://onlinelibrary.wiley.com/termsAndConditions#am http://onlinelibrary.wiley.com/termsAndConditions#vor Numerical Linear Algebra with Applications volume 25, issue 3 ISSN 1070-5325 1099-1506 journal-article 2018 crwiley https://doi.org/10.1002/nla.2161 2024-05-03T11:10:31Z Summary This paper develops two first‐order system—in this context, first‐order refers to the order of the PDE not to the model—least‐squares, fluidity‐based formulations of a nonlinear Stokes flow model for ice sheets that attempt to overcome the difficulties introduced by unbounded viscosity. One commonly used way to define viscosity, Glen's law, allows viscosity to become unbounded as the strain rates approach zero. Often, numerical approaches overcome these singularities by modifying viscosity to limit its maximum. The formulations in this paper, however, reframe the problem to avoid viscosity altogether by defining the system in terms of the inverse of viscosity, which is known as fluidity. This results in a quantity that approaches zero as viscosity approaches infinity. Additionally, a set of equations that represent the curl of the velocity gradient is added to help approximate the solution in a space closer to H 1 , which improves algebraic multigrid convergence. Previous research revealed that the first‐order system least‐squares formulation has difficulties in maintaining optimal discretization convergence on more complex domains. This paper discovers that this problem is linked to how the curl equations are scaled and that stronger scalings result in better solver performance but worse discretization convergence. Determining if there is an optimal scaling that balances performance and convergence is still an open question. Additionally, the fluidity‐based formulations are tested with three 2D benchmark problems. Two of these benchmark problems involve basal sliding and one involves a time‐dependent free surface. The fluidity‐based solutions are consistent with the standard Galerkin method using Taylor–hood elements while better resolving viscosity. Article in Journal/Newspaper Ice Sheet Wiley Online Library Curl ENVELOPE(-63.071,-63.071,-70.797,-70.797) Numerical Linear Algebra with Applications 25 3 |
institution |
Open Polar |
collection |
Wiley Online Library |
op_collection_id |
crwiley |
language |
English |
description |
Summary This paper develops two first‐order system—in this context, first‐order refers to the order of the PDE not to the model—least‐squares, fluidity‐based formulations of a nonlinear Stokes flow model for ice sheets that attempt to overcome the difficulties introduced by unbounded viscosity. One commonly used way to define viscosity, Glen's law, allows viscosity to become unbounded as the strain rates approach zero. Often, numerical approaches overcome these singularities by modifying viscosity to limit its maximum. The formulations in this paper, however, reframe the problem to avoid viscosity altogether by defining the system in terms of the inverse of viscosity, which is known as fluidity. This results in a quantity that approaches zero as viscosity approaches infinity. Additionally, a set of equations that represent the curl of the velocity gradient is added to help approximate the solution in a space closer to H 1 , which improves algebraic multigrid convergence. Previous research revealed that the first‐order system least‐squares formulation has difficulties in maintaining optimal discretization convergence on more complex domains. This paper discovers that this problem is linked to how the curl equations are scaled and that stronger scalings result in better solver performance but worse discretization convergence. Determining if there is an optimal scaling that balances performance and convergence is still an open question. Additionally, the fluidity‐based formulations are tested with three 2D benchmark problems. Two of these benchmark problems involve basal sliding and one involves a time‐dependent free surface. The fluidity‐based solutions are consistent with the standard Galerkin method using Taylor–hood elements while better resolving viscosity. |
author2 |
U.S. Department of Energy Lawrence Livermore National Laboratory |
format |
Article in Journal/Newspaper |
author |
Allen, Jeffery Manteuffel, Tom Rajaram, Harihar |
spellingShingle |
Allen, Jeffery Manteuffel, Tom Rajaram, Harihar Exploring basal sliding with a fluidity‐based, ice‐sheet model using FOSLS |
author_facet |
Allen, Jeffery Manteuffel, Tom Rajaram, Harihar |
author_sort |
Allen, Jeffery |
title |
Exploring basal sliding with a fluidity‐based, ice‐sheet model using FOSLS |
title_short |
Exploring basal sliding with a fluidity‐based, ice‐sheet model using FOSLS |
title_full |
Exploring basal sliding with a fluidity‐based, ice‐sheet model using FOSLS |
title_fullStr |
Exploring basal sliding with a fluidity‐based, ice‐sheet model using FOSLS |
title_full_unstemmed |
Exploring basal sliding with a fluidity‐based, ice‐sheet model using FOSLS |
title_sort |
exploring basal sliding with a fluidity‐based, ice‐sheet model using fosls |
publisher |
Wiley |
publishDate |
2018 |
url |
http://dx.doi.org/10.1002/nla.2161 https://api.wiley.com/onlinelibrary/tdm/v1/articles/10.1002%2Fnla.2161 https://onlinelibrary.wiley.com/doi/pdf/10.1002/nla.2161 https://onlinelibrary.wiley.com/doi/full-xml/10.1002/nla.2161 https://onlinelibrary.wiley.com/doi/am-pdf/10.1002/nla.2161 |
long_lat |
ENVELOPE(-63.071,-63.071,-70.797,-70.797) |
geographic |
Curl |
geographic_facet |
Curl |
genre |
Ice Sheet |
genre_facet |
Ice Sheet |
op_source |
Numerical Linear Algebra with Applications volume 25, issue 3 ISSN 1070-5325 1099-1506 |
op_rights |
http://onlinelibrary.wiley.com/termsAndConditions#am http://onlinelibrary.wiley.com/termsAndConditions#vor |
op_doi |
https://doi.org/10.1002/nla.2161 |
container_title |
Numerical Linear Algebra with Applications |
container_volume |
25 |
container_issue |
3 |
_version_ |
1800753546024452096 |