The main goal of this research is to model the three-dimensional motion of a glacier. Given the shape of the glacier, a nonlinear Stokes problem has to be solved. The shape of the glacier is then solved using a lagrangian Volume of Fluid formulation.
This project is a collaboration with Heinz Blatter Institute for Atmospheric and Climate Science, ETHZ and With Martin Funk Wasserbau/Hydrologie/Glaziologie, ETHZ.
Ice is a very viscous fluid. The simulation of ice flow can therefore be achieved with cfsFlow, in order to determine the evolution of alpine glaciers. Equations of fluid dynamics take into account the dynamics of the ice flow on the bed rock, as well as the sliding of the ice on the mountain rock. The model used in cfsFlow incorporates state-of-the-art models for the ice physics and rheology, and efficient and accurate numerical algorithms for the simulation of alpine glaciers.
The behavior of alpine glaciers greatly depends on the underlying climate model, in particular on radiation data, temperatures, precipitations, and bed rock topology. cfsFlow is able to incorporate in a flexible manner various mass balance models to determine ice melt and snow addition (accumulation and ablation).
Starting with the bed rock topography in terms of altitudes of grid points, and from measurements for the ice surface at a given time, one can construct an approximation of the ice mass. We define the physical laws governing the ice flow and determine if the glacier slides on and sticks to the bedrock. The mass balance definition alows to incorporate many climatic effects in existing mass balance models, through the definition of grids. Mass balance models with various complexity levels can be used.
Example : The Rhone glacier evolution in a median climatic scenario from 2007 to 2100.