Numerical approximation of viscous contact problems applied to glacial sliding

TL;DR

This paper introduces a new numerical method for viscous contact problems, crucial in glaciology, enabling accurate simulation of cavity evolution and stability analysis of sliding laws.

Contribution

A novel mixed formulation with Lagrange multipliers and an upwinding scheme for robustly solving viscous contact problems in glaciology.

Findings

Method accurately reproduces analytical results.

Reveals instability in certain friction law regions.

Simulates cavity evolution under oscillating pressures.

Abstract

Viscous contact problems describe the time evolution of fluid flows in contact with a surface from which they can detach and reattach. These problems are of particular importance in glaciology, where they arise in the study of grounding lines and subglacial cavities. In this work, we propose a novel numerical method for solving viscous contact problems based on a mixed formulation with Lagrange multipliers of a variational inequality involving the Stokes equation. The advection equation for evolving the geometry of the domain occupied by the fluid is then solved via a specially-built upwinding scheme, leading to a robust and accurate algorithm for viscous contact problems. We first verify the method by comparing the numerical results to analytical results obtained by a linearised method. Then, we use this numerical scheme to reconstruct friction laws for glacial sliding with cavitation. Finally, we compute the evolution of cavities from a steady state under oscillating water pressures. The results depend strongly on the location of the initial steady state along the friction law. In particular, we find that if the steady state is located on the downsloping or rate-weakening part of the friction law, the cavity evolves towards the upsloping section, indicating that the downsloping part is unstable.