Numerical methods for a fixed domain formulation of the glacier profile problem with alternative boundary conditions

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:235

Article

Calvo, N.²  Durany, J.²  Toja, R.¹  Vazquez, C.¹ 

[1] Univ A Coruna, Dep Matemat, La Coruna 15071, Spain

[2] Univ Vigo, ETSI Telecomunicac, Dep Matemat Apilcada 2, Vigo 36310, Spain

关键词: Glacier models;    Shallow ice approximation;    Moving boundaries;    Complementarity problems;    Duality methods;    Characteristics scheme;    Finite elements;   

DOI  :  10.1016/j.cam.2010.08.026

来源: Elsevier

【 摘 要 】

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. (C) 2010 Elsevier B.V. All rights reserved.

【 授权许可】

Free