Chris Howland

Motivated by the ablation of vertical ice faces in salt water, we use three-dimensional direct numerical simulations to investigate the heat and salt fluxes in two-scalar vertical convection. For parameters relevant to ice-ocean interfaces in the convection-dominated regime, we observe that the salinity field drives the convection and that heat is essentially transported as a passive scalar. By varying the diffusivity ratio of heat and salt (i.e., the Lewis number \(Le\)), we identify how the different molecular diffusivities affect the scalar fluxes through the system. Away from the walls, we find that the heat transport is determined by a turbulent Prandtl number of \(Pr_t ≈1\) and that double-diffusive effects are practically negligible. However, the difference in molecular diffusivities plays an important role close to the boundaries. In the (unrealistic) case where salt diffused faster than heat, the ratio of salt-to-heat fluxes would scale as \(Le^{1/3}\), consistent with classical nested scalar boundary layers. However, in the realistic case of faster heat diffusion (relative to salt), we observe a transition towards a \(Le^{1/2}\) scaling of the ratio of the fluxes. This coincides with the thermal boundary layer width growing beyond the thickness of the viscous boundary layer. We find that this transition is not determined by a critical Lewis number, but rather by a critical Prandtl number \(Pr ≈10\), slightly below that for cold seawater where \(Pr=14\). We compare our results to similar studies of sheared and double-diffusive flow under ice shelves, and discuss the implications for fluxes in large-scale ice-ocean models. By coupling our results to ice-ocean interface thermodynamics, we describe how the flux ratio impacts the interfacial salinity, and hence the strength of solutal convection and the ablation rate.

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 \(Ra\) and Prandtl number \(Pr\). For the considered parameter range \(1≤Pr≤100\), \(10^6≤Ra≤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_∗\). Finally, we discuss the implications of our results for the parameterisation of heat and salt fluxes at vertical ice-ocean interfaces.

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, \(η(t):=χ/(χ+ \varepsilon)\) (where \(\varepsilon\) and \(χ\) 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.

Droplets splash when they impact dry, flat substrates above a critical velocity that depends on parameters such as droplet size, viscosity, and air pressure. By imaging ethanol drops impacting silicone gels of different stiffnesses, we show that substrate stiffness also affects the splashing threshold. Splashing is reduced or even eliminated: droplets on the softest substrates need over 70% more kinetic energy to splash than they do on rigid substrates. We show that this is due to energy losses caused by deformations of soft substrates during the first few microseconds of impact. We find that solids with Young’s moduli \( ≲100  \textrm{kPa} \) reduce splashing, in agreement with simple scaling arguments. Thus, materials like soft gels and elastomers can be used as simple coatings for effective splash prevention. Soft substrates also serve as a useful system for testing splash-formation theories and sheet-ejection mechanisms, as they allow the characteristics of ejection sheets to be controlled independently of the bulk impact dynamics of droplets.