Search

Abstract: In this paper we develop a set of numerical techniques for the simulation of the profile evolution of a valley glacier in the framework of isothermal shallow ice approximation models. The different mathematical formulations are given in terms of a highly nonlinear parabolic equation. A first nonlinearity comes from the free boundary problem associated with the unknown basal extension of the glacier region. This feature is treated using a fixed domain complementarity formulation which is solved numerically by a duality method. The nonlinear diffusive term is explicitly treated in the time marching scheme. A convection dominated problem arises, so a characteristic scheme is proposed for the time discretization, while piecewise linear finite elements are used for the spatial discretization. The presence of infinite slopes in polar regimes motivates an alternative formulation based on a prescribed flux boundary condition at the head of the glacier instead a homogeneous Dirichlet one. Finally, several numerical examples illustrate the performance of the proposed methods. [ABSTRACT FROM AUTHOR]

This work deals with the numerical solution of a complex mathematical model arising in theoretical glaciology. The global moving boundary problem governs thermomechanical processes jointly with ice sheet hydrodynamics. One major novelty is the inclusion of the ice velocity field computation in the framework of the shallow ice model so that it can be coupled with profile and temperature equations. Moreover, the proposed basal velocity and shear stress laws allow the integration of basal sliding effects in the global model. Both features were not taking into account in a previous paper (Math. Model. Methods Appl. Sci. 12 (2) (2002) 229) and provide more realistic convective terms and more complete Signorini boundary conditions for the thermal problem. In the proposed numerical algorithm, one- and two-dimensional piecewise linear Lagrange finite elements in space and a semi-implicit upwinding scheme in time are combined with duality and Newton''s methods for nonlinearities. A simulation example involving real data issued from Antarctic shows the temperature, profile and velocity qualitative behaviour as well as the free boundaries and basal effects. [Copyright &y& Elsevier]