Numerical approximation of viscous contact problems applied to glacial sliding
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...
| Main Authors: | , , |
|---|---|
| Format: | Journal article |
| Language: | English |
| Published: |
Cambridge University Press
2022
|
| _version_ | 1826308512429375488 |
|---|---|
| author | de Diego, GG Farrell, PE Hewitt, IJ |
| author_facet | de Diego, GG Farrell, PE Hewitt, IJ |
| author_sort | de Diego, GG |
| collection | OXFORD |
| description | 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. |
| first_indexed | 2024-03-07T07:20:35Z |
| format | Journal article |
| id | oxford-uuid:ea3a9fbf-0162-4516-8cb5-bc010719e172 |
| institution | University of Oxford |
| language | English |
| last_indexed | 2024-03-07T07:20:35Z |
| publishDate | 2022 |
| publisher | Cambridge University Press |
| record_format | dspace |
| spelling | oxford-uuid:ea3a9fbf-0162-4516-8cb5-bc010719e1722022-09-15T11:02:59ZNumerical approximation of viscous contact problems applied to glacial slidingJournal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:ea3a9fbf-0162-4516-8cb5-bc010719e172EnglishSymplectic ElementsCambridge University Press2022de Diego, GGFarrell, PEHewitt, IJViscous 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. |
| spellingShingle | de Diego, GG Farrell, PE Hewitt, IJ Numerical approximation of viscous contact problems applied to glacial sliding |
| title | Numerical approximation of viscous contact problems applied to glacial sliding |
| title_full | Numerical approximation of viscous contact problems applied to glacial sliding |
| title_fullStr | Numerical approximation of viscous contact problems applied to glacial sliding |
| title_full_unstemmed | Numerical approximation of viscous contact problems applied to glacial sliding |
| title_short | Numerical approximation of viscous contact problems applied to glacial sliding |
| title_sort | numerical approximation of viscous contact problems applied to glacial sliding |
| work_keys_str_mv | AT dediegogg numericalapproximationofviscouscontactproblemsappliedtoglacialsliding AT farrellpe numericalapproximationofviscouscontactproblemsappliedtoglacialsliding AT hewittij numericalapproximationofviscouscontactproblemsappliedtoglacialsliding |