This study investigated the turbulent wake flow behind a flat-back Ahmed body using a combination of time-resolved tomographic particle image velocimetry measurements and spectral proper orthogonal decomposition (SPOD). The experiments were conducted at a Reynolds number of ReH = 10 000, which is defined as U∞H/ν, where U∞ represents the free-stream velocity, H is the height of the body and ν represents the viscosity. The SPOD analysis revealed four distinct flow motions present in the wake of the Ahmed body, each occupying a specific range of Strouhal number, StH. Here, StH is defined as f × H/U∞, where f is the frequency of the motion. At the lowest resolved StH of 0.007, the system exhibited a bi-stability mode, in which the wake switched between asymmetric states consisting of a tilted toroidal vortex and a streamwise vortex. In the next StH range of 0.014 to 0.123, the flow demonstrated swinging/flapping motions, characterized by small spanwise and vertical movements of the wake barycentre. These movements were attributed to the tilting of the toroidal vortex. The third category included the vortex-shedding motions and consisted of quasi-streamwise vortices. These vortices advected in the downstream direction, causing oblique displacements of the barycentre in the streamwise-vertical plane. The peak energy of the vortex shedding was observed at StH = 0.164. Finally, shear layer instabilities induced small vertical and spanwise velocity fluctuations along the shear layers at a high StH of 1.147.

This paper investigates the transport of drugs delivered by direct injection into the cerebrospinal fluid (CSF) that fills the intrathecal space surrounding the spinal cord. Because of the small drug diffusivity, the dispersion of neutrally buoyant drugs has been shown in previous work to rely mainly on the mean Lagrangian flow associated with the CSF oscillatory motion. Attention is given here to effects of buoyancy, arising when the drug density differs from the CSF density. For the typical density differences found in applications, the associated Richardson number is shown to be of order unity, so that the Lagrangian drift includes a buoyancy-induced component that depends on the spatial distribution of the drug, resulting in a slowly evolving cycle-averaged flow problem that can be analysed with two-time scale methods. The asymptotic analysis leads to a nonlinear integro-differential equation for the spatiotemporal solute evolution that describes accurately drug dispersion at a fraction of the cost involved in direct numerical simulations of the oscillatory flow. The model equation is used to predict drug dispersion of positively and negatively buoyant drugs in an anatomically correct spinal canal, with separate attention given to drug delivery via bolus injection and constant infusion.

We study coherent structures in subsonic turbulent jets subject to a flight stream. A thorough characterisation of the effects of a flight stream on the turbulent field was recently performed by Maia et al. (Phys. Rev. Fluids, vol. 8, 2023, 063902) and fluctuation energy attenuations were observed over a broad range of frequencies and azimuthal wavenumbers. The Kelvin–Helmholtz, Orr and lift-up mechanisms were all shown to be weakened by the flight stream. Here we expand upon that study and model the changes in the dynamics of jets in flight using global resolvent analysis. The resolvent model is found to correctly capture the main effects of the flight stream on the dynamics of coherent structures, which are educed from a large-eddy simulation database using spectral proper orthogonal decomposition. Three modifications of note are: the damping of low-frequency streaky/Orr structures that carry most of the fluctuation energy; a degradation of the low-rank behaviour of the jet in frequencies where modal instability mechanisms are dominant; and a rank decrease at very low Strouhal numbers. The latter effect is underpinned by larger gain separations predicted by the resolvent analysis, due to a reduction in the wavelength of associated flow structures. This leads to a clearer relative dominance of streaky structures generated by the lift-up mechanism, despite the fact that the lift-up mechanism has been weakened with respect to the static jet.

We conducted experiments in a laboratory to study turbulent flow over wind generated water waves. The experiments were performed in a wind-wave-current flume with three free stream wind speeds of Uref = 6.0, 8.0 and 10.0 m s−1, corresponding to 10 m equivalent wind speed of U10 = 10.2, 12.2 and 14.1 m s−1 and the root-mean-square wave height of 0.7, 1.1 and 1.7 cm, respectively, at a fetch of 6.2 m. The instantaneous velocity fields above the waves were obtained by using a particle image velocimetry (PIV) technique. The velocity fields were decomposed into the mean, wave-induced and turbulent velocity components. The tested wind waves were primarily dissipated by capillaries and microscale breaking waves. The Bond number and the shear velocity-fetch based Reynolds number were found to correlate with the wind wave regimes well. The turbulent dissipation rates above the water surface were determined based on resolved spatial gradient of instantaneous velocities, where the time-averaged dissipation rate values were calibrated using those estimated from the one-dimensional velocity spectrum in the temporal space. Subsequently, the turbulent kinetic energy (TKE) budget including its production, dissipation, advection and turbulent transport was presented. In addition, conditional averaging analysis of the TKE budgets over leeward, windward sides and all phases was performed. The results showed a strong dependency with the wave phase in the TKE budget terms except for the dissipation. The production-dissipation ratio increased significantly as the wind speed increased, likely attributed to the increased roughness over the substantial coverage of micro-breaking waves.

The coherent dynamics of bubble clusters are of fundamental and industrial importance, and are elusive due to the complex interactions of disordered bubble oscillations. Here we introduce and demonstrate a method for decomposition of the Lagrangian time series of bubble dynamics data by combining theory and principal component analysis. The decomposition extracts coherent features of bubble oscillations based on their energy, in a way similar to proper orthogonal decomposition of Eulerian flow field data. This method is applied to a dataset of spherical clusters under harmonic excitation at different amplitudes, with various nuclei density and polydispersity parameters. Results indicate that the underlying correlated mode of oscillations is isolated in a single dominant feature in cavitating regimes, independent of the nuclei's parameters. A systematic data analysis procedure further suggests that this feature is globally controlled by the dynamic cloud interaction parameter of Maeda and Colonius (J. Fluid Mech., vol. 862, 2019, pp. 1105–1134) that quantifies the mean-field interactions, regardless of initial polydispersity or nonlinearity. The method provides a simplified and comprehensive representation of complex bubble dynamics as well as a new path to reduced-order modelling of cavitation and nucleation.

Through a combination of laboratory experiments and theoretical models, we investigate the interaction of a mean upwelling through a closed basin with a vertical buoyancy flux. The fluid is mixed by a horizontally oscillating rake, which either traverses the whole basin or which oscillates just near one vertical boundary. We first review the steady state and demonstrate that, in both mixing regimes, the vertical density profile across the basin is controlled by the steady-state balance between the upward advective and diffusive fluxes of salinity as described by the classical model introduced by Munk (Deep-Sea Res., vol. 13, issue 4, 1966, pp. 707–730). However, with boundary mixing, we show that both the upwelling and the buoyancy transport are localised to the mixing zone near the boundary, and the interior fluid is stagnant. We then develop a model to describe the transient evolution of the system if there is either a discrete increase or gradual decrease to the buoyancy flux. In the boundary mixing case, the change in the buoyancy flux at the lower boundary leads to a change in the buoyancy of the fluid in the boundary mixing region, and this induces a transient, buoyancy-driven flow in the boundary region in addition to the steady upwelling. In turn, an equal and opposite vertical flow develops in the interior, and this leads to a change in the density stratification of the interior fluid as the system adjusts to a new equilibrium. However, in our experiments, there is no vertical mixing in the interior and interior fluid may upwell or downwell dependent on the change to the buoyancy forcing. We discuss the implications of our results for the transport and mixing in the deep ocean, and the associated interpretation of field experiments.

This experimental work explores the flow field around a three-dimensional expansion–compression geometry on a slender cone at Mach 8 using high-frequency pressure sensors, high-frame-rate schlieren, temperature-sensitive paint, shear-stress measurements and oil-flow visualizations. The $7^\circ$ cone geometry has a hyperbolic slice which acts as an expansion corner and suppresses the disturbances present in the upstream boundary layer. Downstream of the cone-slice corner, high-frequency boundary-layer disturbances attenuate in all cases. Under laminar conditions, second-mode instabilities from the cone diminish and lower-frequency second-mode waves develop on the slice at a frequency commensurate with the increased boundary layer thickness. For fully turbulent cases, the boundary layer over the slice shows evidence of a two-layered nature with a turbulent outer region and a near-wall region with strong attenuation of high-frequency disturbances and reappearance of lower-frequency instability waves. When a downstream compression ramp is added to the slice, the expanded boundary layer shows enhanced susceptibility to separation such that separation is observed at a $10^{\circ }$ deflection, which is smaller than expected for turbulent conditions. For a $30^{\circ }$ ramp, boundary-layer separation occurs further upstream where the heat flux contours show a decrease in heating that is characteristic of a transitional separation. These results demonstrate the effect of relaminarization caused by an upstream expansion on a subsequent shock-wave/boundary-layer interaction.

The influence of the boundary layer (BL) thickness on the tonal noise generated by subsonic impinging jets is investigated. For that, initially laminar jets at Mach numbers $0.6$ and $0.9$ with BL thicknesses $0.05r_0$, $0.1r_0$ and $0.2r_0$, where $r_0$ is the pipe-nozzle radius, impinging on a plate at $6r_0$ from the nozzle, are simulated. For Mach number $0.9$, acoustic tones due to feedback loops of axisymmetric nature between the nozzle and the plate are generated at frequencies that do not vary with the BL thickness. The two dominant tones are, however, $17$ and $26\ \mathrm {dB}$ stronger for the thickest BL compared with the thinnest one. For Mach number $0.6$, for the thinnest BL no acoustic peaks appear, as observed in the experiments of the literature, but narrow peaks resulting from axisymmetric feedback loops emerge for thicker BLs. Therefore, low subsonic impinging jets can be resonant for specific nozzle-exit conditions. The increase in tone amplitude for Mach number $0.9$, and the establishment of feedback loops for Mach number $0.6$ with increasing BL thickness, are found to result from two changes in the jet flow. The first change is that the shear-layer laminar–turbulent transition occurs farther downstream for a thicker BL, leading to a greater predominance of the axisymmetric aerodynamic fluctuations near the plate. The second change is that the amplification of the flow fluctuations between the nozzle and the plate at the tone frequencies is stronger for thicker BLs.

We extend the resolvent framework to two-phase flows with low-inertia particles. The particle velocities are modelled using the equilibrium Eulerian model. We analyse the turbulent flow in a vertical pipe with Reynolds number of $5300$ (based on diameter and bulk velocity), for Stokes numbers $St^+=0-1$, Froude numbers $Fr_z=-4,-0.4,0.4,4$ and $1/Fr_z = 0$ (gravity omitted). The governing equations are written in input–output form and a singular value decomposition is performed on the resolvent operator. As for single-phase flows, the operator is low rank around the critical layer, and the true response can be approximated using one singular vector. Even with a crude forcing model, the formulation can predict physical phenomena observed in Lagrangian simulations, such as particle clustering and gravitational effects. Increasing the Stokes number shifts the predicted concentration spectra to lower wavelengths; this shift also appears in the direct numerical simulation spectra and is due to particle clustering. When gravity is present, there are two critical layers, one for the concentration field, and one for the velocity field. For upward flow, the peak of concentration fluctuations shifts closer to the wall, in agreement with the literature. We explain this with the aid of the different locations of the two critical layers. Finally, the model correctly predicts the interaction of near-wall vortices with particle clusters. Overall, the resolvent operator provides a useful framework to explain and interpret many features observed in Lagrangian simulations. The application of the resolvent framework to higher $St^+$ flows in combination with Lagrangian simulations is also discussed.

This paper extends the work in Tee et al. (Intl J. Multiphase Flow, vol. 133, 2020, 103462) to investigate the effect of turbulent fluid motions on the translation and rotation of lifting and wall-interacting spheres in boundary layers. Each sphere was released from rest in smooth-wall boundary layers with $Re_\tau =670$ and 1300 ($d^+=56$ and 116, respectively) and allowed to propagate with the incoming fluid. Sphere and surrounding fluid motions were tracked simultaneously via three-dimensional particle tracking velocimetry and stereoscopic particle image velocimetry in streamwise–spanwise planes. Two-point correlations of sphere and fluid streamwise velocities yielded long positive regions associated with long fast- and slow-moving zones that approach and move over the spheres. The related spanwise correlations were shorter due to the shorter coherence length of spanwise fluid structures. In general, spheres lag the surrounding fluid. The less-dense lifting sphere had smaller particle Reynolds numbers varying from near zero up to 300. Its lift-offs coincided with oncoming fast-moving zones and fluid upwash. Wall friction initially retarded the acceleration of the denser sphere. Later, fluid torque associated with approaching high-velocity regions initiated forward rotation. The rotation, which was long-lived, induced sufficient Magnus lift to initiate repeated small lift-offs, reduce wall friction, and accelerate the sphere to higher sustained velocity. Particle Reynolds numbers remained above 200, and vortex shedding was omnipresent such that the spheres clearly altered the fluid motion. Spanwise fluid shear occasionally initiated wall-normal sphere rotation and relatively long-lasting Magnus side lift. Hence the finite sphere size contributed to multiple dynamical effects not present in point-particle models.

This study experimentally investigates the wake structure of a porous square cylinder in terms of permeability over two decades of $Da$ (i.e. $2.4 \times 10^{-5} < Da < 2.9 \times 10^{-3}$). The porous cylinder, featuring a simple cubic lattice structure, was fabricated using an additive manufacturing technique. This unique method, combined with a periodic and scalable lattice structure, effectively isolates permeability from porosity, making it suitable for an in-depth parametric study. The key parameter, permeability, was directly estimated by measuring the pressure drop and superficial velocity for each porous case in an open-loop pipe flow system. The downstream flow fields were obtained using standard planar particle image velocimetry measurements in an open-loop wind tunnel. Based on the experimental data, structural modifications in the near wake were examined in relation to permeability, leading to the identification of four distinct flow regimes depending on $Da$. Additionally, the downstream flow adjustment length ($L_i$) was assessed by introducing a permeability-based source term into the momentum equation, facilitating the development of an analytical model for $L_i$. The present experimental data support this analytical model, and our results further confirmed that $L_i$ plays a crucial role as a characteristic length scale in the near wake.

We study the motion and deformation of a single bubble rising inside a cylindrical container filled with a viscoelastic material. We solve numerically the mass and momentum balances along with the constitutive equation for the viscoelastic stresses, without assuming axial symmetry to allow the growth of three-dimensional disturbances. Hence, we may predict the emergence of the notorious knife-edge shape of the bubble, which is a result of a purely elastic instability triggered in the locality of the trailing edge. Our results compare well with existing experiments. We visualize, for the first time to the best of our knowledge, the flow kinematics and dynamics that arise downstream of the bubble. We propose two quantities, one kinematical and one geometrical, for the determination of the onset of the instability. We demonstrate that extension-rate thinning in the constitutive law is necessary for the emergence of the instability. Moreover, our results indicate that increasing (a) the deformability of the bubble and (b) the initial extension rate hardening of the viscoelastic material, prior to thinning, triggers the instability earlier. These novel findings help us formulate and propose a mechanism that controls the onset of the instability and explain why the knife-edge shape is not encountered as frequently.

Combustion instabilities in annular systems raise fundamental issues that are also of practical importance to aircraft engines and ground-based gas turbine combustors. Recent studies indicate that the injector plays a significant role in the stability of combustors by defining the flame dynamical response and setting the inlet impedance of the system. The present investigation examines the effects of combinations of injectors of two different types ($U$ and $S$) on thermoacoustic instabilities in a laboratory-scale annular combustor and compares different circumferential staging strategies. The combustor operates in a stable fashion when all injection units belong to the $S$-family, but exhibits large amplitude pressure oscillations when all these units are of the $U$-type. When the system comprises a mix of $U$- and $S$-injectors, it is possible to determine the number of $S$-injectors leading to stable operation. For a fixed proportion of $U$- and $S$-injectors, some arrangements give rise to stable operation while others do not. Results also show that introducing symmetry-breaking elements affects the system's modal dynamics. These experimental observations are interpreted in an acoustic energy balance framework used to derive an expression for the growth rate as a function of the describing functions of the flames formed by the different injectors and their respective azimuthal locations. Growth rates are determined for the different configurations and used to explain the various observations, estimate the system damping rate and predict the location of the nodal line when the standing mode prevails.

First-order sensitivities and adjoint analysis are used widely to control the linear stability of unstable flows. Second-order sensitivities have recently helped to increase accuracy. In this paper, a method is presented to calculate arbitrary high-order sensitivities based on Taylor expansions of the incompressible base flow and its eigenproblem around a scalar parameter. For the incompressible Navier–Stokes equations, general expressions for the sensitivities are derived, into which parameter-specific information can be inserted. The computational costs are low since, for all orders, a linear equation system has to be solved, of which the left-hand-side matrix stays constant and thus its preconditioning can be exploited. Two flow scenarios are examined. First, the cylinder flow equations are expanded around the inverse of the Reynolds number, enabling the prediction of the two-dimensional cylinder base flow and its leading eigenvalue as a function of the Reynolds number. This approach computes accurately the base flow and eigenvalue even in the unstable regime, providing, when executed subsequently, a mean to calculate unstable base flows. This case gives a clear introduction into the method and allows us to discuss its constraints regarding convergence behaviour. Second, a small control cylinder is introduced into the domain of the cylinder flow for stabilization. Higher-order sensitivity maps are calculated by modelling the small cylinder with a steady forcing. These maps help to identify stabilizing areas of the flow field for Reynolds numbers within the laminar vortex shedding regime, with the required number of orders increasing as the Reynolds number rises. The results obtained through the proposed method align well with numerically calculated eigenvalues that incorporate the cylinder directly into the grid.

Oscillatory flows have become an indispensable tool in microfluidics, inducing inertial effects for displacing and manipulating fluid-borne objects in a reliable, controllable and label-free fashion. However, the quantitative description of such effects has been confined to limit cases and specialized scenarios. Here we develop an analytical formalism yielding the equation of motion of density-mismatched spherical particles in oscillatory background flows, generalizing previous work. Inertial force terms are systematically derived from the geometry of the flow field together with analytically known Stokes number dependences. Supported by independent, first-principles direct numerical simulations, we find that these forces are important even for nearly density-matched objects such as cells or bacteria, enabling their fast displacement and separation. Our formalism thus consistently incorporates particle inertia into the Maxey–Riley equation, and in doing so provides a generalization of Auton's modification to added mass, as well as recovering the description of acoustic radiation forces on particles as a limiting case.

We derive the spin Euler equation for ideal flows by applying the spherical Clebsch mapping. This equation is based on the spin vector, a unit vector field encoding vortex lines, instead of the velocity. The spin Euler equation enables a feasible Lagrangian study of fluid dynamics, as the isosurface of a spin-vector component is a vortex surface and material surface in ideal flows. We establish a non-blowup criterion for the spin Euler equation, suggesting that the Laplacian of the spin vector must diverge if the solution forms a singularity at some finite time. The direct numerical simulations (DNS) of three ideal flows – the vortex knot, the vortex link and the modified Taylor–Green flow – are conducted by solving the spin Euler equation. The evolution of the Lagrangian vortex surface illustrates that the regions with large vorticity are rapidly stretched into spiral sheets. The DNS result exhibits a pronounced double-exponential growth of the maximum norm of Laplacian of the spin vector, showing no evidence of the finite-time singularity formation if the double-exponential growth holds at later times. Moreover, the present criterion with Lagrangian nature appears to be more sensitive than the Beale–Kato–Majda criterion in detecting the flows that are incapable of producing finite-time singularities.

The probability density function (PDF) for the free surface elevation in an irregular sea has an integral formulation when based on the cumulant generating function. To leading order, the result is Gaussian, whereas nonlinear extensions have long been limited to Gram–Charlier series approximations. As shown recently by Fuhrman et al. (J. Fluid Mech., vol. 970, 2023, A38), however, the second-order integral can be represented exactly in closed form. The present work extends this further, enabling determination of this PDF to even higher orders. Towards this end, a new ordinary differential equation (ODE) governing the PDF is first derived. Asymptotic solutions in the limit of large surface elevation are then found, utilizing the method of dominant balance. These provide new analytical forms for the positive tail of the PDF beyond second order. These likewise clarify how high-order cumulants (involving statistical moments such as the kurtosis) govern the tail, which is shown to get heavier with each successive order. The asymptotic solutions are finally utilized to generate boundary conditions, such that the governing ODE may be solved numerically, enabling novel determination of the PDF at third and higher order. Successful comparisons with challenging data sets confirm accuracy. The methodology thus enables the PDF of the surface elevation to be determined numerically, and the asymptotic tail analytically, to any desired order. Results are worked out explicitly up to fifth order. The theoretical probability of extreme surface elevations (typical of rogue waves) may thus be assessed quantitatively for highly nonlinear irregular seas, requiring only relevant statistical quantities as input.

Air-blast atomizers are extensively used for a variety of purposes. Due to its complexity, the atomization mechanism has not been elucidated. In this study, a mechanistic model is proposed to predict the droplet diameter distribution based on the atomization process of a planar liquid film with co-current gas flows, and its validity is examined by comparing the estimated and measured droplet diameters using high-speed image analysis and laser measurement. As a result, using high-speed imaging, we clarified that the bag film rupture is caused not by the turbulence of the gas flow but by the impact of floating droplets on the liquid film of the expanding bag when the film is thin enough. The average thickness of the liquid film at the bag breakup is of the order of micrometres and varies greatly, resulting in a dispersed distribution of droplet diameters. After the film ruptures, the bag film shrinks towards its transversal and vertical rims due to surface tension, forming large-diameter ligaments. During the contraction process of the bag film, tiny droplets of the order of micrometers are formed at the edge of the perforation. Finally, the remaining ligaments with large diameters fragment into large droplets with submillimetre diameters. The good agreement between the measured and predicted droplet diameter distributions validated the mechanistic model.

It is well known that inertia-free shearing flows of a viscoelastic fluid with curved streamlines, such as the torsional flow between a rotating cone and plate or the flow in a Taylor–Couette geometry, can become unstable to a three-dimensional time-dependent instability at conditions exceeding a critical Weissenberg ($Wi$) number. However, the combined effects of fluid elasticity, shear thinning and finite inertia (as quantified by the Reynolds number $Re$) on the onset of elasto-inertial instabilities are not fully understood. Using a set of cone–plate geometries, we experimentally explore the entire $Wi$–$Re$ phase space for a series of nonlinear viscoelastic fluids (with the dependence on shear rate $\dot{\gamma}$ quantified using a shear-thinning parameter $\beta _P(\dot {\gamma })$). We tune $\beta _P(\dot {\gamma })$ by varying the dissolved polymer concentration in solution. This progressively reduces shear thinning but leads to finite inertial effects before the onset of elastic instability, and thus naturally results in elasto-inertial coupling. Time-resolved rheometric measurements and flow visualization experiments allow us to investigate the effects of flow geometry, and document the combined effects of varying $Wi, Re$ and $\beta _P(\dot {\gamma })$ on the emergence of secondary motions at the onset of instability. The resulting critical state diagram quantitatively depicts the competition between the stabilizing effects of shear thinning and the destabilizing effects of inertia. We extend the curved streamline instability criterion of Pakdel & McKinley (Phys. Rev. Lett., vol. 77, no. 12, 1996, p. 2459) for the onset of purely elastic instability in curvilinear geometries by using scaling arguments to incorporate shear thinning and finite inertial effects. The augmented condition facilitates predictions of the onset of instability over a broader range of flow conditions, and thus bridges the gap between purely elastic and elasto-inertial curved streamline instabilities.

We use Bayesian inference, accelerated by adjoint methods, to construct a quantitatively accurate model of the thermoacoustic behaviour of a conical flame in a duct. We first perform a series of automated experiments on a ducted flame rig. Next, we propose several candidate models of the rig's components and assimilate data into each model to find the most probable parameters for that model. We rank the candidate models based on their marginal likelihood (evidence) and select the most likely model for each component. We begin this process by rigorously characterizing the acoustics of the cold rig. When the flame is introduced, we propose several candidate models for the fluctuating heat release rate. We find that the most likely flame model considers velocity perturbations in both the burner feed tube and the outer duct, even though studies in the literature typically neglect either one of these. Using the most likely model, we infer the flame transfer functions for 24 flames and quantify their uncertainties. We do this with the flames in situ, using only pressure measurements. We find that the inferred flame transfer functions render the model quantitatively accurate, and, where comparable, broadly consistent with direct measurements from several studies in the literature.