We provide scaling relations for the Nusselt number $Nu$ and the friction coefficient $C_{S}$ in sheared Rayleigh–Bénard convection, i.e. in Rayleigh–Bénard flow with Couette- or Poiseuille-type shear forcing, by extending the Grossmann & Lohse (J. Fluid Mech., vol. 407, 2000, pp. 27–56, Phys. Rev. Lett., vol. 86, 2001, pp. 3316–3319, Phys. Rev. E, vol. 66, 2002, 016305, Phys. Fluids, vol. 16, 2004, pp. 4462–4472) theory to sheared thermal convection. The control parameters for these systems are the Rayleigh number $Ra$, the Prandtl number $Pr$ and the Reynolds number $Re_S$ that characterises the strength of the imposed shear. By direct numerical simulations and theoretical considerations, we show that, in turbulent Rayleigh–Bénard convection, the friction coefficients associated with the applied shear and the shear generated by the large-scale convection rolls are both well described by Prandtl's (Ergeb. Aerodyn. Vers. Gött., vol. 4, 1932, pp. 18–29) logarithmic friction law, suggesting some kind of universality between purely shear-driven flows and thermal convection. These scaling relations hold well for $10^6 \leq Ra \leq 10^8$, $0.5 \leq Pr \leq 5.0$, and $0 \leq Re_S \leq 10^4$.

Turbulent shear flows driven by a combination of a pressure gradient and buoyancy forcing are investigated using direct numerical simulations. Specifically, we consider the set-up of a differentially heated vertical channel subject to a Poiseuille-like horizontal pressure gradient. We explore the response of the system to its three control parameters: the Grashof number $Gr$, the Prandtl number $Pr$, and the Reynolds number $Re$ of the pressure-driven flow. From these input parameters, the relative strength of buoyancy driving to the pressure gradient can be quantified by the Richardson number $Ri=Gr/Re^2$. We compare the response of the mixed vertical convection configuration to that of mixed Rayleigh–Bénard convection, and find a nearly identical behaviour, including an increase in wall friction at higher $Gr$, and a drop in the heat flux relative to natural convection for $Ri=O(1)$. This closely matched response is despite vastly different flow structures in the systems. No large-scale organisation is visible in visualisations of mixed vertical convection – an observation that is confirmed quantitatively by spectral analysis. This analysis, combined with a statistical description of the wall heat flux, highlights how moderate shear suppresses the growth of small-scale plumes and reduces the likelihood of extreme events in the local wall heat flux. Vice versa, starting from a pure shear flow, the addition of thermal driving enhances the drag due to the emission of thermal plumes.

We consider the process of convective dissolution in a homogeneous and isotropic porous medium. The flow is unstable due to the presence of a solute that induces a density difference responsible for driving the flow. The mixing dynamics is thus driven by a Rayleigh–Taylor instability at the pore scale. We investigate the flow at the scale of the pores using Hele-Shaw type experiment with bead packs, two-dimensional direct numerical simulations and physical models. Experiments and simulations have been specifically designed to mimic the same flow conditions, namely matching porosities, high Schmidt numbers and linear dependency of fluid density with solute concentration. In addition, the solid obstacles of the medium are impermeable to fluid and solute. We characterise the evolution of the flow via the mixing length, which quantifies the extension of the mixing region and grows linearly in time. The flow structure, analysed via the centreline mean wavelength, is observed to grow in agreement with theoretical predictions. Finally, we analyse the dissolution dynamics of the system, quantified through the mean scalar dissipation, and three mixing regimes are observed. Initially, the evolution is controlled by diffusion, which produces solute mixing across the initial horizontal interface. Then, when the interfacial diffusive layer is sufficiently thick, it becomes unstable, forming finger-like structures and driving the system into a convection-dominated phase. Finally, when the fingers have grown sufficiently to touch the horizontal boundaries of the domain, the mixing reduces dramatically due to the absence of fresh unmixed fluid. With the aid of simple physical models, we explain the physics of the results obtained numerically and experimentally. The solute evolution presents a self-similar behaviour, and it is controlled by different length scales in each stage of the mixing process, namely the length scale of diffusion, the pore size and the domain height.

Iceberg melting is a critical factor for climate change. However, the shape of an iceberg is an often neglected aspect of its melting process. Our study investigates the influence of different ice shapes and ambient flow velocities on melt rates by conducting direct numerical simulations of a simplified system of bluff body flow. Our study focuses on the ellipsoidal shape, with the aspect ratio as the control parameter. We found the shape plays a crucial role in the melting process, resulting in significant variations in the melt rate between different shapes. Without flow, the optimal shape for a minimal melt rate is the disk (two-dimensional) or sphere (three-dimensional), due to the minimal surface area. However, as the ambient flow velocity increases, the optimal shape changes with the aspect ratio. We find that ice with an elliptical shape (when the long axis is aligned with the flow direction) can melt up to 10 % slower than a circular shape when exposed to flowing water. Following the approach considered by Huang et al. (J. Fluid Mech., vol. 765, 2015, R3) for dissolving bodies, we provide a quantitative theoretical explanation for this optimal shape, based on the combined contributions from both surface-area effects and convective-heat-transfer effects. Our findings provide insight into the interplay between phase transitions and ambient flows, contributing to our understanding of the iceberg melting process and highlighting the need to consider the aspect-ratio effect in estimates of iceberg melt rates.

The presence of salt in seawater strongly affects the melt rate and the shape evolution of ice, both of utmost relevance in ice–ocean interactions and thus for the climate. To get a better quantitative understanding of the physical mechanics at play in ice melting in salty water, we numerically investigate the lateral melting of an ice block in stably stratified saline water. The developing ice shape from our numerical results shows good agreement with the experiments and theory from Huppert & Turner (J. Fluid Mech., vol. 100, 1980, pp. 367–384). Furthermore, we find that the melt rate of ice depends non-monotonically on the mean ambient salinity: it first decreases for increasing salt concentration until a local minimum is attained, and then increases again. This non-monotonic behaviour of the ice melt rate is due to the competition among salinity-driven buoyancy, temperature-driven buoyancy and salinity-induced stratification. We develop a theoretical model based on the force balance which gives a prediction of the salt concentration for which the melt rate is minimal, and is consistent with our data. Our findings give insight into the interplay between phase transitions and double-diffusive convective flows.

We numerically study the melting process of a solid layer heated from below such that a liquid melt layer develops underneath. The objective is to quantitatively describe and understand the emerging topography of the structures (characterized by the amplitude and wavelength), which evolve out of an initially smooth surface. By performing both two-dimensional (achieving Rayleigh number up to $Ra=10^{11}$) and three-dimensional (achieving Rayleigh number up to $Ra=10^9$) direct numerical simulations with an advanced finite difference solver coupled to the phase-field method, we show how the interface roughness is spontaneously generated by thermal convection. With increasing height of the melt the convective flow intensifies and eventually even becomes turbulent. As a consequence, the interface becomes rougher but is still coupled to the flow topology. The emerging structure of the interface coincides with the regions of rising hot plumes and descending cold plumes. We find that the roughness amplitude scales with the mean height of the liquid layer. We derive this scaling relation from the Stefan boundary condition and relate it to the non-uniform distribution of heat flux at the interface. For two-dimensional cases, we further quantify the horizontal length scale of the morphology, based on the theoretical upper and lower bounds given for the size of convective cells known from Wang et al. (Phys. Rev. Lett., vol. 125, 2020, 074501). These bounds agree with our simulation results. Our findings highlight the key connection between the morphology of the melting solid and the convective flow structures in the melt beneath.

Many environmental flows arise due to natural convection at a vertical surface, from flows in buildings to dissolving ice faces at marine-terminating glaciers. We use three-dimensional direct numerical simulations of a vertical channel with differentially heated walls to investigate such convective, turbulent boundary layers. Through the implementation of a multiple-resolution technique, we are able to perform simulations at a wide range of Prandtl numbers ${Pr}$. This allows us to distinguish the parameter dependences of the horizontal heat flux and the boundary layer widths in terms of the Rayleigh number $\mbox {{Ra}}$ and Prandtl number ${Pr}$. For the considered parameter range $1\leq {Pr} \leq 100$, $10^{6} \leq \mbox {{Ra}} \leq 10^{9}$, we find the flow to be consistent with a ‘buoyancy-controlled’ regime where the heat flux is independent of the wall separation. For given ${Pr}$, the heat flux is found to scale linearly with the friction velocity $V_\ast$. Finally, we discuss the implications of our results for the parameterisation of heat and salt fluxes at vertical ice–ocean interfaces.

Motivated by observations of turbulence in the strongly stratified ocean thermocline, we use direct numerical simulations to investigate the interaction of a sinusoidal shear flow and a large-amplitude internal gravity wave. Despite strong nonlinearities in the flow and a lack of scale separation, we find that linear ray-tracing theory is qualitatively useful in describing the early development of the flow as the wave is refracted by the shear. Consistent with the linear theory, the energy of the wave accumulates in regions of negative mean shear where we observe evidence of convective and shear instabilities. Streamwise-aligned convective rolls emerge the fastest, but their contribution to irreversible mixing is dwarfed by shear-driven billow structures that develop later. Although the wave strongly distorts the buoyancy field on which these billows develop, the mixing efficiency of the subsequent turbulence is similar to that arising from Kelvin–Helmholtz instability in a stratified shear layer. We run simulations at Reynolds numbers Re of 5000 and 8000, and vary the initial amplitude of the internal gravity wave. For high values of initial wave amplitude, the results are qualitatively independent of $Re$. Smaller initial wave amplitudes delay the onset of the instabilities, and allow for significant laminar diffusion of the internal wave, leading to reduced turbulent activity. We discuss the complex interaction between the mean flow, internal gravity wave and turbulence, and its implications for internal wave-driven mixing in the ocean.

Turbulent mixing exerts a significant influence on many physical processes in the ocean. In a stably stratified Boussinesq fluid, this irreversible mixing describes the conversion of available potential energy (APE) to background potential energy (BPE). In some settings the APE framework is difficult to apply and approximate measures are used to estimate irreversible mixing. For example, numerical simulations of stratified turbulence often use triply periodic domains to increase computational efficiency. In this set-up, however, BPE is not uniquely defined and the method of Winters et al. (J. Fluid Mech., vol. 289, 1995, pp. 115–128) cannot be directly applied to calculate the APE. We propose a new technique to calculate APE in periodic domains with a mean stratification. By defining a control volume bounded by surfaces of constant buoyancy, we can construct an appropriate background buoyancy profile $b_\ast (z,t)$ and accurately quantify diapycnal mixing in such systems. This technique also permits the accurate calculation of a finite-amplitude local APE density in periodic domains. The evolution of APE is analysed in various turbulent stratified flow simulations. We show that the mean dissipation rate of buoyancy variance $\chi$ provides a good approximation to the mean diapycnal mixing rate, even in flows with significant variations in local stratification. When quantifying measures of mixing efficiency in transient flows, we find significant variation depending on whether laminar diffusion of a mean flow is included in the kinetic energy dissipation rate. We discuss how best to interpret these results in the context of quantifying diapycnal diffusivity in real oceanographic flows.

We study direct numerical simulations of turbulence arising from the interaction of an initial background shear, a linear background stratification and an external body force. In each simulation the turbulence produced is spatially intermittent, with dissipation rates varying over orders of magnitude in the vertical. We focus analysis on the statistically quasi-steady states achieved by applying large-scale body forcing to the domain, and compare flows forced by internal gravity waves with those forced by vertically uniform vortical modes. By considering the turbulent energy budgets for each simulation, we find that the injection of potential energy from the wave forcing permits a reversal in the sign of the mean buoyancy flux. This change in the sign of the buoyancy flux is associated with large, convective density overturnings, which in turn lead to more efficient mixing in the wave-forced simulations. The inhomogeneous dissipation in each simulation allows us to investigate localised correlations between the kinetic and potential energy dissipation rates. These correlations lead us to the conclusion that an appropriate definition of an instantaneous mixing efficiency, $\unicode[STIX]{x1D702}(t):=\unicode[STIX]{x1D712}/(\unicode[STIX]{x1D712}+\unicode[STIX]{x1D700})$ (where $\unicode[STIX]{x1D700}$ and $\unicode[STIX]{x1D712}$ are the volume-averaged turbulent viscous dissipation rate and fluctuation density variance destruction rate respectively) in the wave-forced cases is independent of an appropriately defined local turbulent Froude number, consistent with scalings proposed for low Froude number stratified turbulence.