A new dataset of uniform and steady sheet-flow experiments is presented in this paper. An acoustic concentration and velocity profiler (ACVP) is used to measure time-resolved profiles of collocated 2C velocity ($u,w$) and sediment concentration and to measure the time evolution of the bed interface position. Ensemble averaging over 11 similar experiment realisations is done to evaluate the mean profiles of streamwise velocity, concentration, sediment flux and Reynolds shear stress. The repeatability, stationarity and uniformity of the flow are carefully checked for a Shields number ${\it\theta}\approx 0.5$ and a suspension number of $S=1.1$. The mean profile analysis allows to separate the flow into two distinct layers: a suspension layer dominated by turbulence and a bed layer dominated by granular interactions. The bed layer can be further subdivided into a frictional layer capped by a collisional layer. In the suspension layer, the mixing length profile is linear with a strongly reduced von Karman parameter equal to 0.225. The Schmidt number is found to be constant in this region with a mean value of ${\it\sigma}_{s}=0.44$. The present results are then interpreted in terms of existing modelling approaches and the underlying assumptions are discussed. In particular, the well-known Rouse profile is shown to predict the concentration profile adequately in the suspension layer provided that all the required parameters can be evaluated separately. However, the strong intermittency of the flow observed in the bed layer under the impact of turbulent large-scale coherent flow structures suggests the limitations of averaged steady two-phase flow models.

We report on a novel method to generate fast transient microjets and study their characteristics. The simple device consists of two electrodes on a substrate with a hole in between. The side of the substrate with the electrodes is submerged in a liquid. Two separate microjets exit through the tapered hole after an electrical discharge is induced between the electrodes. They are formed during the expansion and collapse of a single cavitation bubble. The cavitation bubble dynamics as well as the jets were studied with high-speed photography at up to 500 000 f.p.s. With increasing jet velocity they become unstable and spray formation is observed. The jet created during expansion (first jet) is in most cases slower than the jet created during bubble collapse, which can reach up to $400~\text{m}~\text{s}^{-1}$. The spray exiting the orifice is at least in part due to the presence of cavitation in the microchannel as observed by high-speed recording. The effect of viscosity was tested using silicone oil of 10, 50 and 100 cSt. Interestingly, for all liquids the transition from a stable to an unstable jet occurs at $We\sim 4600$. We demonstrate that these microjets can penetrate into soft material; thus they can be potentially used as a needleless drug delivery device.

We consider strongly nonlinear waves in fluids in a uniform gravity field, and demonstrate that an incompressible wave motion, in which pressure remains constant in each fluid parcel, is supported by compressible fluids with free and rigid boundaries. We present exact analytic solutions of nonlinear hydrodynamics equations which describe the incompressible wave motion. The solutions provide an extension of the Gerstner wave in an incompressible fluid with a free boundary to waves in compressible three-dimensionally inhomogeneous moving fluids such as oceans and planetary atmospheres.

We present a study of the effect of finite detector integration/exposure time $E$, in relation to interrogation time interval ${\rm\Delta}t$, on analysis of Brownian motion of small particles using numerical simulation of the Langevin equation for both free diffusion and hindered diffusion near a solid wall. The simulation result for free diffusion recovers the known scaling law for the dependence of estimated diffusion coefficient on $E/{\rm\Delta}t$, i.e. for $0\leqslant E/{\rm\Delta}t\leqslant 1$ the estimated diffusion coefficient scales linearly as $1-(E/{\rm\Delta}t)/3$. Extending the analysis to the parameter range $E/{\rm\Delta}t\geqslant 1$, we find a new nonlinear scaling behaviour given by $(E/{\rm\Delta}t)^{-1}[1-((E/{\rm\Delta}t)^{-1})/3]$, for which we also provide an exact analytical solution. The simulation of near-wall diffusion shows that hindered diffusion of particles parallel to a solid wall, when normalized appropriately, follows with a high degree of accuracy the same form of scaling laws given above for free diffusion. Specifically, the scaling laws in this case are well represented by $1-((1+{\it\epsilon})(E/{\rm\Delta}t))/3$, for $E/{\rm\Delta}t\leqslant 1$, and $(E/{\rm\Delta}t)^{-1}[1-((1+{\it\epsilon})(E/{\rm\Delta}t)^{-1})/3]$, for $E/{\rm\Delta}t\geqslant 1$, where the small parameter ${\it\epsilon}$ depends on the size of the near-wall domain used in the estimation of the diffusion coefficient and value of $E$. For the range of parameters reported in the literature, we estimate ${\it\epsilon}<0.03$. The near-wall simulations also show a bias in the estimated diffusion coefficient parallel to the wall even in the limit $E=0$, indicating an overestimation which increases with increasing time delay ${\rm\Delta}t$. This diffusion-induced overestimation is caused by the same underlying mechanism responsible for the previously reported overestimation of mean velocity in near-wall velocimetry.

The assumption of low Reynolds number flow, or Stokes flow, is often applied to the understanding of a broad range of microfluidic devices, including micro-reactors, biomedical devices, and membraneless electrochemical cells. However, recent studies have shown that various inertial effects can play a significant role, even in microfluidic systems. In this work, two- and three-dimensional secondary flows are identified in a generic rectangular flow channel design consisting of a secondary channel feeding fluid into a main channel. We identify a scaling argument which is able to predict the occurrence of these secondary flows as a function of system parameters. The impact of these behaviours on the assumption of fully developed colaminar flow is investigated. This work considers a representative geometry, and identifies a set of conditions where inertial effects can play a key role in a microfluidic device.

Hydrodynamic interactions between two identical elastic dumbbells settling under gravity in a viscous fluid at low Reynolds number are investigated using the point-particle model. The evolution of a benchmark initial configuration is studied, in which the dumbbells are vertical and their centres are aligned horizontally. Rigid dumbbells and pairs of separate beads starting from the same positions tumble periodically while settling. We find that elasticity (which breaks the time-reversal symmetry of the motion) significantly affects the system dynamics. This is remarkable when taking into account that elastic forces are always much smaller than gravity. We observe oscillating motion of the elastic dumbbells, which tumble and change their length non-periodically. Independently of the value of the spring constant, a horizontal hydrodynamic repulsion appears between the dumbbells: their centres of mass move apart from each other horizontally. This motion is fast for moderate values of the spring constant $k$, and slows down when $k$ tends to zero or to infinity; in these limiting cases we recover the periodic dynamics reported in the literature. For moderate values of the spring constant, and different initial configurations, we observe the existence of a universal time-dependent solution to which the system converges after an initial relaxation phase. The tumbling time and the width of the trajectories in the centre-of-mass frame increase with time. In addition to its fundamental significance, the benchmark solution presented here is important to understanding general features of systems with a larger number of elastic particles, in regular and random configurations.

In this paper we describe a general and systematic approach to the centre-manifold reduction and normal form computation of flows undergoing complicated bifurcations. The proposed algorithm is based on the theoretical work of Coullet & Spiegel (SIAM J. Appl. Maths, vol. 43(4), 1983, pp. 776–821) and can be used to approximate centre manifolds of arbitrary dimension for large-scale dynamical systems depending on a scalar parameter. Compared with the classical multiple-scale technique frequently employed in hydrodynamic stability, the proposed method can be coded in a rather general way without any need to resort to the introduction and tuning of additional time scales. The method is applied to the dynamical system described by the incompressible Navier–Stokes equations showing that high-order, weakly nonlinear models of bifurcating flows can be derived automatically, even for multiple codimension bifurcations. We first validate the method on the primary Hopf bifurcation of the flow past a circular cylinder and after we illustrate its application to a codimension-two bifurcation arising in the flow past two side-by-side circular cylinders.

We investigate azimuthal instabilities which exist on the periphery of a non-turbulent liquid jet injected transversely into a gaseous cross-flow. We predict that the temporal growth of such instabilities may lead to the formation of interface corrugations, which are eventually sheared off of the jet surface (known as the jet ‘surface breakup’). In this study we employ temporal linear stability analyses to understand the nature of these instabilities. The analysis is based on a continuous formulation of momentum equations in which the jet and cross-flow are considered to be slightly miscible at the vicinity of the interface. We identify the shear instability as the primary destabilization mechanism in the flow. This inherently inviscid mechanism opposes the previously suggested mechanism of surface breakup (known as ‘boundary-layer stripping’), which is based on a viscous interpretation. The results show that the wavelengths of instabilities increase by moving away from the jet windward stagnation point toward the leeward point. We also investigate the influence of the jet-to-cross-flow density ratio on the flow stability and find that a higher ratio leads to formation of instabilities with higher wavenumbers on the jet surface. The results show that the density may have a non-monotonic stabilizing/destabilizing effect on the flow.

This paper reports on the experimental testing of oxygen-enriched porous fuel injection in a scramjet engine. Fuel was injected via inlet mounted, oxide-based ceramic matrix composite (CMC) injectors on both flow path surfaces that covered a total of 9.2 % of the intake surface area. All experiments were performed at an enthalpy of $3.93{-}4.25\pm 3.2\,\%~\text{MJ}~\text{kg}^{-1}$, flight Mach number 9.2–9.6 and an equivalence ratio of $0.493\pm 3\,\%$. At this condition, the engine was shown to be on the verge of achieving appreciable combustion. Oxygen was then added to the fuel prior to injection such that two distinct enrichment levels were achieved. Combustion was found to increase, by as much as 40 % in terms of combustion-induced pressure rise, over the fuel-only case with increasing oxygen enrichment. Further, the onset of combustion was found to move upstream with increasing levels of oxygen enrichment. Thrust, both uninstalled and specific, and specific impulse were found to be improved with oxygen enrichment. Enhanced fuel–air mixing due to the pre-mixing of oxygen with the fuel together with the porous fuel injection are believed to be the main contributors to the observed enhanced performance of the tested engine.

The interacting vorticity wave formalism for shear flow instabilities is extended here to the magnetohydrodynamic (MHD) setting, to provide a mechanistic description for stabilising and destabilising shear instabilities by the presence of a background magnetic field. The interpretation relies on local vorticity anomalies inducing a non-local velocity field, resulting in action at a distance. It is shown here that the waves supported by the system are able to propagate vorticity via the Lorentz force, and waves may interact. The existence of instability then rests upon whether the choice of basic state allows for phase locking and constructive interference of the vorticity waves via mutual interaction. To substantiate this claim, we solve the instability problem of two representative basic states, one where a background magnetic field stabilises an unstable flow and the other where the field destabilises a stable flow, and perform relevant analyses to show how this mechanism operates in MHD.

We investigate the longitudinal dispersion of a passive tracer by a gravity-driven flow in a porous medium consisting of a series of independent horizontal layers connected to a constant pressure source. We show that in a formation of given vertical extent, the total flux is only weakly dependent on the number of layers, and is very similar to that in a single layer of the same total depth. However, although the flow speed in each layer is approximately uniform, the speed gradually increases with layer depth. As a result, if a pulse of tracer is released in the flow it will migrate more rapidly through the lower layers, leading to longitudinal dispersion of the tracer. Eventually, the location of the tracer in the different layers may become separated in space so that a sufficiently distant observation well would detect a series of discrete pulses of tracer rather than the original coherent input, as would occur in a single permeable layer. For a constant pressure source, at long times, the standard deviation of the longitudinal distribution of tracer asymptotes to a fraction of order 0.1 of the position of the centre of mass, depending on the number of layers and the overpressure of the source.

We consider the behaviour of the flow within a fluid-filled torus when there is a sudden change in the rotation rate of the torus. Experimental work on this problem by Madden & Mullin (J. Fluid Mech., vol. 265, 1994, p. 217) demonstrated a flow with a rich and complex dynamics. In particular, planar (top-down) flow visualisation images show a well-defined laminar band at both the inner and outer bend of the toroidal pipe. Hewitt et al. (J. Fluid Mech., vol. 688, 2011, pp. 88–119) demonstrated the existence of finite-time singularities in the resulting viscous boundary layers, and linked the post-singularity structure to one of the laminar bands identified in experiments (Madden & Mullin J. Fluid Mech., vol. 265, 1994, p. 217; del Pino et al.Phys. Fluids, vol. 20 (12), 2008, 124104). The second band (or laminar front) identified by Madden & Mullin was conjectured by Hewitt et al. to be the result of a centrifugal instability, perhaps generated by small imperfections in the experimental apparatus. Here we explore this conjecture further, demonstrating that a small seam imperfection can generate substantial secondary motion but with considerably different dynamics than the centrifugally driven instability of Hewitt et al.

The linear growth of the spanwise correlation length scale with the distance from the wall in the logarithmic region of wall-bounded turbulent flows has been understood as a reflection of Townsend’s attached eddies. Based on this observation, in the present study, we perform a numerical experiment, which simulates energy-containing motions only at a given spanwise length scale in the logarithmic region, using their self-sustaining nature found recently. The self-sustaining energy-containing motions at each of the spanwise length scales are found to be self-similar with respect to the given spanwise length. Furthermore, their statistical structures are consistent with those of the attached eddies in the original theory, providing direct evidence on the existence of Townsend’s attached eddies. It is shown that a single self-sustaining attached eddy is composed of two distinct elements, one of which is a long streaky motion reaching the near-wall region, and the other is a relatively short vortical structure carrying all the velocity components. For the given spanwise length ${\it\lambda}_{z}$ between ${\it\lambda}_{z}^{+}=100$ and ${\it\lambda}_{z}\simeq 1.5h$, where $h$ is half the height of the channel, the former is found to be self-similar along $y\simeq 0.1{\it\lambda}_{z}$ and ${\it\lambda}_{x}\simeq 10{\it\lambda}_{z}$, while the latter is self-similar along $y\simeq 0.5{\it\lambda}_{z}\sim 0.7{\it\lambda}_{z}$ and ${\it\lambda}_{x}\simeq 2{\it\lambda}_{z}\sim 3{\it\lambda}_{z}$ where $y$ is the wall-normal direction. The scaling suggests that the smallest attached eddy would be a near-wall coherent motion in the form of a streak and quasi-streamwise vortices aligned to that, whereas the largest one would be an outer motion with a very-large-scale motion (VLSM) and large-scale motions (LSMs) aligned to that. The attached eddies in between, the size of which is proportional to their distance from the wall, contribute to the logarithmic region and fill the space caused by the length scale separation. The scaling is also found to yield behaviour consistent with the emergence of $k_{x}^{-1}$ spectra in a number of previous studies. Finally, a further discussion is provided, in particular on Townsend’s inactive motion and several recent theoretical findings.

We examine the dynamics of a binary mixture in a cubic cell subjected to a temperature differential and oscillatory forcing. The Soret effect, which is negative in the present study, provides a coupling mechanism by which a temperature gradient establishes a concentration gradient in a mixture. We present the results of experiments that were performed on the International Space Station (ISS) and compare the observations with the results of direct numerical simulations. The evolution of temperature and concentration fields is investigated by optical digital interferometry. One advantage of the experimental technique is the observation of the fields along two perpendicular directions of the cell, allowing us to restore the three-dimensional field. Experimental evidence disproves speculations that the ISS microgravity environment always affects diffusion-controlled processes. Furthermore, we demonstrate that imposed vibrations with constant frequency and amplitude create slow mean flows and that they do influence the diffusion kinetics. The perturbation of the diffusive fields scales as the square of the vibrational velocity. In addition to calculations of the full three-dimensional Navier–Stokes equations, a two-time-scale computational methodology is used for situations in which the forcing period is very small compared to the natural time scales of the problem. The simulations show excellent agreement with experimental observations.

A perching bird is able to rapidly decelerate while maintaining lift and control, but the underlying aerodynamic mechanism is poorly understood. In this work we perform a study on a simultaneously decelerating and pitching aerofoil section to increase our understanding of the unsteady aerodynamics of perching. We first explore the problem analytically, developing expressions for the added-mass and circulatory forces arising from boundary-layer separation on a flat-plate aerofoil. Next, we study the model problem through a detailed series of experiments at $\mathit{Re}=22\,000$ and two-dimensional simulations at $\mathit{Re}=2000$. Simulated vorticity fields agree with particle image velocimetry measurements, showing the same wake features and vorticity magnitudes. Peak lift and drag forces during rapid perching are measured to be more than 10 times the quasi-steady values. The majority of these forces can be attributed to added-mass energy transfer between the fluid and aerofoil, and to energy lost to the fluid by flow separation at the leading and trailing edges. Thus, despite the large angles of attack and decreasing flow velocity, this simple pitch-up manoeuvre provides a means through which a perching bird can maintain high lift and drag simultaneously while slowing to a controlled stop.

Following Malkus’s (J. Fluid Mech., vol. 1, 1956, pp. 521–539) proposal that turbulent Poiseuille channel flow maximises total viscous dissipation $D$, a variety of variational procedures have been explored involving the maximisation of different flow quantities under different constraints. However, the physical justification for these variational procedures has remained unclear. Here we address more recent claims that mean flow viscous dissipation $D_{m}$ should be maximised on the basis of a statistical stability argument, and that maximising $D_{m}$ yields realistic mean velocity profiles (Malkus, J. Fluid Mech., vol. 489, 2003, pp. 185–198). We clarify the connection between maximising $D_{m}$ and other flow quantities, verify Malkus & Smith’s, (J. Fluid Mech., vol. 208, 1989, pp. 479–507) claim that maximising the ‘efficiency’ yields realistic profiles and show that, in contrast, maximising $D_{m}$ does not yield realistic mean velocity profiles as recently claimed. This leads us to revisit Malkus’s statistical stability argument for maximising $D_{m}$ and to address some of its limitations. We propose an alternative statistical stability argument leading to a principle of minimum kinetic energy for fixed pressure gradient, which suggests a principle of maximum $D$ for fixed Reynolds number under certain conditions. We discuss possible ways to reconcile these conflicting results, focusing on the choice of constraints.

The scope of this work is to present and discuss the results obtained from simulating three-dimensional plunging breaking waves by solving the Navier–Stokes equations, in air and water. Recent progress in computational capabilities has allowed us to run fine three-dimensional simulations, giving us the opportunity to study for the first time fine vortex filaments generated during the early stage of the wave breaking phenomenon. To date, no experimental observations have been made in laboratories, and these structures have only been visualised in rare documentary footage (e.g. BBC 2009 South Pacific. Available on YouTube, 7BOhDaJH0m4). These fine coherent structures are three-dimensional streamwise vortical tubes, like vortex filaments, connecting the splash-up and the main tube of air, elongated in the main flow direction. The first part of the paper is devoted to the presentation of the model and numerical methods. The air entrainment occurring when waves break is then carefully described. Thanks to the high resolution of the grid, these fine elongated structures are simulated and explained.

Atmospheric turbulence is encountered frequently in flight and it creates oncoming flow disturbances of varying direction and magnitude for aircraft passing through turbulent zones. The unique measurement set-up on a motorised glider enables the investigation of the flow processes acting on a laminar wing section in flight through atmospheric turbulence. The expected quasi-steady aerofoil characteristics are deduced from an investigation of boundary-layer transition under calm flight conditions. Spanwise frequency–wavenumber spectra and comparisons with linear stability theory (LST) yield insight into the linear and weakly nonlinear stages of transition. Simultaneous measurement of the oncoming flow, characteristic flow quantities on the wing section and the motion of the aerofoil enables correlations between these quantities and provides insight into the unsteady flight physics. Emphasis is placed on the response of laminar–turbulent transition to moderate free-stream turbulence on both sides of the wing section. On the lower side of the aerofoil significant and rapid upstream fluctuations of transition are observed, which correspond closely to variations in the pressure distribution. Wavelet analysis is applied to gain insight into the composition of these streamwise excursions of the transition front in the time-pseudo-frequency domain. It is shown that they are driven by rapid transient base-flow changes and that transition is initiated by a short growth stage of Tollmien–Schlichting (TS) waves.

Experimental gradient-based optimization is used to maximize the propulsive efficiency of a heaving and pitching flexible panel. Optimum and near-optimum conditions are studied via direct force measurements and particle image velocimetry (PIV). The net thrust and power scale predictably with the frequency and amplitude of the leading edge, but the efficiency shows a complex multimodal response. Optimum pitch and heave motions are found to produce nearly twice the efficiencies of optimum heave-only motions. Efficiency is globally optimized when (i) the Strouhal number is within an optimal range that varies weakly with amplitude and boundary conditions; (ii) the panel is actuated at a resonant frequency of the fluid–panel system; (iii) heave amplitude is tuned such that trailing-edge amplitude is maximized while the flow along the body remains attached; and (iv) the maximum pitch angle and phase lag are chosen so that the effective angle of attack is minimized. The multi-dimensionality and multi-modality of the efficiency response demonstrate that experimental optimization is well-suited for the design of flexible underwater propulsors.

Recent field observations and modelling of breaking surface gravity waves suggest that air-entraining breaking is not sufficiently dissipative of surface gravity waves to balance the dynamics of wind-wave growth and nonlinear interactions with dissipation for the shorter gravity waves of $O(10)$  cm wavelength. Theories of parasitic capillary waves that form at the crest and forward face of shorter steep gravity waves have shown that the dissipative effects of these waves may be one to two orders of magnitude greater than the viscous dissipation of the underlying gravity waves. Thus the parasitic capillaries may provide the required dissipation of the short wind-generated gravity waves. This has been the subject of speculation and conjecture in the literature. Using the nonlinear theory of Fedorov & Melville (J. Fluid Mech., vol. 354, 1998, pp. 1–42), we show that the dissipation due to the parasitic capillaries is sufficient to balance the wind input to the short gravity waves over some range of wave ages and wave slopes. The range of gravity wave lengths on which these parasitic capillary waves are dynamically significant approximately corresponds to the range of short gravity waves that Cox & Munk (J. Mar. Res., vol. 13, 1954, pp. 198–227) found contributed significantly to the mean square slope of the ocean surface, which they measured to be proportional to the wind speed. Here we show that the mean square slope predicted by the theory is proportional to the square of the friction velocity of the wind, ${u_{\ast }}^{2}$, for small wave slopes, and approximately $u_{\ast }$ for larger slopes.