The quadrant analysis of the intense tangential Reynolds stress in plane turbulent channels is generalized to three-dimensional structures (Qs), with special emphasis on the logarithmic and outer layers. Wall-detached Qs are background stress fluctuations. They are small and isotropically oriented, and their contributions to the mean stress cancel. Wall-attached Qs are larger, and carry most of the mean Reynolds stresses. They form a family of roughly self-similar objects that become increasingly complex away from the wall, resembling the vortex clusters in del Álamo et al. (J. Fluid Mech., vol. 561, 2006, pp. 329–358). Individual Qs have fractal dimensions of the order of , slightly fuller than the clusters. They can be described as ‘sponges of flakes’, while vortex clusters are ‘sponges of strings’. The number of attached Qs decays away from the wall, but the fraction of the stress that they carry is independent of their sizes. A substantial fraction of the stress resides in a few large objects extending beyond the centreline, reminiscent of the very large structures of several authors. The predominant logarithmic-layer structure is a side-by-side pair of a sweep (Q4) and an ejection (Q2), with an associated cluster, and shares dimensions and stresses with the conjectured attached eddies of Townsend (J. Fluid Mech., vol. 11, 1961, pp. 97–120). Those attached eddies tend to be aligned streamwise from each other, located near the side walls between the low- and high-velocity large-scale streaks, but that organization does not extend far enough to explain the very long structures in the centre of the channel.

The statistics of a turbulent passive scalar (temperature) and their Reynolds number dependence are studied in decaying grid turbulence for the Taylor-microscale Reynolds number, Rλ, varying from 30 to 731 (21[les ]Peλ[les ]512). A principal objective is, using a single (and simple) flow, to bridge the gap between the existing passive grid-generated low-Péclet-number laboratory experiments and those done at high Péclet number in the atmosphere and oceans. The turbulence is generated by means of an active grid and the passive temperature fluctuations are generated by a mean transverse temperature gradient, formed at the entrance to the wind tunnel plenum chamber by an array of differentially heated elements. A well-defined inertial–convective scaling range for the scalar with a slope, nθ, close to the Obukhov–Corrsin value of 5/3, is observed for all Reynolds numbers. This is in sharp contrast with the velocity field, in which a 5/3 slope is only approached at high Rλ. The Obukhov–Corrsin constant, Cθ, is estimated to be 0.45–0.55. Unlike the velocity spectrum, a bump occurs in the spectrum of the scalar at the dissipation scales, with increasing prominence as the Reynolds number is increased. A scaling range for the heat flux cospectrum was also observed, but with a slope around 2, less than the 7/3 expected from scaling theory. Transverse structure functions of temperature exist at the third and fifth orders, and, as for even-order structure functions, the width of their inertial subranges dilates with Reynolds number in a systematic way. As previously shown for shear flows, the existence of these odd-order structure functions is a violation of local isotropy for the scalar differences, as is the existence of non-zero values of the transverse temperature derivative skewness (of order unity) and hyperskewness (of order 100). The ratio of the temperature derivative standard deviation along and normal to the gradient is 1.2±0.1, and is independent of Reynolds number. The refined similarity hypothesis for the passive scalar was found to hold for all Rλ, which was not the case for the velocity field. The intermittency exponent for the scalar, μθ, was found to be 0.25±0.05 with a possible weak Rλ dependence, unlike the velocity field, where μ was a strong function of Reynolds number. New, higher-Reynolds-number results for the velocity field, which smoothly follow the trends of Mydlarski & Warhaft (1996), are also presented.

The steady distribution function for homogeneous turbulence is studied starting from Liouville's equation, modified by the introduction of an instantaneously fluctuating external force, which acts as a random source of energy. A new technique for solving Liouville's equation is presented giving a systematic development of the concepts of turbulent diffusion and turbulent viscosity. It amounts to a consistent generalization of the random phase approximation. When the rate of input of energy into the kth Fourier component uk has a power form h|k|−α, the functional form of the mean value 〈 uku−k 〉 can be determined exactly in the limit of large Reynolds number; it is $Ah^{\frac{2}{3}}|\bf K|^ {-{\frac {1}{3}}(5+2 \alpha)}$. Liouville's equation proves an inadequate basis for the steady time-dependent mean $\langle u_k(t)u_{-k}(t^ \prime) \rangle $ and a more general equation is derived. The new equation can be solved in a similar way and shows that the time-dependent correlation starts like a Gaussian in time, then passes through an exponentially decaying state, then eventually has a power dependence $|t-t^\prime|- \gamma^ {|k|}$.

In this investigation the creeping flow due to the motion of a liquid drop or a bubble in another immiscible fluid is examined when the interface is partially covered by a stagnant layer of surfactant. The associated boundary-value problem involves mixed boundary conditions at the interface, which lead to a set of dual series equations. An inversion of these equations yields the exact solution to the stagnant cap problem.

Several useful results are obtained in closed form. Among these are the expressions for the drag force, the difference between the maximum and the minimum interfacial tensions, and the amount of adsorbed surfactant. A shifting of the centre of the internal vortex is observed.

A numerical study is made of the temporal eigenvalue spectrum of the Orr-Sommerfeld equation for the Blasius boundary layer. Unlike channel flows, there is no mathematical proof that this flow has an infinite spectrum of discrete eigenvalues. The Orr-Sommerfeld equation is integrated numerically, and the eigenvalues located by tracing out the contour lines in the complex wave velocity (c = cr + ici) plane on which the real and imaginary parts of the secular determinant are zero. This method gives only a finite and small number of discrete eigenvalues for a wide range of Reynolds numbers and wavenumbers. The spectrum of plane Poiseuille flow is used as a guide to study the spectrum of an artificial two wall flow which consists of two Blasius boundary layers. As the upper boundary of this flow moves to infinity, it is found that the portion of the spectrum with an infinite number of eigenvalues moves towards cr = 1 and the spacing between eigenvalues goes to zero. It is concluded, on the basis of this result and the contour method, that the original few eigenvalues found are the only discrete eigenvalues that exist for Blasius flow over a wide portion of the c plane for cr < 1 and cr > 1. It is suggested that the discrete spectrum is supplemented by a continuous spectrum which lies along the cr = 1 axis for ci < −α/R.

Explosive volcanic jets present an unusual dynamic situation of reversing buoyancy. Their initially negative buoyancy with respect to ambient fluid first opposes the motion, but can change sign to drive a convective plume if a sufficient amount of entrainment occurs. The key unknown is the entrainment behaviour for the initial flow regime in which buoyancy acts against the momentum jet. To describe and constrain this regime, we present an experimental study of entrainment into turbulent jets of negative and reversing buoyancy. Using an original technique based on the influence of the injection radius on the threshold between buoyant convection and partial collapse, we show that entrainment is significantly reduced by negative buoyancy. We develop a new theoretical parameterization of entrainment as a function of the local (negative) Richardson number that (i) predicts the observed reduction of entrainment and (ii) introduces a similarity drift in the velocity and buoyancy profiles as a function of distance from source. This similarity drift allows us to reconcile the different estimates found in the literature for entrainment in plumes.

An experimental investigation is conducted into the collapse of granular columns inside rectangular channels. The final shape is documented for slumps inside relatively wide channels, and for collapses inside much narrower slots. In both cases, the collapse is initiated by withdrawing a swinging gate or sliding door, and the flow remains fairly two-dimensional. Four different granular media are used; the properties of the materials vary significantly, notably in their angles of friction for basal sliding and internal deformation. If $H$ is the initial height of the column, $h_{\infty}$ the maximum final height of the column and $a$ the initial aspect ratio, then the data suggest that $H/h_{\infty} \,{\sim} a^{0.6}$ in wide channels and $H/h_{\infty} \,{\sim}\, a^{0.5}$ for narrow slots. For the runout, we find that $(l_{\infty}\,{-}\,L)/L \,{\sim}\, a^{0.9\pm 0.1}$ for wide channels, and $(l_{\infty}\,{-}\,L)/L \,{\sim}\, a^{0.65\pm0.05}$ or $l_\infty/L \,{\sim}\, a^{0.55\pm0.05}$ for narrow slots, where $l_{\infty}$ is the maximum runout of the material and $L$ the initial length of the column along the channel ($a\,{:=}\,H/L$). In all cases, the numerical constant of proportionality in these scaling relations shows clear material dependence. In wide slots, there is no obvious universal scaling behaviour of the final profile, but such a behaviour is evident in narrow slots. The experimental results are compared with theoretical results based on a shallow granular-flow model. The qualitative behaviour of the slump in the wide slot is reproduced by the theoretical model. However, there is qualitative disagreement between theory and the experiments in the narrow slot because of the occurrence of secondary surface avalanching.

We investigate numerically the transition from laminar to chaotic flow of a Boussinesq fluid with Pr = 0.71 in two-dimensional closed, differentially heated, vertical cavities having aspect ratios near unity. The cavities have rigid conducting sidewalls, and rigid insulating top and bottom walls. The physical nature of the resulting flow is a function of the aspect ratio and Rayleigh number.

It is shown that an oscillatory approach to steady-state, oscillatory instabilities, quasi-periodic flow, and chaotic flow exist for the flow regimes investigated. We find that for aspect ratios of approximately three or larger the the first transition from steady-state is due to instability of the sidewall boundary layers, while for small aspect ratios, but larger than ½, it is due to internal waves near the departing corners. For both instabilities we obtain the critical Rayleigh number as a function of aspect ratio and write expressions relating the fundamental frequencies of the oscillatory flow to the Rayleigh number and aspect ratio. When Ra is increased significantly above the first critical value, the flow becomes complex since both types of instabilities can be present. With a further increase in Rayleigh number the flow becomes chaotic and eventually turbulent. The above results are illustrated for different Rayleigh numbers and aspect ratios using time histories, spectral analysis, and streamlines at different values of time.

The time development of the symmetrical standing zone of recirculation, which is formed in the early stage of the flow due to a circular cylinder impulsively set in motion perpendicular to its generators, has been studied using a flow visualization technique. The Reynolds numbers (based upon the diameter) range from 40 to 104. Some new phenomena indicated in the flow patterns are revealed, and several different regimes are differentiated by a detailed analysis of the evolution of the main flow characteristics. A correlation with some theoretical results is established.

It is well known that the drag in a turbulent flow of a polymer solution is significantly reduced compared to Newtonian flow. Here we consider this phenomenon by means of a direct numerical simulation of a turbulent channel flow. The polymers are modelled as elastic dumbbells using the FENE-P model. In the computations the polymer model is solved simultaneously with the flow equations, i.e. the polymers are deformed by the flow and in their turn influence the flow structures by exerting a polymer stress. We have studied the results of varying the polymer parameters, such as the maximum extension, the elasticity and the concentration. For the case of highly extensible polymers the results of our simulations are very close to the maximum drag reduction or Virk (1975) asymptote. Our simulation results show that at approximately maximum drag reduction the slope of the mean velocity profile is increased compared to the standard logarithmic profile in turbulent wall flows. For the r.m.s. of the streamwise velocity fluctuations we find initially an increase in magnitude which near maximum drag reduction changes to a decrease. For the velocity fluctuations in the spanwise and wall-normal directions we find a continuous decrease as a function of drag reduction. The Reynolds shear stress is strongly reduced, especially near the wall, and this is compensated by a polymer stress, which at maximum drag reduction amounts to about 40% of the total stress. These results have been compared with LDV experiments of Ptasinski et al. (2001) and the agreement, both qualitatively and quantitatively, is in most cases very good. In addition we have performed an analysis of the turbulent kinetic energy budgets. The main result is a reduction of energy transfer from the streamwise direction, where the production of turbulent kinetic energy takes place, to the other directions. A substantial part of the energy production by the mean flow is transferred directly into elastic energy of the polymers. The turbulent velocity fluctuations also contribute energy to the polymers. The elastic energy of the polymers is subsequently dissipated by polymer relaxation. We have also computed the various contributions to the pressure fluctuations and identified how these change as a function of drag reduction. Finally, we discuss some cross-correlations and various length scales. These simulation results are explained here by two mechanisms. First, as suggested by Lumley (1969) the polymers damp the cross-stream or wall-normal velocity fluctuations and suppress the bursting in the buffer layer. Secondly, the ‘shear sheltering’ mechanism acts to amplify the streamwise fluctuations in the thickened buffer layer, while reducing and decoupling the motions within and above this layer. The expression for the substantial reduction in the wall drag derived by considering the long time scales of the nonlinear fluctuations of this damped shear layer, is shown to be consistent with the experimental data of Virk et al. (1967) and Virk (1975).

The Koopman operator provides a powerful way of analysing nonlinear flow dynamics using linear techniques. The operator defines how observables evolve in time along a nonlinear flow trajectory. In this paper, we perform a Koopman analysis of the first Hopf bifurcation of the flow past a circular cylinder. First, we decompose the flow into a sequence of Koopman modes, where each mode evolves in time with one single frequency/growth rate and amplitude/phase, corresponding to the complex eigenvalues and eigenfunctions of the Koopman operator, respectively. The analytical construction of these modes shows how the amplitudes and phases of nonlinear global modes oscillating with the vortex shedding frequency or its harmonics evolve as the flow develops and later sustains self-excited oscillations. Second, we compute the dynamic modes using the dynamic mode decomposition (DMD) algorithm, which fits a linear combination of exponential terms to a sequence of snapshots spaced equally in time. It is shown that under certain conditions the DMD algorithm approximates Koopman modes, and hence provides a viable method to decompose the flow into saturated and transient oscillatory modes. Finally, the relevance of the analysis to frequency selection, global modes and shift modes is discussed.

Pressure fluctuations in a turbulent channel flow are investigated by analysing a database obtained from a direct numerical simulation. Detailed statistics associated with the pressure fluctuations are presented. Characteristics associated with the rapid (linear) and slow (nonlinear) pressure are discussed. It is found that the slow pressure fluctuations are larger than the rapid pressure fluctuations throughout the channel except very near the wall, where they are about the same magnitude. This is contrary to the common belief that the nonlinear source terms are negligible compared with the linear source terms. Probability density distributions, power spectra, and two-point correlations are examined to reveal the characteristics of the pressure fluctuations. The global dependence of the pressure fluctuations and pressure–strain correlations are also examined by evaluating the integral associated with Green-function representations of them. In the wall region where the pressure–strain terms are large, most contributions to the pressure–strain terms are from the wall region (i.e. local), whereas away from the wall where the pressure–strain terms are small, contributions are global. Structures of instantaneous pressure and pressure gradients at the wall and the corresponding vorticity field are examined.

Four alternative theoretical treatments of ‘displacement thickness’, and, generally, of the influence of boundary layers and wakes on the flow outside them, are set out, first for two-dimensional, and then for three-dimensional, laminar or turbulent, incompressible flow. They may be called the methods of ‘flow reduction’, ‘equivalent sources’, ‘velocity comparison’ and ‘mean vorticity’.

The principal expression obtained for the displacement thickness δ1 in three-dimensional flow may be written $\delta_1 = \delta_x - \frac{1}{Uh_y} \frac {\partial}{\partial y}\int^x_0 \delta_y dx,$ if, as orthogonal coordinates (x, y) specifying position on the surface, we choose x as the velocity potential of the external flow, and y as a coordinate, constant along the external-flow streamlines, such that hydy is the distance between (x, y) and (x, y + dy); and if also δx and δy are the streamwise and transverse ‘volume-flow thicknesses’ $\delta_x = \frac {1}{U}\int _0^\infty (U - u)\;dz,\;\;\;\;\;\;\delta_y = \frac {1}{U}\int ^\infty_0 v\; dz,$z is the distance from the surface, u and v are the x and y components of velocity, and u takes the value U just outside the boundary layer.

In a turbulent two-dimensional flow enstrophy systematically cascades to very small scales, at which it is dissipated. The kinetic energy, on the other hand, remains at large scales and the total kinetic energy is constant. Above random topography an initially turbulent flow tends to a steady state with streamlines parallel to contours of constant depth, anticyclonic around a bump. A numerical experiment verifies this prediction. In a closed basin on a beta-plane the solution with minimum enstrophy implies a westward flow in the interior, returning in narrow boundary layers to the north and south. This result is interpreted using a parameterization of the effects of the eddies on the large-scale flow. The numerical solution is in qualitative agreement, but corresponds to a minimum of a more complex measure of the total enstrophy than the usual quadratic integral.

Foam in porous media exhibits an unusually high apparent viscosity, making it useful in many industrial processes. The rheology of foam, however, is complex and not well understood. Previous pore-level models of foam are based primarily on studies of bubble flow in circular capillaries. A circular capillary, however, lacks the corners that characterize the geometry of the pores. We study the pressure–velocity relation of bubble flow in polygonal capillaries. A long bubble in a polygonal capillary acts as a leaky piston. The ‘piston’ is reluctant to move because of a large drag exerted by the capillary sidewalls. The liquid in the capillary therefore bypasses the bubble through the leaky corners at a speed an order higher than that of the bubble. Consequently, the pressure work is dissipated predominantly by the motion of the fluid and not by the motion of the bubble. This is opposite to the conclusion based on bubble flow in circular capillaries. The discovery of this new flow regime reconciles two groups of contradictory foam-flow experiments.

Part 1 of this work studies the fluid films deposited on capillary walls in the limit Ca → 0 (Ca ≡ μU/σ, where μ is the fluid viscosity, U the bubble velocity, and σ the surface tension). Part 2 (Wong et al. 1995) uses the film profile at the back end to calculate the drag of the bubble. Since the bubble length is arbitrary, the film profile is determined here as a general function of the dimensionless downstream distance x. For 1 [Lt ] x [Lt ] Ca−1, the film profile is frozen with a thickness of order Ca2/3 at the centre and order Ca at the sides. For x ∼ Ca−1, surface tension rearranges the film at the centre into a parabolic shape while the film at the sides thins to order Ca4/3. For x [Gt ] Ca−1, the film is still parabolic, but the height decreases as film fluid leaks through the side constrictions. For x ∼ Ca−5/3, the height of the parabola is order Ca2/3. Finally, for x [Gt ] Ca−5/3, the height decreases as Ca1/4x−1/4.

Interaction of isotropic quasi-incompressible turbulence with a weak shock wave was studied by direct numerical simulations. The effects of the fluctuation Mach number Mt of the upstream turbulence and the shock strength M21 — 1 on the turbulence statistics were investigated. The ranges investigated were 0.0567 ≤ Mt ≤ 0.110 and 1.05 ≤ M1 ≤ 1.20. A linear analysis of the interaction of isotropic turbulence with a normal shock wave was adopted for comparisons with the simulations.

Both numerical simulations and the linear analysis of the interaction show that turbulence is enhanced during the interaction with a shock wave. Turbulent kinetic energy and transverse vorticity components are amplified, and turbulent lengthscales are decreased. The predictions of the linear analysis compare favourably with simulation results for flows with M2t < a(M21 — 1) with a ≈ 0.1, which suggests that the amplification mechanism is primarily linear. Simulations also showed a rapid evolution of turbulent kinetic energy just downstream of the shock, a behaviour not reproduced by the linear analysis. Investigation of the budget of the turbulent kinetic energy transport equation shows that this behaviour can be attributed to the pressure transport term.

Shock waves were found to be distorted by the upstream turbulence, but still had a well-defined shock front for M2t < a(M21— 1) with a ≈ 0.1). In this regime, the statistics of shock front distortions compare favourably with the linear analysis predictions. For flows with M2t > a(M21— 1 with a ≈ 0.1, shock waves no longer had well-defined fronts: shock wave thickness and strength varied widely along the transverse directions. Multiple compression peaks were found along the mean streamlines at locations where the local shock thickness had increased significantly.

Stability characteristics of a single-phase free convection loop are presented. In the experiments, water was placed inside a toroidal glass loop oriented in a vertical plane. The lower half of the loop was heated and the upper half was cooled. At low heat-transfer rates and also a t high heat-transfer rates the free convection flow was observed to be steady. For the intermediate range, however, the flow was found to be highly oscillatory. Stability predictions are also developed. The comparison between theory and experiment yields favourable agreement.

Observations of unstable behaviour have been reported previously for single-phase fluids in the vicinity of the thermodynamic critical point. In these situations it has been assumed that the unusual behaviour of the fluid properties in the near-critical region necessarily constitutes the underlying cause of such instabilities. In contrast t o this view, analyses by Keller (1966) and Welander (1967) indicate that instabilities can occur for ordinary fluids as well. Results of the present study confirm this contention, since instabilities were clearly observed for water at atmospheric pressure and moderate temperatures.

We study the velocity fields of unforced, high Reynolds number, subsonic jets, issuing from round nozzles with turbulent boundary layers. The objective of the study is to educe wavepackets in such flows and to explore their relationship with the radiated sound. The velocity field is measured using a hot-wire anemometer and a stereoscopic, time-resolved PIV system. The field can be decomposed into frequency and azimuthal Fourier modes. The low-angle sound radiation is measured synchronously with a microphone ring array. Consistent with previous observations, the azimuthal wavenumber spectra of the velocity and acoustic pressure fields are distinct. The velocity spectrum of the initial mixing layer exhibits a peak at azimuthal wavenumbers $m$ ranging from 4 to 11, and the peak is found to scale with the local momentum thickness of the mixing layer. The acoustic pressure field is, on the other hand, predominantly axisymmetric, suggesting an increased relative acoustic efficiency of the axisymmetric mode of the velocity field, a characteristic that can be shown theoretically to be caused by the radial compactness of the sound source. This is confirmed by significant correlations, as high as 10 %, between the axisymmetric modes of the velocity and acoustic pressure fields, these values being significantly higher than those reported for two-point flow–acoustic correlations in subsonic jets. The axisymmetric and first helical modes of the velocity field are then compared with solutions of linear parabolized stability equations (PSE) to ascertain if these modes correspond to linear wavepackets. For all but the lowest frequencies close agreement is obtained for the spatial amplification, up to the end of the potential core. The radial shapes of the linear PSE solutions also agree with the experimental results over the same region. The results suggests that, despite the broadband character of the turbulence, the evolution of Strouhal numbers $0. 3\leq St\leq 0. 9$ and azimuthal modes 0 and 1 can be modelled as linear wavepackets, and these are associated with the sound radiated to low polar angles.

The inviscid instability of columnar vortex flows in unbounded domains to three-dimensional perturbations is considered. The undisturbed flows may have axial and swirl velocity components with a general dependence on distance from the swirl axis. The equation governing the disturbance is found to simplify when the azimuthal wavenumber n is large. This permits us to develop the solution in an asymptotic expansion and reveals a class of unstable modes. The asymptotic results are confirmed by comparisons with numerical solutions of the full problem for a specific flow modelling the trailing vortex. It is found that the asymptotic theory predicts the most-unstable wave with reasonable accuracy for values of n as low as 3, and improves rapidly in accuracy as n increases. This study enables us to formulate a sufficient condition for the instability of columnar vortices as follows. Let the vortex have axial velocity W(r), azimuthal velocity V(r), where r is distance from the axis, let Ω be the angular velocity V/r, and let Γ be the circulation rV. Then the flow is unstable if $ V\frac{d\Omega}{dr}\left[ \frac{d\Omega}{dr}\frac{d\Gamma}{dr} + \left(\frac{dW}{dr}\right)^2\right] < 0.$

Current macroscopic theories of two-phase flow in porous media are based on the extended Darcy's law and an algebraic relationship between capillary pressure and saturation. Both of these equations have been challenged in recent years, primarily based on theoretical works using a thermodynamic approach, which have led to new governing equations for two-phase flow in porous media. In these equations, new terms appear related to the fluid–fluid interfacial area and non-equilibrium capillarity effects. Although there has been a growing number of experimental works aimed at investigating the new equations, a full study of their significance has been difficult as some quantities are hard to measure and experiments are costly and time-consuming. In this regard, pore-scale computational tools can play a valuable role. In this paper, we develop a new dynamic pore-network simulator for two-phase flow in porous media, called DYPOSIT. Using this tool, we investigate macroscopic relationships among average capillary pressure, average phase pressures, saturation and specific interfacial area. We provide evidence that at macroscale, average capillary pressure–saturation–interfacial area points fall on a single surface regardless of flow conditions and fluid properties. We demonstrate that the traditional capillary pressure–saturation relationship is not valid under dynamic conditions, as predicted by the theory. Instead, one has to employ the non-equilibrium capillary theory, according to which the fluids pressure difference is a function of the time rate of saturation change. We study the behaviour of non-equilibrium capillarity coefficient, specific interfacial area, and its production rate versus saturation and viscosity ratio.

A major feature of our pore-network model is a new computational algorithm, which considers capillary diffusion. Pressure field is calculated for each fluid separately, and saturation is computed in a semi-implicit way. This provides more numerical stability, compared with previous models, especially for unfavourable viscosity ratios and small capillary number values.