Application of MPRGP algorithm for ice sheet melting problems

Abstract

Melting of glaciers and ice sheets poses several ecosystem threats, making the accurate modelling of their thermodynamic behaviour crucial. When simulating these dynamics, boundary conditions often involve inequality constraints that can be formulated as quadratic programming (QP) problems. The Modified Proportioning with Reduced Gradient Projections algorithm (MPRGP) is highly effective for solving QP problems.

The thesis focuses on applying MPRGP to problems arising in glaciology, specifically on solving the heat equation with a simple temperature upper bound inequality constraint. The main goal is to develop an interface between the PERMON toolbox, utilising its QP solver methods, and the Elmer/Ice software used for glacier modelling. The proposed solution is then benchmarked on a thermo-mechanically coupled steady state simulation for the Midtre Lovénbreen glacier in Svalbard.

Tání ledovců a ledových příkrovů představuje významnou hrozbu pro ekosystémy, a proto je přesné modelování jejich termodynamického chování zcela zásadní. Při simulacích často vznikají okrajové podmínky s nerovnostími omezeními, které lze formulovat jako úlohy kvadratického programování (QP). K řešení těchto úloh je velmi efektivní algoritmus Modified Proportioning with Reduced Gradient Projections (MPRGP).

Tato práce se zaměřuje na aplikaci algoritmu MPRGP na problémy v glaciologii, konkrétně na řešení rovnice vedení tepla s jednoduchým horním nerovnostním omezením teploty. Hlavním cílem je vývoj rozhraní mezi knihovnou PERMON a softwarem Elmer/Ice (určeným pro modelování ledovců) a využití metod PERMON knihovny pro řešení QP úloh. Navržené řešení je následně testováno na termomechanicky sdružené simulaci ustáleného stavu ledovce Midtre Lovénbreen na Špicberkách.