The transient dynamics of stirred tanks whose impeller speed undergoes smooth or step changes is investigated. First, a low-order model is developed, linking the impeller torque with the ‘extent’ of the solid-body rotation in the tank, derived from an angular momentum balance in a control volume around the impeller. Utilisation of this model enables the prediction of the torque ‘spike’ appearing after an impulsive change of the shaft speed, and of the torque evolution during a quasi-steady transition. For the case of a small impulsive change in the shaft speed, a characteristic spin-up time is also proposed. Torque measurements performed in an unbaffled stirred tank show considerable agreement with the theoretical predictions.

Many transitional wall-bounded shear flows are characterised by the coexistence in state space of laminar and turbulent regimes. Probing the edge boundary between the two attractors has led in the last decade to the numerical discovery of new (unstable) solutions to the incompressible Navier–Stokes equations. However, the iterative bisection method used to compute edge states can become prohibitively costly for large systems. Here we suggest a simple feedback control strategy to stabilise edge states, hence accelerating their numerical identification by several orders of magnitude. The method is illustrated for several configurations of cylindrical pipe flow. Travelling waves solutions are identified as edge states, and can be isolated rapidly in only one short numerical run. A new branch of solutions is also identified. When the edge state is a periodic orbit or chaotic state, the feedback control does not converge precisely to solutions of the uncontrolled system, but nevertheless brings the dynamics very close to the original edge manifold in a single run. We discuss the opportunities offered by the speed and simplicity of this new method to probe the structure of both state space and parameter space.

The behaviour of clean and contaminated bubbles in solid-body rotating flows is compared in terms of drag and lift forces. Both spherical and deformed bubbles are considered. For that comparison, we have completed the data published in Rastello et al. (J. Fluid Mech., vol. 624, 2009, pp. 159–178; J. Fluid Mech., vol. 682, 2011, pp. 434–459) by a new series of measurements. When they are contaminated, bubbles are subject to an additional lift force due to the spinning of their surfaces, while the clean ones are not. A detailed description of this spinning motion is presented and an expression for the Magnus-like lift it induces is given in the light of the new information. The component of the lift induced by flow rotation depends on the Rossby number $Ro$, contrary to the case of clean bubbles. Including the ‘spin’ induced lift component in the dynamical equations provides a better prediction of the bubble’s trajectory in contaminated fluid. The presence of contaminants immobilizes the rear part of the bubble and reduces significantly the deformation. The laws of deformation according to the nature of the surface are presented. The way deformation influences the drag and lift coefficients in pure and contaminated fluids is quantified and discussed. Expressions for these various coefficients are proposed.

This study focuses on the response of premixed flames to a transient hydrodynamic perturbation in an intermediate situation between laminar stretched flames and turbulent flames: an axisymmetric vortex interacting with a flame. The reasons motivating this choice are discussed in the framework of turbulent combustion models and flame response to the stretch rate. We experimentally quantify the dependence of the flame kinematic properties (displacement and consumption speeds) to geometrical scalars (stretch rate and curvature) in flames characterized by different effective Lewis numbers. Whilst the displacement speed can be readily measured using particle image velocimetry and tomographic diagnostics, providing a reliable estimate of the consumption speed from experiments remains particularly challenging. In the present work, a method based on a budget of fuel on a well chosen domain is proposed and validated both experimentally and numerically using two-dimensional direct numerical simulations of flame/vortex interactions. It is demonstrated that the Lewis number impact neither the geometrical nor the kinematic features of the flames, these quantities being much more influenced by the vortex intensity. While interacting with the vortex, the flame displacement (at an isotherm close to the leading edge) and consumption speeds are found to increase almost independently of the type of fuel. We show that the total stretch rate is not the only scalar quantity impacting the flame displacement and consumption speeds and that curvature has a significant influence. Experimental data are interpreted in the light of asymptotic theories revealing the existence of two distinct Markstein numbers, one characterizing the dependence of flame speed to curvature, the other to the total stretch rate. This theory appears to be well suited for representing the evolution of the displacement speed with respect to either the total stretch rate, curvature or strain rate. It also explains the limited dependence of the flame displacement speed to Lewis number and the strong correlation with curvature observed in the experiments. An explicit relationship between displacement and consumption speeds is also given, indicating that the fuel consumption rate is likely to be altered by both the total stretch rate and curvature.

Phytoplankton are the foundation of aquatic food webs. Through photosynthesis, phytoplankton draw down $\text{CO}_{2}$ at magnitudes equivalent to forests and other terrestrial plants and convert it to organic material that is then consumed by other planktonic organisms in higher trophic levels. Mechanisms that affect local concentrations and velocities are of primary significance to many encounter-based processes in the plankton, including prey–predator interactions, fertilization and aggregate formation. We report results from simulations of sinking phytoplankton, considered as elongated spheroids, in homogenous isotropic turbulence to answer the question of whether trajectories and velocities of sinking phytoplankton are altered by turbulence. We show in particular that settling spheroids with physical characteristics similar to those of diatoms weakly cluster and preferentially sample regions of downwelling flow, corresponding to an increase of the mean settling speed with respect to the mean settling speed in quiescent fluid. We explain how different parameters can affect the settling speed and what underlying mechanisms might be involved. Interestingly, we observe that the increase in the aspect ratio of the prolate spheroids can affect the clustering and the average settling speed of particles by two mechanisms: first is the effect of aspect ratio on the rotation rate of the particles, which saturates faster than the second mechanism of increasing drag anisotropy.

Numerical simulations of the turbulent flow over the flat backed Ahmed model at Reynolds number $Re\simeq 4\times 10^{5}$ are conducted using a lattice Boltzmann solver to clarify the mean topology of the static symmetry-breaking mode of the wake. It is shown that the recirculation region is occupied by a skewed low pressure torus, whose part closest to the body is responsible for an extra low pressure imprint on the base. Shedding of one-sided vortex loops is also reported, indicating global quasi-periodic dynamics in conformity with the seminal work of Grandemange et al. (J. Fluid Mech., vol. 722, 2013, pp. 51–84). Despite the limited low frequency resolution of the simulation, power spectra of the lateral velocity fluctuations at different locations corroborate the presence of this quasi-periodic mode at a Strouhal number of $St=0.16\pm 0.03$. A shallow base cavity of $5\,\%$ of the body height reduces the drag coefficient by $3\,\%$ but keeps the recirculating torus and its interaction with the base mostly unchanged. The drag reduction lies in a global constant positive shift of the base pressure distribution. For a deep base cavity of $33\,\%$ of the body height, a drag reduction of $9.5\,\%$ is obtained. It is accompanied by a large elongation of the recirculation inside the cavity that considerably attenuates the low pressure sources therein together with a symmetrization of the low pressure torus. The global quasi-periodic mode is found to be inhibited by the cavity.

This article describes unexplored details of the intriguing spectral manifestation of the small-amplitude waves at the surfaces of a bubble-laden drop. Its natural frequencies of interfacial pulsation reveal a non-trivial variation with the position of the cavity inside the liquid. This configurational dependence of spectra is calculated for arbitrary location of the void by using a novel approach under low capillary number and low Bond number limits. The analysis is based on expansion in two sets of basis functions where their mutual transformations are utilized to enforce interfacial boundary conditions. The obtained results quantify a few important features which have both scientific and technological significance. For a concentric geometry, the inherent azimuthal degeneracy makes the frequencies for a number of vibrational modes exactly the same. For an eccentric position of the bubble, however, this degeneracy disappears, creating small deviations in the spectral values corresponding to different azimuthal modes. Such behaviour is akin to fine-structure split in an atomic system, where different quantum numbers ensure small deviation in energy levels of the states. The formulated mathematical procedure can determine the individual frequency values for the interfacial oscillation even if these are grouped closely together in bands. The paper shows how the number of fine structures inside a band and their specific values can be exploited to predict the size and position of the cavity in an opaque drop without any direct visualization of its interior.

This paper presents a probability density function (PDF) form of the population balance equation (PBE) for polysized and polyshaped droplets and solid particles in turbulent flows. A key contribution of this paper lies in the inclusion of an explicit consideration of the inertial effects and the shape of particles in the PDF-PBE formulation. The number density is taken as a function of droplet or particle size (volume) and shape as well as space and time. Potentially, other particle properties could also be included in the formulation. Inertial effects are quantified through the Stokes number, leading to accurate modelling of the different trajectories that are followed by droplets and/or particles with different sizes and shapes. To treat these effects, a new affordable approach is proposed and referred to as the method of Stokes binning. Here, the inertial dispersed elements are accelerated due to fluid dynamic forces associated with an averaged Stokes number in each bin. The model is validated against two data sets. The first data set includes a series of numerical test cases involving the injection of polyshaped droplets ranging in size from 1 to 50 $\unicode[STIX]{x03BC}\text{m}$ into a turbulent jet resulting in inlet Stokes numbers ranging from 0.03 to 75.2. The second data set consists of an experimental case focusing on the dispersion of 60 and 90 $\unicode[STIX]{x03BC}\text{m}$ spherical droplets in a turbulent round jet, resulting in inlet Stokes numbers of 53 and 122, respectively. The results confirm the ability of the approach to accurately model the polysized and polyshaped droplet dispersion using as few as eight Stokes bins. This approach has the potential to greatly reduce the computational cost of modelling the evolution of inertial droplets and particles in turbulent flows.

Tomographic particle image velocimetry (TPIV) measurement with six high-resolution charge-coupled device (CCD) cameras is conducted to investigate flow structures over a finite circular cylinder with an aspect ratio of 2 ($h/d=2$). This short wall-mounted cylinder is fully immersed in a thick turbulent boundary layer ($\unicode[STIX]{x1D6FF}/h=1.025$). Focus is placed on the three-dimensional instantaneous vortex structures and their dynamic characteristics in the wake flow fields. Based on the present results, a refined topological model of the mean wake field behind the finite circular cylinder is proposed, where the spatial locations of the typical vortex structures and their interactions are described in more detail. Among the reported typical vortex structures (i.e. the horseshoe, tip, base, trailing and arch vortex), emphasis is laid on discussion of the tip and arch vortex. The instantaneous 3D M-shape arch vortex and an alternating large-scale streamwise vortex are first found in the present experiment, and their developments are also discussed. Therefore, it is suggested that the instantaneous finite-cylinder wake is dominated by the arch vortex system and the large-scale streamwise vortices. Moreover, in the instantaneous volumetric flow fields, both the antisymmetric and the symmetric wake behaviours are observed. With proper orthogonal decomposition (POD) analysis, the dynamic characteristics of the wake field are clarified. Different from the flow around an infinite cylinder without control, the third and fourth POD modes are characterized by low-frequency symmetric shedding. The low-frequency feature shown in the second mode pair is observed and associated with the occurrence of instantaneous symmetric 3D wake behaviour triggered by the low-aspect-ratio effect and the extension of the separated shear layer. The low frequency seems be attributed to the flapping phenomenon, i.e. oscillation of the recirculation in the backward-facing step flow. It is found that the flapping motion has a modulating effect on the occurrence of the antisymmetric shedding vortex and thus the large-scale streamwise vortex.

The effect of initial conditions on transition to turbulence is studied in a variable-density shock-driven flow. Richtmyer–Meshkov instability (RMI) evolution of fluid interfaces with two different imposed initial perturbations is observed before and after interaction with a second shock reflected from the end wall of a shock tube (reshock). The first perturbation is a predominantly single-mode long-wavelength interface which is formed by inclining the entire tube to 80$^{\circ }$ relative to the horizontal, yielding an amplitude-to-wavelength ratio, $\unicode[STIX]{x1D702}/\unicode[STIX]{x1D706}=0.088$, and thus can be considered as half the wavelength of a triangular wave. The second interface is multi-mode, and contains additional shorter-wavelength perturbations due to the imposition of shear and buoyancy on the inclined perturbation of the first case. In both cases, the interface consists of a nitrogen-acetone mixture as the light gas over carbon dioxide as the heavy gas (Atwood number, $A\sim 0.22$) and the shock Mach number is $M\approx 1.55$. The initial condition was characterized through Proper Orthogonal Decomposition and density energy spectra from a large set of initial condition images. The evolving density and velocity fields are measured simultaneously using planar laser-induced fluorescence (PLIF) and particle image velocimetry (PIV) techniques. Density, velocity, and density–velocity cross-statistics are calculated using ensemble averaging to investigate the effects of additional modes on the mixing and turbulence quantities. The density and velocity data show that a distinct memory of the initial conditions is maintained in the flow before interaction with reshock. After reshock, the influence of the long-wavelength inclined perturbation present in both initial conditions is still apparent, but the distinction between the two cases becomes less evident as smaller scales are present even in the single-mode case. Several methods are used to calculate the Reynolds number and turbulence length scales, which indicate a transition to a more turbulent state after reshock. Further evidence of transition to turbulence after reshock is observed in the velocity and density fluctuation spectra, where a scaling close to $k^{-5/3}$ is observed for almost one decade, and in the enstrophy fluctuation spectra, where a scaling close to $k^{1/3}$ is observed for a similar range. Also, based on normalized cross correlation spectra, local isotropy is reached at lower wave numbers in the multi-mode case compared with the single-mode case before reshock. By breakdown of large scales to small scales after reshock, rapid decay can be observed in cross-correlation spectra in both cases.

The response of the free liquid surface in a Hele-Shaw cell subjected to a horizontal oscillation is investigated. We study the low-oscillation-amplitude regime and we show, by varying the fluid viscosity, $\unicode[STIX]{x1D708}$, and the forcing frequency, $\unicode[STIX]{x1D714}$, that the ratio between the Stokes viscous length, $\sqrt{2\unicode[STIX]{x1D708}/\unicode[STIX]{x1D714}}$, and the cell thickness greatly affects the amplitude and phase lag of the gravity waves. In particular, the sloshing system undergoes an underdamped/overdamped transition for sufficiently large viscosities. A consistent theoretical model, based on a modification of Darcy’s law to include unsteadiness, is then introduced to rationalize the experimental observations. Contrary to traditional sloshing wave theory, the viscous flow dissipation comes at leading order in the analysis, rather than as a higher-order asymptotic correction to the inviscid sloshing dynamics. The analytical expression for the resonance curves agrees well with experimental results without tunable parameters.

Most commonly used models for turbulent mixing in the ocean rely on a background stratification against which turbulence must work to stir the fluid. While this background stratification is typically well defined in idealized numerical models, it is more difficult to capture in observations. Here, a potential discrepancy in ocean mixing estimates due to the chosen calculation of the background stratification is explored using direct numerical simulation data of breaking internal waves on slopes. Two different methods for computing the buoyancy frequency $N$, one based on a three-dimensionally sorted density field (often used in numerical models) and the other based on locally sorted vertical density profiles (often used in the field), are used to quantify the effect of $N$ on turbulence quantities. It is shown that how $N$ is calculated changes not only the flux Richardson number $R_{f}$, which is often used to parameterize turbulent mixing, but also the turbulence activity number or the Gibson number $Gi$, leading to potential errors in estimates of the mixing efficiency using $Gi$-based parameterizations.

There is a long-standing debate in the literature of stratified flows over topography concerning the correct dimensionless number to refer to as a Froude number. Common definitions using external quantities of the flow include $U/(ND)$, $U/(Nh_{0})$, and $Uk/N$, where $U$ and $N$ are, respectively, scales for the background velocity and buoyancy frequency, $D$ is the depth, and $h_{0}$ and $k^{-1}$ are, respectively, height and width scales of the topography. It is also possible to define an internal Froude number $Fr_{\unicode[STIX]{x1D6FF}}=u_{0}/\sqrt{g^{\prime }\unicode[STIX]{x1D6FF}}$, where $u_{0}$, $g^{\prime }$, and $\unicode[STIX]{x1D6FF}$ are, respectively, the characteristic velocity, reduced gravity, and vertical length scale of the perturbation above the topography. For the case of hydrostatic lee waves in a deep ocean, both $U/(ND)$ and $Uk/N$ are insignificantly small, rendering the dimensionless number $Nh_{0}/U$ the only relevant dynamical parameter. However, although it appears to be an inverse Froude number, such an interpretation is incorrect. By non-dimensionalizing the stratified Euler equations describing the flow of an infinitely deep fluid over topography, we show that $Nh_{0}/U$ is in fact the square of the internal Froude number because it can identically be written in terms of the inner variables, $Fr_{\unicode[STIX]{x1D6FF}}^{2}=Nh_{0}/U=u_{0}^{2}/(g^{\prime }\unicode[STIX]{x1D6FF})$. Our scaling also identifies $Nh_{0}/U$ as the ratio of the vertical velocity scale within the lee wave to the group velocity of the lee wave, which we term the vertical Froude number, $Fr_{vert}=Nh_{0}/U=w_{0}/c_{g}$. To encapsulate such behaviour, we suggest referring to $Nh_{0}/U$ as the lee-wave Froude number, $Fr_{lee}$.