Improved convergence and stability properties in a three-dimensional higher-order ice sheet model
We present a finite difference implementation of a three-dimensional higher-order ice sheet model. In comparison to a conventional centred difference discretisation it enhances both numerical stability and convergence. In order to achieve these benefits the discretisation of the governing force bala...
Published in: | Geoscientific Model Development |
---|---|
Main Authors: | , , , , , |
Format: | Other/Unknown Material |
Language: | English |
Published: |
2018
|
Subjects: | |
Online Access: | https://doi.org/10.5194/gmd-4-1133-2011 https://gmd.copernicus.org/articles/4/1133/2011/ |
id |
ftcopernicus:oai:publications.copernicus.org:gmd11914 |
---|---|
record_format |
openpolar |
spelling |
ftcopernicus:oai:publications.copernicus.org:gmd11914 2023-05-15T16:40:35+02:00 Improved convergence and stability properties in a three-dimensional higher-order ice sheet model Fürst, J. J. Rybak, O. Goelzer, H. Smedt, B. Groen, P. Huybrechts, P. 2018-09-27 info:eu-repo/semantics/application/pdf https://doi.org/10.5194/gmd-4-1133-2011 https://gmd.copernicus.org/articles/4/1133/2011/ eng eng info:eu-repo/grantAgreement/EC/FP7/226375 doi:10.5194/gmd-4-1133-2011 https://gmd.copernicus.org/articles/4/1133/2011/ info:eu-repo/semantics/openAccess eISSN: 1991-9603 info:eu-repo/semantics/Text 2018 ftcopernicus https://doi.org/10.5194/gmd-4-1133-2011 2020-07-20T16:25:57Z We present a finite difference implementation of a three-dimensional higher-order ice sheet model. In comparison to a conventional centred difference discretisation it enhances both numerical stability and convergence. In order to achieve these benefits the discretisation of the governing force balance equation makes extensive use of information on staggered grid points. Using the same iterative solver, a centred difference discretisation that operates exclusively on the regular grid serves as a reference. The reprise of the ISMIP-HOM experiments indicates that both discretisations are capable of reproducing the higher-order model inter-comparison results. This setup allows a direct comparison of the two numerical implementations also with respect to their convergence behaviour. First and foremost, the new finite difference scheme facilitates convergence by a factor of up to 7 and 2.6 in average. In addition to this decrease in computational costs, the accuracy for the resultant velocity field can be chosen higher in the novel finite difference implementation. Changing the discretisation also prevents build-up of local field irregularites that occasionally cause divergence of the solution for the reference discretisation. The improved behaviour makes the new discretisation more reliable for extensive application to real ice geometries. Higher accuracy and robust numerics are crucial in time dependent applications since numerical oscillations in the velocity field of subsequent time steps are attenuated and divergence of the solution is prevented. Other/Unknown Material Ice Sheet Copernicus Publications: E-Journals Geoscientific Model Development 4 4 1133 1149 |
institution |
Open Polar |
collection |
Copernicus Publications: E-Journals |
op_collection_id |
ftcopernicus |
language |
English |
description |
We present a finite difference implementation of a three-dimensional higher-order ice sheet model. In comparison to a conventional centred difference discretisation it enhances both numerical stability and convergence. In order to achieve these benefits the discretisation of the governing force balance equation makes extensive use of information on staggered grid points. Using the same iterative solver, a centred difference discretisation that operates exclusively on the regular grid serves as a reference. The reprise of the ISMIP-HOM experiments indicates that both discretisations are capable of reproducing the higher-order model inter-comparison results. This setup allows a direct comparison of the two numerical implementations also with respect to their convergence behaviour. First and foremost, the new finite difference scheme facilitates convergence by a factor of up to 7 and 2.6 in average. In addition to this decrease in computational costs, the accuracy for the resultant velocity field can be chosen higher in the novel finite difference implementation. Changing the discretisation also prevents build-up of local field irregularites that occasionally cause divergence of the solution for the reference discretisation. The improved behaviour makes the new discretisation more reliable for extensive application to real ice geometries. Higher accuracy and robust numerics are crucial in time dependent applications since numerical oscillations in the velocity field of subsequent time steps are attenuated and divergence of the solution is prevented. |
format |
Other/Unknown Material |
author |
Fürst, J. J. Rybak, O. Goelzer, H. Smedt, B. Groen, P. Huybrechts, P. |
spellingShingle |
Fürst, J. J. Rybak, O. Goelzer, H. Smedt, B. Groen, P. Huybrechts, P. Improved convergence and stability properties in a three-dimensional higher-order ice sheet model |
author_facet |
Fürst, J. J. Rybak, O. Goelzer, H. Smedt, B. Groen, P. Huybrechts, P. |
author_sort |
Fürst, J. J. |
title |
Improved convergence and stability properties in a three-dimensional higher-order ice sheet model |
title_short |
Improved convergence and stability properties in a three-dimensional higher-order ice sheet model |
title_full |
Improved convergence and stability properties in a three-dimensional higher-order ice sheet model |
title_fullStr |
Improved convergence and stability properties in a three-dimensional higher-order ice sheet model |
title_full_unstemmed |
Improved convergence and stability properties in a three-dimensional higher-order ice sheet model |
title_sort |
improved convergence and stability properties in a three-dimensional higher-order ice sheet model |
publishDate |
2018 |
url |
https://doi.org/10.5194/gmd-4-1133-2011 https://gmd.copernicus.org/articles/4/1133/2011/ |
genre |
Ice Sheet |
genre_facet |
Ice Sheet |
op_source |
eISSN: 1991-9603 |
op_relation |
info:eu-repo/grantAgreement/EC/FP7/226375 doi:10.5194/gmd-4-1133-2011 https://gmd.copernicus.org/articles/4/1133/2011/ |
op_rights |
info:eu-repo/semantics/openAccess |
op_doi |
https://doi.org/10.5194/gmd-4-1133-2011 |
container_title |
Geoscientific Model Development |
container_volume |
4 |
container_issue |
4 |
container_start_page |
1133 |
op_container_end_page |
1149 |
_version_ |
1766030995020054528 |