Anisotropic Radial Basis Function Methods for Continental Size Ice Sheet Simulations

• Published in: Journal of Computational Physics
• Main Authors: Cheng, Gong, Shcherbakov, Victor
• Format: Text
• Language: unknown
• Published: 2017
• Subjects: Mathematics - Numerical Analysis; Physics - Fluid Dynamics; 65N15; 65N35; 76D07; Ice Sheet
• Online Access: http://arxiv.org/abs/1711.09947; https://doi.org/10.1016/j.jcp.2018.06.020

In this paper we develop and implement anisotropic radial basis function methods for simulating the dynamics of ice sheets and glaciers. We test the methods on two problems: the well-known benchmark ISMIP-HOM B that corresponds to a glacier size ice and a synthetic ice sheet whose geometry is inspired by the EISMINT benchmark that corresponds to a continental size ice sheet. We illustrate the advantages of the radial basis function methods over a standard finite element method. We also show how the use of anisotropic radial basis functions allows for accurate simulation of the velocities on a large ice sheet, which was not possible with standard isotropic radial basis function methods due to a large aspect ratio between the ice length and the ice thickness. Additionally, we implement a partition of unity method in order to improve the computational efficiency of the radial basis function methods. Comment: The authors contributed equally to this work