Ice and Fluid Dynamics

Abstract:

Thermodynamic ablation of ice in contact with the ocean is an essential element of ice sheet and ocean interactions but is challenging to model and quantify. Building on earlier observations of sea ice ablation, a variety of recent theoretical, experimental and observational studies have considered ice ablation in contrasting geometries, from vertical to near-horizontal ice faces, and reveal different scaling behaviour for predicted ablation rates in different dynamical regimes. However, uncertainties remain about when the contrasting results should be applied, as existing model parameterisations do not capture all relevant regimes of ice–ocean ablation. To progress towards improved models of ice–ocean​ interaction, we synthesise current understanding into a classification of ablation types. We examine the effect of the classification on the parameterisation of turbulent fluxes from the ocean towards the ice, and identify the dominant processes next to ice interfaces of different orientation. Four ablation types are defined: melting and dissolving based on ocean temperatures, and shear-controlled and buoyancy-controlled regimes based on the dynamics of the near-ice molecular sublayer. We describe existing observational and modelling studies of sea ice, ice shelves, and glacier termini, as well as laboratory studies, to show how they fit into this classification. Two sets of observations from the Ross and Ronne Ice Shelf cavities suggest that both the buoyancy-controlled and shear-controlled regimes may be relevant under different oceanographic conditions. Overall, buoyancy-controlled dynamics are more likely when the molecular sublayer has lower Reynolds number, and shear for higher Reynolds number, although the observations suggest some variability about this trend.

Abstract:

Marine-terminating glaciers, such as those along the coastline of Greenland, often release meltwater into the ocean in the form of subglacial discharge plumes. Though these plumes can dramatically alter the mass loss along the front of a glacier, the conditions surrounding their genesis remain poorly constrained. In particular, little is known about the geometry of subglacial outlets and the extent to which seawater may intrude into them. Here, the latter is addressed by exploring the dynamics of an arrested salt wedge – a steady-state, two-layer flow system where salty water partially intrudes a channel carrying fresh water. Building on existing theory, we formulate a model that predicts the length of a non-entraining salt wedge as a function of the Froude number, the slope of the channel and coefficients for interfacial and wall drag. In conjunction, a series of laboratory experiments were conducted to observe a salt wedge within a rectangular channel. For experiments conducted with laminar flow (Reynolds number Re < 800), good agreement with theoretical predictions are obtained when the drag coefficients are modelled as being inversely proportional to Re. However, for fully turbulent flows on geophysical scales, these drag coefficients are expected to asymptote toward finite values. Adopting reasonable drag coefficient estimates for this flow regime, our theoretical model suggests that typical subglacial channels may permit seawater intrusions of the order of several kilometres. While crude, these results indicate that the ocean has a strong tendency to penetrate subglacial channels and potentially undercut the face of marine-terminating glaciers.

Abstract:

© 2019 The solidification of a binary alloy results in the formation of a porous mushy layer, within which spontaneous localisation of fluid flow can lead to the emergence of features over a range of spatial scales. We describe a finite volume method for simulating binary alloy solidification in two dimensions with local mesh refinement in space and time. The coupled heat, solute, and mass transport is described using an enthalpy method with flow described by a Darcy-Brinkman equation for flow across porous and liquid regions. The resulting equations are solved on a hierarchy of block-structured adaptive grids. A projection method is used to compute the fluid velocity, whilst the viscous and nonlinear diffusive terms are calculated using a semi-implicit scheme. A series of synchronization steps ensure that the scheme is flux-conservative and correct for errors that arise at the boundaries between different levels of refinement. We also develop a corresponding method using Darcy's law for flow in a porous medium/narrow Hele-Shaw cell. We demonstrate the accuracy and efficiency of our method using established benchmarks for solidification without flow and convection in a fixed porous medium, along with convergence tests for the fully coupled code. Finally, we demonstrate the ability of our method to simulate transient mushy layer growth with narrow liquid channels which evolve over time.