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.

A complete understanding of jet dynamics is greatly enabled by the accurate separation of acoustically efficient wavepackets from their higher-energy convecting turbulent counterparts. Momentum potential theory (MPT) has proven highly effective in filtering the desired acoustic component irrespective of operating conditions or nozzle complexity. However, MPT is a data-intensive method predicated on the knowledge of fluctuation quantities in the entire flow field; as such, it has to date been applied only to numerically obtained data. This work develops an approach to extend its application to extract coherent wavepacket data from high-speed schlieren images. Pixel intensities from the schlieren image are mapped to a scaled surrogate for the line-of-sight integrated density gradient. The linear relation between the irrotational scalar MPT potential and time derivatives of density fluctuations is then exploited to perform the filtering. The effectiveness of the procedure is demonstrated using experimental as well as numerical schlieren images representing a wide range of imperfectly expanded free and impinging jet configurations. When combined with spectral proper orthogonal decomposition (SPOD), the method yields modes that accurately capture (i) the Mach wave radiation from a military-style jet, (ii) the mode shapes of the feedback tones in an impinging jet and (iii) the screech signature in twin rectangular jets, without recourse to user adjusted parameters. This technique can greatly enhance the use of high-speed diagnostics for real-time monitoring of the near-field acoustic content for potential feedback control. Additionally, the general nature of the approach allows a straightforward application to other flows, such as cavity and airfoil flow-acoustic interactions.

In many geophysical and astrophysical flows, suppression of fluctuations along one direction of the flow drives a quasi-two-dimensional upscale flux of kinetic energy, leading to the formation of strong vortex condensates at the largest scales. Recent studies have shown that the transition towards this condensate state is hysteretic, giving rise to a limited bistable range in which both the condensate state as well as the regular three-dimensional state can exist at the same parameter values. In this work, we use direct numerical simulations of thin-layer flow to investigate whether this bistable range survives as the domain size and turbulence intensity are increased. By studying the time scales at which rare transitions occur from one state into the other, we find that the bistable range grows as the box size and/or Reynolds number $Re$ are increased, showing that the bistability is neither a finite-size nor a finite-$Re$ effect. We furthermore predict a cross-over from a bimodal regime at low box size, low $Re$ to a regime of pure hysteresis at high box size, high $Re$, in which any transition from one state to the other is prohibited at any finite time scale.

The interplay between convective, rotational and magnetic forces defines the dynamics within the electrically conducting regions of planets and stars. Yet their triadic effects are separated from one another in most studies, arguably due to the richness of each subset. In a single laboratory experiment, we apply a fixed heat flux, two different magnetic field strengths and one rotation rate, allowing us to chart a continuous path through Rayleigh–Bénard convection (RBC), two regimes of magnetoconvection, rotating convection and two regimes of rotating magnetoconvection, before finishing back at RBC. Dynamically rapid transitions are determined to exist between jump rope vortex states, thermoelectrically driven magnetoprecessional modes, mixed wall- and oscillatory-mode rotating convection and a novel magnetostrophic wall mode. Thus, our laboratory ‘pub crawl’ provides a coherent intercomparison of the broadly varying responses arising as a function of the magnetorotational forces imposed on a liquid-metal convection system.

The effect of centrifugal force in turbulent rotating Rayleigh–Bénard convection (RRBC) is studied experimentally in an aspect-ratio $\varGamma =1$ cylindrical convection cell and in the ranges of the Froude number $0.004\leq Fr \leq 0.363$ and the Rayleigh number $2.8\times 10^8 \leq Ra \leq 9.5\times 10^9$, and with the Prandtl number fixed at $Pr=4.34$. We use the bulk temperature anomaly to determine the onset Froude number $Fr_c$, beyond which the centrifugal effects cannot be regarded as insignificant. It is found that $Fr_c$ depends on $Ra$ as $Fr_c\sim Ra^{0.53}$, which may be understood qualitatively by the idea of local force balance. For $Fr>Fr_c$, the centrifugal effect is more pronounced for smaller $Ra$, which is also found for larger constant $1/Ro$. This implies that the response of the system to the centrifugal force depends on the flow states, which, in RRBC, is mainly determined by the competition between the buoyancy and Coriolis forces. Detailed analysis of the sidewall temperature signal shows results consistent with those obtained from the bulk temperature. Based on the above results, we propose a different division of the $1/Ro$–$Fr$ phase space than previously suggested. For the heat transport, the results under fixed $1/Ro$ show well-defined $Nu$–$Ra$ scalings, which can provide a better prediction for the heat transport when extrapolating to the unexplored regions in the phase space.

We develop direct inversion methods for inferring the rheology of a fluid from observations of its shallow flow. First, the evolution equation for the free-surface flow of an inertia-less current with general constitutive law is derived. The relationship between the volume flux of fluid and the basal stress, $\tau _b$, is encapsulated by a single function $F(\tau _b)$, which depends only on the constitutive law. The inversion method consists of (i) determining the flux and basal stress from the free-surface evolution, (ii) comparing the flux with the basal stress to constrain $F$ and (iii) inferring the constitutive law from $F$. Examples are presented for both steady and transient free-surface flows demonstrating that a wide range of constitutive laws can be directly obtained. For flows in which the free-surface velocity is known, we derive a different method, which circumvents the need to calculate the flux.

Recently, an unusual scaling law has been observed in circular hydraulic jumps and has been attributed to a supposed missing term in the local energy balance of the flow (Bhagat et al., J. Fluid Mech., vol. 851, 2018, R5). In this paper, we show that – though the experimental observation is valuable and interesting – this interpretation is presumably not the right one. When transposed to the case of an axial sheet formed by two impinging liquid jets, the assumed principle leads in fact to a velocity distribution in contradiction with the present knowledge for this kind of flow. We show here how to correct this approach by maintaining consistency with surface tension thermodynamics: for Savart–Taylor sheets, when adequately corrected, we recover the well-known $1/r$ liquid thickness with a constant and uniform velocity dictated by Bernoulli's principle. In the case of circular hydraulic jumps, we propose here a simple approach based on Watson's description of the flow in the central region (Watson, J. Fluid Mech., vol. 20, 1964, pp. 481–499), combined with appropriate boundary conditions on the circular front formed. Depending on the specific condition, we find in turn the new scaling by Bhagat et al. (2018) and the more conventional scaling law found long ago by Bohr et al. (J. Fluid Mech., vol. 254, 1993, pp. 635–648). We clarify here a few situations in which one should hold rather than the other, hoping to reconcile the observations of Bhagat et al. with the present knowledge of circular hydraulic jump modelling. However, the question of a possible critical Froude number imposed at the jump exit and dictating logarithmic corrections to scaling remains an open and unsolved question.

Many recent studies have diagnosed opposite-signed diffusion eigenvalues to be a prevalent feature of the transfer tensor for diffusive tracer transport by oceanic mesoscale eddies. This diagnosed tensor, which we refer to as the diffusion tensor, therefore accounts for tracer filamentation effects. The preferential orientation of this filamentation is quantified by the principal axis of the diffusion tensor, namely the diffusion axis. Parameterisations of eddy diffusion commonly invoke a diffusion tensor, typically one with non-negative eigenvalues to avoid numerical issues. Motivated by the need to parameterise tracer filamentation, in this study we examine diffusion of a Gaussian tracer patch with imposed opposite-signed diffusion eigenvalues, and in particular we focus on the time scale for the onset of instability. For a fixed diffusion axis, numerical instability is an inevitable consequence of persistent up-gradient fluxes associated with the negative eigenvalue. For typical oceanic scales and diffusion magnitudes, this time scale is of the order of $100$ days, but is shorter for larger negative eigenvalues or for finer grid resolutions. We show that imposing a time-dependent diffusion axis can lead to simulations with no onset of instability after 100 000 days of tracer evolution. Although motivated by oceanographic fluid dynamics, our results have much broader applications since diffusive processes are present in a wide range of fluid flows.

We develop a model to predict the fragmentation limit of drops colliding off-centre. The prediction is excellent over a wide range of liquid properties and it can be used without adjusting any parameter. The so-called stretching separation is attributed to the extension of the merged drop above a critical aspect ratio of 3.25. The evolution of this aspect ratio is influenced by the liquid viscosity and can be interpreted via an energy balance. This approach is then adapted to drop–jet collisions, which we model as consecutive drop–drop collisions. The fragmentation criterion is similar to that observed for drop–drop collisions, while the evolution of the stretched jet aspect ratio is modified to account for the different flow fields and geometry.

The unsteady response of nozzles with steady heat transfer forced by acoustic and/or entropy waves is modelled. The approach is based on the quasi-one-dimensional linearised Euler equations. The equations are cast in terms of three variables, namely the dimensionless mass, stagnation temperature and entropy fluctuations, which are invariants of the system at zero frequency and with no heat transfer. The resulting first-order system of differential equations is then solved using the Magnus expansion method, where the perturbation parameters are the normalised frequency and the volumetric heat transfer. In this work, a measure of the flow non-isentropicity (in this case the steady heat transfer) is used for the first time as an expansion parameter. The solution method was applied to a converging–diverging nozzle with constant heat transfer for both subcritical and supercritical flow cases, showing good agreement with numerical predictions. It was observed that the acoustic and entropy transfer functions of the nozzle strongly depend on the frequency and heat transfer.

A scaling law to predict the onset of the primary Hopf instability in the steady flow past two-dimensional symmetric bluff bodies is proposed. It uses a measure of the spatial extent of the separation bubble as the length scale and the largest reverse-flow speed within it as the velocity scale. The ensuing Reynolds number, evaluated at the onset of the primary Hopf bifurcation, collapses quite nicely for bodies of different shape and aspect ratio even when a small angle of attack perturbs the symmetry; its relative variation is one order of magnitude smaller than that of the usual Reynolds number defined with the free-stream velocity and the cross-stream body size. With the new scaling, it can be roughly assessed whether the steady flow past a two-dimensional bluff body is absolutely and globally unstable to two-dimensional perturbations without a computationally expensive stability analysis: only the inspection of the base flow is required. More importantly, the scaling provides an insight into the flow mechanism that produces the instability.

We present a non-perturbative mesh-free spectral method for the axisymmetric scenario of a radially highly deformed sphere in Stokes flow. Spectra of harmonic and biharmonic Stokes flow modes are used to provide a general algebraic spectral solution to axisymmetric Stokes flow. A solution for the flow field around a sphere with prescribed surface velocity field is presented, and is used to obtain the velocity field around a deformed sphere. The method is demonstrated on two problems: hydrodynamic radii of radially deformed spheres for a range of strength and angular dependence of the deformation, with results in good agreement with the boundary element method; and self-phoretic velocity of spheroids with surface flux of the driving field in a source/sink or source/inert configuration.

Using simulations and a virtual-spring-based approach, we measure the segregation force, $F_{seg},$ in size-bidisperse sphere mixtures over a range of concentrations, particle-size ratios and shear rates to develop a semiempirical model for $F_{seg}$ that extends its applicability from the well-studied non-interacting intruders regime to finite-concentration mixtures where cooperative phenomena occur. The model predicts the concentration below which the single-intruder assumption applies and provides an accurate description of the pressure partitioning between species.

Results from seven direct and large-eddy simulations of gravity currents on slopes ranging from 0.14$^{\circ }$ to 2.86$^{\circ }$ that span from the subcritical to the supercritical regime are studied. By considering a long domain, attention is focused on the near-self-similar state approached by these currents far downstream. In the self-similar limit, the various shape factors, local Richardson number, entrainment coefficient, velocity scale and basal drag coefficient reach a constant value, while the current height, volume and momentum fluxes continue to increase linearly. Their dependence on slope is presented.

Multiphase flow in porous materials is conventionally described by an empirical extension to Darcy's law, which assumes that the pressure gradient is proportional to the flow rate. Through a series of two-phase flow experiments, we demonstrate that even when capillary forces are dominant at the pore scale, there is a nonlinear intermittent flow regime with a power-law dependence between pressure gradient and flow rate. Energy balance is used to predict accurately the start of the intermittent regime in hydrophobic porous media. The pore-scale explanation of the behaviour based on the periodic filling of critical flow pathways is confirmed through 3D micron-resolution X-ray imaging.

The thermal resistance associated with the interface between a solid and a liquid is analysed from an atomistic point of view. Partial evaluation of the associated Green–Kubo integral elucidates the various factors governing heat transport across the interface and leads to a quantitative model for the thermal resistance in terms of atomistic-level system parameters. The model is validated using molecular dynamics simulations.

The interplay between chemical reaction and substrate deformation is discussed by adapting Ranz's formulation for scalar mixing to the case of a reactive mixture between segregated reactants, initially separated by an interface whose thickness may not be vanishingly small. Experiments in a simple shear flow demonstrate the existence of three regimes depending on the Damköhler number $Da=t_s/t_c$ where $t_s$ is the mixing time of the interface width and $t_c$ is the chemical time. Instead of treating explicitly the chemical cross-term, we rationalize these different regimes by globalizing it as a production term involving a flux which depends on the rate at which the reaction zone is fed by the reactants, a formulation valid for $Da>1$. For $Da<1$, the reactants interpenetrate before they react, giving rise to a ‘diffusio-chemical’ regime where chemical production occurs within a substrate whose width is controlled by molecular diffusion.

The central open question about Rayleigh–Bénard convection – buoyancy-driven flow in a fluid layer heated from below and cooled from above – is how vertical heat flux depends on the imposed temperature gradient in the strongly nonlinear regime where the flows are typically turbulent. The quantitative challenge is to determine how the Nusselt number $Nu$ depends on the Rayleigh number $Ra$ in the $Ra\to \infty$ limit for fluids of fixed finite Prandtl number $Pr$ in fixed spatial domains. Laboratory experiments, numerical simulations and analysis of Rayleigh's mathematical model have yet to rule out either of the proposed ‘classical’ $Nu \sim Ra^{1/3}$ or ‘ultimate’ $Nu \sim Ra^{1/2}$ asymptotic scaling theories. Among the many solutions of the equations of motion at high $Ra$ are steady convection rolls that are dynamically unstable but share features of the turbulent attractor. We have computed these steady solutions for $Ra$ up to $10^{14}$ with $Pr=1$ and various horizontal periods. By choosing the horizontal period of these rolls at each $Ra$ to maximize $Nu$, we find that steady convection rolls achieve classical asymptotic scaling. Moreover, they transport more heat than turbulent convection in experiments or simulations at comparable parameters. If heat transport in turbulent convection continues to be dominated by heat transport in steady rolls as $Ra\to \infty$, it cannot achieve the ultimate scaling.