Anisotropic Radial Basis Function Methods for Continental Size Ice Sheet Simulations

• Main Authors: Cheng, Gong, Shcherbakov, Victor
• Format: Text
• Language: unknown
• Published: arXiv 2017
• Subjects: Numerical Analysis math.NA; Fluid Dynamics physics.flu-dyn; FOS Mathematics; FOS Physical sciences; 65N15, 65N35, 76D07; Ice Sheet
• Online Access: https://dx.doi.org/10.48550/arxiv.1711.09947; https://arxiv.org/abs/1711.09947

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. : The authors contributed equally to this work