Hostname: page-component-586b7cd67f-g8jcs Total loading time: 0 Render date: 2024-11-27T03:31:48.007Z Has data issue: false hasContentIssue false

Numerical approximation of viscous contact problems applied to glacial sliding

Published online by Cambridge University Press:  15 March 2022

Gonzalo G. de Diego*
Affiliation:
Mathematical Institute, University of Oxford, Oxford OX2 6GG, UK
Patrick E. Farrell
Affiliation:
Mathematical Institute, University of Oxford, Oxford OX2 6GG, UK
Ian J. Hewitt
Affiliation:
Mathematical Institute, University of Oxford, Oxford OX2 6GG, UK
*
Email address for correspondence: gonzalezdedi@maths.ox.ac.uk

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 equations. 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, then the cavity evolves towards the upsloping section, indicating that the downsloping part is unstable.

Type
JFM Papers
Copyright
© The Author(s), 2022. Published by Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

REFERENCES

Adams, R.A. & Fournier, J.J.F. 2003 Sobolev Spaces, Pure and Applied Mathematics, vol. 140. Academic Press, Elsevier.Google Scholar
Andrews, L.C., Catania, G.A., Hoffman, M.J., Gulley, J.D., Lüthi, M.P., Ryser, C., Hawley, R.L. & Neumann, T.A. 2014 Direct observations of evolving subglacial drainage beneath the Greenland Ice Sheet. Nature 514 (7520), 8083.CrossRefGoogle ScholarPubMed
Boffi, D., Brezzi, F. & Fortin, M. 2013 Mixed Finite Element Methods and Applications. Springer.CrossRefGoogle Scholar
de Diego, G.G. 2021 Finite element solver for subglacial cavities. https://doi.org/10.5281/zenodo.5887302.CrossRefGoogle Scholar
de Diego, G.G., Farrell, P.E. & Hewitt, I.J. 2021 On the finite element approximation of a semicoercive Stokes variational inequality arising in glaciology. arXiv:2108.00046.Google Scholar
Fowler, A.C. 1981 A theoretical treatment of the sliding of glaciers in the absence of cavitation. Phil. Trans. R. Soc. A 298 (1445), 637681.Google Scholar
Fowler, A.C. 1986 A sliding law for glaciers of constant viscosity in the presence of subglacial cavitation. Proc. R. Soc. A 407 (1832), 147170.Google Scholar
Fowler, A.C. 1987 Sliding with cavity formation. J. Glaciol. 33 (115), 255267.CrossRefGoogle Scholar
Gagliardini, O., Cohen, D., Raback, P. & Zwinger, T. 2007 Finite-element modeling of subglacial cavities and related friction law. J. Geophys. Res. 112 (F2), F02027.Google Scholar
Gagliardini, O. & Zwinger, T. 2008 The ISMIP-HOM benchmark experiments performed using the finite-element code Elmer. Cryosphere 2 (1), 6776.CrossRefGoogle Scholar
Glen, J.W. 1958 The flow law of ice: a discussion of the assumptions made in glacier theory, their experimental foundation and consequences. In Physics of the Movement of Ice: Symposium at Chamonix 1958, pp. 171–183. International Association of Hydrological Sciences.Google Scholar
Glowinski, R., Lions, J.L. & Trémoliéres, R. 1981 Numerical Analysis of Variational Inequalities. North Holland.Google Scholar
Gudmundsson, G.H. 1997 Basal-flow characteristics of a non-linear flow sliding frictionless over strongly undulating bedrock. J. Glaciol. 43 (143), 8089.CrossRefGoogle Scholar
Helanow, C., Iverson, N.R., Woodard, J.B. & Zoet, L.K. 2021 A slip law for hard-bedded glaciers derived from observed bed topography. Sci. Adv. 7 (20), eabe7798.CrossRefGoogle ScholarPubMed
Helanow, C., Iverson, N.R., Zoet, L.K. & Gagliardini, O. 2020 Sliding relations for glacier slip with cavities over three-dimensional beds. Geophys. Res. Lett. 47 (3), e2019GL084924.CrossRefGoogle Scholar
Hintermüller, M., Ito, K. & Kunisch, K. 2002 The primal-dual active set strategy as a semismooth Newton method. SIAM J. Optim. 13 (3), 865888.CrossRefGoogle Scholar
Hoffman, M.J., Andrews, L.C., Price, S.F., Catania, G.A., Neumann, T.A., Lüthi, M.P., Gulley, J., Ryser, C., Hawley, R.L. & Morriss, B. 2016 Greenland subglacial drainage evolution regulated by weakly connected regions of the bed. Nat. Commun. 7 (1), 13903.CrossRefGoogle Scholar
Iken, A. 1981 The effect of the subglacial water pressure on the sliding velocity of a glacier in an idealized numerical model. J. Glaciol. 27 (97), 407421.CrossRefGoogle Scholar
Iken, A. & Bindschadler, R.A. 1986 Combined measurements of subglacial water pressure and surface velocity of Findelengletscher, Switzerland: conclusions about drainage system and sliding mechanism. J. Glaciol. 32 (110), 101119.CrossRefGoogle Scholar
Kamb, B. 1970 Sliding motion of glaciers: theory and observation. Rev. Geophys. 8 (4), 673728.CrossRefGoogle Scholar
LeVeque, R.J. 2007 Finite Difference Methods for Ordinary and Partial Differential Equations: Steady-State and Time-Dependent Problems. SIAM.CrossRefGoogle Scholar
Mair, D.W.F., Sharp, M.J. & Willis, I.C. 2002 Evidence for basal cavity opening from analysis of surface uplift during a high-velocity event: Haut Glacier d'Arolla, Switzerland. J. Glaciol. 48 (161), 208216.CrossRefGoogle Scholar
Rathgeber, F., Ham, D.A., Mitchell, L., Lange, M., Luporini, F., McRae, A.T.T., Bercea, G.-T., Markall, G.R. & Kelly, P.H.J. 2016 Firedrake: automating the finite element method by composing abstractions. ACM Trans. Math. Softw. 43 (3), 127.CrossRefGoogle Scholar
Ritz, C., Edwards, T.L., Durand, G., Payne, A.J., Peyaud, V. & Hindmarsh, R.C.A. 2015 Potential sea-level rise from Antarctic ice-sheet instability constrained by observations. Nature 528 (7580), 115118.CrossRefGoogle ScholarPubMed
Schoof, C. 2005 The effect of cavitation on glacier sliding. Proc. R. Soc. A 461 (2055), 609627.CrossRefGoogle Scholar
Schoof, C. 2007 Ice sheet grounding line dynamics: steady states, stability, and hysteresis. J Geophys. Res. 112 (F3), F03S28.Google Scholar
Schoof, C. 2011 Marine ice sheet dynamics. Part 2. A Stokes flow contact problem. J. Fluid Mech. 679, 122155.CrossRefGoogle Scholar
Schoof, C. & Hewitt, I.J. 2013 Ice-sheet dynamics. Annu. Rev. Fluid Mech. 45 (1), 217239.CrossRefGoogle Scholar
Stadler, G. 2007 Path-following and augmented Lagrangian methods for contact problems in linear elasticity. J. Comput. Appl. Maths 203 (2), 533547.CrossRefGoogle Scholar
Stubblefield, A.G., Spiegelman, M. & Creyts, T.T. 2021 Variational formulation of marine ice-sheet and subglacial-lake grounding-line dynamics. J. Fluid Mech. 919, A23.CrossRefGoogle Scholar
Sugiyama, S. & Gudmundsson, G.H. 2004 Short-term variations in glacier flow controlled by subglacial water pressure at Lauteraargletscher, Bernese Alps, Switzerland. J. Glaciol. 50 (170), 353362.CrossRefGoogle Scholar
Weertman, J. 1957 On the sliding of glaciers. J. Glaciol. 3 (21), 3338.CrossRefGoogle Scholar
zenodo/Firedrake-20211103.0 2021 Software used in ‘Numerical approximation of viscous contact problems applied to glacial sliding’. https://doi.org/10.5281/zenodo.5643646.CrossRefGoogle Scholar
Zoet, L.K. & Iverson, N.R. 2015 Experimental determination of a double-valued drag relationship for glacier sliding. J. Glaciol. 61 (225), 17.CrossRefGoogle Scholar