期刊
COMMUNICATIONS IN NONLINEAR SCIENCE AND NUMERICAL SIMULATION
卷 22, 期 1-3, 页码 526-544出版社
ELSEVIER
DOI: 10.1016/j.cnsns.2014.08.007
关键词
Ice sheet models; Free boundaries; Enthalpy formulation; Finite elements; Duality methods
类别
资金
- Research Project of MCII [MTM2010-21135-C02]
- Research Project of MINECO [MTM2013-47800-C2]
In this paper we propose and numerically solve an original enthalpy formulation for the problem governing the thermal behaviour of polythermal ice sheets. Although the modelling follows some ideas introduced in Aschwanden and Blatter (2009), nonlinear basal boundary conditions in both cold and temperate regions are also considered, thus including the sliding effects in the frame of a fully coupled shallow ice approximation (SIA) model. One of the main novelties of this work comes from the introduction of the Heaviside multivalued operator to take into account the discontinuity of the thermal diffusion function at the cold-temperate transition surface (CTS) free boundary. Moreover, we propose a duality method for maximal monotone operators to solve simultaneously the nonlinear diffusive term and the free boundary. Some numerical simulation examples with real data from Antarctica are presented and illustrate the small differences between the computed results from the enthalpy formulation here proposed and the alternative formulation in terms of the temperature (Calvo et al., 2001). (C) 2014 Elsevier B.V. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据