A strongly nonlinear problem arising in glaciology
The computation of glacier movements leads to a system of nonlinear partial differential equations. The existence and uniqueness of a weak solution is established by using the calculus of variations. A discretization by the finite element method is done. The solution of the discrete problem is proved to be convergent to the exact solution. A first simple numerical algorithm is proposed and its convergence numerically studied.
The main goal of this article is to establish a priori and a posteriori error estimates for the numerical approximation of some non linear elliptic problems arising in glaciology. The stationary motion of a glacier is given by a non-newtonian fluid flow model which becomes, in a first two-dimensional approximation, the so-called infinite parallel sided slab model. The approximation of this model is made by a finite element method with piecewise polynomial functions of degree 1. Numerical results...
The main goal of this article is to establish a priori and a posteriori error estimates for the numerical approximation of some non linear elliptic problems arising in glaciology. The stationary motion of a glacier is given by a non-Newtonian fluid flow model which becomes, in a first two-dimensional approximation, the so-called infinite parallel sided slab model. The approximation of this model is made by a finite element method with piecewise polynomial functions of degree 1. Numerical results...
In this paper we study a model problem describing the movement of a glacier under Glen’s flow law and investigated by Colinge and Rappaz [Colinge and Rappaz, ESAIM: M2AN 33 (1999) 395–406]. We establish error estimates for finite element approximation using the results of Chow [Chow, SIAM J. Numer. Analysis 29 (1992) 769–780] and Liu and Barrett [Liu and Barrett, SIAM J. Numer. Analysis 33 (1996) 98–106] and give an analysis of the convergence of the successive approximations used in [Colinge and...
In this paper we study a model problem describing the movement of a glacier under Glen's flow law and investigated by Colinge and Rappaz [Colinge and Rappaz, ESAIM: M2AN33 (1999) 395–406]. We establish error estimates for finite element approximation using the results of Chow [Chow, SIAM J. Numer. Analysis29 (1992) 769–780] and Liu and Barrett [Liu and Barrett, SIAM J. Numer. Analysis33 (1996) 98–106] and give an analysis of the convergence of the successive approximations used in [Colinge and...
En estas notas presentamos algunos modelos físicos que han sido propuestos recientemente para tratar el problema de los movimientos repentinos y casi periódicos del hielo, así como la aparición de corrientes de hielo rápidas en los grandes mantos glaciares que se deslizan sobre lechos blandos y deformables. Estos fenómenos están relacionados con la transición de un régimen de flujo lento a uno rápido y pueden aparecer debido a una modificación del sistema de drenaje del glaciar. Los fenómenos en...
The motion of a three-dimensional glacier is considered. Ice is modeled as an incompressible non-Newtonian fluid. At each time step, given the shape of the glacier, a nonlinear elliptic system has to be solved in order to obtain the two components of the horizontal velocity field. Then, the shape of the glacier is updated by solving a transport equation. Finite element techniques are used to compute the velocity field and to solve the transport equation. Numerical results are compared to experiments...
This paper discusses finite element discretization and preconditioning strategies for the iterative solution of nonsymmetric indefinite linear algebraic systems of equations arising in modelling of glacial rebound processes. Some numerical experiments for the purely elastic model setting are provided. Comparisons of the performance of the iterative solution method with a direct solution method are included as well.
The linear heat equation predicts that the variations of temperature along a cold ice sheet {i.e. at a temperature less than is freezing point) due to a sudden increase in air temperature, are very very slow. Based on this we represent the nonlinear evolution of an ice sheet as a sequence of steady states. As a first fundamental indication that this model is correct well posedness with respect to the variations of initial and boundary data is proved. Further an estimate of the error made in evaluating...