The present investigation considers the vortex wakes behind finite-length yawed cylinders that are stationary or vibrating transversly in a uniform flow. In the Reynolds-number range 160–103 a number of three-dimensional and nominally two-dimensional wake flows are observed and attributed to the dominance of end conditions over rather large spans of the cylinder (aspect ratios up to 100).

In very large wind farms, the vertical interaction with the atmospheric boundary layer plays an important role, i.e. the total energy extraction is governed by the vertical transport of kinetic energy from higher regions in the boundary layer towards the turbine level. In the current study, we investigate optimal control of wind-farm boundary layers, considering the individual wind turbines as flow actuators, whose energy extraction can be dynamically regulated in time so as to optimally influence the flow field and the vertical energy transport. To this end, we use large-eddy simulations of a fully developed pressure-driven wind-farm boundary layer in a receding-horizon optimal control framework. For the optimization of the wind-turbine controls, a conjugate-gradient optimization method is used in combination with adjoint large-eddy simulations for the determination of the gradients of the cost functional. In a first control study, wind-farm energy extraction is optimized in an aligned wind farm. Results are accumulated over one hour of operation. We find that the energy extraction is increased by 16 % compared to the uncontrolled reference. This is directly related to an increase of the vertical fluxes of energy towards the wind turbines, and vertical shear stresses increase considerably. A further analysis, decomposing the total stresses into dispersive and Reynolds stresses, shows that the dispersive stresses increase drastically, and that the Reynolds stresses decrease on average, but increase in the wake region, leading to better wake recovery. We further observe also that turbulent dissipation levels in the boundary layer increase, and overall the outer layer of the boundary layer enters into a transient decelerating regime, while the inner layer and the turbine region attain a new statistically steady equilibrium within approximately one wind-farm through-flow time. Two additional optimal control cases study penalization of turbulent dissipation. For the current wind-farm geometry, it is found that the ratio between wind-farm energy extraction and turbulent boundary-layer dissipation remains roughly around 70 %, but can be slightly increased by a few per cent by penalizing the dissipation in the optimization objective. For a pressure-driven boundary layer in equilibrium, we estimate that such a shift can lead to an increase in wind-farm energy extraction of 6 %.

Direct numerical simulations of decaying high-Reynolds-number turbulence are presented at resolutions up to 8002 for general periodic flows and 20482 for periodic flows with large-scale symmetries. For turbulence initially excited at large scales, we characterize a transition of the inertial energy-spectrum exponent from n ≈ − 4 at early times to n ≈ − 3 when the turbulence becomes more mature. In physical space, the first regime is associated with isolated vorticity-gradient sheets, as predicted by Saffman (1971). The second regime, which is essentially statistical, corresponds to an enstrophy cascade (Kraichnan 1967; Batchelor 1969) and reflects the formation of layers resulting from the packing of vorticity-gradient sheets. In addition to these small-scale structures, the simulation displays vorticity macro-eddies which will survive long after the vorticity-gradient layers have been dissipated (McWilliams 1984). We validate the linear description of two-dimensional turbulence suggested by Weiss (1981), which predicts that coherent vortices will survive in regions where vorticity dominates strain, while vorticity-gradient sheets will be formed in regions where strain dominates. We show that this analysis remains valid even after vorticity-gradient sheets have been formed.

Digital particle image velocimetry (DPIV) measurements of the velocity field under breaking waves in the laboratory are presented. The region of turbulent fluid directly generated by breaking is too large to be imaged in one video frame and so an ensemble-averaged representation of the flow is built up from a mosaic of image frames. It is found that breaking generates at least one coherent vortex that slowly propagates downstream at a speed consistent with the velocity induced by its image in the free surface. Both the kinetic energy of the flow and the vorticity decay approximately as t−1. The Reynolds stress of the turbulence also decays as t−1 and is, within the accuracy of the measurements, everywhere negative, consistent with downward transport of streamwise momentum. Estimates of the mometum flux from waves to currents based on the measurements of the Reynolds stress are consistent with earlier estimates. The implications of the measurements for breaking in the field are discussed. Based on geometrical optics and wave action conservation, we suggest that the presence of the breaking-induced vortex provides an explanation for the suppression of short waves by breaking. Finally, in Appendices, estimates of the majority of the terms in the turbulent kinetic energy budget are presented at an early stage in the evolution of the turbulence, and comparisons with independent acoustical measurements of breaking are presented.

A planar liquid layer is bounded below by a rigid plate and above by an interface with a passive gas. A steady shear flow is set up by imposing a temperature gradient along the layer and driving the motion by thermocapillarity. This dynamic state is susceptible to surface-wave instabilities that couple the interfacial deflection to the underlying shear flow. These instabilities are found to be directly related to the two-dimensional waves on an isothermal layer subject to wind shear as described by Miles and by Smith & Davis. Hence the surface-tension gradients are important only in that they drive the basic shear flow. The surface-wave stability characteristics for liquid layers with and without return-flow profiles are presented, and special attention is paid to long-wave instabilities. Comparisons are made with available experimental observations.

This paper presents a comparison between theory and experiment for the in-line forces on cylinders of general cross-section in planar oscillatory flows of small amplitude. The theoretical analysis evaluates corrections to the standard inviscid inertial force at low Keulegan-Carpenter numbers which arise from the presence of viscous laminar boundary layers and from the development of vortex shedding. The boundary-layer contribution due to both skin friction and displacement effects is calculated to first order in the Stokes parameter β−½. The contribution to the in-line force from separation and vortex shedding, for which the results presented only apply to sharp-edged bodies, is taken from previous work on vortex shedding from isolated edges using the discrete vortex modelling technique. The resulting force has components both in phase with the fluid acceleration (inertia) and in phase with the velocity (drag).

The theoretical results are compared to measurements taken in a [xcup ]-tube water channel on a number of cylinders of different cross-section including circular cylinders and sharp-edged sections. The comparisons suggest that the theory is valid for Keulegan–Carpenter numbers below about 3 and for moderately high values of the β parameter.

Known structures and self-sustaining mechanisms of wall turbulence are reviewed and explored in the context of the scale interactions implied by the nonlinear advective term in the Navier–Stokes equations. The viewpoint is shaped by the systems approach provided by the resolvent framework for wall turbulence proposed by McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382), in which the nonlinearity is interpreted as providing the forcing to the linear Navier–Stokes operator (the resolvent). Elements of the structure of wall turbulence that can be uncovered as the treatment of the nonlinearity ranges from data-informed approximation to analysis of exact solutions of the Navier–Stokes equations (so-called exact coherent states) are discussed. The article concludes with an outline of the feasibility of extending this kind of approach to high-Reynolds-number wall turbulence in canonical flows and beyond.

We investigate the flow induced by a localized insoluble surfactant on a thin film. This problem is intended to model the behaviour of the lung's thin-film lining after an aerosol droplet lands on its surface. The surfactant-induced surface-tension gradients drive convection (Marangoni convection) within the film, disrupting the film surface and causing the surfactant to spread. The surfactant may also spread on the film's surface by surface diffusion without inducing convection. Gravity provides a restoring force that decreases film disturbances.

Lubrication theory is employed to derive equations that describe the evolution of the film thickness and surfactant concentration. A nonlinear surface-tension equation of state describes the relationship between the surfactant concentration and the surface tension. Solutions of the evolution equations are found numerically using the method of lines and analytically under limiting cases of small and large surface diffusivity. The results elucidate the behaviour of the thin-film/surfactant system.

We find that surface-tension-induced convection creates film disturbances that increase the film thickness near the surfactant's leading edge, and thins the film in the central region. Surface diffusion causes more rapid spreading of the surfactant, and decreases the film disturbances. Gravity decreases the film disturbances by creating bi-directional flow in the form of a ring vortex. This behaviour may have implications for the delivery of medications or toxins by aerosol inhalation.

A better understanding of the transition process in open flows can be obtained through identification of the possible asymptotic response states in the flow. In the present work, the asymptotic states in laminar wakes behind circular cylinders at low supercritical Reynolds numbers are investigated. Direct numerical simulation of the flow is performed, using spectral-element techniques. Naturally produced wakes, and periodically forced wakes are considered separately.

It is shown that, in the absence of external forcing, a periodic state is obtained, the frequency of which is selected by the absolute instability of the time-average flow. The non-dimensional frequency of the vortex street (Strouhal number) is a continuous function of the Reynolds number. In periodically forced wakes, however, non-periodic states are also possible, resulting from the bifurcation of the natural periodic state. The response of forced wakes can be characterized as: (i) lock-in, if the dominant frequency in the wake equals the excitation frequency, or (ii) non-lock-in, when the dominant frequency in the wake equals the Strouhal frequency. Both types of response can be periodic or quasi-periodic, depending on the combination of the amplitude and frequency of the forcing. At the boundary separating the two types of response transitional states develop, which are found to exhibit a low-order chaotic behaviour. Finally, all states resulting from the bifurcation of the natural state can be represented in a two-parameter space inside ‘resonant horn’ type of regions.

It is shown that the approximation of a thin thermal boundary layer gives an accurate description of the growth of spherical vapour bubbles in a superheated liquid except for very small superheats. If the further approximations of a linear variation of vapour pressure with temperature and of constant physical properties are made, then scaled variables can be introduced which describe the growth under any conditions. This scaled description is not valid during the early, surface-tension dominated, portion of the growth. The rate of bubble growth for large superheats is somewhat overestimated in the intermediate stage in which both inertial and thermal effects play a role. This overestimate does not lead to a serious error in the radius-time behaviour for ranges of practical interest. The asymptotic, or thermally controlled, stage of growth is accurately described by the scaled formulation.

An experimental investigation into the structure of a supersonic jet impinging on a large plate is presented. Digital particle image velocimetry (DPIV), shadowgraph photography and acoustic measurements are used to understand the relationship between the unsteady jet structure and the production of tones for nozzle-to-plate spacings between 1 and 5 nozzle exit diameters at a nozzle–pressure ratio equal to 4. Results indicate that the instability of the jet depends on the location of the plate in the shock cell structure of the corresponding free jet and the strength of the standoff shock wave, rather than on the occurrence of recirculation zones in the impingement region. Phase-locked studies show streamwise displacements of the stand-off shock wave, a moving recirculation zone in the subsonic flow in front of the plate, and significant oscillations of both the compression and expansion regions in the peripheral supersonic flow when tones are produced. Sound is shown to be generated by periodic pulsing of the wall jet boundary resulting from periodic motion of the flow in the impingement and near-wall regions of the flow.

Recent work on single-bubble sonoluminescence (SBSL) has shown that many features of this phenomenon, especially the dependence of SBSL intensity and stability on experimental parameters, can be explained within a hydrodynamic approach. More specifically, many important properties can be derived from an analysis of bubble wall dynamics. This dynamics is conveniently described by the Rayleigh–Plesset (RP) equation. Here we derive analytical approximations for RP dynamics and subsequent analytical laws for parameter dependences. These results include (i) an expression for the onset threshold of SL, (ii) an analytical explanation of the transition from diffusively unstable to stable equilibria for the bubble ambient radius (unstable and stable sonoluminescence), and (iii) a detailed understanding of the resonance structure of the RP equation. It is found that the threshold for SL emission is shifted to larger bubble radii and larger driving pressures if surface tension is increased, whereas even a considerable change in liquid viscosity leaves this threshold virtually unaltered. As an enhanced viscosity stabilizes the bubbles to surface oscillations, we conclude that the ideal liquid for violently collapsing, surface-stable SL bubbles should have small surface tension and large viscosity, although too large viscosity (ηl[ges ]40ηwater) will again preclude collapses.

The flow at the leading edge and in the interior of a rectangular porous obstruction is described through experiments and scaling. The porous obstruction consists of an emergent, rectangular array of cylinders in shallow flow, a configuration that mimics aquatic vegetation. The main features of the flow depend upon the non-dimensional canopy flow-blockage, which is a function of the obstruction width and porosity. For the ranges of canopy flow-blockage tested in this paper, the fluid decelerates upstream of the obstruction over a length scale proportional to the array width. For high flow-blockage, the interior adjustment length within the porous obstruction is set by the array width. For low flow-blockage, the array's frontal area per unit volume sets the interior adjustment length. Downstream of the adjustment regions, the interior velocity is governed by a balance between the lateral divergence of the turbulent stress and canopy drag, or by a balance between the pressure gradient and canopy drag, depending on the lateral penetration into the array of Kelvin–Helmholtz (KH) vortices, which is set by the non-dimensional canopy flow-blockage. For a porous obstruction with two stream-parallel edges, the KH vortex streets along the two edges are in communication across the width of the array: a phenomenon that results in cross-array vortex organization, which significantly enhances the vortex strength and creates significant lateral transport within the porous obstruction.

A quantitative study of some nonlinear aspects of drop-shape oscillations in a liquid–liquid system has been completed. The results suggest a soft nonlinearity in the fundamental resonant mode frequency as the oscillation amplitude is increased. Indications of an increase in the rate of decay have also been obtained. A study of the internal flow fields has revealed patterns of circulation not present at low amplitude.

In order to predict response and wake modes for elastically mounted circular cylinders in a fluid flow, we employ controlled-vibration experiments, comprised of prescribed transverse vibration of a cylinder in the flow, over a wide regime of amplitude and frequency. A key to this study is the compilation of high-resolution contour plots of fluid force, in the plane of normalized amplitude and wavelength. With such resolution, we are able to discover discontinuities in the force and phase contours, which enable us to clearly identify boundaries separating different fluid-forcing regimes. These appear remarkably similar to boundaries separating different vortex-formation modes in the map of regimes by Williamson & Roshko (J. Fluids Struct., vol. 2, 1988, pp. 355–381). Vorticity measurements exhibit the 2S, 2P and P + S vortex modes, as well as a regime in which the vortex formation is not synchronized with the body vibration. By employing such fine-resolution data, we discover a high-amplitude regime in which two vortex-formation modes overlap. Associated with this overlap regime, we identify a new distinct mode of vortex formation comprised of two pairs of vortices formed per cycle, where the secondary vortex in each pair is much weaker than the primary vortex. This vortex mode, which we define as the 2POVERLAP mode (2PO), is significant because it is responsible for generating the peak resonant response of the body. We find that the wake can switch intermittently between the 2P and 2PO modes, even as the cylinder is vibrating with constant amplitude and frequency. By examining the energy transfer from fluid to body motion, we predict a free-vibration response which agrees closely with measurements for an elastically mounted cylinder. In this work, we introduce the concept of an ‘energy portrait’, which is a plot of the energy transfer into the body motion and the energy dissipated by damping, as a function of normalized amplitude. Such a plot allows us to identify stable and unstable amplitude-response solutions, dependent on the rate of change of net energy transfer with amplitude (the sign of dE*/dA*). Our energy portraits show how the vibration system may exhibit a hysteretic mode transition or intermittent mode switching, both of which correspond with such phenomena measured from free vibration. Finally, we define the complete regime in the amplitude–wavelength plane in which free vibration may exist, which requires not only a periodic component of positive excitation but also stability of the equilibrium solutions.

The design of fluids management processes in the low-gravity environment of space requires an accurate description of capillarity-controlled flow in containers. Here we consider the spontaneous redistribution of fluid along an interior corner of a container due to capillary forces. The analytical portion of the work presents an asymptotic formulation in the limit of a slender fluid column, slight surface curvature along the flow direction z, small inertia, and low gravity. The scaling introduced explicitly accounts for much of the variation of flow resistance due to geometry and so the effects of corner geometry can be distinguished from those of surface curvature. For the special cases of a constant height boundary condition and a constant flow condition, the similarity solutions yield that the length of the fluid column increases as t1/2 and t3/5, respectively. In the experimental portion of the work, measurements from a 2.2 s drop tower are reported. An extensive data set, collected over a previously unexplored range of flow parameters, includes estimates of repeatability and accuracy, the role of inertia and column slenderness, and the effects of corner angle, container geometry, and fluid properties. At short times, the fluid is governed by inertia (t[lsim ]tLc). Afterwards, an intermediate regime (tLc[lsim ]t[lsim ] tH) can be shown to be modelled by a constant-flow-like similarity solution. For t[ges ]tH it is found that there exists a location zH at which the interface height remains constant at a value h(zH, t)=H which can be shown to be well predicted. Comprehensive comparison is made between the analysis and measurements using the constant height boundary condition. As time increases, it is found that the constant height similarity solution describes the flow over a lengthening interval which extends from the origin to the invariant tip solution. For t[Gt ]tH, the constant height solution describes the entire flow domain. A formulation applicable throughout the container (not just in corners) is presented in the limit of long times.

This paper uses numerical simulation to analyse the effects of uniform rotation on homogeneous turbulence. Both large-eddy and full simulations were made. The results indicate that the predominant effect of rotation is to decrease the rate of dissipation of the turbulence and increase the lengthscales, especially those along the axis of rotation. These effects are a consequence of the reduction, due to the generation of inertial waves, of the net energy transfer from large eddies to small ones. Experiments are also influenced by a more complicated interaction between the rotation and the wakes of the turbulence-generating grid which modifies the nominal initial conditions in the experiment. The latter effect is accounted for in simulations by modifying the initial conditions. Finally, a two-equation model is proposed that accounts for the effects of rotation and is able to reproduce the experimental decay of the turbulent kinetic energy.

Results of laboratory experiments are used to determine the morphology and the ascent rate of growing buoyant plumes in a homogeneous, viscous fluid. The plumes were formed by injecting a glucose solution through a small orifice into another glucose solution of different density. Two classes of creeping (low-Reynolds-number) plumes are investigated: (i) diapiric plumes, for which the plume viscosity is approximately equal to the ambient-fluid viscosity, and (ii) cavity plumes, for which the plume fluid is much less viscous than the ambient fluid. Fully developed diapirs consist of a tapered cylindrical stem capped by a mushroom-shaped vortex at its leading edge. Fully developed cavity plumes consist of a nearly spherical chamber connected to the source by a narrow umbilical conduit. It is observed that the ascent velocity of cavity plumes increases with time as t⅖. The ascent velocity of diapirs is found to be proportional to the terminal velocity of a cylinder moving parallel to its axis. The presence of pre-existing conduits alters the morphology of cavity plumes and greatly increases their ascent rate. Fossil conduits act as plume guides by offering low-resistance ascent paths. Finally, a series of experiments have been made on the interaction between cavity plumes and a large-scale background circulation. A low-viscosity plume generated by a source towed steadily through a highly viscous fluid breaks into a chain of regularly spaced, individual cavities, as first demonstrated by Skilbeck & Whitehead. The cavities ascend as an inclined linear array of Stokes droplets. Dimensional analysis is used to derive scaling laws for the cavity volumes and their replication rate in terms of the source parameters and the tow speed. The qualitative results from these experiments generally lend support to the hypothesis that buoyant plumes in the Earth's mantle are the source of hot-spot volcanism. In particular the experiments suggest an explanation for the observation that hot spots remain nearly fixed in the presence of mantle convection.

The manner by which external vortical disturbances penetrate the laminar boundary layer and induce transition is explored. Linear theory suggests that the well-known Klebanoff mode precursor to transition can be understood as a superposition of Squire continuous modes. Shear sheltering influences the ability of free-stream disturbances to generate a packet of Squire modes. A coupling coefficient between continuous spectrum spectrum Orr–Sommerfeld and Squire modes is used to characterize the interaction. Full numerical simulations with prescribed modes at the inlet substantiate this approach. With two weakly coupled modes at the inlet, the boundary layer is little perturbed; with two strongly coupled modes, Klebanoff modes are produced; with one strongly coupled and one weakly coupled high-frequency mode, the complete transition process is simulated.

Inspired by the correlation between the propulsion efficiency of a flapping foil propeller and stability of the wake behind it (which leads to the optimal Strouhal number for propulsion), we numerically simulated a heaving/pitching foil in energy harvesting regime, and investigated the relation between wake stability and the energy harvesting efficiency. The base flow is computed using a Navier–Stokes algorithm and the stability analysis is performed via the Orr–Sommerfeld equation. The wake is found to be convectively unstable and the frequency of the most unstable mode fw is determined. The case when fw ~ f coincides with maximum energy harvesting efficiency of the system (f is the frequency of foil oscillation), suggesting that flow energy extraction is closely related to efficient evolution of the wake. This occurs at a frequency of f ~ 0.15 (f is normalized by the chord length and the flow speed), under the constraint that there is significant vortex shedding from the leading edge at sufficiently large effective angles of attack. Indeed, this ‘foil–wake resonance’ is usually associated with multi-vortex shedding from the leading edge. Furthermore, detailed examination of energy extractions from the heaving and the pitching motions indicates that near the optimal performance point the average energy extraction from the pitching motion is close to zero. This suggests the feasibility of achieving high-efficient energy harvesting through a simple fully passive system we proposed earlier in which no activation is needed.