FEniCS-full-Stokes ...

• Main Author: Schmidt, Niko
• Format: Article in Journal/Newspaper
• Language: unknown
• Published: Zenodo 2023
• Subjects: Ice Sheet
• Online Access: https://dx.doi.org/10.5281/zenodo.8154332; https://zenodo.org/record/8154332

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: - Removed commented out code - Added post-processing Python scripts - Added dataset from simulation ...