Thermally driven vertical convection (VC) – the flow in a box heated on one side and cooled on the other side, is investigated using direct numerical simulations with Rayleigh numbers over the wide range of $10^7\le Ra\le 10^{14}$ and a fixed Prandtl number $Pr=10$ in a two-dimensional convection cell with unit aspect ratio. It is found that the dependence of the mean vertical centre temperature gradient $S$ on $Ra$ shows three different regimes: in regime I ($Ra \lesssim 5\times 10^{10}$), $S$ is almost independent of $Ra$; in the newly identified regime II ($5\times 10^{10} \lesssim Ra \lesssim 10^{13}$), $S$ first increases with increasing $Ra$ (regime $\textrm {{II}}_a$), reaches its maximum and then decreases again (regime $\textrm {{II}}_b$); and in regime III ($Ra\gtrsim 10^{13}$), $S$ again becomes only weakly dependent on $Ra$, being slightly smaller than in regime I. The transition from regime I to regime II is related to the onset of unsteady flows arising from the ejection of plumes from the sidewall boundary layers. The maximum of $S$ occurs when these plumes are ejected over approximately half of the area (downstream) of the sidewalls. The onset of regime III is characterized by the appearance of layered structures near the top and bottom horizontal walls. The flow in regime III is characterized by a well-mixed bulk region owing to continuous ejection of plumes over large fractions of the sidewalls, and, as a result of the efficient mixing, the mean temperature gradient in the centre $S$ is smaller than that of regime I. In the three different regimes, significantly different flow organizations are identified: in regime I and regime $\textrm {{II}}_a$, the location of the maximal horizontal velocity is close to the top and bottom walls; however, in regime $\textrm {{II}}_b$ and regime III, banded zonal flow structures develop and the maximal horizontal velocity now is in the bulk region. The different flow organizations in the three regimes are also reflected in the scaling exponents in the effective power law scalings $Nu\sim Ra^\beta$ and $Re\sim Ra^\gamma$. Here, $Nu$ is the Nusselt number and $Re$ is the Reynolds number based on maximal vertical velocity (averaged over vertical direction). In regime I, the fitted scaling exponents ($\beta \approx 0.26$ and $\gamma \approx 0.51$) are in excellent agreement with the theoretical predictions of $\beta =1/4$ and $\gamma =1/2$ for laminar VC (Shishkina, Phys. Rev. E., vol. 93, 2016, 051102). However, in regimes II and III, $\beta$ increases to a value close to 1/3 and $\gamma$ decreases to a value close to 4/9. The stronger $Ra$ dependence of $Nu$ is related to the ejection of plumes and the larger local heat flux at the walls. The mean kinetic dissipation rate also shows different scaling relations with $Ra$ in the different regimes.

Misaligned wind turbine rotors redirect the wake, and manipulate the wake shape by introducing a counter-rotating vortex pair. This mechanism is of great interest for improving wind farm power output through static or dynamic misalignment. In this study, cross-plane stereo-particle image velocimetry measurements are used to characterize the wake evolution for tilt misalignment and verify differences with yaw misalignment. Blockage from the ground, shear in the velocity profile, turbulence levels, hub-vortices and tip-vortices are found to strongly affect wake evolution for a tilted wind turbine resulting in a non-symmetric behaviour for upwards deflecting or downwards deflecting tilt. The downwards deflection of a negatively tilted wind turbine is found to result in the most benefits for wake recovery and power availability downstream through increased wake-curling, faster wake-recovery, and downdraft of high-momentum flow.

In this experimental study, multiscale rough surfaces with regular (cuboid) elements are used to examine the effects of roughness-scale hierarchy on turbulent boundary layers. Three iterations have been used with a first iteration of large-scale cuboids onto which subsequent smaller cuboids are uniformly added, with their size decreasing with a power-law as the number increases. The drag is directly measured through a floating-element drag balance, while particle image velocimetry allowed the assessment of the flow field. The drag measurements revealed the smallest roughness iteration can contribute to nearly 7 $\%$ of the overall drag of a full surface, while the intermediate iterations are responsible for over $12\,\%$ (at the highest Reynolds number tested). It is shown that the aerodynamic roughness length scale between subsequent iterations varies linearly, and can be described with a geometrical parameter proportional to the frontal solidity. Mean and turbulent statistics are evaluated using the drag information, and highlighted substantial changes within the canopy region as well as in the outer flow, with modifications to the inertial sublayer (ISL) and the wake region. These changes are shown to be caused by the presence of large-scale secondary motions in the cross-plane, which itself is believed to be a consequence of the largest multiscale roughness phase (spacing between largest cuboids), shown to be of the same order of magnitude as the boundary-layer thickness. Implications on the classical similarity laws are additionally discussed.

Compound particles are a class of composite systems in which solid particles encapsulated in a fluid droplet are suspended in another fluid. They are encountered in various natural and biological processes, for e.g. nucleated cells, hydrogels, microcapsules etc. Generation and transportation of such multiphase structures in microfluidic devices is associated with several challenges because of the poor understanding of their structural stability in a background flow and the rheological characteristics of their dispersions. Hence, in this work, we analyse the flow in and around a concentric compound particle and investigate the deformation dynamics of the confining drop and its stability against breakup in imposed linear flows. In the inertia-less limit (Reynolds number, $Re \ll 1$) and assuming that the surface tension force dominates the viscous forces (low capillary number, $Ca$, limit), we obtain analytical expressions for the velocity and pressure fields up to ${O}(Ca)$ for a compound particle subjected to a linear flow using a domain perturbation technique. Simultaneously, we determine the deformed shape of the confining drop correct up to ${O}(Ca^2)$, facilitating the following. (i) Since ${O}(Ca^2)$ calculations account for the rotation of the anisotropically deformed interface, the reorientation dynamics of the deformed compound particles is determined. (ii) Calculations involving the ${O}(Ca^2)$ shape of the confining interface are found to be important for compound particles as ${O}(Ca)$ calculations make qualitatively different predictions in generalised extensional flows. (iii) An ${O}(Ca)$ constitutive equation for the volume-averaged stress for a dilute dispersion of compound particles was developed to study both shear and extensional rheology in a unified framework. Our analysis shows that the presence of an encapsulated particle always enhances all the measured rheological quantities such as the effective shear viscosity, extensional viscosity and normal stress differences. (iv) Moreover, linear viscoelastic behaviour of a dilute dispersion of compound particles is characterised in terms of complex modulus by subjecting the dilute dispersion to a small-amplitude oscillatory shear (SAOS) flow. (v) Various expressions pertaining to a suspension of particles, drops, and particles coated with a fluid film are also derived as limiting cases of compound particles.

The natural frequencies for sloshing without coupling with lateral tank motions differ from the natural sloshing frequencies with coupling. A consequence is that the nonlinear multimodal sloshing theory for prescribed tank motion should be revised when studying the liquid sloshing dynamics in connection with marine structures in ocean waves. The needed revisions are done for a rectangular rigid tank with finite liquid depth. Steady-state resonant solutions of the constructed nonlinear modal equations are derived to analytically describe the coupled resonant sloshing and sway of a floating rigid body in regular incident deep water waves with two-dimensional flow conditions at the lowest coupled sloshing–sway natural frequency. The steady-state theoretical results are validated by comparing them with the model tests by Rognebakke & Faltinsen (J. Ship Res., vol. 47, issue 3, 2003, pp. 208–221). The occurrence of an instability frequency range is theoretically justified.

Particles of differing sizes are notoriously prone to segregation in shear driven flows under the action of gravity. This has important implications in many industrial processes, where particle-size segregation can lead to flow problems and reduced product quality, as well as longer product development and start-up times. Particle-size segregation also readily occurs in many hazardous geophysical mass flows (such as snow avalanches, debris flows and volcanic pyroclastic flows) and can lead to the formation of destructive bouldery flow fronts and significantly longer runouts. Although general theories exist to model particle-size segregation, the detailed functional dependence of the segregation flux on the shear rate, gravity, pressure, particle concentration, grain size and grain-size ratio is still not known. This paper describes refractive-index matched oscillatory shear-cell experiments that shed light on the segregation velocity in the two extreme cases of (i) a single large intruder rising up through a matrix of smaller grains, and (ii) a single small intruder percolating down through a matrix of large particles. Despite the sometimes markedly different time scales for segregation in these two situations, a unifying scaling law has been found that is able to collapse all the experimental data over a wide range of shear rates and grain-size ratios in the range $[1.17,4.17]$. The resulting functional form is easily generalizable to intermediate concentrations and can quantitatively capture laboratory experiments and numerical simulations with a $50\ {:}\ 50$ mix of large and small grains.

The evolution of the interfacial instabilities on a smooth and laminar liquid sheet formed by the oblique impingement of two liquid jets is investigated in the presence of external acoustic forcing. The forced liquid sheet exhibits different instability patterns, such as flapping flag, sinuous breakup and vibrating membrane structure, depending on the forcing frequency. The transition frequencies, where the instability pattern changes, are obtained using the linear theory of sheet instability and are shown to be increasing functions of the jet Weber number. These different instability patterns show distinct growth rate behaviours, which cannot be predicted using a unique theory such as the conventional linear theory of sheet breakup. The dynamics of liquid sheet breakup in the presence of forcing is further assessed with the help of established instability theories. Our measurements show that the sheet thickness distribution plays a major role in dictating the relative influence of different mechanisms on the instability characteristics. All the results are finally consolidated to arrive at the instability regime map on the Weber number versus frequency ($We$–$\tilde {\omega }$) space to serve as a guideline for adopting a suitable modelling approach.

Formation of a microlayer in nucleate boiling results in substantial changes to the heat-transfer dynamics of the process, and establishing a criterion to distinguish between the microlayer and contact-line regimes encountered is necessary to obtain a quantitative description of the nucleate-boiling phenomenon. In this work, we describe microlayer formation as a dewetting transition in the presence of phase change. We derive a specific criterion for modelling the transition based on a synthesis of existing theoretical, experimental and numerical data. As a result, a new formulation of the transition criterion is developed, and applied to (i) reference data from a dewetting experiment of a volatile liquid, and (ii) results derived from a high-resolution direct numerical simulation of nucleate boiling. In both cases, very good agreement is found for Jakob numbers $\lesssim$75, a range which covers many important boiling situations.

Spontaneous motion of liquid droplets can occur on hydrophobic, micro-structured, solid surfaces comprising a structural gradient. In this study, we examine such motion experimentally and explain our observations by invoking variable droplet–surface interactions (both actuation and resistance forces) arising from the structural gradient. The oscillatory motion of the droplet constitutes an integral aspect of the behaviour and this is incorporated into the overall modelling. The theoretical model features a truncated spheroid for the drop shape (flattened in the region of solid contact) coupled with the oscillatory and alternate leading and trailing motion of the contact line. Results from the model and experiments provide good qualitative and quantitative agreement. The component of the vertical oscillation is found to help overcome wetting hysteresis and actuate the motion, this being a key element for the completeness of the model.

We study the reactive displacement of two miscible fluids in channel flows and establish a quantitative link between fluid stretching and chemical reactivity. At the mixing interface, the two fluids react according to the instantaneous irreversible bimolecular reaction $A + B \rightarrow C$. We simulate the advection–diffusion–reaction problem using a random walk based reactive particle method that is free of numerical dispersion. The relative contributions of stretching and diffusion to mixing-limited reaction is controlled by changing the Péclet number, and the channel roughness is also systematically varied. We observe optimal ranges of fluid stretching that maximize reactivity, which are captured by a Lagrangian stretching measure based on an effective time period that honours the stretching history. We show that the optimality originates from the competition between the enhanced mixing by fluid stretching and the mass depletion of the reactants. We analytically derive the spatial distribution of reaction products using a lamellar formulation and successfully predict the optimal ranges of fluid stretching, which are consistent across different levels of channel roughness. These findings provide a mechanistic understanding of how the interplay between fluid stretching, diffusion and channel roughness controls mixing-limited reactions in rough channel flows, and show how reaction hot spots can be predicted from the concept of optimal fluid stretching.

When two converging flows of unequal velocities and different flow directions come into contact, large-scale coherent structures are generated. For shallow conditions and small angles between the incoming streams, the mixing layer (ML) forming downstream of the confluence apex contains quasi-two-dimensional, Kelvin–Helmholtz (KH) vortices. Due to flow shallowness (e.g. stabilizing effect of bed friction), these vortices gradually lose their coherence at large distances from the ML origin. For large angles between the incoming streams, the spatial development of the shallow ML is more complex as strongly coherent, streamwise-oriented-vortical (SOV) cells form in the vicinity of the shallow ML and helical cells of secondary flow are generated due to curvature effects. The present paper uses three-dimensional eddy-resolving numerical simulations to study the effects of varying the angles between the two incoming channels and the downstream channel, the velocity ratio (VR) of the incoming flows and the flow depth on flow, turbulence structure and sediment entrainment mechanisms inside the ML and its surroundings. The simulations are performed for highly idealized conditions in which the ML develops in a wide channel (no interactions with the channel banks), over a flat bed and the flow depth is constant. Simulation results show that the SOV cells play an important role in the redistribution of the streamwise momentum. Some of the SOV cells are subject to bimodal oscillations in the lateral direction which induce strong interactions between the SOV cells and the ML vortices and sharply increase mixing. As for the case of a shallow ML developing between parallel streams, vortex pairing ceases some distance from the ML origin, the KH vortices start losing their coherence and the ML assumes an undulatory shape. The paper describes the effects of VR, angle between the incoming streams (α = 0° and 60°), planform geometry (symmetric vs. asymmetric confluences with α = 60°) and flow shallowness on the transverse shift of the ML centreline, ML width, dynamics of the KH vortices and on the formation, position and circulation of the SOV cells. In most of the α = 60° cases, the largest bed shear stresses are induced beneath the SOV cells rather than beneath the high-speed stream and the SOV cells play a major role in enhancing mixing between the two streams.

Direct numerical simulations are performed to investigate a stratified shear layer at high Reynolds number ($Re$) in a study where the Richardson number ($Ri$) is varied among cases. Unlike previous work on a two-layer configuration in which the shear layer resides between two layers with constant density, an unbounded fluid with uniform stratification is considered here. The evolution of the shear layer includes a primary Kelvin–Helmholtz shear instability followed by a wide range of secondary shear and convective instabilities, similar to the two-layer configuration. During transition to turbulence, the shear layers at low $Ri$ exhibit a period of thickness contraction (not observed at lower $Re$) when the momentum and buoyancy fluxes are counter-gradient. The behaviour in the turbulent regime is significantly different from the case with a two-layer density profile. The transition layers, which are zones with elevated shear and stratification that form at the shear-layer edges, are stronger and also able to support a significant internal wave flux. After the shear layer becomes turbulent, mixing in the transition layers is shown to be more efficient than that which develops in the centre of the shear layer. Overall, the cumulative mixing efficiency ($E^C$) is larger than the often assumed value of 1/6. Also, $E^C$ is found to be smaller than that in the two-layer configuration at moderate Ri. It is relatively less sensitive to background stratification, exhibiting little variation for $0.08 \leqslant Ri \leqslant 0.2$. The dependence of mixing efficiency on buoyancy Reynolds number during the turbulence phase is qualitatively similar to homogeneous sheared turbulence.

Here we provide a unified theoretical description of two different physical situations in which liquid jets are expelled out of the bulk of a liquid as a consequence of the capillary collapse of a void. We demonstrate that the velocity field giving rise to the emergence of these jets can be calculated as the flow generated by a line of sinks with a length and an intensity that can be expressed in terms of the initial cavity radius and the wavelength and velocity of the capillary waves propagating along the cavity walls. The predicted jet speeds, which are expressed through algebraic equations, are in good quantitative agreement with those obtained from experiments and from the simulations of bubbles bursting on a free surface or after the implosion of the crater formed when a drop impacts a liquid pool.

Oscillatory flow features are common in the unstart of hypersonic mixed-compression intakes and can be classified as low-amplitude or high-amplitude oscillatory unstarted flows. The low-amplitude oscillatory unstarted flow is driven by the shear layer from shock interactions ahead of the cowl, while the high-amplitude oscillatory unstarted flow is driven by the separation caused by shock–boundary-layer interaction on the ramp. While previous studies have observed these flow features and reported their associated frequency, there is no simple criterion available for predicting which mode will occur, and there is a lack of consensus on the appropriate frequency scaling parameter. We study a mixed-compression hypersonic intake in a hypersonic wind tunnel by varying the internal contraction ratio and the throttling ratio to observe various kinds of unstart regimes. Two significant conclusions emerge from considering the results for high-throttling-ratio conditions $(TR > 0.55)$ from the current as well as previous studies. Firstly, the actual shock-on-lip condition at the cowl corresponding to the unthrottled condition, as observed from schlieren images, demarcates the boundary between the two modes of oscillatory unstart flows upon throttling. Secondly, a suitable length scale $(l^*)$, defined as the extent of the subsonic region in the unstarted flow (as observed from the experimental schlieren images), gives the appropriate frequency scaling parameter ($f^* = a_0/4l^*$ where $a_0$ is the stagnation acoustic speed).

The paper reports on three direct numerical simulations of sediment transport with different prototypical non-spherical particle shapes, ranging from prolate to oblate to a triaxial ellipsoidal shape, and a fourth simulation with spherical particles. All physical and numerical parameters of these simulations are identical, including bulk Reynolds number, equivalent particle diameter and sphericity of the non-spherical particles, so that the shape is the only difference. Numerous statistical quantities assessing the computed solutions are presented, such as mean particle velocity, sediment-bed porosity, equivalent roughness height, mean fluid velocity and Reynolds stresses. These allow a detailed analysis and reveal the considerable influence of the particle shape on the entire flow: fluid as well as particles. It is found that oblate spheroids yield the highest porosity among the four cases and spherical particles the lowest mean porosity. Prolate particles exhibit the strongest tendency to form spanwise clusters and generate the highest total bed friction. Triaxial ellipsoids representing natural medium-size sand particles yield results in between these cases, still exhibiting substantial differences with respect to spherical particles. The mean sediment transport rate is determined in all four cases showing that this quantity cannot be correlated to the mean shear stress of the fluid alone. This emphasises the importance of the particle shape in modelling the transport of natural sediment, and the results provide well-controlled data points for this task.

The paper addresses the fluid–structure interaction of submerged aquatic canopies, with particular focus on the complex interplay between coherent flow structures and the motion of vegetation elements. New insights into the underlying mechanisms are gained from a large eddy simulation of a submerged model canopy flow. The model canopy is made up of 800 highly flexible blades, each individually resolved by an immersed boundary method. The obtained high-resolution flow data reveal well-known coherent turbulent structures, including velocity streaks, Kelvin–Helmholtz (KH) vortices in the mixing layer as well as hairpin (HP) vortices in the outer flow region. The present results show that the interaction of these prototypical structures plays a key role creating unique turbulent features such as composite KH/HP vortices located between a high-speed and low-speed streak. Under the influence of these pronounced eddies, groups of blades respond by a strong local reconfiguration. Due to the convection of the coherent structures by the mean flow this causes an apparent wave-like motion of the canopy in streamwise direction, known as monami. A frequency analysis of this phenomenon shows that the vegetation responds almost passively, merely reflecting local flow conditions.

The motion of glaciers over their bedrock or drops of fluid along a solid surface can vary dramatically when these substrates are lubricated. We investigate the coupled flow of a gravity current (GC) of a strain-rate softening fluid that is lubricated by a denser, lower-viscosity Newtonian fluid. We present a set of experiments in which such GCs are discharged axisymmetrically and at constant flux over a flat surface. Using imaging techniques, we follow the evolution of the front and thickness field of each fluid. We find that, unlike purely Newtonian lubricated GCs, the fronts of the non-Newtonian and lubricating Newtonian fluids have a power-law time evolution with different exponents. Each of these exponents is similar to that of a non-lubricated GC of the same fluid. Nevertheless, the fronts of our lubricated GCs evolve faster than those of the corresponding non-lubricated GCs owing to larger intercepts. In addition, in contrast with the monotonically declining thickness of non-lubricated GCs, the thickness of the lubricated, non-Newtonian fluid is nearly uniform, and that of the lubricating fluid is non-monotonic with localised spikes. Despite these complex thickness patterns, lubricated GCs remain axisymmetric as long as the flux of the lubricating fluid is sufficiently smaller than that of the non-Newtonian fluid.

We investigate the removal of a dense bottom layer by a gravity current, via Navier–Stokes simulations in the Boussinesq limit. The problem is governed by a dimensionless thickness parameter for the bottom layer, and by the ratio of the density differences between bottom layer, gravity current and ambient fluids. A quasisteady gravity current forms that propagates along the interface and displaces some of the dense bottom fluid, which accumulates ahead of the gravity current and forms an undular bore or a series of internal gravity waves. Depending on the ratio of the gravity current front velocity to the linear shallow-water wave velocity, we observe the existence of different regimes, characterized by small-amplitude waves or by a train of steep, nonlinear internal waves. We develop a semiempirical model that provides reasonable estimates of several important flow properties. We also formulate a more sophisticated, self-contained model based on the conservation principles for mass and vorticity that does not require empirical closure assumptions. This model is able to predict such quantities as the gravity current height and the internal wave or bore velocity as a function of the governing dimensionless parameters, generally to within approximately a 10 $\%$ accuracy. An energy budget analysis provides information on the rates at which potential energy is converted into kinetic energy and then dissipated, and on the processes by which energy is transferred from the gravity current fluid to the dense and ambient fluids. We observe that the energy content of thicker and denser bottom layers grows more rapidly.

A rigorous bound on the energy dissipation in the flow around a sphere is derived analytically. The bound is allowed to depend on the viscosity of the fluid, the radius of the sphere, the free-stream velocity and the amplitude of the difference between the actual flow field and the irrotational flow around the sphere. The bound becomes asymptotically independent of viscosity as the Reynolds number tends to infinity. Provided that the velocity fluctuations are of the same order as the free-stream velocity, the scaling of the bound is consistent with a Reynolds-number-independent drag coefficient at high Reynolds numbers.

We study numerically the melting of a horizontal layer of a pure solid above a convecting layer of its fluid rotating about the vertical axis. In the rotating regime studied here, with Rayleigh numbers of order $10^{7}$, convection takes the form of columnar vortices, the number and size of which depend upon the Ekman and Prandtl numbers, as well as the geometry – periodic or confined. As the Ekman and Rayleigh numbers vary, the number and average area of vortices vary in inverse proportion, becoming thinner and more numerous with decreasing Ekman number. The vortices transport heat to the phase boundary, thereby controlling its morphology characterized by the number and size of the voids formed in the solid, and the overall melt rate, which increases when the lower boundary is governed by a no-slip rather than a stress-free velocity boundary condition. Moreover, the number and size of voids formed are relatively insensitive to the Stefan number, here inversely proportional to the latent heat of fusion. For small values of the Stefan number, the convection in the fluid reaches a slowly evolving geostrophic state wherein columnar vortices transport nearly all the heat from the lower boundary to melt the solid at an approximately constant rate. In this quasi-steady state, we find that the Nusselt number, characterizing the heat flux, co-varies with the interfacial roughness, for all the flow parameters and Stefan numbers considered here. This confluence of processes should influence the treatment of moving boundary problems, particularly those in astrophysical and geophysical problems where rotational effects are important.