The properties of turbulence subgrid-scale stresses are studied using experimental data in the far field of a round jet, at a Reynolds number of Rλ ≈ 310. Measurements are performed using two-dimensional particle displacement velocimetry. Three elements of the subgrid-scale stress tensor are calculated using planar filtering of the data. Using a priori testing, eddy-viscosity closures are shown to display very little correlation with the real stresses, in accord with earlier findings based on direct numerical simulations at lower Reynolds numbers. Detailed analysis of subgrid energy fluxes and of the velocity field decomposed into logarithmic bands leads to a new similarity subgrid-scale model. It is based on the ‘resolved stress’ tensor Lij, which is obtained by filtering products of resolved velocities at a scale equal to twice the grid scale. The correlation coefficient of this model with the real stress is shown to be substantially higher than that of the eddy-viscosity closures. It is shown that mixed models display similar levels of correlation. During the a priori test, care is taken to only employ resolved data in a fashion that is consistent with the information that would be available during large-eddy simulation. The influence of the filter shape on the correlation is documented in detail, and the model is compared to the original similarity model of Bardina et al. (1980). A relationship between Lij and a nonlinear subgrid-scale model is established. In order to control the amount of kinetic energy backscatter, which could potentially lead to numerical instability, an ad hoc weighting function that depends on the alignment between Lij and the strain-rate tensor, is introduced. A ‘dynamic’ version of the model is shown, based on the data, to allow a self-consistent determination of the coefficient. In addition, all tensor elements of the model are shown to display the correct scaling with normal distance near a solid boundary.

Although there are a great many papers dedicated to the problem of a cylinder vibrating transverse to a fluid flow ($Y$-motion), there are almost no papers on the more practical case of vortex-induced vibration in two degrees of freedom ($X,Y$ motion) where the mass and natural frequencies are precisely the same in both $X$- and $Y$-directions. We have designed the present pendulum apparatus to achieve both of these criteria. Even down to the low mass ratios, where $m^*\,{=}\,6$, it is remarkable that the freedom to oscillate in-line with the flow affects the transverse vibration surprisingly little. The same response branches, peak amplitudes, and vortex wake modes are found for both $Y$-only and $X,Y$ motion. There is, however, a dramatic change in the fluid–structure interactions when mass ratios are reduced below $m^*\,{=}\,6$. A new amplitude response branch with significant streamwise motion appears, in what we call the ‘super-upper’ branch, yielding massive amplitudes of 3 diameters peak-to-peak ($A^*_Y \,{\sim}\, 1.5$). We discover a corresponding periodic vortex wake mode, comprising a triplet of vortices being formed in each half-cycle, in what we define as a ‘2T’ mode. We qualitatively interpret the principal vortex dynamics and vortex forces which yield a positive rate of energy transfer ($\dot{e}_V$) causing the body vibration, using the following simple equation: \[\dot{e}_V = 2 {\Gamma}^* U^*_V \skew3\dot{Y}\] where $\Gamma^*$ is vortex strength, $U^*_V$ is the speed downstream of the dominant near-wake vorticity, and $\skew3\dot{Y}$ is the transverse velocity of the body. This simple approach suggests that the massive amplitude of vibration for the 2T mode is principally attributed to the energy transfer from the ‘third’ vortex of each triplet, which is not present in the lower-amplitude 2P mode. We also find two low-speed streamwise vibration modes, which is not unexpected, since they correspond to the first and second excitation modes of vibration for flexible cantilevers. By considering equations of motion for the two degrees of freedom, we find a critical mass, $m^*_{\hbox{\scriptsize\it crit}} \,{=}\, 0.52$, similar to recent $Y$-only studies, below which the large-amplitude vibrations persist to infinite flow velocity. We show that the critical mass $m^*_{\hbox{\scriptsize\it crit}}$ is the same for the $X$- and $Y$-directions, which ensures that the shapes of $X, Y$ trajectories can retain their form as the velocity becomes large. The extensive studies of vortex-induced vibration for $Y$-only body motions, built up over the last 35 years, remain of strong relevance to the case of two degrees of freedom, for $m^* \,{>}\, 6$. It is only for ‘small’ mass ratios, $m^* \,{<}\, 6$, that one observes a rather dramatic departure from previous results, which would suggest a possible modification to offshore design codes.

Thermal instabilities of a contained fluid are investigated for a fairly general class of problems in which the dynamics are dominated by the effects of rotation. In systems of constant depth in the direction of the axis of rotation the instability develops when the buoyancy forces suffice to overcome the stabilizing effects of thermal conduction and of viscous dissipation in the Ekman boundary layers. Owing to the Taylor–Proudman theorem, a slight gradient in depth exerts a strongly stabilizing influence. The theory is applied to describe the instability of the ‘lower symmetric régime’ in the rotating annulus experiments at high rotation rates. An example of geophysical relevance is the instability of a self-gravitating, internally heated, rotating fluid sphere. The results of the perturbation theory for this problem agree reasonably well with the results of an extension of the analysis by Roberts (1968).

Two rigid spheres of radii a and b are immersed in infinite fluid whose velocity at infinity is a linear function of position. No external force or couple acts on the spheres, and the effect of inertia forces on the motion of the fluid and the spheres is neglected. The purpose of the paper is to provide a systematic and explicit description of those aspects of the interaction between the two spheres that are relevant in a calculation of the mean stress in a suspension of spherical particles subjected to bulk deformation. The most relevant aspects are the relative velocity of the two sphere centres (V) and the force dipole strengths of the two spheres (S′ij, S″ij), as functions of the vector r separating the two centres.

It is shown that V, S′ij and S″ij depend linearly on the rate of strain at infinity and can be represented in terms of several scalar parameters which are functions of r/a and b/a alone. These scalar functions provide a framework for the expression of the many results previously obtained for particular linear ambient flows or for particular values of r/a or of b/a. Some new results are established for the asymptotic forms of the functions both for r/(a + b) [Gt ] 1 and for values of r − (a + b) small compared with a and b. A reasonably complete numerical description of the interaction of two rigid spheres of equal size is assembled, the main deficiency being accurate values of the scalar functions describing the force dipole strength of a sphere in the intermediate range of sphere separations.

In the case of steady simple shearing motion at infinity, some of the trajectories of one sphere centre relative to another are closed, a fact which has consequences for the rheological problem. These closed forms are described analytically, and also numerically in the case b/a = 1.

Plate tectonics provides a remarkably accurate kinematic description of the motion of the earth's crust but a fully dynamical theory requires an understanding of convection in the mantle. Thus the properties of plates and of the mantle must be related to a systematic study of convection. This paper reviews both the geophysical information and the fluid dynamics of convection in a Boussinesq fluid of infinite Prandtl number. Numerical experiments have been carried out on several simple two-dimensional models, in which convection is driven by imposed horizontal temperature gradients or else by heating either internally or from below. The results are presented and analysed in terms of simple physical models. Although the computations are highly idealized and omit variation of viscosity and other major features of mantle convection, they can be related to geophysical measurements. In particular, the external gravity field depends on changes in surface elevation; this suggests an observational means of investigating convection in the upper mantle.

The level-set approach is applied to a regime of premixed turbulent combustion where the Kolmogorov scale is smaller than the flame thickness. This regime is called the thin reaction zones regime. It is characterized by the condition that small eddies can penetrate into the preheat zone, but not into the reaction zone.

By considering the iso-scalar surface of the deficient-species mass fraction Y immediately ahead of the reaction zone a field equation for the scalar quantity G(x, t) is derived, which describes the location of the thin reaction zone. It resembles the level-set equation used in the corrugated flamelet regime, but the resulting propagation velocity s*L normal to the front is a fluctuating quantity and the curvature term is multiplied by the diffusivity of the deficient species rather than the Markstein diffusivity. It is shown that in the thin reaction zones regime diffusive effects are dominant and the contribution of s*L to the solution of the level-set equation is small.

In order to model turbulent premixed combustion an equation is used that contains only the leading-order terms of both regimes, the previously analysed corrugated flamelets regime and the thin reaction zones regime. That equation accounts for non-constant density but not for gas expansion effects within the flame front which are important in the corrugated flamelets regime. By splitting G into a mean and a fluctuation, equations for the Favre mean [Gtilde]and the variance [Gtilde]″2 are derived. These quantities describe the mean flame position and the turbulent flame brush thickness, respectively. The equation for [Gtilde]″2 is closed by considering two-point statistics. Scaling arguments are then used to derive a model equation for the flame surface area ratio [rhotilde]. The balance between production, kinematic restoration and dissipation in this equation leads to a quadratic equation for the turbulent burning velocity. Its solution shows the ‘bending’ behaviour of the turbulent to laminar burning velocity ratio sT/sL, plotted as a function of v′/sL. It is shown that the bending results from the transition from the corrugated amelets to the thin reaction zones regimes. This is equivalent to a transition from Damköhler's large-scale to his small-scale turbulence regime.

Early treatments of flames as gasdynamic discontinuities in a fluid flow are based on several hypotheses and/or on phenomenological assumptions. The simplest and earliest of such analyses, by Landau and by Darrieus prescribed the flame speed to be constant. Thus, in their analysis they ignored the structure of the flame, i.e. the details of chemical reactions, and transport processes. Employing this model to study the stability of a plane flame, they concluded that plane flames are unconditionally unstable. Yet plane flames are observed in the laboratory. To overcome this difficulty, others have attempted to improve on this model, generally through phenomenological assumptions to replace the assumption of constant velocity. In the present work we take flame structure into account and derive an equation for the propagation of the discontinuity surface for arbitrary flame shapes in general fluid flows. The structure of the flame is considered to consist of a boundary layer in which the chemical reactions occur, located inside another boundary layer in which transport processes dominate. We employ the method of matched asymptotic expansions to obtain an equation for the evolution of the shape and location of the flame front. Matching the boundary-layer solutions to the outer gasdynamic flow, we derive the appropriate jump conditions across the front. We also derive an equation for the vorticity produced in the flame, and briefly discuss the stability of a plane flame, obtaining corrections to the formula of Landau and Darrieus.

This paper attempts to emulate the great study by Goldstein (1929) ‘On the vortex wake of a screw propeller’, by looking for a dynamical theory of how another type of propulsion system has evolved towards ever higher performance. An ‘undulatory’ mode of animal propulsion in water is rather common among invertebrates, and this paper offers a preliminary quantitative analysis of how a series of modifications of that basic undulatory mode, found in the vertebrates (and especially in the fishes), tends to improve speed and hydromechanical efficiency.

Posterior lateral compression is the most important of these. It is studied first in ‘pure anguilliform’ (eel-like) motion of fishes whose posterior cross-sections are laterally compressed, although maintaining their depth (while the body tapers) by means of long continuous dorsal and ventral fins all the way to a vertical ‘trailing edge’. Lateral motion of such a cross-section produces a large and immediate exchange of momentum with a considerable ‘virtual mass’ of water near it.

In § 2, ‘elongated-body theory’ (an extended version of inviscid slender-body theory) is developed in detail for pure anguilliform motion and subjected to several careful checks and critical studies. Provided that longitudinal variation of cross-sectional properties is slow on a scale of the cross-sectional depth s (say, if the wavelength of significant harmonic components of that variation exceeds 5s), the basic approach is applicable and lateral water momentum per unit length is closely proportional to the square of the local cross-section depth.

The vertical trailing edge can be thought of as acting with a lateral force on the wake through lateral water momentum shed as the fish moves on. The fish's mean rate of working is the mean product of this lateral force with the lateral component of trailing-edge movement, and is enhanced by the virtual-mass effect, which makes for good correlation between lateral movement and local water momentum. The mean rate of shedding of energy of lateral water motions into the vortex wake represents the wasted element in this mean rate of working, and it is from the difference of these two rates that thrust and efficiency can best be calculated.

Section 3, still from the standpoint of inviscid theory, studies the effect of any development of discrete dorsal and ventral fins, through calculations on vortex sheets shed by fins. A multiplicity of discrete dorsal (or ventral) fins might be thought to destroy the slow variation of cross-sectional properties on which elongated-body theory depends, but the vortex sheets filling the gaps between them are shown to maintain continuity rather effectively, avoiding thrust reduction and permitting a slight decrease in drag.

Further advantage may accrue from a modification of such a system in which (while essentially anguilliform movement is retained) the anterior dorsal and ventral fins become the only prominent ones. Vortex sheets in the gaps between them and the caudal fin may largely be reabsorbed into the caudal-fin boundary layer, without any significant increase in wasted wake energy. The mean rate of working can be improved, however, because the trailing edges of the dorsal and ventral fins do work that is not cancelled at the caudal fin's leading edge, as phase shifts destroy the correlation of that edge's lateral movement with the vortex-sheet momentum reabsorbed there.

Tentative improvements to elongated-body theory through taking into account lateral forces of viscous origin are made in §4. These add to both the momentumandenergyof the water's lateral motions, but mayreduce the efficiencyof anguilliform motion because the extra momentum at the trailing edge, resulting from forces exerted by anterior sections, is badly correlated with that edge's lateral movements. Adoption of the ‘carangiform’ mode, in which the amplitude of the basic undulation grows steeply from almost zero over the first half or even two-thirds of a fish's length to a large value at the caudal fin, avoids this difficulty.

Any movement which a fish attempts to make, however, is liable to be accompanied by ‘recoil’, that is, by extra movements of pure translation and rotation required for overall conservation of momentum and angular momentum. These recoil movements, a potentially serious source of thrust and efficiency loss in carangiform motion, are calculated in § 4, which shows how they are minimized with the right distribution of total inertia (the sum of fish mass and the water's virtual mass). It seems to be no coincidence that carangiform motion goes always with a long anterior region of high depth (possessing a substantial moment of total inertia) and a region of greatly reduced depth just before the caudal fin.

The theory suggests (§5) that reduction of caudal-fin area in relation to depth by development of a caudal fin into a herring-like ‘pair of highly sweptback wings’ should reduce drag without significant loss of thrust. The same effect can be expected (although elongated-body theory ceases to be applicable) from widening of the wing pair (sweepback reduction). That line of development of the carangiform mode in many of the Percomorphi leads towards the lunate tail, a culminating point in the enhancement of speed and propulsive efficiency which has been reached also along some quite different lines of evolution.

A beginning in the analysis of its advantages is made here using a ‘twodimensional’ linearized theory. Movements of any horizontal section of caudal fin, with yaw angle fluctuating in phase with its velocity of lateral translation, are studied for different positions of the yawing axis. The wasted energy in the wake has a sharp minimum when that axis is at the ‘three-quarter-chord point’, but rate of working increases somewhat for axis positions distal to that. Something like an optimum regarding efficiency, thrust and the proportion of thrust derived from suction at the section's rounded leading edge is found when the yawing axis is along the trailing edge.

This leads on the present over-simplified theory to the suggestion that a hydromechanically advantageous configuration has the leading edge bowed forward but the trailing edge straight. Finally, there is a brief discussion of possible future work, taking three-dimensional and non-linear effects into account, that might throw light on the commonness of a trailing edge that is itself slightly bowed forward among the fastest marine animals.

The structure of round jets in cross-flow was studied using flow visualization techniques and flying-hot-wire measurements. The study was restricted to jet to freestream velocity ratios ranging from 2.0 to 6.0 and Reynolds numbers based on the jet diameter and free-stream velocity in the range of 440 to 6200.

Flow visualization studies, together with time-averaged flying-hot-wire measurements in both vertical and horizontal sectional planes, have allowed the mean topological features of the jet in cross-flow to be identified using critical point theory. These features include the horseshoe (or necklace) vortex system originating just upstream of the jet, a separation region inside the pipe upstream of the pipe exit, the roll-up of the jet shear layer which initiates the counter-rotating vortex pair and the separation of the flat-wall boundary layer leading to the formation of the wake vortex system beneath the downstream side of the jet.

The topology of the vortex ring roll-up of the jet shear layer was studied in detail using phase-averaged flying-hot-wire measurements of the velocity field when the roll-up was forced. From these data it is possible to examine the evolution of the shear layer topology. These results are supported by the flow visualization studies which also aid in their interpretation.

The study also shows that, for velocity ratios ranging from 4.0 to 6.0, the unsteady upright vortices in the wake may form by different mechanisms, depending on the Reynolds number. It is found that at high Reynolds numbers, the upright vortex orientation in the wake may change intermittently from one configuration of vortex street to another. Three mechanisms are proposed to explain these observations.

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 two types of thermal-convective instabilities: (i) stationary longitudinal rolls, which involve the classical Marangoni instability studied by Pearson; and (ii) unsteady hydrothermal waves, which involve a new mechanism of instability deriving its energy from the horizontal temperature gradients. Thermal stability characteristics for liquid layers with and without return-flow profiles are presented as functions of the Prandtl number of the liquid and the Biot number of the interface. Comparisons are made with available experimental observations.

The turbulent energy equation is converted into a differential equation for the turbulent shear stress by defining three empirical functions relating the turbulent intensity, diffusion and dissipation to the shear stress profile. This equation, the mean momentum equation and the mean continuity equation form a hyperbolic system. Numerical integrations by the method of characteristics with preliminary choices of the three empirical functions compare favourably with the results of conventional calculation methods over a wide range of pressure gradients. Nearly all the empirical information required has been derived solely from the boundary layer in zero pressure gradient.

Measurements of the mean velocity profile and pressure drop were performed in a fully developed, smooth pipe flow for Reynolds numbers from 31×103 to 35×106. Analysis of the mean velocity profiles indicates two overlap regions: a power law for 60<y+<500 or y+<0.15R+, the outer limit depending on whether the Kármán number R+ is greater or less than 9×103; and a log law for 600<y+<0.07R+. The log law is only evident if the Reynolds number is greater than approximately 400×103 (R+>9×103). Von Kármán's constant was shown to be 0.436 which is consistent with the friction factor data and the mean velocity profiles for 600<y+<0.07R+, and the additive constant was shown to be 6.15 when the log law is expressed in inner scaling variables.

A new theory is developed to explain the scaling in both overlap regions. This theory requires a velocity scale for the outer region such that the ratio of the outer velocity scale to the inner velocity scale (the friction velocity) is a function of Reynolds number at low Reynolds numbers, and approaches a constant value at high Reynolds numbers. A reasonable candidate for the outer velocity scale is the velocity deficit in the pipe, UCL−Ū, which is a true outer velocity scale, in contrast to the friction velocity which is a velocity scale associated with the near-wall region which is ‘impressed’ on the outer region. The proposed velocity scale was used to normalize the velocity profiles in the outer region and was found to give significantly better agreement between different Reynolds numbers than the friction velocity.

The friction factor data at high Reynolds numbers were found to be significantly larger (>5%) than those predicted by Prandtl's relation. A new friction factor relation is proposed which is within ±1.2% of the data for Reynolds numbers between 10×103 and 35×106, and includes a term to account for the near-wall velocity profile.

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 963 grid using Fourier collocation methods. Extensive comparisons between the direct and modelled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence – conducted on a coarse 323 grid – is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

In recent work it has been shown that there can be substantial transient growth in the energy of small perturbations to plane Poiseuille and Couette flows if the Reynolds number is below the critical value predicted by linear stability analysis. This growth, which may be as large as O(1000), occurs in the absence of nonlinear effects and can be explained by the non-normality of the governing linear operator - that is, the non-orthogonality of the associated eigenfunctions. In this paper we study various aspects of this energy growth for two- and three-dimensional Poiseuille and Couette flows using energy methods, linear stability analysis, and a direct numerical procedure for computing the transient growth. We examine conditions for no energy growth, the dependence of the growth on the streamwise and spanwise wavenumbers, the time dependence of the growth, and the effects of degenerate eigenvalues. We show that the maximum transient growth behaves like O(R2), where R is the Reynolds number. We derive conditions for no energy growth by applying the Hille–Yosida theorem to the governing linear operator and show that these conditions yield the same results as those derived by energy methods, which can be applied to perturbations of arbitrary amplitude. These results emphasize the fact that subcritical transition can occur for Poiseuille and Couette flows because the governing linear operator is non-normal.

An exact and very general Lagrangian-mean description of the back effect of oscillatory disturbances upon the mean state is given. The basic formalism applies to any problem whose governing equations are given in the usual Eulerian form, and irrespective of whether spatial, temporal, ensemble, or ‘two-timing’ averages are appropriate. The generalized Lagrangian-mean velocity cannot be defined exactly as the ‘mean following a single fluid particle’, but in cases where spatial averages are taken can easily be visualized, for instance, as the motion of the centre of mass of a tube of fluid particles which lay along the direction of averaging in a hypothetical initial state of no disturbance.

The equations for the Lagrangian-mean flow are more useful than their Eulerian-mean counterparts in significant respects, for instance in explicitly representing the effect upon mean-flow evolution of wave dissipation or forcing. Applications to irrotational acoustic or water waves, and to astrogeophysical problems of waves on axisymmetric mean flows are discussed. In the latter context the equations embody generalizations of the Eliassen-Palm and Charney-Drazin theorems showing the effects on the mean flow of departures from steady, conservative waves, for arbitrary, finite-amplitude disturbances to a stratified, rotating fluid, with allowance for self-gravitation as well as for an external gravitational field.

The equations show generally how the pseudomomentum (or wave ‘momentum’) enters problems of mean-flow evolution. They also indicate the extent to which the net effect of the waves on the mean flow can be described by a ‘radiation stress’, and provide a general framework for explaining the asymmetry of radiation-stress tensors along the lines proposed by Jones (1973).

Direct numerical simulations (DNS) and experiments are carried out to study fully developed turbulent pipe flow at Reynolds number Rec ≈ 7000 based on centreline velocity and pipe diameter. The agreement between numerical and experimental results is excellent for the lower-order statistics (mean flow and turbulence intensities) and reasonably good for the higher-order statistics (skewness and flatness factors). To investigate the differences between fully developed turbulent flow in an axisymmetric pipe and a plane channel geometry, the present DNS results are compared to those obtained from a channel flow simulation. Beside the mean flow properties and turbulence statistics up to fourth order, the energy budgets of the Reynolds-stress components are computed and compared. The present results show that the mean velocity profile in the pipe fails to conform to the accepted law of the wall, in contrast to the channel flow. This confirms earlier observations reported in the literature. The statistics on fluctuating velocities, including the energy budgets of the Reynolds stresses, appear to be less affected by the axisymmetric pipe geometry. Only the skewness factor of the normal-to-the-wall velocity fluctuations differs in the pipe flow compared to the channel flow. The energy budgets illustrate that the normal-to-the-wall velocity fluctuations in the pipe are altered owing to a different ‘impingement’ or ‘splatting’ mechanism close to the curved wall.

In a previous paper, Germano, et al. (1991) proposed a method for computing coefficients of subgrid-scale eddy viscosity models as a function of space and time. This procedure has the distinct advantage of being self-calibrating and requires no a priori specification of model coefficients or the use of wall damping functions. However, the original formulation contained some mathematical inconsistencies that limited the utility of the model. In particular, the applicability of the model was restricted to flows that are statistically homogeneous in at least one direction. These inconsistencies and limitations are discussed and a new formulation that rectifies them is proposed. The new formulation leads to an integral equation whose solution yields the model coefficient as a function of position and time. The method can be applied to general inhomogeneous flows and does not suffer from the mathematical inconsistencies inherent in the previous formulation. The model has been tested in isotropic turbulence and in the flow over a backward-facing step.

We describe the results of a numerical investigation of the dynamics of breakup of streams of immiscible fluids in the confined geometry of a microfluidic T-junction. We identify three distinct regimes of formation of droplets: squeezing, dripping and jetting, providing a unifying picture of emulsification processes typical for microfluidic systems. The squeezing mechanism of breakup is particular to microfluidic systems, since the physical confinement of the fluids has pronounced effects on the interfacial dynamics. In this regime, the breakup process is driven chiefly by the buildup of pressure upstream of an emerging droplet and both the dynamics of breakup and the scaling of the sizes of droplets are influenced only very weakly by the value of the capillary number. The dripping regime, while apparently homologous to the unbounded case, is also significantly influenced by the constrained geometry; these effects modify the scaling law for the size of the droplets derived from the balance of interfacial and viscous stresses. Finally, the jetting regime sets in only at very high flow rates, or with low interfacial tension, i.e. higher values of the capillary number, similar to the unbounded case.

In a large range of Reynolds numbers, 6 × 104 < Re < 5 × 106, the flow around single cylinders with smooth surfaces has been investigated. The high values of the Reynolds numbers were obtained in a test channel which could be pressurized up to 40 bar of static pressure. New experiments were performed to measure the local pressure and skin friction distribution around the cylinder. From these results the total drag, the pressure drag and the friction drag were calculated. By means of the skin friction distribution the position of the separation points, separation bubbles or transition points can be localized. These data allow one to define three states of the flow: the subcritical flow, where the boundary layer separates laminarly; the critical flow, in which a separation bubble, followed by a turbulent reattachment, occurs; and the supercritical flow, where an immediate transition from the laminar to the turbulent boundary layer is observed at a critical distance from the stagnation point. According to the total drag coefficient the values found in this paper connect the subcritical region represented by the measurements of Wieselsberger (1923) and Fage & Warsap (1930) with the supercritical range in which Roshko (1961) carried out his experiments.

A solid long slender body is considered placed in a fluid undergoing a given undisturbed flow. Under conditions in which fluid inertia is negligible, the force per unit length on the body is obtained as an asymptotic expansion in terms of the ratio of the cross-sectional radius to body length. Specific examples are given for the resistance to translation of long slender bodies for cases in which the body centre-line is curved as well as for those for which the centre-line is straight.