The collision of two liquid drops is both an applied question (in rain formation or combustion, for example) and a beautiful basic situation, where impact involves liquid phases, making this problem worth studying in its own right. In a stimulating paper, Planchette, Lorenceau & Brenn (J. Fluid Mech., this issue, vol. 702, 2012, pp. 5–25) consider collisions between oil and water, which often lead to water drops protected by a shell of oil. By looking at the deformations during impact, they characterize the dynamical conditions leading to single encapsulation, and derive a criterion for avoiding fragmentation.

Binary collisions of drops of immiscible liquids are investigated experimentally at well-defined conditions of impact. In the experiments we vary all relevant properties of an aqueous and an oil phase, the impact parameter, the drop size and the relative velocity. The drops observed after the collisions exhibit three main phenomena: full encapsulation, head-on fragmentation, and off-centre fragmentation. The regimes characterized by these phenomena replace the ones observed in binary collisions of drops of the same liquid: coalescence, reflexive separation, and stretching separation. Our aim is a universal description of the two fragmentation thresholds of such collisions. Based on the capillary instability and an energy balance, we establish for head-on collisions a scaling law for the evolution of the threshold impact velocity with the properties of the liquids and the droplet size. The fragmentation threshold for off-centre collisions is compared to established models from the literature, which appear unsatisfactory. Introducing an effective impact parameter, which accounts empirically for the deformation and rotation of the drops upon impact, we describe this fragmentation threshold in a universal way. For both fragmentation thresholds, the agreement between experimental data and their theoretical representation is very good. Our work yields new insight into binary collisions of drops and proposes a perspective to develop a more general description with implications for binary collisions of drops of a single liquid as well.

Control of amplifier flows poses a great challenge, since the influence of environmental noise sources and measurement contamination is a crucial component in the design of models and the subsequent performance of the controller. A model-based approach that makes a priori assumptions on the noise characteristics often yields unsatisfactory results when the true noise environment is different from the assumed one. An alternative approach is proposed that consists of a data-based system-identification technique for modelling the flow; it avoids the model-based shortcomings by directly incorporating noise influences into an auto-regressive (ARMAX) design. This technique is applied to flow over a backward-facing step, a typical example of a noise-amplifier flow. Physical insight into the specifics of the flow is used to interpret and tailor the various terms of the auto-regressive model. The designed compensator shows an impressive performance as well as a remarkable robustness to increased noise levels and to off-design operating conditions. Owing to its reliance on only time-sequences of observable data, the proposed technique should be attractive in the design of control strategies directly from experimental data and should result in effective compensators that maintain performance in a realistic disturbance environment.

Large-amplitude internal solitary waves generate shear flows that intensify from the wings of the waves to their maxima. Upstream perturbations of the hydrostatic equilibrium in the form of wave packets along the path of wave propagation are expected to trigger shear instability and ultimately generate Kelvin–Helmholtz roll-ups. In contrast, as shown here with accurate simulations of incompressible stratified Euler equations, large internal waves can act as suppressors of perturbations. The precise understanding of the mechanisms leading to different outcomes, including whether instability is excited, is the focus of this work. Under the action of shear flows, small-amplitude wave packets undergo stretching and filamentation, which lead to significant absorption of perturbation energy into the background shear. It is found that this typical behaviour is present in the self-induced shear by internal waves, regardless of whether the shear is stable or unstable, and can leave a quieter state in the wave’s wake for a wide range of perturbation parameters. In the unstable case, even once perturbations are selected to excite the instability, our results show that this absorption can act to reduce growth in the strong-shear region, effectively making roll-up development observable only downstream of the wave crest. Our approach is both analytical and numerical; a model valid for relatively thin pycnoclines and suitable for local spectral analysis is devised and used. Energy diagnostics on the simulations are implemented to validate the numerics and illustrate the energy exchanges between background wave flow and its shear. A link between the absorption mechanism and the clustering of local eigenvalues along the wave is proposed. This promotes an energetic coupling among neutral modes stronger than what may be expected to occur in slowly varying flows, and gives rise to multi-modal transient dynamics of the kind often referred to as non-normality effects. For those cases in which the wave-induced shear meets the conditions for local instability, it is found that the growth of disturbances is selective with respect to the sign of the mode excited upstream. Elements of this phenomenon are interpreted by asymptotic analysis for spatial growth in time-independent slowly varying media.

Turbulent Couette–Poiseuille and Couette flows at different mean strain rates (velocity ratio of Couette wall to bulk flow, –), in a square duct at a bulk Reynolds number 10 000 are investigated by large eddy simulation. The numerical framework consists of a finite-volume method with a staggered-grid arrangement of dependent variables. Spatial derivatives are approximated using second-order centred differencing, and a fractional-step method is used for temporal integration. Simulations are conducted with grids. Secondary flow near the Couette wall shows a gradual change of vortex size and position as the moving wall velocity increased, where the two clockwise rotating vortices gradually merge in tandem with speed of the moving wall and form a large clockwise vortex. A linear relation is observed to exist between the angle of the two vortices and the parameter , and the angle saturates beyond . Also, at , together with a small counter-clockwise corner vortex, this vortex pattern is similar to that observed in the corner region of the duct flow with a free surface. The change of the vortex patterns also influences the dominant transport terms in the streamwise vorticity transport equation. Near the moving wall due to the reduction of the streamwise velocity fluctuation at the moving wall, turbulence structure gradually moves towards a rod-like axisymmetric turbulence, and as increases beyond , turbulence reverts to the disc-like structure.

This paper aims at understanding the influence of three-dimensional effects in hovering flapping flight. Numerical simulations at a Reynolds number of 1000 are performed to compare two types of flapping kinematics whose plunging phase is characterized by either a rectilinear translation or a revolving motion. In this way, we are able to isolate the three-dimensional effects induced by the free end condition from that induced by the spanwise incident velocity gradient (and the associated implicit Coriolis and centrifugal effects). In the rectilinear translation case, the analysis of the wake and of the aerodynamic loads reveals that the wingspan can be compartmented into three distinct regions whether it is predominantly subjected to an unstable two-dimensional flow, a stable three-dimensional flow or both two-dimensional and three-dimensional effects. It is found that this partitioning exhibits common features for three different aspect ratios of the wing. In conjunction with the previous results of Ringuette, Milano & Gharib (J. Fluid Mech., vol. 581, 2007, pp. 453–468), this suggests that the influence of the tip vortex over the wingspan is driven by a characteristic length scale. In addition, this length scale matches the position of the connecting point between leading and tip vortices observed in the revolving case, providing insight into the connecting process. In both translating and revolving cases, leading edge vortex attachment and strong spanwise velocities are found to be strongly correlated phenomena. Spanwise velocities (that mostly confine at the periphery of the vortices), together with downward velocities, do not only affect the leading edge vortex but also act as an inhibitor for the trailing edge vortex growth. As a consequence, cross-wake interactions between leading and trailing edge vortices are locally limited, hence contributing to flow stabilization.

We present a continuum model for melt water drainage through a spatially distributed system of connected subglacial cavities, and consider in this context the complications introduced when effective pressure or water pressure drops to zero. Instead of unphysically allowing water pressure to become negative, we model the formation of a partially vapour- or air-filled space between ice and bed. Likewise, instead of allowing sustained negative effective pressures, we allow ice to separate from the bed at zero effective pressure. The resulting model is a free boundary problem in which an elliptic obstacle problem determines hydraulic potential, and therefore also determines regions of zero effective pressure and zero water pressure. This is coupled with a transport problem for stored water, and the coupled system bears some similarities with Hele-Shaw and squeeze-film models. We present a numerical method for computing time-dependent solutions, and find close agreement with semi-analytical travelling wave and steady-state solutions. As may be expected, we find that ice–bed separation is favoured by high fluxes and low ice surface slopes and low bed slopes, while partially filled cavities are favoured by low fluxes and high slopes. At the boundaries of regions with zero water or effective pressure, discontinuities in water level are frequently present, either in the form of propagating shocks or as stationary hydraulic jumps accompanied by discontinuities in potential gradient.

We present a new model of subglacial drainage incorporating flow in a network of channels and a porous sheet, with water exchange between the two determined by pressure gradients. The sheet represents the average effect of many linked cavities, whilst the channels emerge from individual cavities that enlarge due to dissipation-induced melting. The model distinguishes cases when the water pressure drops to zero, in which case it allows for the drainage space to be only partially filled with water (free surface flow), and when the pressure reaches the ice overburden pressure, in which case it allows for uplift of the ice to whatever extent is needed to accommodate the water (flotation). Numerical solutions are found for a one-dimensional flow-line version of the model. The results capture typically observed or inferred features of subglacial drainage systems, including open channel flow at the ice margin, seasonal channel evolution, and high water pressures and uplift of the ice surface driven by rapid changes in water supply.

It has been pointed out that the anomalous enstrophy dissipation in non-smooth weak solutions of the two-dimensional Euler equations has a clue to the emergence of the inertial range in the energy density spectrum of two-dimensional turbulence corresponding to the enstrophy cascade as the viscosity coefficient tends to zero. However, it is uncertain how non-smooth weak solutions can dissipate the enstrophy. In the present paper, we construct a weak solution of the two-dimensional Euler equations from that of the Euler- equations proposed by Holm, Marsden & Ratiu (Phys. Rev. Lett., vol. 80, 1998, pp. 4173–4176) by taking the limit of . To accomplish this task, we introduce the -point-vortex () system, whose evolution corresponds to a unique global weak solution of the two-dimensional Euler- equations in the sense of distributions (Oliver & Shkoller, Commun. Part. Diff. Equ., vol. 26, 2001, pp. 295–314). Since the system is a formal regularization of the point-vortex system and it is known that, under a certain special condition, three point vortices collapse self-similarly in finite time (Kimura, J. Phys. Soc. Japan, vol. 56, 1987, pp. 2024–2030), we expect that the evolution of three -point vortices for the same condition converges to a singular weak solution of the Euler- equations that is close to the triple collapse as , which is examined in the paper. As a result, we find that the three -point vortices collapse to a point and then expand to infinity self-similarly beyond the critical time in the limit. We also show that the Hamiltonian energy and a kinematic energy acquire a finite jump discontinuity at the critical time, but the energy dissipation rate converges to zero in the sense of distributions. On the other hand, an enstrophy variation converges to the measure with a negative mass, which indicates that the enstrophy dissipates in the distributional sense via the self-similar triple collapse. Moreover, even if the special condition is perturbed, we can confirm numerically the convergence to the singular self-similar evolution with the enstrophy dissipation. This indicates that the self-similar triple collapse is a robust mechanism of the anomalous enstrophy dissipation in the sense that it is observed for a certain range of the parameter region.

High-speed particle image velocimetry (PIV) measurements of bypass transition reveal the breakdown of the ubiquitous streaks into turbulent spots. Individual streak velocity profiles are examined and contrasted with the root mean square profiles typically reported. An estimation of streak amplitude based on the modulation of the instantaneous boundary layer thickness is proposed. Examination of the PIV velocity fields shows how turbulent spot precursors, identified with concurrent hot-film recordings, consist of streamwise arrangements of positive and negative streaks. As secondary instability progresses, the interface between these streaks is observed to result in turbulent structures. In an attempt to further elucidate the role of the free stream turbulence, correlation maps are generated to determine the extent of the wall-normal fluctuations. Significant damping of the free stream is found within the boundary layer for all Reynolds numbers prior to the onset of spot precursors.

A reactive fluid dissolving the surface of a uniform fracture will trigger an instability in the dissolution front, leading to spontaneous formation of pronounced well-spaced channels in the surrounding rock matrix. Although the underlying mechanism is similar to the wormhole instability in porous rocks there are significant differences in the physics, due to the absence of a steadily propagating reaction front. In previous work we have described the geophysical implications of this instability in regard to the formation of long conduits in soluble rocks. Here we describe a more general linear stability analysis, including axial diffusion, transport-limited dissolution, nonlinear kinetics, and a finite-length system.

We investigate the mechanisms that affect the formation of columnar vortices for spin-up in a cylinder where the temperatures at the horizontal walls are prescribed. Numerical results from the three-dimensional Navier–Stokes equations show that a short-lived instability, suppressed by the combined effect of rotation and stratification, generates small temperature variations in the azimuthal direction. Temperature-gradient anomalies produce vorticity, and these vortices stir the fluid at the interface of the central vortex core thus reinforcing the temperature gradients. For sufficiently strong temperature gradients, the central vortex core breaks up into several columnar vortices. It is found, in particular, that small aspect ratios (height over radius of the cylindrical fluid layer) tend to inhibit the instability, while larger ones, , have the opposite effect. The main source of instability is the baroclinic vorticity production and not the presence of a solid sidewall since, counter-intuitively, the flow is more unstable for a free-slip boundary than for a no-slip one. Finally the effect of the temperature boundary conditions (isothermal versus adiabatic) on the horizontal boundaries has been investigated. The adiabatic boundaries help to preserve for longer times the sloping density interfaces that feed, with their potential energy, the baroclinic vorticity production; this results in more unstable flows.

Small planktonic organisms ubiquitously display unsteady or impulsive motion to attack a prey or escape a predator in natural environments. Despite this, the role of unsteady forces such as history and added mass forces on the low-Reynolds-number propulsion of small organisms, e.g. Paramecium, is poorly understood. In this paper, we derive the fundamental equation of motion for an organism swimming by means of the surface distortion in a non-uniform background flow field at a low-Reynolds-number regime. We show that the history and added mass forces are important as the product of Reynolds number and Strouhal number increases above unity. Our results for an unsteady squirmer show that unsteady inertial effects can lead to a non-zero mean velocity for the cases with zero streaming parameters, which have zero mean velocity in the absence of inertia.

Hydrodynamic effects of the relationship between the roll and pitch oscillations in low-aspect-ratio fins, with a laminar section and a rounded leading edge, flapping at transitional to moderately high Reynolds numbers, are considered. The fin is hinged at one end and its roll amplitude is large. Also examined is how this relationship is affected by spanwise twist, which alters the pitch oscillation amplitude and its phase relative to the roll motion. Force, efficiency and surface hot-film-anemometry measurements, and flow visualization are carried out in a tow tank. A fin of an abstracted penguin-wing planform and a NACA 0012 cross-section is used, and the chord Reynolds number varies from 3558 to 150 000 based on total speed. The fin is forced near the natural shedding frequency. Strouhal number and pitch amplitude are directly related when thrust is produced, and efficiency is maximized in narrow combinations of Strouhal number and pitch amplitude when oscillation of the leading-edge stagnation point is minimal. Twist makes the angle of attack uniform along the span and enhances thrust by up to 24 %, while maintaining high efficiency. Only 5 % of the power required to roll is spent to pitch, and yet roll and pitch are directly related. During hovering, dye visualization shows that a diffused leading-edge vortex is produced in rigid fins, which enlarges along the span; however, twist makes the vortex more uniform and the fin in turn requires less power to roll. Low-order phase maps of the measurements of force oscillation versus its derivative are modelled as due to van der Pol oscillators; the higher-order maps show trends in the sub-regimes of the transitional Reynolds number. Fin oscillation imparts a chordwise fluid motion, yielding a Stokes wave in the near-wall vorticity layer. When the roll and pitch oscillations are directly related, the wave is optimized: causing vorticity lift-up as the fin is decelerated at the roll extremity; the potential energy at the stagnation point is converted into kinetic energy; a vortex is produced as the lifted vorticity is wrapped around the leading edge; and free-stream reattachment keeps the vortex trapped. When the twist oscillation is phased along the span, this vortex becomes self-preserving at all amplitudes of twist, indicating the most stable (low-bandwidth) tuned nature.

The anisotropy of a low-Reynolds-number grid turbulence is investigated through direct numerical simulations based on the lattice Boltzmann method. The focus is on the anisotropy of the Reynolds-stress () and Reynolds-stress dissipation-rate () tensors and the approach taken is that using the invariant analysis introduced by Lumley & Newman (J. Fluid Mech., vol. 82, 1977, pp. 161–178). The grid is made up of thin square floating elements in an aligned configuration.

The anisotropy is initially high behind the grid and decays quickly as the downstream distance increases. The anisotropy invariant map (AIM) analysis shows that the return-to-isotropic trend of both and is fast and follows a perfectly axisymmetic expansion, although just behind the grid there is an initial tendency toward a one-component state. It is found that the linear relation with is satisfied during the return-to-isotropy phase of the turbulence decay, although close to the grid a form , where is a nonlinear function of , is more appropriate. For large downstream distances, becomes almost independent of , suggesting that despite the absence of an inertial range, the (dissipative) small scales present a high degree of isotropy. It is argued that (i) the very small values of the mean strain rate and (ii) the weak anisotropy of the large scales are in fact responsible for this result.

Surface stimulation of any physical origin (electrohydrodynamic, thermocapillary, etc.) has the goal of generating localized perturbations on the free surface or the velocity field of a capillary jet. Among these perturbations, only the axisymmetric ones are determinant for the jet breakup. Often, the stimulation is weak enough for a linear model to be applicable. Then, the stimulation can be described by means of the Green functions for stresses, both normal and tangential to the interface, the calculations of which are, in addition, uncoupled from the hydrodynamic variables. If a harmonic forcing is applied, these Green functions are combinations of the spatial modes whose associated poles lie inside the appropriate integration contour of the complex wavenumber plane. This is the motivation for a comprehensive enumeration and description of the spatial modes, which has not been done up to now. Modes familiar from a temporal analysis, the dominant and subdominant capillary modes and the hydrodynamic modes, are present, along with modes specific to a spatial analysis. Most of the latter have already been mentioned in the literature for inviscid jets, but not analysed. A mode not previously found is reported. In addition, a description of the velocity field associated with each mode is provided, as a tool to understand their physical origin and behaviour. The relative importance of each mode in both normal- and tangential-stress stimulations is discussed. Finally, the well-known merging of poles below a critical jet velocity, leading to absolute instability, is analysed in the light of the modal description.

The two-dimensional, incompressible flow in a plane sudden expansion is investigated numerically for a systematic variation of the geometry, covering expansion ratios (steps-to-outlet heights) from to . By means of a three-dimensional linear stability analysis global temporal modes are scrutinized. In a symmetric expansion the primary bifurcation is stationary and two-dimensional, breaking the mirror symmetry with respect to the mid-plane. The secondary asymmetric flow experiences a secondary instability to different three-dimensional modes, depending on the expansion ratio. For a moderately asymmetric expansion only one of the two secondary flows (the connected branch) is realized at low Reynolds numbers. Since the perturbed secondary flow does not deviate much from the symmetric secondary flow, both secondary stability boundaries are very close to each other. For very small and very large expansion ratios an asymptotic behaviour is found for suitably scaled critical Reynolds numbers and wavenumbers. Representative instabilities are analysed in detail using an a posteriori energy transfer analysis to reveal the physical nature of the instabilities. Depending on the geometry, pure centrifugal and elliptical amplification processes are identified. We also find that the basic flow can become unstable due to the effects of flow deceleration, streamline convergence and high shear stresses, respectively.

Linearized Navier–Stokes equations are solved to investigate the impact on the growth of near-wall turbulent streaks that arises from streamwise-travelling waves of spanwise wall velocity. The percentage change in streak amplification due to the travelling waves, over a range of wave parameters, is compared to published direct numerical simulation (DNS) predictions of turbulent skin-friction reduction; a clear correlation between the two is observed. Linearized simulations at a much higher Reynolds number, more relevant to aerospace applications, produce results that show no marked differences to those obtained at low Reynolds number. It is also observed that there is a close correlation between DNS data of drag reduction and a very simple characteristic of the ‘generalized’ Stokes layer generated by the streamwise-travelling waves.

We propose a general strategy for determining the minimal finite amplitude disturbance that triggers transition to turbulence in shear flows. This involves constructing a variational problem that searches over all disturbances of fixed initial amplitude which respect the boundary conditions, incompressibility and the Navier–Stokes equations, to maximize a chosen functional over an asymptotically long time period. The functional must be selected such that it identifies turbulent velocity fields by taking significantly enhanced values compared to those for laminar fields. We illustrate this approach using the ratio of the final to initial perturbation kinetic energies (energy growth) as the functional and the energy norm to measure amplitudes in the context of pipe flow. Our results indicate that the variational problem yields a smooth converged solution provided that the initial amplitude is below the threshold for transition. This optimal is the nonlinear analogue of the well-studied (linear) transient growth optimal. At the critical threshold, the optimization seeks out a disturbance that is on the ‘edge’ of turbulence during the period. Above this threshold, when disturbances trigger turbulence by the end of the period, convergence is then practically impossible. The first disturbance found to trigger turbulence as the amplitude is increased identifies the ‘minimal seed’ for the given geometry and forcing (Reynolds number). We conjecture that it may be possible to select a functional such that the converged optimal below threshold smoothly converges to the minimal seed at threshold. Our choice of the energy growth functional is shown to come close to this for the pipe flow geometry investigated here.

The dynamic pressure distribution on the bottom of a wave flume, due to the interaction of water waves with a submerged structure, is investigated experimentally and analytically, for both first- and second-order gravity waves of finite amplitude. The dynamic pressure excess is found to be very important, even for incoming waves propagating in deep water conditions. In this depth condition, a high pressure zone, thirty times larger than the dynamic pressure excess expected in the absence of the obstacle, is found in its vicinity. On the other hand, a low pressure zone is observed in the vicinity of the submerged obstacle for incoming waves propagating in smaller depth conditions. In any case, pressure gradients remain important. The second-order disturbance is found to be larger than first order in deep water conditions, for some specific conditions and locations. This result is interpreted in terms of nonlinear coupling of first-order components, including local modes.