A model problem is analysed to study the microscopic flow near the surface of two-dimensional porous media. In the idealized problem we consider axial flow through infinite and semi-infinite lattices of cylindrical inclusions. The effect of lattice geometry and inclusion shape on the permeability and surface flow are examined. Calculations show that the definition of a slip coefficient for a porous medium is meaningful only for extremely dilute systems. Brinkman's equation gives reasonable predictions for the rate of decay of the mean velocity for certain simple geometries, but fails for to predict the correct behaviour for media anisotropic in the plane normal to the flow direction.

Based on energy principles, we propose a statistical model to describe the bubble size probability density function of the daughter bubbles resulting from the shattering of a mother bubble of size D0 immersed in a fully developed turbulent water flow. The model shows that the bubble size p.d.f. depends not only on D0, but also on the value of the dissipation rate of turbulent kinetic energy of the underlying turbulence of the water, ε. The phenomenological model is simple, yet it predicts detailed experimental measurements of the transient bubble size p.d.f.s performed over a range of bubble sizes and dissipation rates ε in a very consistent manner. The agreement between the model and the experiments is particularly good for low and moderate bubble turbulent Weber numbers, Wet = ρΔu2(D0)D0/σ where the assumption of the binary breakup is shown to be consistent with the experimental observations. At larger values of Wet, it was found that the most probable number of daughter bubbles increases and the assumption of tertiary breakup is shown to lead to a better fit of the experimental measurements.

Direct numerical simulation (DNS) of the near field of a three-dimensional spatially developing turbulent lifted hydrogen jet flame in heated coflow is performed with a detailed mechanism to determine the stabilization mechanism and the flame structure. The DNS was performed at a jet Reynolds number of 11,000 with over 940 million grid points. The results show that auto-ignition in a fuel-lean mixture at the flame base is the main source of stabilization of the lifted jet flame. A chemical flux analysis shows the occurrence of near-isothermal chemical chain branching preceding thermal runaway upstream of the stabilization point, indicative of hydrogen auto-ignition in the second limit. The Damköhler number and key intermediate-species behaviour near the leading edge of the lifted flame also verify that auto-ignition occurs at the flame base. At the lifted-flame base, it is found that heat release occurs predominantly through ignition in which the gradients of reactants are opposed. Downstream of the flame base, both rich-premixed and non-premixed flames develop and coexist with auto-ignition. In addition to auto-ignition, Lagrangian tracking of the flame base reveals the passage of large-scale flow structures and their correlation with the fluctuations of the flame base. In particular, the relative position of the flame base and the coherent flow structure induces a cyclic motion of the flame base in the transverse and axial directions about a mean lift-off height. This is confirmed by Lagrangian tracking of key scalars, heat release rate and velocity at the stabilization point.

Combined thermocapillary–buoyancy convection in a thin rectangular geometry is investigated experimentally, with an emphasis on the generation of hydrothermal-wave instabilities. For sufficiently thin layers, pure hydrothermal waves are observed, and are found to be oblique as predicted by a previous linear-stability analysis (Smith & Davis 1983). For thicker layers, both a steady multicell state and an oscillatory state are found to exist, but the latter is not in the form of a pure hydrothermal wave.

The instant of transition to cellularity of centrally ignited, outwardly propagating spherical flames in a reactive environment of fuelx–oxidizer mixture, at atmospheric and elevated pressures, was experimentally determined using high-speed schlieren imaging and subsequently interpreted on the basis of hydrodynamic and diffusional–thermal instabilities. Experimental results show that the transition Péclet number, Pec = RcℓL, assumes an almost constant value for the near-equidiffusive acetylene flames with wide ranges in the mixture stoichiometry, oxygen concentration and pressure, where Rc is the flame radius at transition and ℓL the laminar flame thickness. However, for the non-equidiffusive hydrogen and propane flames, Pec respectively increases and decreases somewhat linearly with the mixture equivalence ratio. Evaluation of Pec using previous theory shows complete qualitative agreement and satisfactory quantitative agreement, demonstrating the insensitivity of Pec to all system parameters for equidiffusive mixtures, and the dominance of the Markstein number, Ze(Le – 1), in destabilization for non-equidiffusive mixtures, where Ze is the Zel'dovich number and Le the Lewis number. The importance of using locally evaluated values of ℓL, Ze and Le, extracted from either computationally simulated one-dimensional flame structure with detailed chemistry and transport, or experimentally determined response of stretched flames, in the evaluation of Pec is emphasized.

We present experimental results for the acoustic field of jets with Mach numbers between 0.35 and 0.6. An azimuthal ring array of six microphones, whose polar angle, , was progressively varied, allows the decomposition of the acoustic pressure into azimuthal Fourier modes. In agreement with past observations, the sound field for low polar angles (measured with respect to the jet axis) is found to be dominated by the axisymmetric mode, particularly at the peak Strouhal number. The axisymmetric mode of the acoustic field can be clearly associated with an axially non-compact source, in the form of a wavepacket: the sound pressure level for peak frequencies is found be superdirective for all Mach numbers considered, with exponential decay as a function of , where is the Mach number based on the phase velocity of the convected wave. While the mode spectrum scales with Strouhal number, suggesting that its energy content is associated with turbulence scales, the axisymmetric mode scales with Helmholtz number – the ratio between source length scale and acoustic wavelength. The axisymmetric radiation has a stronger velocity dependence than the higher-order azimuthal modes, again in agreement with predictions of wavepacket models. We estimate the axial extent of the source of the axisymmetric component of the sound field to be of the order of six to eight jet diameters. This estimate is obtained in two different ways, using, respectively, the directivity shape and the velocity exponent of the sound radiation. The analysis furthermore shows that compressibility plays a significant role in the wavepacket dynamics, even at this low Mach number. Velocity fluctuations on the jet centreline are reduced as the Mach number is increased, an effect that must be accounted for in order to obtain a correct estimation of the velocity dependence of sound radiation. Finally, the higher-order azimuthal modes of the sound field are considered, and a model for the low-angle sound radiation by helical wavepackets is developed. The measured sound for azimuthal modes 1 and 2 at low Strouhal numbers is seen to correspond closely to the predicted directivity shapes.

Conditionally sampled measurements with particle image velocimetry (PIV) of a turbulent round submerged liquid jet in a laboratory have been taken at Re = 2 × 103 between 60 and 100 nozzle diameters from the nozzle in order to investigate the dynamics and transport processes at the continuous and well-defined bounding interface between the turbulent and non-turbulent regions of flow. The jet carries a fluorescent dye measured with planar laser-induced fluorescence (LIF), and the surface discontinuity in the scalar concentration is identified as the fluctuating turbulent jet interface. Thence the mean outward ‘boundary entrainment’ velocity is derived and shown to be a constant fraction (about 0.07) of the the mean jet velocity on the centreline. Profiles of the conditional mean velocity, mean scalar and momentum flux show that at the interface there are clear discontinuities in the mean axial velocity and mean scalar and a tendency towards a singularity in mean vorticity. These actual or asymptotic discontinuities are consistent with the conditional mean momentum and scalar transport equations integrated across the interface. Measurements of the fluxes of turbulent kinetic energy and enstrophy are consistent with computations by Mathew & Basu (Phys. Fluids, vol. 14, 2002, pp. 2065–2072) in showing that for a jet flow (without forcing) the entrainment process is dominated by small-scale eddying at the highly sheared interface (‘nibbling’), with large-scale engulfing making a small (less than 10%) contribution consistent with concentration measurements showing that the interior of the jet is well mixed. (Turbulent jets differ greatly from the free shear layer in this respect.) To explain the difference between velocity and scalar profiles, their conditional mean gradients are defined in terms of a local eddy viscosity and eddy diffusivity and the momentum and scalar fluxes inside the interface. Since the eddy diffusivity is larger than the eddy viscosity, the scalar profile is flatter inside the interface so that the scalar discontinuity is relatively greater than the mean velocity discontinuity. Theoretical arguments, following Hunt, Eames & Westerweel (in Proc. of the IUTAM Symp. on Computational Physics and New Perspectives in Turbulence, ed. Y. Kaneda, vol. 4, 2008, pp. 331–338, Springer), are proposed for how the vortex sheet develops, how the internal structure of the interface layer relates to the inhomogeneous rotational and irrotational motions on each side and why the dominant entrainment process of jets and wakes differs from that of free shear layers.

Burst structures in the near wall region of turbulent flows are associated with a large portion of the turbulent momentum transport from the wall. However, quantitative measures of the timescales associated with the burst event are not well defined, largely due to ambiguities associated with the methods used to detect a burst.

In the present study, Eulerian burst-detection schemes were developed through extensions of the uv quadrant 2, VITA, and u-level techniques. Each of the basic techniques detects ejections. One or more ejections are contained in each burst and hence the key idea is to identify and to group those ejections from a single burst into a single-burst detection. When the ejection detections were grouped appropriately into burst detections, all of the extended techniques yielded the same average time between bursts as deduced from flow visualization for fully-developed channel flow in the range 8700 [les ] Reh [les ] 17 800. The present results show that inner variables (wall shear stress and kinematic viscosity) are the best candidates for the proper scaling of the average time between bursts. Conditional velocity sampling during burst and ejection detections shows that these burst events are closely correlated with slower-than-average moving fluid, moving both away from the wall and toward the wall.

This paper presents an infinite-series solution to the creeping viscous motion of a fluid through low- and moderate-aspect-ratio pores. The flow field is divided into two simply bounded regions: a cylindrical volume bounded by the walls of the pore and the entrance and exit planes, and an infinite half-space outside the pore. Analytic solutions are first obtained in each region for unknown functions representing arbitrary axial and radial velocity profiles at the pore entrance (exit). These unknown functions are then determined by matching the normal and tangential stress at the pore opening.

The results indicate that the velocity profile approaches to within 1·5 per cent of a Poiseuille profile after a short entrance distance of half the pore radius. In the far field the solution matches exactly the streamline pattern for a flow through an orifice of zero thickness obtained by Sampson (1891). The pressure drop across the pore exhibits linear dependence on the aspect ratio and is closely approximated (less than one per cent error) by a simple algebraic expression.

Many experimenters have used oscillating grids to produce turbulence for various laboratory purposes, especially in studies of mixing, but there have been few direct measurements of the properties of the turbulence itself. In the present paper we report experiments which attempt to relate the turbulent velocity and length scales to the external parameters, the frequency and amplitude, for three forms of grid oscillated in a tank of water. Turbulent velocities have been measured in the absence of a mean flow by using a hot film moved through the fluid to provide its own mean velocity. The output is stored and analysed in a small computer, which rapidly evaluates velocity and length scale statistics from an ensemble of records. The spatial variation of these quantities with distance from the stirrer is of special interest. It agrees with results suggested by an inertial-decay theory, and with previous measurements made by Bouvard & Dumas (1967) using a different form of stirrer. A particular purpose of the work has been to ‘calibrate’ the entrainment experiments of Turner (1968), by providing absolute scales of velocity and length in the fluid near a mixing interface, for the same grid as was used in the earlier experiments. Evidence is presented which suggests that other forms of grid may not be calibrated simply by extrapolating these results.

The high-Reynolds-number behaviour of the canonical incompressible turbulent channel flow is investigated through large-scale direct numerical simulation (DNS). A Reynolds number is achieved ($Re_{\tau } = h/\delta _v \approx 4000$, where $h$ is the channel half-height, and $\delta _v$ is the viscous length scale) at which theory predicts the onset of phenomena typical of the asymptotic Reynolds number regime, namely a sensible layer with logarithmic variation of the mean velocity profile, and Kolmogorov scaling of the velocity spectra. Although higher Reynolds numbers can be achieved in experiments, the main advantage of the present DNS study is access to the full three-dimensional flow field. Consistent with refined overlap arguments (Afzal & Yajnik, J. Fluid Mech. vol. 61, 1973, pp. 23–31; Jiménez & Moser, Phil. Trans. R. Soc. Lond. A, vol. 365, 2007, pp. 715–732), our results suggest that the mean velocity profile never achieves a truly logarithmic profile, and the logarithmic diagnostic function instead exhibits a linear variation in the outer layer whose slope decreases with the Reynolds number. The extrapolated value of the von Kármán constant is $k \approx 0.41$. A near logarithmic layer is observed in the spanwise velocity variance, as predicted by Townsend’s attached eddy hypothesis, whereas the streamwise variance seems to exhibit a shoulder, perhaps being still affected by low-Reynolds-number effects. Comparison with previous DNS data at lower Reynolds number suggests enhancement of the imprinting effect of outer-layer eddies onto the near-wall region. This mechanisms is associated with excess turbulence kinetic energy production in the outer layer, and it reflects in flow visualizations and in the streamwise velocity spectra, which exhibit sharp peaks in the outer layer. Associated with the outer energy production site, we find evidence of a Kolmogorov-like inertial range, limited to the spanwise spectral density of $u$, whereas power laws with different exponents are found for the other spectra. Finally, arguments are given to explain the ‘odd’ scaling of the streamwise velocity variances, based on the analysis of the kinetic energy production term.

In addition to the random processes described as horizontal turbulence, there are certain more regular processes which may contribute to horizontal mixing. One of these occurs in a shearing current, where the vertical gradient of velocity combined with vertical turbulent mixing leads to an effective diffusion in the horizontal direction. It is shown that this effect occurs in an alternating flow, such as a tidal current, as well as in a steady flow. In estuaries and coastal waters horizontal mixing by the shear effect may be associated with tidal currents, density currents or wind-driven currents. In each case the effective coefficient of horizontal diffusion, Kx or Ky, is inversely proportional to the coefficient of vertical eddy diffusion Kz. The occurrence of a stable gradient of density therefore increases the effective horizontal mixing very considerably. Results obtained from observations in the Mersey estuary and Irish Sea are compared with theoretical estimates of these effects.

This experimental study describes the turbulent wake behind a two-dimensional porous obstruction, consisting of a circular array of cylinders. The cylinders extend from the channel bed through the water surface, mimicking a patch of emergent vegetation. Three patch diameters () and seven solid volume fractions () are tested. Because flow can pass through the patch, directly downstream there is a region of steady, non-zero, streamwise velocity, , called the steady wake. For the patch diameters and solid volume fractions considered here, is a function of only. The length of the steady wake () increases as decreases and can be predicted from the growth of a plane shear layer. The formation of the von-Kármán vortex street is delayed until the end of the steady wake. There are two regions of elevated transverse velocity fluctuation (): directly behind the patch, associated with the wake turbulence of individual cylinders; and at the distance from the patch, associated with the formation of large-scale wake oscillation. Velocity along the centreline of the wake starts to increase only after the patch-scale vortex street is formed, and it approaches the free-stream velocity over a distance . The dimensionless length of the entire wake, , increases with patch porosity.

It is shown that the fundamental features of both thermal instabilities and the corresponding nonlinear convection in rapidly rotating spherical systems (in the range of the Taylor number 109 < T < 1012) are determined by the fluid properties characterized by the size of the Prandtl number. Coefficients of the asymptotic power law for the onset of convection at large Taylor number are estimated in the range of the Prandtl number 0.1 ≤ Pr ≤ 100. For fluids of moderately small Prandtl number, a new type of convective instability in the form of prograde spiralling drifting columnar rolls is discovered. The linear columnar rolls extend spirally from near latitude 60° to the equatorial region, and each spans azimuthally approximately five wavelengths with the inclination angle between a spirally elongated roll and the radial direction exceeding 45°. As a consequence, the radial lengthscale of the linear roll becomes comparable with the azimuthal lengthscale. A particularly significant finding is the connection between the new instability and the predominantly axisymmetric convection. Though non-axisymmetric motions are preferred at the onset of convection, the nonlinear convection (at the Rayleigh number of the order of (R—Rc)/Rc = O(0.1)) bifurcating supercritically from the spiralling mode is primarily dominated by the component of the axisymmetric zonal flow, which contains nearly 90% of the total kinetic energy. For fluids of moderately large Prandtl numbers, thermal instabilities at the onset of convection are concentrated in a cylindrical annulus coaxial with the axis of rotation; the position of the convection cylinder is strongly dependent on the size of the Prandtl number. The associated nonlinear convection consists of predominantly non-axisymmetric columnar rolls together with a superimposed weak mean flow that contains less than 10% of the total kinetic energy at (R—Rc)/Rc = O(0.1). A double-layer structure of the temperature field (with respect to the basic state) forms as a result of strong nonlinear interactions between the nonlinear flow and the temperature field. It is also demonstrated that the aspect ratio of the spherical shell does not substantially influence the fundamental properties of convection.

The low-frequency dynamics of the shock-induced separation region in a Mach 2 compression ramp interaction is investigated by performing high-speed particle image velocimetry (HSPIV) measurements, at a rate of 6kHz, in a streamwise–spanwise plane. The HSPIV measurements made in the upstream turbulent boundary layer indicate the presence of spanwise strips of elongated regions of uniform streamwise velocity that extend to lengths greater than 30δ, validating previous results based on planar laser scattering measurements obtained by Ganapathisubramani, Clemens & Dolling (J. Fluid Mech., vol. 585, 2007, p. 369). At a wall normal-location of y/δ=0.2, a surrogate for separation based on a velocity threshold is found to fluctuate over a streamwise range of ±1.2δ, consistent with previous studies. The amplitude of unsteadiness has contributions from at least two sources that are related to the incoming boundary layer. First, the velocity threshold based surrogate separation line exhibits large-scale undulations along the spanwise direction that conform to the passage of elongated low- and high-speed regions in the upstream boundary layer. This motion is classified as the local influence of the upstream boundary layer. Second, the spanwise-averaged surrogate separation is found to respond to the overall change in streamwise velocity in the incoming boundary layer and is classified as the global influence of the upstream boundary layer. However, this global influence includes the contributions from the elongated low- and high-speed regions. Preliminary findings based on statistical analysis suggest that the local influence contributes nearly 50% more than the global influence. Regardless, the low-frequency unsteadiness of the separation-region can be attributed to the local and global influences of the incoming boundary layer.

A study has been made of the wake of a cylinder vibrating in line with an incident steady flow. The Reynolds number for the experiments was 190, and the vortex shedding was at all times synchronized with the vibrations of the cylinder, which were in a range of frequencies near twice the Strouhal shedding frequency for the stationary cylinder. Two distinct vortex wake patterns were encountered. The first is a complex regime in which two vortices are shed during each cycle of the vibration and form an alternating pattern of vortex pairs downstream. The second pattern is an alternating street which results from the shedding of a single vortex during each cycle of the cylinder's motion. The street geometry in the latter case shares many basic characteristics with the wake of a cylinder vibrating in cross-flow. These include the effects of vibration amplitude and frequency on the longitudinal and transverse spacing of the vortices. The results obtained from these experiments in air are in agreement with previous findings from free- and forced-vibration experiments in water at both higher and lower Reynolds numbers.

Perturbations of channel geometry (like variations of channel curvature or channel width) in meandering rivers give rise to morphodynamic effects which display themselves through the development of large-scale perturbations of bottom topography in the form of stationary bars developing in the longitudinal direction. The latter may then drive the lateral migration of the channel by enhancing bank erosion at bar pools: through this mechanism local perturbations of channel geometry may affect the planimetric development of meandering rivers on large timescales. The problem tackled herein is whether such morphodynamic influence is invariably felt downstream as the commonly employed model of river meandering would suggest.

In order to solve this problem, we derive the exact solution of the linearized form of the mathematical problem of river morphodynamics. Linear analysis had pointed out the existence of a resonance phenomenon: in a linear (hence ideal) context, resonance occurs when the meander wavenumber and the width ratio of the channel take values (λR and βR, respectively) such as to force free spatial modes of the system consisting of free bars which neither grow nor decay either in time or in space. Channels characterized by values of the width ratio β larger (smaller) than βR are called super- (sub-)resonant. The present solution, which applies to channels with constant width and arbitrary curvature distribution, shows that two distinct scenarios may occur: downstream influence is associated with sub-resonant channels and vice versa dominant upstream influence occurs in super-resonant channels. Small-amplitude waves of bottom topography are shown to migrate downstream in the former case and may migrate upstream in the latter, as resonance also defines the threshold conditions below (above) which small-amplitude alternate bar perturbations (may) migrate downstream (upstream).

These results have several implications. In the present paper we examine the overdeepening phenomenon whereby abrupt variations of channel curvature, as in sequences of straight and constant curvature reaches, lead to sequences of stationary alternate bars with amplitude decaying in the longitudinal direction. We show that, along with downstream overdeepening, an upstream overdeepening scenario is predicted in the super-resonant regime.

Implications of the upstream influence on planimetric development of meandering rivers are investigated in Part 2.

Nonlinear self-excited thermoacoustic oscillations appear in systems involving confined combustion in the form of coupled acoustic pressure oscillations and unsteady heat release rate. In this paper, we investigate the nonlinear transition undergone by thermoacoustic oscillations to flame blowout via intermittency, in response to variation in the location of the combustion source with respect to the acoustic field of the confinement. A ducted laminar premixed conical flame, stabilized on a circular jet exit with a fully developed exit velocity profile, was investigated. Transition to limit cycle oscillations from a non-oscillatory state was observed to occur via a subcritical Hopf bifurcation. Limit cycle oscillations underwent a further bifurcation to quasi-periodic oscillations characterized by the repeated formation of elongated necks in the flame that pinch off as pockets of unburned fuel–air mixture. The quasi-periodic state loses stability, resulting in an intermittent state identified as type II through recurrence analysis of phase space trajectories reconstructed from the acoustic pressure time trace. In this state, the flame undergoes repeated lift-off and reattachment. Instantaneous flame images suggest that the intermittent flame behaviour is influenced by jet flow dynamics.

We present a stochastic analysis of a data set consisting of 1.25 × 107 samples of the local velocity measured in the turbulent region of a round free jet. We find evidence that the statistics of the longitudinal velocity increment v(r) can be described as a Markov process. This new approach to characterize small-scale turbulence leads to a Fokker–Planck equation for the r-evolution of the probability density function (p.d.f.) of v(r). This equation for p(v, r) is completely determined by two coefficients D1(v, r) and D2(v, r) (drift and diffusion coefficient, respectively). It is shown how these coefficients can be estimated directly from the experimental data without using any assumptions or models for the underlying stochastic process. The solutions of the resulting Fokker–Planck equation are compared with experimentally determined probability density functions. It is shown that the Fokker–Planck equation describes the measured p.d.f.(s) correctly, including intermittency effects. Furthermore, knowledge of the Fokker–Planck equation also allows the joint probability density of N increments on N different scales p(v1, r1, …, vN, rN) to be determined.

The sparse identification of nonlinear dynamics (SINDy) is a recently proposed data-driven modelling framework that uses sparse regression techniques to identify nonlinear low-order models. With the goal of low-order models of a fluid flow, we combine this approach with dimensionality reduction techniques (e.g. proper orthogonal decomposition) and extend it to enforce physical constraints in the regression, e.g. energy-preserving quadratic nonlinearities. The resulting models, hereafter referred to as Galerkin regression models, incorporate many beneficial aspects of Galerkin projection, but without the need for a high-fidelity solver to project the Navier–Stokes equations. Instead, the most parsimonious nonlinear model is determined that is consistent with observed measurement data and satisfies necessary constraints. Galerkin regression models also readily generalize to include higher-order nonlinear terms that model the effect of truncated modes. The effectiveness of such an approach is demonstrated on two canonical flow configurations: the two-dimensional flow past a circular cylinder and the shear-driven cavity flow. For both cases, the accuracy of the identified models compare favourably against reduced-order models obtained from a standard Galerkin projection procedure. Finally, the entire code base for our constrained sparse Galerkin regression algorithm is freely available online.