Measurements have been made of the mean and fluctuating pressure distributions on long circular cylinders, having smooth and rough surfaces, at Reynolds numbers of 1·11 × 105 and 2·35 × 105 in both uniform and turbulent streams. The presence of free-stream turbulence a t these Reynolds numbers was found to suppress coherent vortex shedding on the smooth cylinder and give rise to a complex pressure field in which the mean pressure distribution was almost independent of Reynolds number over the small range of Reynolds numbers tested. The pressure distributions on the rough cylinder were found to be completely different in uniform and turbulent streams; the presence of turbulence gave rise to an increase in the level of vortex shedding energy, and produced mean pressure distributions similar to those obtained on smooth cylinders at Reynolds numbers of the order of 107.

We establish through numerical simulation conditions for optimal undulatory propulsion for a single fish, and for a pair of hydrodynamically interacting fish, accounting for linear and angular recoil. We first employ systematic two-dimensional (2-D) simulations to identify conditions for minimal propulsive power of a self-propelled fish, and continue with targeted 3-D simulations for a danio-like fish; all at Reynolds number 5000. We find that the Strouhal number, phase angle between heave and pitch at the trailing edge, and angle of attack are principal parameters. For 2-D simulations, imposing a deformation based on measured displacement for carangiform swimming provides, at best, efficiency of 35 %, which increases to 50 % for an optimized motion; for a 3-D fish, the efficiency increases from 22 % to 34 %. Indeed, angular recoil has significant impact on efficiency, and optimized body bending requires maximum bending amplitude upstream of the trailing edge. Next, we turn to 2-D simulation of two hydrodynamically interacting fish. We find that the upstream fish benefits energetically only for small distances. In contrast, the downstream fish can benefit at any position that allows interaction with the upstream wake, provided its body motion is timed appropriately with respect to the oncoming vortices. For an in-line configuration, one body length apart, the efficiency of the downstream fish can increase from 50 % to 60 %; for an offset arrangement it can reach 80 %. This proves that in groups of fish, energy savings can be achieved for downstream fish through interaction with oncoming vortices, even when the downstream fish lies directly inside the jet-like flow of an upstream fish.

The interaction of two laser-induced bubbles in bulk water is investigated. The strength and direction of the emerging liquid jets can be controlled by adjusting the relative bubble positions, the time difference between bubble generation, and the laser pulse energies determining the bubble sizes. Experimental and numerical studies are performed for millimetre-sized bubble pairs. Taking bubbles of equal energy, a maximum jet velocity is found for close anti-phase bubbles, i.e. when the second bubble is produced at the maximum volume of the first one and the bubble walls are almost touching and not merging. Under these conditions, one bubble produces a fast jet with a peak velocity of about $150~\text{m}~\text{s}^{-1}$ that reaches a distance into the surrounding liquid of at least three times the maximum bubble radius. Collapse of the other bubble results in a slow jet of large mass that rapidly converts into a ring vortex. Correspondingly, the interaction with adjacent structures is dominated either by localized jet impact or by shear stresses extending over a larger area. Furthermore, interactions between micrometre-sized bubble pairs are investigated numerically to understand and predict how the effects of the physical parameters on bubble dynamics would change when the bubbles become smaller. The results are discussed with respect to micropumping and opto-injection.

The mass transfer mechanism across a sheared air–water interface without bubble entrainment due to wave breaking was experimentally investigated in terms of the turbulence structure of the organized motions in the interfacial region in a wind-wave tank. The transfer velocity of the carbon dioxide (CO2) on the water side was measured through reaeration experiments of CO2, and the fluid velocities in the air and water flows were measured using both a hot-wire anemometer and a laser-Doppler velocimeter. The results show that the mass transfer across a sheared air–water interface is more intensively promoted in wind shear, compared to an unsheared interface. However, the effect of the wind shear on the mass transfer tends to saturate in the high-shear region in the present wind-wave tank, where the increasing rate of mass transfer velocity with the wind shear decreases rapidly. The effect of the wind shear on the mass transfer can be well explained on the basis of the turbulence structure near the air–water interface. That is, surface-renewal eddies are induced on the water side through the high wind shear on the air–water interface by the strong organized motion generated in the air flow above the interface, and the renewal eddies control the mass transfer across a sheared interface. The mass transfer velocity is correlated with the frequency of the appearance of the surface-renewal eddies, as it is in open-channel flows with unsheared interfaces, and it increases approximately in proportion to the root of the surface-renewal frequency. The surface-renewal frequency increases with increasing the wind shear, but for high shear the rate of increase slows. This results in the saturated effect of the wind shear on the mass transfer in the high-shear region in the present wind-wave tank. The mass transfer velocity can be well estimated by the surface-renewal eddy-cell model based on the concept of the time fraction when the surface renewal occurs.

The viscous, creeping flow through a cylindrical tube of a liquid, which contains rigid, spherical particles, is investigated analytically. The spheres are located on the axis of the cylinder and are equally spaced. Solutions are derived for particles in motion and fixed, with and without fluid discharge. Numerical results are presented for the drag on each sphere and the mean pressure drop for a wide range of sizes and spacings of the spheres. The study is motivated by possible application to blood flow in capillaries, where red blood cells represent particles of the same order of magnitude as the diameter of the capillary itself. The results may also be of interest in other applications, such as sedimentation and fluidized beds. It is shown that there is little interaction between particles if the spacing is more than one tube diameter, and that the additional pressure drop over that for Poiseuille flow is less than 50% if the sphere diameter is less than 0·8 of the tube diameter.

The flow around a finite circular cylinder with a height-to-diameter ratio of 6 was studied using large eddy simulation (LES). The cylinder was mounted on a ground plane where a thin boundary layer with a thickness of only 7% of the cylinder's diameter was formed. The flow resulting from the present LES was used to present a detailed picture of both the instantaneous and the time-averaged flow. Among the results of the present study is the explanation for the formation of two time-averaged tip vortices from two pairs of instantaneous vortices. The main recirculation region on the cylinder's top was found to be one vortex with both legs attached to the top surface, in agreement with some previous investigations. The mean flow in the near wake region was found to contain two horseshoe vortices, one large arch vortex stretching from the cylinder's top and another vortex resulting from the flow moving down below the rear stagnation point on the cylinder. The instantaneous horseshoe vortex near the ground was found to be unsteady and changed shape over time. The time-averaged flow showed three complete horseshoe vortices, in agreement with previous knowledge. The downwash process above the free end into the near wake occurs at a rather constant angle with respect to the plane normal to the streamwise direction. The same angle was preserved in the inclination of the arch vortex in the time-averaged flow.

We study the Reynolds-number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier–Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behaviour with Reynolds number in the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it is shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. Under the practical assumption that the mean velocity consists of a logarithmic region, we identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_{\tau } \approx 1{0}^{3} {\unicode{x2013}} 1{0}^{10} $). Results from this low-rank model of the Navier–Stokes equations compare favourably with experimental results in the literature.

The two-dimensional flow around a rotating circular cylinder is studied at Re = 100. The instability mechanisms for the first and second shedding modes are analysed. The region in the flow with a role of ‘wavemaker’ in the excitation of the global instability is identified by considering the structural sensitivity of the unstable mode. This approach is compared with the analysis of the perturbation kinetic energy production, a classic approach in linear stability analysis. Multiple steady-state solutions are found at high rotation rates, explaining the quenching of the second shedding mode. Turning points in phase space are associated with the movement of the flow stagnation point. In addition, a method to examine which structural variation of the base flow has the largest impact on the instability features is proposed. This has relevant implications for the passive control of instabilities. Finally, numerical simulations of the flow are performed to verify that the structural sensitivity analysis is able to provide correct indications on where to position passive control devices, e.g. small obstacles, in order to suppress the shedding modes.

Micro-cavitation bubbles generated by ultrasound have wide and important applications in medical ultrasonics and sonochemistry. An approximate theory is developed for nonlinear and non-spherical bubbles in a compressible liquid by using the method of matched asymptotic expansions. The perturbation is performed to the second order in terms of a small parameter, the bubble-wall Mach number. The inner flow near the bubble can be approximated as incompressible at the first and second orders, leading to the use of Laplace's equation, whereas the outer flow far away from the bubble can be described by the linear wave equation, also for the first and second orders. Matching between the two expansions provides the model for the non-spherical bubble behaviour in a compressible fluid. A numerical model using the mixed Eulerian–Lagrangian method and a modified boundary integral method is used to obtain the evolving bubble shapes. The primary advantage of this method is its computational efficiency over using the wave equation throughout the fluid domain. The numerical model is validated against the Keller–Herring equation for spherical bubbles in weakly compressible liquids with excellent agreement being obtained for the bubble radius evolution up to the fourth oscillation. Numerical analyses are further performed for non-spherical oscillating acoustic bubbles. Bubble evolution and jet formation are simulated. Outputs also include the bubble volume, bubble displacement, Kelvin impulse and liquid jet tip velocity. Bubble behaviour is studied in terms of the wave frequency and amplitude. Particular attention is paid to the conditions if/when the bubble jet is formed and when the bubble becomes multiply connected, often forming a toroidal bubble. When subjected to a weak acoustic wave, bubble jets may develop at the two poles of the bubble surface after several cycles of oscillations. A resonant phenomenon occurs when the wave frequency is equal to the natural oscillation frequency of the bubble. When subjected to a strong acoustic wave, a vigorous liquid jet develops along the direction of wave propagation in only a few cycles of the acoustic wave.

We investigate experimentally the force generated by the unsteady vortex formation of low-aspect-ratio normal flat plates with one end free. The objective of this study is to determine the role of the free end, or tip, vortex. Understanding this simple case provides insight into flapping-wing propulsion, which involves the unsteady motion of low-aspect-ratio appendages. As a simple model of a propulsive half-stroke, we consider a rectangular normal flat plate undergoing a translating start-up motion in a towing tank. Digital particle image velocimetry is used to measure multiple perpendicular sections of the flow velocity and vorticity, in order to correlate vortex circulation with the measured plate force. The three-dimensional wake structure is captured using flow visualization. We show that the tip vortex produces a significant maximum in the plate force. Suppressing its formation results in a force minimum. Comparing plates of aspect ratio six and two, the flow is similar in terms of absolute distance from the tip, but evolves faster for aspect ratio two. The plate drag coefficient increases with decreasing aspect ratio.

Experiments have been carried out to study the three-dimensional characteristics of vortex shedding from a half-ellipse shape with a blunt trailing edge. In order to control the occurrence of vortex dislocations, the trailing edges of the models used were constructed with a series of periodic waves across their spans. Flow visualization was carried out in a water tunnel at a Reynolds number of 2500, based on trailing-edge thickness. A number of shedding modes were observed and the sequence of mode transitions recorded. Quantitative data were obtained from wind tunnel measurements performed at a Reynolds number of 40000. Two shedding frequencies were recorded with the higher frequency occurring at spanwise positions coinciding with minima in the chord. At these same positions the base pressure was lowest and the vortex formation length longest. Arguments are put forward to explain these observations. It is shown that the concept of a universal Strouhal number holds, even when the flow is three-dimensional. The spanwise variation in time-average base pressure is predicted using the estimated amount of time the flow spends at the two shedding frequencies.

The break-up mechanism of axisymmetric liquid sheets formed by the impingement of two co-axial jets has been examined. Three break-up regimes in the Weber number range from 100 to 3 × 104 are reported. In the first break-up regime, droplets are formed through successive mergings of liquid beads along the nearly circular periphery of the sheet. The formation of beads is caused by Rayleigh instability. In the transition regime, Taylor's cardioid wave pattern prevails in the first half of this regime, while the sheet begins to flap in the second half.

In the second break-up regime, antisymmetric waves on the sheet grow radially. A semi-empirical equation has been deduced to predict the break-up radius of the sheet. The motion of an axisymmetric vibrating membrane with radially decreasing thickness has been studied to include Helmholtz instability as an analogue of the wave motion of the expanding circular sheet. A distorted progressive wave equation has been solved by the WKBJ method to indicate the effect of cylindrical geometry. The calculated wave speed agrees fairly well with experimental data at low Weber numbers.

It is now well known that the turbulence structure of thin shear layers can be strongly affected by the application of extra rates of strain in addition to the shear velocity gradient. Examples of such extra strain rates include lateral divergence or convergence, and streamline curvature in the plane of the mean shear. The changes in Reynolds stress are an order of magnitude larger than would be expected from the explicit extra terms in the Reynolds-stress transport equations, and therefore an order of magnitude larger than predicted by conventional calculation methods. In the present paper, one of a series on ‘complex’ turbulent flows, we show that bulk compression or dilatation (i.e. an extra strain rate div U) also appears to affect turbulent shear layers, typical values of Reynolds stress being increased by compression and decreased by dilatation. The fractional change in Reynolds stress is an order of magnitude larger than the fractional change in volume of a fluid element. The physical mechanism is probably analogous to that responsible for the large effects of divergence or convergence in incompressible flow. Because the phenomenon seems to be of great practical importance we discuss it in the context of engineering calculation methods. An empirical correction formula, analogous to those used to allow for divergence or curvature effects, greatly reduces the large discrepancies found between recent experiments on supersonic boundary layers and calculations by conventional extensions of successful incompressible-flow methods.

The diffraction of regular waves by a vertical circular cylinder in finite depth water is considered, within the frame of potential theory. The wave slope kA is assumed to be small so that successive boundary value problems at orders kA, k2A2, and k3A3 can be formulated. Here we focus on the third-order (k3A3) problem but restrict ourselves to the triple-frequency component of the diffraction potential. The method of resolution is based on eigenfunction expansions and on the integral equation technique with the classical Green function expressed in cylindrical coordinates. Third-order (triple-frequency) loads are calculated and compared with experimental measurements and approximate methods based on long-wave theories.

As part of a general investigation of complex turbulent flows, extensive one-point measurements have been made of the turbulence structure of the mixing layer bounding a normally impinging plane jet with an irrotational core. The ratio of shear-layer thickness to streamline radius of curvature reaches a maximum of about 0.2, the sense of the curvature being stabilizing. Downstream of the impingement region the shear layer returns asymptotically to being a classical plane mixing layer. The most striking feature of the results is that the return is not monotonic: after decreasing in the region of stabilizing curvature, the Reynolds stresses, triple products, energy dissipation rate and other turbulence quantities overshoot the plane-layer values before finally decreasing. Some conclusions are drawn about the nature of the turbulent transport of Reynolds stress, and about the representation of this and other processes in calculation methods for complex turbulent flows. An incidental result of the work is a comprehensive set of measurements in a plane mixing layer.

High-frequency fluctuations of temperature and of longitudinal and vertical velocity components have been measured with high-resolution probes in order to test the local-isotropy assumption. The simultaneous measurements of u’, w’, θ’ and the measurements in two space points with various separations in either the longitudinal or transverse directions were made in the large boundary layer (Rλ = 616) of the I.M.S.T. Air-Sea Interaction Simulation Tunnel. There is consistent evidence that the local-isotropy assumption is satisfied by the velocity field at all scales smaller than twenty times the Kolmogorov microscale (η ≈ 0.27 × 10−3 m), i.e. in the dissipative range of scales but not in the expected inertial subrange. The direct comparisons of the lateral and longitudinal temperature autocorrelation and structure functions show that the temperature field does not verify the isotropy assumption at all scales greater thanor equal to 37 and presumably at even smaller scales. This is confirmed by the study of the temperature-increment skewness and flatness factors. The spectral distribution -of the non-zero derivative skewness (S(θ) = +0.9) shows that it is essentially contributed by those scales for which the dynamic field satisfies isotropy.

We develop a time-dependent numerical algorithm, using a boundary-integral approach, to investigate fingering in Hele-Shaw cells. Starting from a sinusoidal variation in the initial interface, stable fingers quickly form for a wide range of the dimensionless surface-tension parameter. For very low values of the parameter, the incipient finger bifurcates. The stable fingers are clearly the same as those obtained by McLean & Saffman (1981) using a steady-state algorithm. These steady-state solutions were found to be linearly unstable. We resolve this apparent discrepancy regarding stability by tracing the fate of small disturbances placed on and about the finger tip. We show that some small disturbances do, indeed, grow initially; however, they reach a maximum amplitude and decay as they convect backward from the tip of the finger to regions where stabilizing surface tension is the major physical force. Relatively large imposed disturbances, on the other hand, cause a finger to bifurcate; the critical disturbance amplitude decreases as the surface tension is reduced

Several conditional-sampling techniques are applied to data bases generated by large-eddy and direct numerical simulations. It is shown that the bursting process is associated with well-organized vortical structures inclined at about 45° to the wall. These vortical structures are identified by examining the vortex lines of three-dimensional instantaneous and ensemble-averaged vorticity fields. Two distinct horseshoe-shaped vortical structures corresponding to the sweep and ejection events are detected. These vortical structures are associated with high Reynolds shear stress and hence make a significant contribution to turbulent-energy production.

We present results from an experimental and numerical investigation into the structure of vortex patterns and the dynamics of their interactions for the incompressible flow in the near field of a round coaxial jet issuing into a quiescent ambient fluid. A two-colour planar laser-induced-fluorescence technique is used to document the flow field via still photographs and ciné sequences over a limited range of parameters. We examine the effects of varying the velocity ratio as well as the absolute velocities of the two coaxial streams for equal densities and for a single area ratio. Results show that a variety of widely differing near-field vortex patterns can arise, with very different interaction dynamics, which can depend both on the velocity ratio and on the absolute velocities of the two streams. The observed vortex structures and their dynamics are interpreted in terms of the instability of the initially cylindrical and concentric vorticity layers separating each of the fluid streams, and their subsequent rollup to form wake-like or shear-layer-like vortices. Our results show that in addition to the velocity jump across each of these vorticity layers, an accounting of the layer thicknesses and the wake defect within each layer can be essential to understanding the resulting near-field structure that occurs. Ensuing dynamical interactions between the vortices formed from each layer can produce a strong coupling between the development of the two layers. These resulting vortex structures and interaction dynamics are also seen to produce widely differing mixing patterns in the jet near field.

This paper describes an experimental study of turbulent boundary layers over two-dimensional spanwise groove and three-dimensional sandgrain roughnesses in the ‘transition regime’ between hydraulically smooth and fully rough conditions. Mean-flow measurements show that a state of kinematic near-self-preservation is also reached by sandgrain roughness and not just by d-type grooved roughness alone as commonly believed; sandgrain roughness simply requires an order-of-magnitude-longer length to reach such a state. The two roughness Reynolds numbers demarcating the boundaries of the transition regime of k-type roughnesses are found to decrease with increasing roughness-element spanwise aspect ratio (span/height). A more important role of the upper-Reynolds-number limit of the transition regime in the drag behaviour is indicated. The two Reynolds-number limits of the transition regime correlate with the two critical Reynolds numbers that describe the stability of the vortex-shedding process existing behind a similar but isolated roughness element lying submerged in an otherwise laminar boundary layer. The results provide a guideline for reducing k-type rough-wall drag by lowering the spanwise aspect ratio of the roughness elements. The vortex-shedding process in rough-wall turbulent boundary layers is described by the stability parameter $U\tau (\overline{T}/\nu)^{\frac{1}{2}}$ whose value is the same for all roughnesses examined herein; here Uτ is the friction velocity, $\overline{T}$ is the mean time period of vortex shedding and v is the kinematic viscosity of the fluid.