Turbulent flows behind two side-by-side square cylinders with three different gap ratios, namely, $L_{d}/T_{0}=4,$ 6 and 8 ($L_{d}$ is the separation distance between two cylinders and $T_{0}$ is the cylinder thickness) are investigated by using direct numerical simulations. Depending on the strength of the gap flow, the three cases can generally be characterized into two regimes, one being the weak gap flow regime and the other being the robust gap flow regime. The wake-interaction length scale can only be applied to characterize the spatial evolution of the dual-wake flow in the robust gap flow regime. And only in this regime can the so-called ‘extreme events’ (i.e. non-Gaussian velocity fluctuations with large flatness) be identified. For the case with $L_{d}/T_{0}=6$, two downstream locations, i.e. $X/T_{0}=6$ and 26, at which the turbulent flows are highly non-Gaussian distributed and approximately Gaussian distributed, respectively, are analysed in detail. A well-defined $-5/3$ energy spectrum can be found in the near-field region (i.e. $X/T_{0}=6$), where the turbulent flow is still developing and highly intermittent and Kolmogorov’s universal equilibrium hypothesis does not hold. We confirm that the approximate $-5/3$ power law in the high-frequency range is closely related to the occurrences of the extreme events. As the downstream distance increases, the velocity fluctuations gradually adopt a Gaussian distribution, corresponding to a decrease in the strength of the extreme events. Consequently, the range of the $-5/3$ power law narrows. In the upstream region (i.e. $X/T_{0}=6$), the second-order structure function exhibits a power-law exponent close to $1$, whereas in the far downstream region (i.e. $X/T_{0}=26$) the expected $2/3$ power-law exponent appears. The larger exponent at $X/T_{0}=6$ is related to the fact that fluid motions in the intermediate range can directly ‘feel’ the large-scale vortex shedding.

Multiphase flows in porous media represent fluid dynamics problems of great complexity involving a wide range of physical phenomena. These flows have attracted the attention of an impressive group of renowned researchers and have spawned a number of classic problems in fluid dynamics. These multiphase flows are perhaps best known for their importance in oil recovery from petroleum reservoirs, but they also find application in novel areas such as hydrofracturing for natural gas recovery. In a recent article, Zinchenko & Davis (J. Fluid Mech. 2013, vol. 725, pp. 611–663) present computational simulations that break new ground in the study of emulsions flowing through porous media. These simulations provide sufficient scale to capture the disordered motion and complex break-up patterns of individual droplets while providing sufficient statistical samples for estimating meaningful macroscopic properties of technical interest.

An experimental analysis of the vortex-induced vibrations of a hanging string with variable tension along its length is presented in this paper. It is shown that standing waves develop along the hanging string. First, the evolution of the Strouhal number $\mathit{St}$ with the Reynolds number $\mathit{Re}$ follows a trend similar to what is observed for a circular cylinder in a flow for relatively low Reynolds numbers ($32\lt \mathit{Re}\lt 700$). Second, the extracted mode shapes are self-similar: a rescaling of the spanwise coordinate by a self-similarity coefficient allows all of them to collapse onto a unique function. The self-similar behaviour of the spatial distribution of the vibrations along the hanging string is then explained theoretically by performing a linear stability analysis of an adapted wake-oscillator model. This linear stability analysis finally provides an accurate description of the mode shapes and of the evolution of the self-similarity coefficient with the flow speed.

In this, the third, part we present a complete asymptotic analysis of the distribution of velocity and electric potential in an electrically conducting liquid between two circular electrodes of finite diameter, 2b, when a current is passed between them. The electrodes are set opposite to each other in insulating planes, a distance 2a apart, and a magnetic field is applied perpendicular to these planes.

The asymptotic solution is obtained under the restriction that the Hartmann number M satisfies both the conditions M [Gt ] 1 and M½l [Gt ] 1, where l = b/a. It enables us to calculate the distribution of velocity and electrical potential throughout the flow field, and provides an expansion for the resistance R between the electrodes in descending powers of M½l, which is correct provided terms of order M−1 and (lM½)−3 are neglected. Comparison of the theoretical results with experiment shows good agreement in the measurements of R and in the direct probe measurements of electrical potential within the fluid. This is one of the first experiments in which direct probe measurements of an MHD flow, as well as external measurements, have provided such a satisfactory confirmation of the theory. Direct measurements of the velocity, by means of a Pitot tube as in part 2, or by means of a hot-film anemometer undertaken by Malcolm (1968), agree less well with the theory.

The exact equations governing three-dimensional motion of an inviscid non-diffusive incompressible fluid stratified in density or of an inviscid non-diffusive gas stratified in entropy are given and briefly discussed.

This paper develops a stability analysis for the onset of thermoacoustic oscillations in a gas-filled looped tube with a stack inserted, subject to a temperature gradient. Analysis is carried out based on approximate theories for a thermoviscous diffusion layer derived from the thermoacoustic-wave equation taking account of the temperature dependence of the viscosity and the heat conductivity. Assuming that the stack consists of many pores axially and that the thickness of the diffusion layer is much thicker than the pore radius, the diffusion wave equation with higher-order terms included is applied for the gas in the pores of the stack. For the gas outside of the pores, the theory of a thin diffusion layer is applied. In a section called the buffer tube over which the temperature relaxes from that at the hot end of the stack to room temperature, the effects of the temperature gradient are taken into account. With plausible temperature distributions specified on the walls of the stack and the buffer tube, the solutions to the equations in both theories are obtained and a frequency equation is finally derived analytically by matching the conditions at the junctions between the various sections. Seeking a real solution to the frequency equation, marginal conditions of instability are obtained numerically not only for the one-wave mode but also for the two-wave mode, where the tube length corresponds to one wavelength and two wavelengths, respectively. It is revealed that the marginal conditions depend not only on the thickness of the diffusion layer but also on the porosity of the stack. Although the toroidal geometry allows waves to be propagated in both senses along the tube, it is found that the wave propagating in the sense from the cold to the hot end through the stack is always greater, so that a travelling wave in this sense emerges as a whole. The spatial and temporal variations of excess pressure and mean axial velocity averaged over the cross-section of a flow passage are displayed for the two modes of oscillations at the marginal state. The spatial distribution of mean acoustic energy flux (acoustic intensity) over one period is also shown. It is unveiled that the energy flux is generated only in the stack, and it decays slowly in the other sections by lossy effects due to a boundary layer. Mechanisms for the generation of the acoustic energy flux are also discussed.

In a low Morton number ($\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}M$) regime, the stability of a single drop rising in an immiscible viscous liquid is experimentally and computationally examined for varying viscosity ratio $\eta $ (the viscosity of the drop divided by that of the suspending fluid) and varying Eötvös number ($\mathit{Eo}$). Three-dimensional computations, rather than three-dimensional axisymmetric computations, are necessary since non-axisymmetric unstable drop behaviour is studied. The computations are performed using the sharp-interface coupled level-set and volume-of-fluid (CLSVOF) method in order to capture the deforming drop boundary. In the lower $\eta $ regimes, $\eta = 0.02 $ or 0.1, and when $\mathit{Eo}$ exceeds a critical threshold, it is observed that a rising drop exhibits nonlinear lateral/tilting motion. In the higher $\eta $ regimes, $\eta = 0.1$, 1.94, 10 or 100, and when $\mathit{Eo}$ exceeds another critical threshold, it is found that a rising drop becomes unstable and breaks up into multiple drops. The type of breakup, either ‘dumbbell’, ‘intermediate’ or ‘toroidal’, depends intimately on $\eta $ and $\mathit{Eo}$.

Local linear stability analysis is applied to the flow inside a cyclone separator to investigate the unsteady precession of the vortex core. The results of the stability analysis are compared with experimental measurements of the vortex oscillations using high-speed photography with particle seeding and hot-wire anemometry. The experiments reveal distinct spatial variation in the oscillation behaviour within the cyclone. The unsteady motion is focused at each end of the device, at both the narrow cone tip and just below the exhaust duct at the top of the cone, which is known as a vortex finder. The local stability analysis shows that an absolute instability is present throughout the flow for some non-zero azimuthal wavenumbers. The unsteady flow is observed to be driven by coupling between the shear layer and inertial waves confined within the vortex core. Comparison of the stability analysis with experiments shows the same frequency and mode shape behaviour and suggests that the local analysis accurately predicts the unstable modes of the system. The precessing vortex core is responsible for a narrow-band acoustic noise. Comparisons are also drawn with acoustic measurements made on cyclones in which the system is defined by key non-dimensional parameters, such as the swirl number and outlet diameter ratio. The results in this study demonstrate the applicability of local stability analysis to a complex swirling system and yield credible details about the underlying mechanisms of the unstable flow inside the cyclone.

This contribution is a remake of the one published under the same title by Reynolds (Nord. Hydrol., vol. 7, 1976, pp. 161–183). As in that paper, attention is given to the latest developments in the field of river morphodynamics and, in particular, on the formation of bed patterns. Stimulated by the work by Kidanemariam and Uhlmann (J. Fluid Mech., vol. 750, 2014, R2) the opportunity is taken to examine the most recent developments in terms of sediment transport models.

Experiments have shown that a liquid droplet on an inclined plane can be made to move uphill by sufficiently strong, vertical oscillations (Brunet et al., Phys. Rev. Lett., vol. 99, 2007, 144501). In this paper, we study a two-dimensional, inviscid, irrotational model of this flow, with the velocity of the contact lines a function of contact angle. We use asymptotic analysis to show that, for forcing of sufficiently small amplitude, the motion of the droplet can be separated into an odd and an even mode, and that the weakly nonlinear interaction between these modes determines whether the droplet climbs up or slides down the plane, consistent with earlier work in the limit of small contact angles (Benilov and Billingham, J. Fluid Mech. vol. 674, 2011, pp. 93–119). In this weakly nonlinear limit, we find that, as the static contact angle approaches $\unicode[STIX]{x03C0}$ (the non-wetting limit), the rise velocity of the droplet (specifically the velocity of the droplet averaged over one period of the motion) becomes a highly oscillatory function of static contact angle due to a high frequency mode that is excited by the forcing. We also solve the full nonlinear moving boundary problem numerically using a boundary integral method. We use this to study the effect of contact angle hysteresis, which we find can increase the rise velocity of the droplet, provided that it is not so large as to completely fix the contact lines. We also study a time-dependent modification of the contact line law in an attempt to reproduce the unsteady contact line dynamics observed in experiments, where the apparent contact angle is not a single-valued function of contact line velocity. After adding lag into the contact line model, we find that the rise velocity of the droplet is significantly affected, and that larger rise velocities are possible.

The influence of the large scale organisation of free-stream turbulence on a turbulent boundary layer is investigated experimentally in a wind tunnel through hot-wire measurements. An active grid is used to generate high-intensity free-stream turbulence with turbulence intensities and local turbulent Reynolds numbers in the ranges $7.2\,\%\leqslant u_{\infty }^{\prime }/U_{\infty }\leqslant 13.0\,\%$ and $302\leqslant Re_{\unicode[STIX]{x1D706},\infty }\leqslant 760$, respectively. In particular, several cases are produced with fixed $u_{\infty }^{\prime }/U_{\infty }$ and $Re_{\unicode[STIX]{x1D706},\infty }$, but up to a 65 % change in the free-stream integral scale $L_{u,\infty }/\unicode[STIX]{x1D6FF}$. It is shown that, while qualitatively the spectra at various wall-normal positions in the boundary layer look similar, there are quantifiable differences at the large wavelengths all the way to the wall. Nonetheless, profiles of the longitudinal statistics up to fourth order are well collapsed between cases at the same $u_{\infty }^{\prime }/U_{\infty }$. It is argued that a larger separation of the integral scale would not yield a different result, nor would it be physically realisable. Comparing cases across the wide range of turbulence intensities and free-stream Reynolds numbers tested, it is demonstrated that the near-wall spectral peak is independent of the free-stream turbulence, and seemingly universal. The outer peak was also found to be described by a set of global scaling laws, and hence both the near-wall and outer spectral peaks can be predicted a priori with only knowledge of the free-stream spectrum, the boundary layer thickness ($\unicode[STIX]{x1D6FF}$) and the friction velocity ($U_{\unicode[STIX]{x1D70F}}$). Finally, a conceptual model is suggested that attributes the increase in $U_{\unicode[STIX]{x1D70F}}$ as $u_{\infty }^{\prime }/U_{\infty }$ increases to the build-up of energy at large wavelengths near the wall because that energy cannot be transferred to the universal near-wall spectral peak.

Simulations of strongly stratified turbulence often exhibit coherent large-scale structures called vertically sheared horizontal flows (VSHFs). VSHFs emerge in both two-dimensional (2D) and three-dimensional (3D) stratified turbulence with similar vertical structure. The mechanism responsible for VSHF formation is not fully understood. In this work, the formation and equilibration of VSHFs in a 2D Boussinesq model of stratified turbulence is studied using statistical state dynamics (SSD). In SSD, equations of motion are expressed directly in the statistical variables of the turbulent state. Restriction to 2D turbulence facilitates application of an analytically and computationally attractive implementation of SSD referred to as S3T, in which the SSD is expressed by coupling the equation for the horizontal mean structure with the equation for the ensemble mean perturbation covariance. This second-order SSD produces accurate statistics, through second order, when compared with fully nonlinear simulations. In particular, S3T captures the spontaneous emergence of the VSHF and associated density layers seen in simulations of turbulence maintained by homogeneous large-scale stochastic excitation. An advantage of the S3T system is that the VSHF formation mechanism, which is wave–mean flow interaction between the emergent VSHF and the stochastically excited large-scale gravity waves, is analytically understood in the S3T system. Comparison with fully nonlinear simulations verifies that S3T solutions accurately predict the scale selection, dependence on stochastic excitation strength, and nonlinear equilibrium structure of the VSHF. These results constitute a theory for VSHF formation applicable to interpreting simulations and observations of geophysical examples of turbulent jets such as the ocean’s equatorial deep jets.

We investigate the growth of gas bubbles in a water solution at rest with a supersaturation level that is generally associated with diffusive mass transfer. For $\text{CO}_{2}$ bubbles, it has been previously observed that, after some time of growing in a diffusive regime, a density-driven convective flow enhances the mass transfer rate into the bubble. This is due to the lower density of the gas-depleted liquid which surrounds the bubble. In this work, we report on experiments with different supersaturation values, measuring the time $t_{conv}$ it takes for convection to dominate over the diffusion-driven growth. We demonstrate that by considering buoyancy and drag forces on the depleted liquid around the bubble, we can satisfactorily predict the transition time. In fact, our analysis shows that this onset does not only depend on the supersaturation, but also on the absolute pressure, which we corroborate in experiments. Subsequently, we study how the depletion caused by the growth of successive single bubbles influences the onset of convection. Finally, we study the convection onset around diffusively growing nitrogen $\text{N}_{2}$ bubbles. As $\text{N}_{2}$ is much less soluble in water, the growth takes much longer. However, after waiting long enough and consistent with our theory, convection still occurs as for any gas–liquid combination, provided that the density of the solution sufficiently changes with the gas concentration.