The Laplace–Young equation is solved for axisymmetric menisci, analytically in terms of elliptic integrals for all possible types of pendular rings and liquid bridges when the effect of gravity is negligible, numerically for selected other cases in order to assess gravity's effect. Meniscus shapes, mean curvatures, areas and enclosed volumes are reported, as are capillary forces. It is shown that capillary attraction may become capillary repulsion when wetting is imperfect. The special configurations of vanishing capillary force and of zero mean curvature are treated. The range of utility of the convenient ‘circle approximation’ is evaluated.

Experimental studies of the flow field near the wall in a turbulent boundary layer using hot-wire probes are reported. Measurements of the product uv are studied using the technique of conditional sampling with a large digital computer to single out special events (bursting) when large contributions to turbulent energy and Reynolds stress occur. The criterion used to determine when the product uv is sampled is that the streamwise velocity at the edge of the sublayer should have attained a certain value. With this simple criterion we find that 60% of the contribution to $\overline{uv}$ is produced when the sublayer velocity is lower than the mean. This result is true at both low, Rθ = 4230, and high, Rθ = 38 000, Reynolds numbers. With a more strict sampling criterion, that the filtered sublayer velocity at two side-by-side points should be simultaneously low and decreasing, individual contributions to $\overline{uv}$ as large as 62 $\overline{uv}$ have been identified. Additional measurements using correlations between truncated u and v signals reveal that the largest contributions to the Reynolds stress and turbulent energy occur when u < 0, v > 0 during an intense bursting process and the remainder of the contributions occur during a less intense recovery process. Thus, contributions to the turbulent energy production and Reynolds stress at a point near the wall are of relatively large magnitude, short duration and occur intermittently. A rough measure of the intermittency factor for uv at a point near the wall is 0·55 since 99% of the contribution to $\overline{uv}$ is made during only 55% of the total time.

Transition in the boundary layer on a flat plate is examined from the point of view of intermittent production of turbulent spots. On the hypothesis of localized laminar breakdown, for which there is some expermental evidence, Emmons’ probability calculations can be extended to explain the observed statistical similarity of transition regions. Application of these ideas allows detailed calculations of the boundary layer parameters including mean velocity profiles and skin friction during transition. The mean velocity profiles belong to a universal one-parameter family with the intermittency factor as the parameter. From an examination of experimental data the probable existence of a relation between the transition Reynolds number and the rate of production of the turbulent spots is deduced. A simple new technique for the measurement of the intermittency factor by a Pitot tube is reported.

The Rayleigh stability equation of inviscid linearized stability theory was integrated numerically for amplified disturbances of the hyperbolic-tangent velocity profile. The evaluation of the eigenvalues and eigenfunctions is followed by a discussion of the streamline pattern of the disturbed flow. Here no qualitative distinction is found between an amplified and the neutral disturbance. But considering the vorticity distribution of the disturbed flow it is shown that in the case of amplified disturbances two concentrations of vorticity occur within a disturbance wavelength, while in the neutral case only one maximum of vorticity exists. The results are discussed with respect to the instability mechanism of free boundary-layer flow.

In this paper the dimensional-analysis approach to wall turbulence of Perry & Abell (1977) has been extended in a number of directions. Further recent developments of the attached-eddy hypothesis of Townsend (1976) and the model of Perry & Chong (1982) are given, for example, the incorporation of a Kolmogoroff (1941) spectral region. These previous analyses were applicable only to the ‘wall region’ and are extended here to include the whole turbulent region of the flow. The dimensional-analysis approach and the detailed physical modelling are consistent with each other and with new experimental data presented here.

We use machine learning to perform super-resolution analysis of grossly under-resolved turbulent flow field data to reconstruct the high-resolution flow field. Two machine learning models are developed, namely, the convolutional neural network (CNN) and the hybrid downsampled skip-connection/multi-scale (DSC/MS) models. These machine learning models are applied to a two-dimensional cylinder wake as a preliminary test and show remarkable ability to reconstruct laminar flow from low-resolution flow field data. We further assess the performance of these models for two-dimensional homogeneous turbulence. The CNN and DSC/MS models are found to reconstruct turbulent flows from extremely coarse flow field images with remarkable accuracy. For the turbulent flow problem, the machine-leaning-based super-resolution analysis can greatly enhance the spatial resolution with as little as 50 training snapshot data, holding great potential to reveal subgrid-scale physics of complex turbulent flows. With the growing availability of flow field data from high-fidelity simulations and experiments, the present approach motivates the development of effective super-resolution models for a variety of fluid flows.

The phenomenon examined is the abrupt structural change which can occur at some station along the axis of a swirling flow, notably the leading-edge vortex formed above a delta wing at incidence. Contrary to previous attempts at an explanation, the belief demonstrated herein is that vortex breakdown is not a manifestation of instability or of any other effect indicated by study of infinitesimal disturbances alone. It is instead a finite transition between two dynamically conjugate states of axisymmetric flow, analogous to the hydraulic jump in openchannel flow. A set of properties essential to such a transition, corresponding to a set shown to provide a complete explanation for the hydraulic jump, is demonstrated with wide generality for axisymmetric vortex flows; and the interpretation covers both the case of mild transitions, where an undular structure is developed without the need arising for significant energy dissipation, and the case of strong ones where a region of vigorous turbulence is generated. An important part of the theory depends on the calculus of variations; and the comprehensiveness with which certain properties of conjugate flow pairs are demonstrable by this analytical means suggests that present ideas may be useful in various other problems.

When a layer of fluid is heated uniformly from below and cooled from above, a cellular regime of steady convection is set up at values of the Rayleigh number exceeding a critical value. A method is presented here to determine the form and amplitude of this convection. The non-linear equations describing the fields of motion and temperature are expanded in a sequence of inhomogeneous linear equations dependent upon the solutions of the linear stability problem. We find that there are an infinite number of steady-state finite amplitude solutions (having different horizontal plan-forms) which formally satisfy these equations. A criterion for ‘relative stability’ is deduced which selects as the realized solution that one which has the maximum mean-square temperature gradient. Particular conclusions are that for a large Prandtl number the amplitude of the convection is determined primarily by the distortion of the distribution of mean temperature and only secondarily by the self-distortion of the disturbance, and that when the Prandtl number is less than unity self-distortion plays the dominant role in amplitude determination. The initial heat transport due to convection depends linearly on the Rayleigh number; the heat transport at higher Rayleigh numbers departs only slightly from this linear dependence. Square horizontal plan-forms are preferred to hexagonal plan-forms in ordinary fluids with symmetric boundary conditions. The proposed finite amplitude method is applicable to any model of shear flow or convection with a soluble stability problem.

Force measurements were conducted in a pressurized wind tunnel from subcritical up to transcritical Reynolds numbers 2.3 × 104 [les ] Re [les ] 7.1 × 106 without changing the experimental arrangement. The steady and unsteady forces were measured by means of a piezobalance, which features a high natural frequency, low interferences and a large dynamic range. In the critical Reynolds-number range, two discontinuous transitions were observed, which can be interpreted as bifurcations at two critical Reynolds numbers. In both cases, these transitions are accompanied by critical fluctuations, symmetry breaking (the occurrence of a steady lift) and hysteresis. In addition, both transitions were coupled with a drop of the CD value and a jump of the Strouhal number. Similar phenomena were observed in the upper transitional region between the super- and the transcritical Reynolds-number ranges. The transcritical range begins at about Re ≈ 5 × 106, where a narrow-band spectrum is formed with Sr(Re = 7.1 × 106) = 0.29.

The frequency and amplification rates for a disturbance growing with respect to time are compared with those of a spatially-growing wave having the same wave-number. For small rates of amplification it is shown that the frequencies are equal to a high order of approximation, and that the spatial growth is related to the time growth by the group velocity.

The statistical density function is derived for a variable (such as the surface elevation in a random sea) that is ‘weakly non-linear’. In the first approximation the distribution is Gaussian, as is well known. In higher approximations it is shown that the distribution is given by successive sums of a Gram–Charlier series; not quite in the form that has sometimes been used as an empirical fit for observed distributions, but in a modified form due to Edgeworth.

It is shown that the cumulants of the distribution are much simpler to calculate than the corresponding moments; and the approximate distributions are in fact derived by inversion of the cumulant-generating function.

The theory is applied to random surface waves on water. The third cumulant and hence the skewness of the distribution of surface elevation is evaluated explicity in terms of the directional energy spectrum. It is shown that the skewness λ3 is generally positive, and positive upper and lower bounds for λ3 are derived. The theoretical results are compared with some measurements made by Kinsman (1960).

It is found that for free, undamped surface waves the skewness of the distribution of surface slopes is of a higher order than the skewness of the surface elevation. Hence the observed skewness of the slopes may be a sensitive indicator of energy transfer and dissipation within the water.

A powerful, easy-to-use analytic technique for nonlinear problems, the homotopy analysis method, is employed to give analytic solutions of magnetohydrodynamic viscous flows of non-Newtonian fluids over a stretching sheet. For the so-called second-order and third-order power-law fluids, the explicit analytic solutions are given by recursive formulas with constant coefficients. Also, for real power-law index and magnetic field parameter in a quite large range, an analytic approach is proposed. All of our analytic results agree well with numerical ones. In particular, a simple analytic formula of the dimensionless velocity gradient at the wall is found, which is accurate for all real power-law indices and magnetic field parameters. This analytic formula can give sufficiently accurate results for the skin friction on the moving sheet that it would find wide application in industries. Physically, they indicate that the magnetic field tends to increase the skin friction, and that this effect is more pronounced for shear-thinning than for shear-thickening fluids.

Approximate layer-averaged equations describing the mechanics of turbid underflows are derived. Closure of the equations describing the balance of fluid mass, sediment mass, and mean flow momentum provides for the delineation of a three-equation model. A description of sediment exchange with the bed allows for the possibility of a self-accelerating turbidity current in which sediment entrainment from the bed is linked to flow velocity. A consideration of the balance of the mean energy of the turbulence yields a constraint on physically realistic solutions to the three-equation model. It is shown that the self-acceleration predicted by the three-equation model is so strong that the energy constraint fails to be satisfied. In particular, the turbulent energy consumed in entraining new bed sediment exceeds the supply of energy to the turbulence, so that the turbulence, and thus the turbidity current, must die. The problem is rectified by the formulation of a four-equation model, in which an explicit accounting is made of the mean energy of the turbulence. Sediment entrainment from the bed is linked to the level of turbulence in the four-equation model. Self-acceleration is again predicted, although it is somewhat subdued compared with that predicted by the three-equation model. The predictions of both models are summarized over a wide range of conditions.

A numerical method for studying water entry of a two-dimensional body of arbitrary cross-section is presented. It is a nonlinear boundary element method with a jet flow approximation. The method has been verified by comparisons with new similarity solution results for wedges with deadrise angles varying from 4° to 81°. A simple asymptotic solution for small deadrise angles α based on Wagner (1932) agrees with the similarity solution for small α.

The effect of periodic two-dimensional excitation on the development of a turbulent mixing region was studied experimentally. Controlled oscillations of variable ampli- tude and frequency were applied at the initiation of mixing between two parallel air streams. The frequency of forcing was at least an order of magnitude lower than the initial instability frequency of the flow in order to test its effect far downstream. The effect of the velocity difference between the streams was also investigated in this experiment. A typical Reynolds number based on the velocity difference and the momentum thickness of the shear layer was l04.

It was determined that the spreading rate of the mixing layer is sensitive to periodic surging even if the latter is so small that it does not contribute to the initial energy of the fluctuations. Oscillations at very small amplitudes tend to increase the spreading rate of the flow by enhancing the amalgamation of neighbouring eddies, but at higher amplitudes the flow resonates with the imposed oscillation. The resonance region can extend over a significant fraction of the test section depending on the Strouhal number and a dimensionless velocity-difference parameter. The flow in the resonance region consists of a single array of large, quasi-two-dimensional vortex lumps, which do not interact with one another. The exponential shape of the mean-velocity distribution is not affected in this region, but the spreading rate of the flow with increasing distance downstream is inhibited. The Reynolds stress in this region changes sign, indicating that energy is extracted from the turbulence to the mean motion; the intensity of the spanwise fluctuations is also reduced, suggesting that the flow tends to become more two-dimensional.

Amalgamation of large coherent eddies is resumed beyond the resonance region, but the flow is not universally similar. There are many indications suggesting that the large eddies in the turbulent mixing layer at fairly large Re are governed by an inviscid instability.

The flow of water in straight open channels with prismatic complex cross-sections is considered. Lateral distributions of depth-mean velocity and boundary shear stress are derived theoretically for channels of any shape, provided that the boundary geometry can be discretized into linear elements. The analytical model includes the effects of bed-generated turbulence, lateral shear turbulence and secondary flows. Experimental data from the Science and Engineering Research Council (SERC) Flood Channel Facility are used to illustrate the relative importance of these three effects on internal shear stresses. New experimental evidence concerning the spatial distribution of Reynolds stresses τyx and τzx is presented for the particular case of compound or two-stage channels. In such channels the vertical distributions of τzx are shown to be highly nonlinear in the regions of strongest lateral shear and the depth-averaged values of τyx are shown to be significantly different from the depth mean apparent shear stresses. The importance of secondary flows in the lateral shear layer region is therefore established. The influence of both Reynolds stresses and secondary flows on eddy viscosity values is quantified. A numerical study is undertaken of the lateral distributions of local friction factor and dimensionless eddy viscosity. The results of this study are then used in the analytical model to reproduce lateral distributions of depth-mean velocity and boundary shear stress in a two stage channel. The work will be of interest to engineers engaged in flood channel hydraulics and overbank flow in particular.

A characteristic feature of a steady trailing line vortex from one side of a wing, and of other types of line vortex, is the existence of strong axial currents near the axis of symmetry. The purpose of this paper is to account in general terms for this axial flow in trailing line vortices. the link between the azimuthal and axial components of motion in a steady line vortex is provided by the pressure; the radial pressure gradient balances the centrifugal force, and any change in the azimuthal motion with distance x downstream produces an axial pressure gradient and consequently axial acceleration.

It is suggested, in a discussion of the evolution of an axisymmetric line vortex out of the vortex sheet shed from one side of a wing, that the two processes of rolling-up of the sheet and of concentration of the vorticity into a smaller cross-section should be distinguished; the former always occurs, whereas the latter seems not to be inevitable.

In § 4 there is given a similarity solution for the flow in a trailing vortex far downstream where the departure of the axial velocity from the free stream speed is small. The continual slowing-down of the azimuthal motion by viscosity leads to a positive axial pressure gradient and consequently to continual loss of axial momentum, the asymptotic variation of the axial velocity defect at the centre being as x−1 log x.

The concept of the drag associated with the core of a trailing vortex is introduced, and the drag is expressed as an integral over a transverse plane which is independent of x. This drag is related to the arbitrary constant appearing in the above similarity solution.

Measurements of air and solid-particle velocities were made in a vertical pipe two-phase flow by the use of a laser-Doppler velocimeter (LDV). Five kinds of plastic particles, diameters of which ranged from about 3 mm to 200 μm, were transported in a vertical pipe of 30 mm inner diameter. It was found that, the smaller the particle size, the flatter was the mean air velocity distribution for the same mass flow ratio of solids to air. Large particles increased air turbulence throughout the pipe section, while small particles reduced it. Both effects of promotion and suppression of turbulence were observed at the same time in the presence of particles of medium size, that is, the turbulence was increased around the pipe centre and reduced near the wall. The frequency spectrum of air turbulence normalized by the turbulence intensity was not changed by the large particles. In the presence of the small particles, the higher-frequency parts of the spectrum increased.

Coherent structures (CS) near the wall (i.e. y + ≤ 60) in a numerically simulated turbulent channel flow are educed using a conditional sampling scheme which extracts the entire extent of dominant vortical structures. Such structures are detected from the instantaneous flow field using our newly developed vortex definition (Jeong & Hussain 1995) - a region of negative λ2, the second largest eigenvalue of the tensor SikSkj + ΩikΩkj - which accurately captures the structure details (unlike velocity-, vorticity- or pressure-based eduction). Extensive testing has shown that λ2 correctly captures vortical structures, even in the presence of the strong shear occurring near the wall of a boundary layer. We have shown that the dominant near-wall educed (i.e. ensemble averaged after proper alignment) CS are highly elongated quasi-streamwise vortices; the CS are inclined 9° in the vertical (x, y)-plane and tilted ±4° in the horizontal (x, z)-plane. The vortices of alternating sign overlap in x as a staggered array; there is no indication near the wall of hairpin vortices, not only in the educed data but also in instantaneous fields. Our model of the CS array reproduces nearly all experimentally observed events reported in the literature, such as VITA, Reynolds stress distribution, wall pressure variation, elongated low-speed streaks, spanwise shear, etc. In particular, a phase difference (in space) between streamwise and normal velocity fluctuations created by CS advection causes Q4 ('sweep’) events to dominate Q2 ('ejection’) and also creates counter-gradient Reynolds stresses (such as Ql and Q3 events) above and below the CS. We also show that these effects are adequately modelled by half of a Batchelor's dipole embedded in (and decoupled from) a background shear U(y). The CS tilting (in the (x, z)-plane) is found to be responsible for sustaining CS through redistribution of streamwise turbulent kinetic energy to normal and spanwise components via coherent pressure-strain effects.

This work is dedicated to systematically studying and predicting the wake characteristics of a yawed wind turbine immersed in a turbulent boundary layer. To achieve this goal, wind tunnel experiments were performed to characterize the wake of a horizontal-axis wind turbine model. A high-resolution stereoscopic particle image velocimetry system was used to measure the three velocity components in the turbine wake under different yaw angles and tip-speed ratios. Moreover, power and thrust measurements were carried out to analyse the performance of the wind turbine. These detailed wind tunnel measurements were then used to perform a budget study of the continuity and Reynolds-averaged Navier–Stokes equations for the wake of a yawed turbine. This theoretical analysis revealed some notable features of the wakes of yawed turbines, such as the asymmetric distribution of the wake skew angle with respect to the wake centre. Under highly yawed conditions, the formation of a counter-rotating vortex pair in the wake cross-section as well as the vertical displacement of the wake centre were shown and analysed. Finally, this study enabled us to develop general governing equations upon which a simple and computationally inexpensive analytical model was built. The proposed model aims at predicting the wake deflection and the far-wake velocity distribution for yawed turbines. Comparisons of model predictions with the wind tunnel measurements show that this simple model can acceptably predict the velocity distribution in the far wake of a yawed turbine. Apart from the ability of the model to predict wake flows in yawed conditions, it can provide valuable physical insight on the behaviour of turbine wakes in this complex situation.