Several problems of great importance in the study of glaciers and ice sheets involve processes of attachment and reattachment of the ice from the bedrock. Consider, for example, an ice sheet sliding from the continent into the ocean, where it goes afloat. Another example is that of subglacial cavitation, a fundamental mechanism in glacial sliding where the ice detaches from the bedrock along the downstream area of an obstacle. Such problems are generally modelled as a viscous Stokes flow with a free boundary and contact boundary conditions. In this talk, I will present a framework for solving such problems numerically. I will start by introducing the mathematical formulation of these viscous contact problems and the challenges that arise when trying to approximate them numerically. I will then show how, given a stable scheme for the free boundary equation, one can build a penalty formulation for the viscous contact problem in such a way that the resulting algorithm remains stable and robust.
