Search

Abstract: In this paper the numerical solution of a highly nonlinear model for the thermomechanical behavior of polythermal glaciers is presented. The modeling follows the shallow ice approximation (SIA) for glaciers introduced in Fowler (1997) . The model has been extended to incorporate additional moving boundaries and other nonlinear features. Moreover, a fixed domain formulation is proposed to avoid the computational drawbacks of a time-dependent domain in the numerical simulation with front tracking methods. In this setting, the coupled problem is decomposed into different nonlinear problems which allow one to obtain sequentially the profile evolution, the velocity field, the glacier surface and atmospheric temperatures, basal magnitudes and the temperature distribution inside the ice mass. A fixed point iteration algorithm converges to the solution of the nonlinear coupled problem. Among different numerical methods involved in the solution of the subproblems, characteristic schemes for time discretization, finite elements for spatial discretization, duality methods for the nonlinearities associated to maximal monotone operators and a Newton scheme for the nonlinear viscous term are proposed. Several numerical simulation examples illustrate the performance of the numerical methods and the behavior of the involved physical magnitudes. [Copyright &y& Elsevier]

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]

In this work we propose a new numerical method for solving a thermomechanical coupled problem which arises as a mathematical model for the evolution of ice sheets profiles and the corresponding temperatures. From the accumulation-ablation ratio, the atmospheric temperature and the geothermic flux, the ice sheet profile is the solution of a moving boundary problem governed by a nonlinear convection-diffusion equation. Moreover, the mathematical model which governs the temperature distribution involves three nonlinearities: a reaction term due to the viscous dissipation, a Signorini boundary condition associated to the geothermic flux and an enthalpy term issued from the two-phase Stefan formulation of the polythermal regime. The problems are discretized in time with upwind characteristics schemes and in space with finite elements. The nonlinearities are solved either by fixed point methods or by duality techniques. Several numerical simulation examples involving real data sets issued from the Antartic ice sheet are shown. [ABSTRACT FROM AUTHOR]