Implementation of finite difference method based on explicit and Crank-Nicolson method to the thermomechanical ice sheet model

Thermomechanical ice sheet model is applied to simulate ice sheet behavior in the Antarctic region. The parameters involved are ice thickness, ice temperature and ice velocity. This model is proposed by Bueler (2012)[1]. It corresponding roughly to the two-dimensional European Ice Sheet Modeling Ini...

Full description

Bibliographic Details
Main Authors: Kasmin-Bajuri, N. K. B., Alias, N.
Format: Conference Object
Language:unknown
Published: Institute of Electrical and Electronics Engineers Inc. 2016
Subjects:
Online Access:http://eprints.utm.my/73015/
https://www.scopus.com/inward/record.uri?eid=2-s2.0-84995592275&doi=10.1109%2fISMSC.2015.7594091&partnerID=40&md5=7433a4369e460c60958837eda7608805
Description
Summary:Thermomechanical ice sheet model is applied to simulate ice sheet behavior in the Antarctic region. The parameters involved are ice thickness, ice temperature and ice velocity. This model is proposed by Bueler (2012)[1]. It corresponding roughly to the two-dimensional European Ice Sheet Modeling Initiative (EISMENT) experiment and the thermomechanical coupling equation. Thus, the expression of thermomechanical coupling equations is quite simplified. The numerical discretization scheme for the model is explicit and Crank-Nicholson methods. The efficiency of the finite difference methods such as implicit, explicit and Crank-Nicholson methods for solving the partial differential equation (PDE) of thermomechanical ice sheet modeling is also determined. MATLAB has been choosing as the development platform for the implementations since it is well suited for the kind of computations required.