The effect of damping the longest streaks in wall-bounded turbulence is explored using numerical experiments. It is found that long streaks are not required for the self-sustenance of the bursting process, which is relatively little affected by their absence. In particular, there are turbulence states in which the fluctuations of the streamwise velocity have approximately the same length as the bursts, and are thus presumably associated with the bursts themselves, while the burst structure is essentially indistinguishable from flows in which longer velocity fluctuations are present. This suggests that the long streaks found in unmodified flows may be by-products, rather than active parts of the process.

Compressible mixing layer instabilities are of importance to a wide range of environmental and industrial applications. Past studies have focused on either ideal-gas or real-fluid thermodynamic regimes of single-species mixing layers. However, mixing layers of binary mixtures at supercritical conditions, commonly encountered in fuel injection systems, introduce additional complexities due to the added compositional degree-of-freedom. Moreover, the effect of strong variations in thermodynamic response functions across the Widom line on the binary mixing layer stability remains poorly understood. Thus, the objective of this study is to examine the coupling between the hydrodynamic instability and the real-fluid thermodynamics across the Widom line and its effects on the overall binary mixing layer dynamics. To this end, we develop a linear stability analysis of the full binary-species compressible transport equations coupled with the PC-SAFT equation of state. Analysis shows the existence of a novel instability mechanism that arises from juxtapositioning of the Widom-line transition and the hydrodynamic inflection point. This novel thermodynamically induced instability mechanism has the net effect of destabilizing the binary mixing layer at lowering supercritical conditions towards the critical pressure point. This is in contrast to previous stability analyses of supercritical single-species mixing layers, where increasing pressure destabilizes the flow due to its effect on reducing the density stratification. The discovered thermodynamically induced instability mechanism of binary mixing flows highlights the need for an extension of classical instability criteria to incorporate the effect of strong variations in the thermodynamic response functions across the Widom line on mixing layer instability.

Viscous streaming is an efficient mechanism to exploit inertia at the microscale for flow control. While streaming from rigid features has been thoroughly investigated, when body compliance is involved, as in biological settings, little is known. Here, we investigate body elasticity effects on streaming in the minimal case of an immersed soft cylinder. Our study reveals an additional streaming process, available even in Stokes flows. Paving the way for advanced forms of flow manipulation, we illustrate how gained insights may translate to complex geometries beyond circular cylinders.

Motivated by the variation of local shear produced by internal waves in the ocean, we use direct numerical simulations to investigate the effect of a time-dependent shear forcing on the evolution and mixing of turbulence produced by Kelvin–Helmholtz instability (KHI) at high Reynolds number. The forcing is implemented using a tilting coordinate system which causes the background shear to accelerate and decelerate periodically. We demonstrate that, with suitable timing between development of instability and the shear oscillation cycle, turbulence produced by KHI with a decelerating shear mixes in a distinctly different way from the flow with constant background shear, specifically with the energy for turbulent motions extracted from alternative sources. As a result, the total amount of mixing as measured by the change in background potential energy can in fact be significantly larger for flows in which the shear is decelerated, despite the fact that the total kinetic energy in the flow is significantly smaller. The mixing has characteristics more in common with convectively driven rather than shear-driven flows, supporting the argument for an underlying change in the mechanisms triggering the turbulence.

A simple model for the motion of shape-changing swimmers in Poiseuille flow was recently proposed and numerically explored by Omori et al. (J. Fluid Mech., vol. 930, 2022, A30). These explorations hinted that a small number of interacting mechanics can drive long-time behaviours in this model, cast in the context of the well-studied alga Chlamydomonas and its rheotactic behaviours in such flows. Here, we explore this model analytically via a multiple-scale asymptotic analysis, seeking to formally identify the causal factors that shape the behaviour of these swimmers in Poiseuille flow. By capturing the evolution of a Hamiltonian-like quantity, we reveal the origins of the long-term drift in a single swimmer-dependent constant, whose sign determines the eventual behaviour of the swimmer. This constant captures the nonlinear interaction between the oscillatory speed and effective hydrodynamic shape of deforming swimmers, driving drift either towards or away from rheotaxis.

Drop-on-demand printing applications involve a drop connected to a fluid reservoir between which volume can be exchanged, a situation that can be idealized as a sessile drop with prescribed volume flux across the drop/reservoir boundary. Here we compute the frequency spectrum for these pressure disturbances, as it depends upon the static contact-angle $\alpha$ (CA) and an empirical constant $\chi$ relating the reservoir pressure to volume exchanged, for either (i) pinned or (ii) free contact-lines. Mode shapes are characterized by the mode number pair $[k,\ell ]$ with property $k+\ell =\mathbb {Z}^{+}_{odd}$ that can be associated with the symmetry properties of the Rayleigh drop modes for the free sphere. We report instabilities to the axisymmetric $[1,0]$ and non-axisymmetric rocking $[2,1]$ modes that are related to centre-of-mass motions, and show how the spectral degeneracy of the Rayleigh drop modes breaks with the model parameters.

The inner–outer interaction model (Marusic et al., Science, vol. 329, 2010, pp. 193–196) and the attached-eddy model (Townsend, Cambridge University Press, 1976) are two fundamental models describing the multiscale turbulence interactions and the organization of energy-containing motions in the logarithmic region of high-Reynolds-number wall-bounded turbulence, respectively. In this paper, by coupling the additive description with the attached-eddy model, the generation process of streamwise wall-shear fluctuations, resulting from wall-attached eddies, is portrayed. Then, by resorting to the inner–outer interaction model, the streamwise wall-shear stress fluctuations generated by attached eddies in a turbulent channel flow are isolated. Direct comparison between the statistics from these two models demonstrates that they are consistent with and complement each other. Meanwhile, we further show that the superpositions of attached eddies follow an additive process strictly by verifying the validity of the strong and extended self-similarity. Moreover, we propose a Gaussian model to characterize the instantaneous distribution of streamwise wall-shear stress, resulting from the attached-eddy superpositions. These findings are important for developing an advanced reduced-order wall model.

We report non-monotonic wettability effects on displacement efficiency in heterogeneous porous structures at the post-breakthrough stage, in contrast to the monotonic ones in homogeneous porous structures. Experiments on designed microfluidic chips show that there exists a critical wettability to attain the highest efficiency of displacement in the porous matrix structure combined with a preferential flow pathway, while a stronger wettability of the displacing fluid leads to a higher displacement efficiency on the same matrix structure only. The porous structure with or without a preferential flow pathway results in totally different topological characteristics of phase distribution during displacement. Pore-scale mechanisms are identified to elucidate the formation of this non-monotonic wettability rule: cooperative pore filling under weakly water-wet conditions yields the best displacement; corner flow under strongly water-wet conditions and Haines events under strongly oil-wet conditions decrease the displacement efficiency. The pore-scale findings may provide unique insights into the joint effects of both wettability and flow heterogeneity on fluid displacement in porous media.

We revisit the theoretical solutions of the problem of the reflected flow from a vertical obstacle encountered by a current generated by dambreak of a Cartesian long reservoir over a dry bottom and with tailwaters. Previous investigations derived accurate solutions by quite complex manipulation of the balances along the characteristics, supported by numerical solutions, of the shallow-water equations (a significant extension is the recent paper Hogg & Skevington, Q. J. Mech. Appl. Maths, vol. 74, no. 4, 2021, pp. 441–465). Here we developed a simplified model, based on the major assumption that the fluid between the obstacle (wall) and the reflected jump is stagnant. This allows the solution of the problem by a straightforward numerical integration of one initial-value ordinary differential equation. The model provides the position, height and speed of the jump as functions of time. For long times a simple analytical approximation is also available. The model points out clear-cut effects of the presence of the tailwaters. The model has been validated by comparisons with exact solutions of the shallow-water equations of the recent study of Hogg & Skevington (Q. J. Mech. Appl. Maths, vol. 74, no. 4, 2021, pp. 441–465) (obtained by significantly more complex hodograph-plane methods). In all the tested cases the agreement is good, for long periods of time. This model provides reliable insights and fast quantitative predictions recommended for use in research and engineering problems where a fair approximation is sufficient for the application.

Dynamic stall on airfoils negatively impacts their aerodynamic performance and can lead to structural damage. Accurate prediction and modelling of the dynamic stall loads are crucial for a more robust design of wings and blades that operate under unsteady conditions susceptible to dynamic stall and for widening the range of operation of these lifting surfaces. Many dynamic stall models rely on empirical parameters that need to be obtained from experimental or numerical data which limits their generalisability. Here, we introduce physically derived time scales to replace the empirical parameters in the Goman–Khrabrov dynamic stall model. The physics-based time constants correspond to the dynamic stall delay and the decay of post-stall load fluctuations. The dynamic stall delay is largely independent of the type of motion, the Reynolds number and the airfoil geometry, and is described as a function of a normalised instantaneous pitch rate. The post-stall decay is independent of the motion kinematics and is related to the Strouhal number of the post-stall vortex shedding. The general validity of our physics-based time constants is demonstrated using three sets of experimental dynamic stall data covering various airfoil profiles, Reynolds numbers varying from 75 000 to 1 000 000, and sinusoidal and ramp-up pitching motions. The use of physics-based time constants generalises the Goman–Khrabrov dynamic stall model and extends its range of application.

We study the kinematics, dynamics and flow fields generated by an oscillating, compliant membrane hydrofoil extracting energy from a uniform water stream at a chord-based Reynolds number $Re \approx 3 \times 10^4$. Hydrodynamic forces during the foil's motion cause the membrane to dynamically morph its shape, effectively increasing the camber during the oscillation cycle. The membrane's deflection is modelled using the Young–Laplace equation, with pressure term approximated from thin-airfoil theory. Simultaneous tracking of the membrane deformation and the surrounding flow field using laser profiling and particle image velocimetry, respectively, reveals the role of dynamic cambering in stabilizing the leading-edge vortices on the membrane. In this regime of operation, we obtain up to 160 % higher power extraction when compared to a rigid, symmetric hydrofoil. The present work provides a demonstration of how passive compliance of soft materials interacting with fluids may be exploited in tidal and fluvial energy extraction.

The disturbance flow field for a Mach $6$ flat-plate boundary layer flow with a wall-to-free-stream temperature ratio of $0.5$ is studied using direct numerical simulation (DNS), linear stability theory (LST) and biorthogonal decomposition. In the second-mode instability region, the DNS flow field can be reconstructed using a single LST mode, namely $F^+_1$. However, when the supersonic mode emerges, none of the discrete modes nor the continuous branches alone can precisely match the DNS data. A superposition of the pair of discrete supersonic modes $F^\pm _1$ and the slow acoustic continuous spectrum is required to reproduce the disturbance amplitude distributions and the amplification rates observed in the DNS. Another finding is that the supersonic mode $F^-_1$, which should be decaying following LST, in actuality is amplified based on the projected DNS data.

The rise of a single bubble confined between two vertical plates is investigated over a wide range of Reynolds numbers. In particular, we focus on the evolution of the bubble speed, aspect ratio and drag coefficient during the transition from the viscous to the inertial regime. For sufficiently large bubbles, a simple model based on power balance captures the transition for the bubble velocity and matches all the experimental data despite strong time variations of bubble aspect ratio at large Reynolds numbers. Surprisingly, bubbles in the viscous regime systematically exhibit an ellipse elongated along its direction of motion while bubbles in the inertia-dominated regime are always flattened perpendicularly to it.

Flow control for turbulent skin-friction drag reduction is applied to a transonic airfoil to improve its aerodynamic performance. The study is based on direct numerical simulations (with up to 1.8 billion cells) of the compressible turbulent flow around a supercritical airfoil, at Reynolds and Mach numbers of $Re_\infty = 3 \times 10^{5}$ and $M_\infty =0.7$. Control via spanwise forcing is applied over a fraction of the suction side of the airfoil. Besides locally reducing friction, the control modifies the shock wave and significantly improves the aerodynamic efficiency of the airfoil by increasing lift and decreasing drag. Hence, the airfoil can achieve the required lift at a lower angle of attack and with a lower drag. Estimates at the aircraft level indicate that substantial savings are possible; when control is active, its energy cost becomes negligible thanks to the small application area. We suggest that skin-friction drag reduction should be considered not only as a goal, but also as a tool to improve the global aerodynamics of complex flows.

We report a new interferometric method to study the interactions between two gas bubbles undergoing small-amplitude oscillations in a liquid, based on the extension of a previously developed one-bubble set-up. Nanometric oscillations of millimetre-sized supported bubbles are excited acoustically; the response of each bubble is recorded interferometrically, as a function of the mutual distance (from quasi-contact to greater than the bubbles radii). The interferometric nature of the technique and the resonant nature of the vibration modes enable the accurate measurement of the amplitude (with sub-nanometric sensitivity), frequency and mutual phase of oscillation, whose variations over the bubble–bubble distance range allow the interactions to be probed. The bubbles oscillate at the same frequencies, exhibiting a low-frequency, in-phase and a high-frequency, out-of-phase resonance peak, whose separation is a function of distance, in good agreement with the theory for free interacting bubbles. The technique, here demonstrated for the volume modes of air bubbles in water, can be extended to other gas–liquid and liquid–liquid interfaces, bare or adsorbate-covered, as well as to shape oscillations.

We present a decomposition of the streamwise fluid force for in-line vortex-induced vibration (VIV) to provide insight into how the wake drag acts as a driving force in fluid–structure interaction. This force decomposition is an extension of that proposed in the recent work of Konstantinidis et al. (J. Fluid Mech., vol. 907, 2021, p. A34), and is applied to and validated by our experiments examining a circular cylinder freely vibrating in line with the free stream. It is revealed from the decomposition and linear analysis that two regimes of significant vibration are in phase synchronisation, while they are separated by a desynchronised regime marked by competition between non-stationary frequency responses of the cylinder vibration and the vortex shedding. Of interest, such a near-resonance desynchronisation regime is not seen in the transverse vibration case.

We examine the response of an inertialess viscoelastic channel flow to localized perturbations. We thus performed an experiment in which we perturbed the flow using a localized velocity pulse and probed the perturbed fluid packet downstream from the perturbation location. While for low Weissenberg numbers the perturbed fluid reaches the measurement location as a single velocity pulse, for sufficiently high Weissenberg numbers and perturbation strengths, a random number of pulses arrive at the measurement location. The average number of pulses observed increases with the Weissenberg number. This observation constitutes a transition to a novel elastic pulse-splitting regime. Our results suggest a possible new direction for studying the elastic instability of viscoelastic channel flows at high elasticity through the growth of localized perturbations.

The propulsion of a pitching flexible plate in a uniform flow is investigated numerically. The effects of bending stiffness ($K$), pitching amplitude ($A_L$) and frequency ($St$) on the wake patterns, thrust generations and propulsive performances of the fluid–plate system are analysed. Four typical wake patterns, i.e. von Kármán, reversed von Kármán, deflected and chaotic wakes, emerge from various kinematics, and the $St-A_L$ wake maps are given for various $K$. The drag-to-thrust transitions (DTT) and the wake transitions (WT) between the von Kármán and reversed von Kármán wakes are examined. Results indicate that the WT and DTT boundaries can be scaled by the chord-averaged distance of travel, $\mathcal {L}$, which leads to $\mathcal {L}\times St \approx 1$ and $\mathcal {L}\times St \approx 1.2$, respectively. Further, the resonance mechanism for the performance enhancement is revealed and confirmed in a wide range of parameters. The dimensionless average speed of plate, $\mathcal {U^*}\left (=\mathcal {L}\times St\right )$, is adopted merely to characterize the propulsive performances. For the first time, the $\mathcal {U^*}$-based scaling laws for the thrust and power are revealed in pitching rigid and flexible plates for various $A_L$ and $St$. This study may deepen our understanding of biological swimming and flying, and provide a guide for bionic design.

Wind tunnel experiments were performed to investigate turbulent flow over an array of heterogeneous roughness elements using stereoscopic particle image velocimetry. Nine streamwise planes, covering one periodic cell of a multi-scale roughness element that is arranged in a staggered pattern, are combined to quantify mean flow features and Reynolds stresses. Dispersive stresses, arising from spatial variations in the temporally averaged mean velocity, are also presented. The results highlight that the roughness elements create a large deficit pathway along the surfaces. Outer scaling of the time-averaged streamwise velocity presents features which are nearly independent of the roughness element type, with parameters of the flow revealing values near those observed in smooth wall boundary layers, such as the wake strength parameter, opposing the earlier work containing aligned patterns. The strength of the secondary motion is most accentuated at the ridges of the roughness, showing that the formation of structures is sensitive to the location investigated in the spanwise direction.

The non-dimensional energy of starting vortex rings typically converges to values around 0.33 when they are created by a piston-cylinder or a bluff body translating at a constant speed. To explore the limits of the universality of this value and to analyse the variations that occur outside of those limits, we present an alternative approach to the slug-flow model to predict the non-dimensional energy of a vortex ring. Our approach is based on the self-similar vortex sheet roll-up described by Pullin (J. Fluid Mech., vol. 88, 1978, pp. 401–430). We derive the vorticity distribution for the vortex core resulting from a spiralling shear layer roll-up and compute the associated non-dimensional energy. To demonstrate the validity of our model for vortex rings generated through circular nozzles and in the wake of disks, we consider different velocity profiles of the vortex generator that follow a power law with a variable time exponent $m$. Higher values of $m$ indicate a more uniform vorticity distribution. For a constant velocity ($m=0$), our model yields a non-dimensional energy of ${{E}^{{*}}}=0.33$. For a constant acceleration ($m=1$), we find ${{E}^{{*}}}=0.19$. The limiting value $m \rightarrow \infty$ corresponds to a uniform vorticity distribution and leads to ${{E}^{{*}}}=0.16$, which is close to values found in the literature for Hill's spherical vortex. The radial diffusion of the vorticity within the vortex core results in the decrease of the non-dimensional energy. For a constant velocity, we obtain realistic vorticity distributions by radially diffusing the vorticity distribution of the Pullin spiral and predict a decrease of the non-dimensional energy from 0.33 to 0.28, in accordance with experimental results. Our proposed model offers a practical alternative to the existing slug flow model to predict the minimum non-dimensional energy of a vortex ring. The model is applicable to piston-generated and wake vortex rings and only requires the kinematics of the vortex generator as input.