Ice melting in a turbulent stratified shear flow

<p>In this talk I will present preliminary results of direct numerical simulations of ice melting in a turbulent stratified shear flow. The model solves the evolution of the turbulent fluid phase and of the diffusive solid ice phase, due to melting and freezing, in a fully coupled way. This is done by combining a Direct Numerical Simulation (DNS) code with a novel formulation of the equations for the solid and liquid phases of water based on the phase-field method. DNS enables turbulent motions to be simulated without approximation, i.e. solving Navier Stokes equations, while the phase-field method allows the ice-ocean interface to be rough and evolve in response to melting. I will present results on the turbulent boundary layer and on the self-generated basal topography at the ice-water interface. The ultimate goal of this work is to propose a new DNS-based parameterization of the melting process at rough ice-ocean boundaries that takes into account the effects of temperature and salt stratification, and flow velocities.</p>

Stratified shear instability in a field of pre-existing turbulence

Turbulent mixing of heat and momentum in the stably-stratified ocean interior occurs in discrete events driven by vertical variations of the horizontal velocity. Typically, these events have been modelled assuming an initially laminar stratified shear flow which develops wavelike instabilities, becomes fully turbulent, and then relaminarizes into a stable state. However, in the real ocean there is always some level of turbulence left over from previous events. Using direct numerical simulations, we show that the evolution of a stably-stratified shear layer may be significantly modified by pre-existing turbulence. The classical billow structure associated with Kelvin–Helmholtz instability is suppressed and eventually eliminated as the strength of the initial turbulence is increased. A corresponding energetics analysis shows that potential energy changes and dissipation of kinetic energy depend non-monotonically on initial turbulence strength, with the largest effects when initial turbulence is present but insufficient to prevent billow formation. The mixing efficiency decreases with increasing initial turbulence amplitude as the development of the Kelvin–Helmholtz billow, with its large pre-turbulent mixing efficiency, is arrested.

Instability of sheared density interfaces

Of the canonical flow instabilities (Kelvin–Helmholtz, Holmboe-wave and Taylor–Caulfield) of stratified shear flow, the Taylor–Caulfield instability (TCI) has received relatively little attention, and forms the focus of the current study. First, a diagnostic of the linear instability dynamics is developed that exploits the net pseudomomentum to distinguish TCI from the other two instabilities for any given flow profile. Second, the nonlinear dynamics of TCI is studied across its range of unstable horizontal wavenumbers and bulk Richardson numbers using numerical simulation. At small bulk Richardson numbers, a cascade of billow structures of sequentially smaller size may form. For large bulk Richardson numbers, the primary nonlinear travelling waves formed by the linear instability break down via a small-scale, Kelvin–Helmholtz-like roll-up mechanism with an associated large amount of mixing. In all cases, secondary parasitic nonlinear Holmboe waves appear at late times for high Prandtl number. Third, a nonlinear diagnostic is proposed to distinguish between the saturated states of the three canonical instabilities based on their distinctive density–streamfunction and generalised vorticity–streamfunction relations.

The structure and origin of confined Holmboe waves

Finite-amplitude manifestations of stratified shear flow instabilities and their spatio-temporal coherent structures are believed to play an important role in turbulent geophysical flows. Such shear flows commonly have layers separated by sharp density interfaces, and are therefore susceptible to the so-called Holmboe instability, and its finite-amplitude manifestation, the Holmboe wave. In this paper, we describe and elucidate the origin of an apparently previously unreported long-lived coherent structure in a sustained stratified shear flow generated in the laboratory by exchange flow through an inclined square duct connecting two reservoirs filled with fluids of different densities. Using a novel measurement technique allowing for time-resolved, near-instantaneous measurements of the three-component velocity and density fields simultaneously over a three-dimensional volume, we describe the three-dimensional geometry and spatio-temporal dynamics of this structure. We identify it as a finite-amplitude, nonlinear, asymmetric confined Holmboe wave (CHW), and highlight the importance of its spanwise (lateral) confinement by the duct boundaries. We pay particular attention to the spanwise vorticity, which exhibits a travelling, near-periodic structure of sheared, distorted, prolate spheroids with a wide ‘body’ and a narrower ‘head’. Using temporal linear stability analysis on the two-dimensional streamwise-averaged experimental flow, we solve for three-dimensional perturbations having two-dimensional, cross-sectionally confined eigenfunctions and a streamwise normal mode. We show that the dispersion relation and the three-dimensional spatial structure of the fastest-growing confined Holmboe instability are in good agreement with those of the observed confined Holmboe wave. We also compare those results with a classical linear analysis of two-dimensional perturbations (i.e. with no spanwise dependence) on a one-dimensional base flow. We conclude that the lateral confinement is an important ingredient of the confined Holmboe instability, which gives rise to the CHW, with implications for many inherently confined geophysical flows such as in valleys, estuaries, straits or deep ocean trenches. Our results suggest that the CHW is an example of an experimentally observed, inherently nonlinear, robust, long-lived coherent structure which has developed from a linear instability. We conjecture that the CHW is a promising candidate for a class of exact coherent states underpinning the dynamics of more disordered, yet continually forced stratified shear flows.