Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Numerical ice sheet modelling using FEniCS
Uppsala University, Disciplinary Domain of Science and Technology, Mathematics and Computer Science, Department of Information Technology.
2019 (English)Independent thesis Advanced level (degree of Master (Two Years)), 20 credits / 30 HE creditsStudent thesis
Abstract [en]

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.

Place, publisher, year, edition, pages
2019. , p. 50
Series
IT ; 19039
National Category
Engineering and Technology
Identifiers
URN: urn:nbn:se:uu:diva-396571OAI: oai:DiVA.org:uu-396571DiVA, id: diva2:1368307
Educational program
Master Programme in Computational Science
Supervisors
Examiners
Available from: 2019-11-06 Created: 2019-11-06 Last updated: 2019-11-06Bibliographically approved

Open Access in DiVA

fulltext(1747 kB)11 downloads
File information
File name FULLTEXT01.pdfFile size 1747 kBChecksum SHA-512
eeb82f19fe799a212b12a694ce238094b064149667d659867f792ea37b4bfd249a547b8dbc9556d50d217ce1abd1f800d8701533f8326c6d1336bec44e5e1535
Type fulltextMimetype application/pdf

By organisation
Department of Information Technology
Engineering and Technology

Search outside of DiVA

GoogleGoogle Scholar
Total: 11 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

urn-nbn

Altmetric score

urn-nbn
Total: 31 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf