Numerical ice sheet modelling using FEniCS
The increase in computational power in the last 20 years has made it feasible to solve the full Stokes equation for ice flow, using finite element methods (FEM). However, the numerical properties of these equations remains largely unknown. This is due to their nonlinear nature, which makes them hard...
Main Author: | |
---|---|
Format: | Bachelor Thesis |
Language: | English |
Published: |
Uppsala universitet, Institutionen för informationsteknologi
2019
|
Subjects: | |
Online Access: | http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-396571 |
id |
ftuppsalauniv:oai:DiVA.org:uu-396571 |
---|---|
record_format |
openpolar |
spelling |
ftuppsalauniv:oai:DiVA.org:uu-396571 2023-05-15T16:40:39+02:00 Numerical ice sheet modelling using FEniCS Löfgren, André 2019 application/pdf http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-396571 eng eng Uppsala universitet, Institutionen för informationsteknologi IT 19039 http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-396571 info:eu-repo/semantics/openAccess Engineering and Technology Teknik och teknologier Student thesis info:eu-repo/semantics/bachelorThesis text 2019 ftuppsalauniv 2023-02-23T21:51:03Z The increase in computational power in the last 20 years has made it feasible to solve the full Stokes equation for ice flow, using finite element methods (FEM). However, the numerical properties of these equations remains largely unknown. This is due to their nonlinear nature, which makes them hard to analyze mathematically. For this reason convergence rate and stability has to be established by performing simulations. In this thesis a 2D ice solver has been developed. The solver was then tested on the ISMIP-HOM benchmarks, in order to assert convergence rates and stability. The solver was developed using the FEM software FEniCS. It was verified that the solver could obtain the same convergence rate as in the case of linear Stokes; however, for curved boundaries the convergence rate in velocity dropped by one order. It was also found that discontinuity in the coefficient of basal traction reduced convergence rate to linear in both pressure and velocity. Time dependence was added to the solver by coupling the Stoke's system to the kinematic free surface equation. Time dependence was then tested using one of the EISMINT benchmark. It was found that for a mesh with a mesh size parameter of 10 km, instabilities arised after 2500 years when using a time step of 25 years, resulting in spurious numerical oscillation. However, by introducing a diffusive term into the surface equation, it was possible to attenuate these oscillation without affecting the overall shape of the ice sheet. Bachelor Thesis Ice Sheet Uppsala University: Publications (DiVA) |
institution |
Open Polar |
collection |
Uppsala University: Publications (DiVA) |
op_collection_id |
ftuppsalauniv |
language |
English |
topic |
Engineering and Technology Teknik och teknologier |
spellingShingle |
Engineering and Technology Teknik och teknologier Löfgren, André Numerical ice sheet modelling using FEniCS |
topic_facet |
Engineering and Technology Teknik och teknologier |
description |
The increase in computational power in the last 20 years has made it feasible to solve the full Stokes equation for ice flow, using finite element methods (FEM). However, the numerical properties of these equations remains largely unknown. This is due to their nonlinear nature, which makes them hard to analyze mathematically. For this reason convergence rate and stability has to be established by performing simulations. In this thesis a 2D ice solver has been developed. The solver was then tested on the ISMIP-HOM benchmarks, in order to assert convergence rates and stability. The solver was developed using the FEM software FEniCS. It was verified that the solver could obtain the same convergence rate as in the case of linear Stokes; however, for curved boundaries the convergence rate in velocity dropped by one order. It was also found that discontinuity in the coefficient of basal traction reduced convergence rate to linear in both pressure and velocity. Time dependence was added to the solver by coupling the Stoke's system to the kinematic free surface equation. Time dependence was then tested using one of the EISMINT benchmark. It was found that for a mesh with a mesh size parameter of 10 km, instabilities arised after 2500 years when using a time step of 25 years, resulting in spurious numerical oscillation. However, by introducing a diffusive term into the surface equation, it was possible to attenuate these oscillation without affecting the overall shape of the ice sheet. |
format |
Bachelor Thesis |
author |
Löfgren, André |
author_facet |
Löfgren, André |
author_sort |
Löfgren, André |
title |
Numerical ice sheet modelling using FEniCS |
title_short |
Numerical ice sheet modelling using FEniCS |
title_full |
Numerical ice sheet modelling using FEniCS |
title_fullStr |
Numerical ice sheet modelling using FEniCS |
title_full_unstemmed |
Numerical ice sheet modelling using FEniCS |
title_sort |
numerical ice sheet modelling using fenics |
publisher |
Uppsala universitet, Institutionen för informationsteknologi |
publishDate |
2019 |
url |
http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-396571 |
genre |
Ice Sheet |
genre_facet |
Ice Sheet |
op_relation |
IT 19039 http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-396571 |
op_rights |
info:eu-repo/semantics/openAccess |
_version_ |
1766031051930468352 |