FEniCS-full-Stokes ...
This Python 3 source code solves the full-Stokes equations with different solvers and step size controls. We solve the full-Stokes equations with different algorithms and test these with the experiment ISMIP-HOM B, see (Pattyn et al.; Benchmark experiments for higher-order and full-Stokes ice sheet...
Main Author: | |
---|---|
Format: | Article in Journal/Newspaper |
Language: | unknown |
Published: |
Zenodo
2024
|
Subjects: | |
Online Access: | https://dx.doi.org/10.5281/zenodo.7805478 https://zenodo.org/doi/10.5281/zenodo.7805478 |
id |
ftdatacite:10.5281/zenodo.7805478 |
---|---|
record_format |
openpolar |
spelling |
ftdatacite:10.5281/zenodo.7805478 2024-03-31T07:53:21+00:00 FEniCS-full-Stokes ... Schmidt, Niko 2024 https://dx.doi.org/10.5281/zenodo.7805478 https://zenodo.org/doi/10.5281/zenodo.7805478 unknown Zenodo https://dx.doi.org/10.5281/zenodo.8279222 https://dx.doi.org/10.5281/zenodo.7805479 https://dx.doi.org/10.5281/zenodo.8154332 https://dx.doi.org/10.5281/zenodo.10618661 https://dx.doi.org/10.5281/zenodo.10470996 Creative Commons Attribution 4.0 International https://creativecommons.org/licenses/by/4.0/legalcode cc-by-4.0 SoftwareSourceCode article Software 2024 ftdatacite https://doi.org/10.5281/zenodo.780547810.5281/zenodo.827922210.5281/zenodo.780547910.5281/zenodo.815433210.5281/zenodo.1061866110.5281/zenodo.10470996 2024-03-04T12:06:06Z This Python 3 source code solves the full-Stokes equations with different solvers and step size controls. We solve the full-Stokes equations with different algorithms and test these with the experiment ISMIP-HOM B, see (Pattyn et al.; Benchmark experiments for higher-order and full-Stokes ice sheet models (ISMIP-HOM; 2008; The Cryosphere). The program allows one to choose between Armijo step sizes, exact step sizes, and constant step sizes; the Picard iteration and the Newton iteration; the functional and the residual norm as a minimization term; and an experiment in two and three dimensions. A combination of exact step sizes with the residual norm as the minimization term is not possible as the exact step sizes rely on a convex function as the minimization term. An executable example with comments is in examples/run_fullStokes.py. In this example are written comments on how to switch between algorithms. This source code relies on FEniCS https://fenicsproject.org/download/archive/ version 2019.1.0. FEniCS ... : New in this version: For time-dependent problems, we fixed a mistake for velocities vx<0 in the discretization scheme. Note that the commented out lines 17, 56-58 in picard.py and 16, 50-52 in newton.py are used for the stopping criteria that the residual norm reduced enough compared to the initial guess. They are not used for the experiments in which the number of iterations is fixed. ... Article in Journal/Newspaper Ice Sheet DataCite Metadata Store (German National Library of Science and Technology) |
institution |
Open Polar |
collection |
DataCite Metadata Store (German National Library of Science and Technology) |
op_collection_id |
ftdatacite |
language |
unknown |
description |
This Python 3 source code solves the full-Stokes equations with different solvers and step size controls. We solve the full-Stokes equations with different algorithms and test these with the experiment ISMIP-HOM B, see (Pattyn et al.; Benchmark experiments for higher-order and full-Stokes ice sheet models (ISMIP-HOM; 2008; The Cryosphere). The program allows one to choose between Armijo step sizes, exact step sizes, and constant step sizes; the Picard iteration and the Newton iteration; the functional and the residual norm as a minimization term; and an experiment in two and three dimensions. A combination of exact step sizes with the residual norm as the minimization term is not possible as the exact step sizes rely on a convex function as the minimization term. An executable example with comments is in examples/run_fullStokes.py. In this example are written comments on how to switch between algorithms. This source code relies on FEniCS https://fenicsproject.org/download/archive/ version 2019.1.0. FEniCS ... : New in this version: For time-dependent problems, we fixed a mistake for velocities vx<0 in the discretization scheme. Note that the commented out lines 17, 56-58 in picard.py and 16, 50-52 in newton.py are used for the stopping criteria that the residual norm reduced enough compared to the initial guess. They are not used for the experiments in which the number of iterations is fixed. ... |
format |
Article in Journal/Newspaper |
author |
Schmidt, Niko |
spellingShingle |
Schmidt, Niko FEniCS-full-Stokes ... |
author_facet |
Schmidt, Niko |
author_sort |
Schmidt, Niko |
title |
FEniCS-full-Stokes ... |
title_short |
FEniCS-full-Stokes ... |
title_full |
FEniCS-full-Stokes ... |
title_fullStr |
FEniCS-full-Stokes ... |
title_full_unstemmed |
FEniCS-full-Stokes ... |
title_sort |
fenics-full-stokes ... |
publisher |
Zenodo |
publishDate |
2024 |
url |
https://dx.doi.org/10.5281/zenodo.7805478 https://zenodo.org/doi/10.5281/zenodo.7805478 |
genre |
Ice Sheet |
genre_facet |
Ice Sheet |
op_relation |
https://dx.doi.org/10.5281/zenodo.8279222 https://dx.doi.org/10.5281/zenodo.7805479 https://dx.doi.org/10.5281/zenodo.8154332 https://dx.doi.org/10.5281/zenodo.10618661 https://dx.doi.org/10.5281/zenodo.10470996 |
op_rights |
Creative Commons Attribution 4.0 International https://creativecommons.org/licenses/by/4.0/legalcode cc-by-4.0 |
op_doi |
https://doi.org/10.5281/zenodo.780547810.5281/zenodo.827922210.5281/zenodo.780547910.5281/zenodo.815433210.5281/zenodo.1061866110.5281/zenodo.10470996 |
_version_ |
1795032964866768896 |