We present new experiments of particle-driven turbulent plumes issuing from a constant source of dense particle-laden fluid, with buoyancy flux, $B$, in a uniform horizontal current, $u$. Experiments show that a turbulent, well-mixed plume develops, in which the downward vertical speed $w$ decreases with depth $z$ according to $w = 0.76 (B/uz)^{1/2}$ while the horizontal speed rapidly asymptotes to the current speed $u$, provided that the Stokes settling speed of the particles $v<0.92 w$. For $v > 0.92 w$, the particles separate from the plume fluid, and their depth $z$ increases according to the simple sedimentation trajectory $\textrm {d}z/{\textrm {d}\kern0.7pt x} = v/u$. As the particles sediment, they form clusters of particles, which lead to fluctuations in the particle load with position, but do not appear to change the time-average sedimentation speed. We explore the impact of these results for deep-sea mining, in which the fate of the plume water as well as the particles is key for assessing potential environmental impacts.

The behaviour of turbulent, buoyant, planar plumes is fundamentally coupled to the environment within which they develop. The effect of a background stratification directly influences a plumes buoyancy and has been the subject of numerous studies. Conversely, the effect of an ambient co-flow, which directly influences the vertical momentum of a plume, has not previously been the subject of theoretical investigation. The governing conservation equations for the case of a uniform co-flow are derived and the local dynamical behaviour of the plume is shown to be characterised by the scaled source Richardson number and the relative magnitude of the co-flow and plume source velocities. For forced, pure and lazy plume release conditions the co-flow acts to narrow the plume and reduce both the dilution and the asymptotic Richardson number relative to the classic zero co-flow case. Analytical solutions are developed for pure plumes from line sources, and for highly forced and highly lazy releases from sources of finite width in a weak co-flow. Contrary to releases in quiescent surroundings, our solutions show that all classes of release can exhibit plume contraction and the associated necking. For entraining plumes, a dynamical invariance spatially only occurs for pure and forced releases and we derive the co-flow strengths that lead to this invariance.

A detailed analysis of steady coating flow of a thin film of a viscous fluid on the outside of a uniformly rotating horizontal circular cylinder in the absence of surface-tension effects but in the presence of a non-uniform pressure distribution due to an irrotational airflow with circulation shows that the presence of the airflow can result in qualitatively different behaviour of the fluid film from that in classical coating flow. Full-film solutions corresponding to a continuous film of fluid covering the entire cylinder are possible only when the flux and mass of fluid do not exceed critical values, which are determined in terms of the non-dimensional parameters $F$ and $K$ representing the speed of the far-field airflow and the circulation of the airflow, respectively. The qualitative changes in the behaviour of the film thickness as $F$ and $K$ are varied are described. In particular, the film thickness can have as many as four stationary points and, in general, has neither top-to-bottom nor right-to-left symmetry. In addition, when the circulation of the airflow is in the same direction as the rotation of the cylinder the maximum mass of fluid that can be supported on the cylinder is always less than that in classical coating flow, whereas when the circulation is in the opposite direction the maximum mass of fluid can be greater than that in classical coating flow.

Fish typically swim by periodic bending of their bodies. Bending seems to follow a universal rule; it occurs at about one-third from the posterior end of the fish body with a maximum bending angle of about $30^{\circ }$. However, the hydrodynamic mechanisms that shaped this convergent design and its potential benefit to fish in terms of swimming speed and efficiency are not well understood. It is also unclear to what extent this bending is active or follows passively from the interaction of a flexible posterior with the fluid environment. Here, we use a self-propelled two-link model, with fluid–structure interactions described in the context of the vortex sheet method, to analyse the effects of both active and passive body bending on the swimming performance. We find that passive bending is more efficient but could reduce swimming speed compared with rigid flapping, but the addition of active bending could enhance both speed and efficiency. Importantly, we find that the phase difference between the posterior and anterior sections of the body is an important kinematic factor that influences performance, and that active antiphase flexion, consistent with the passive flexion phase, can simultaneously enhance speed and efficiency in a region of the design space that overlaps with biological observations. Our results are consistent with the hypothesis that fish that actively bend their bodies in a fashion that exploits passive hydrodynamics can at once improve speed and efficiency.

Precession driven flows are of great interest for both, industrial and geophysical applications. While cylindrical, spherical and spheroidal geometries have been investigated in great detail, the numerically and theoretically more challenging case of a non-axisymmetric cavity has received less attention. We report experimental results on the flows in a precessing triaxial ellipsoid, with a focus on the base flow of uniform vorticity, which we show to be in good agreement with existing theoretical models. As predicted, the uniform vorticity component exhibits two branches of solutions leading to a hysteresis cycle as a function of the Poincaré number. The first branch is observed at low forcing and characterized by large amplitude of the total fluid rotation and a moderate tilt angle of the fluid rotation axis. In contrast, the second branch displays only a moderate fluid rotation and a large tilt angle of the fluid rotation axis, which tends to align with the precession axis. In addition, we observe the occurrence of parametric instabilities early in the first branch, which saturate in the second branch, where we observe the same order of the kinetic energy in the base flow and instabilities.

We study the topology and the temporal dynamics of turbulent Rayleigh–Bénard convection in a liquid metal with a Prandtl number of 0.03 located inside a box with a square base area and an aspect ratio of $\varGamma = 5$. Experiments and numerical simulations are focused on the Rayleigh number range $6.7 \times 10^{4} \leqslant Ra \leqslant 3.5 \times 10^{5}$, where a new cellular flow regime has been reported previously (Akashi et al., Phys. Rev. Fluids, vol. 4, 2019, 033501). This flow structure shows symmetries with respect to the vertical planes crossing at the centre of the container. The dynamic behaviour is dominated by strong three-dimensional oscillations with a period length that corresponds to the turnover time. Our analysis reveals that the flow structure in the $\varGamma = 5$ box corresponds in key features to the jump rope vortex structure, which has recently been discovered in a $\varGamma = 2$ cylinder (Vogt et al., Proc. Natl Acad. Sci. USA, vol. 115, 2018, pp. 12674–12679). While in the $\varGamma = 2$ cylinder a single jump rope vortex occurs, the coexistence of four recirculating swirls is detected in this study. Their approach to the lid or the bottom of the convection box causes a temporal deceleration of both the horizontal velocity at the respective boundary and the vertical velocity in the bulk, which in turn is reflected in Nusselt number oscillations. The cellular flow regime shows remarkable similarities to properties commonly attributed to turbulent superstructures.

We address a challenge of active flow control: the optimization of many actuation parameters guaranteeing fast convergence and avoiding suboptimal local minima. This challenge is addressed by a new optimizer, called the explorative gradient method (EGM). EGM alternatively performs one exploitive downhill simplex step and an explorative Latin hypercube sampling iteration. Thus, the convergence rate of a gradient based method is guaranteed while, at the same time, better minima are explored. For an analytical multi-modal test function, EGM is shown to significantly outperform the downhill simplex method, the random restart simplex, Latin hypercube sampling, Monte Carlo sampling and the genetic algorithm. EGM is applied to minimize the net drag power of the two-dimensional fluidic pinball benchmark with three cylinder rotations as actuation parameters. The net drag power is reduced by 29 % employing direct numerical simulations at a Reynolds number of $100$ based on the cylinder diameter. This optimal actuation leads to 52 % drag reduction employing Coanda forcing for boat tailing and partial stabilization of vortex shedding. The price is an actuation energy corresponding to 23 % of the unforced parasitic drag power. EGM is also used to minimize drag of the $35^\circ$ slanted Ahmed body employing distributed steady blowing with 10 inputs. 17 % drag reduction are achieved using Reynolds-averaged Navier–Stokes simulations at the Reynolds number $Re_H=1.9 \times 10^5$ based on the height of the Ahmed body. The wake is controlled with seven local jet-slot actuators at all trailing edges. Symmetric operation corresponds to five independent actuator groups at top, middle, bottom, top sides and bottom sides. Each slot actuator produces a uniform jet with the velocity and angle as free parameters, yielding 10 actuation parameters as free inputs. The optimal actuation emulates boat tailing by inward-directed blowing with velocities which are comparable to the oncoming velocity. We expect that EGM will be employed as efficient optimizer in many future active flow control plants as alternative or augmentation to pure gradient search or explorative methods.

Classical scaling relationships for rheological quantities such as the $\mu (J)$-rheology have become increasingly popular for closures of two-phase flow modelling. However, these frameworks have been derived for monodisperse particles. We aim to extend these considerations to sediment transport modelling by using a more realistic sediment composition. We investigate the rheological behaviour of sheared sediment beds composed of polydisperse spherical particles in a laminar Couette-type shear flow. The sediment beds consist of particles with a diameter size ratio of up to 10, which corresponds to grains ranging from fine to coarse sand. The data was generated using fully coupled, grain resolved direct numerical simulations using a combined lattice Boltzmann–discrete element method. These highly resolved data yield detailed depth-resolved profiles of the relevant physical quantities that determine the rheology, i.e. the local shear rate of the fluid, particle volume fraction, total shear and granular pressure. A comparison against experimental data shows excellent agreement for the monodisperse case. We improve upon the parameterization of the $\mu (J)$-rheology by expressing its empirically derived parameters as a function of the maximum particle volume fraction. Furthermore, we extend these considerations by exploring the creeping regime for viscous numbers much lower than used by previous studies to calibrate these correlations. Considering the low viscous numbers of our data, we found that the friction coefficient governing the quasi-static state in the creeping regime tends to a finite value for vanishing shear, which decreases the critical friction coefficient by a factor of three for all cases investigated.

Indoor ventilation is essential for a healthy and comfortable living environment. A key issue is to discharge anthropogenic air contamination such as CO$_2$ gas or, of potentially more direct consequence, airborne respiratory droplets. Here, by employing direct numerical simulations, we study mechanical displacement ventilation with a wide range of ventilation rates $Q$ from 0.01 to 0.1 m$^3$ s$^{-1}$ person$^{-1}$. For this ventilation scheme, a cool lower zone is established beneath a warm upper zone with interface height $h$, which depends on $Q$. For weak ventilation, we find the scaling relation $h\sim Q^{3/5}$, as suggested by Hunt & Linden (Build. Environ., vol. 34, 1999, pp. 707–720). Also, the CO$_{2}$ concentration decreases with $Q$ within this regime. However, for too strong ventilation, the interface height $h$ becomes insensitive to $Q$, and the ambient averaged CO$_2$ concentration decreases towards the ambient value. At these values of $Q$, the concentrations of pollutants are very low and so further dilution has little effect. We suggest that such scenarios arise when the vertical kinetic energy associated with the ventilation flow is significant compared with the potential energy of the thermal stratification.

We experimentally, numerically and theoretically investigate the nonlinear interaction between a cavitation bubble and the interface of two immiscible fluids (oil and water) on multiple time scales. The underwater electric discharge method is utilized to generate a cavitation bubble near or at the interface. Both the bubble dynamics on a short time scale and the interface evolution on a much longer time scale are recorded via high-speed photography. Two mechanisms are found to contribute to the fluid mixing in our system. First, when a bubble is initiated in the oil phase or at the interface, an inertia-dominated high-speed liquid jet generated from the collapsing bubble penetrates the water–oil interface, and consequently transports fine oil droplets into the water. The critical standoff parameter for jet penetration is found to be highly dependent on the density ratio of the two fluids. Furthermore, the pinch-off of an interface jet produced long after the bubble dynamics stage is reckoned as the second mechanism, carrying water droplets into the oil bulk. The dependence of the bubble jetting behaviours and interface jet dynamics on the governing parameters is systematically studied via experiments and boundary integral simulations. Particularly, we quantitatively demonstrate the respective roles of surface tension and viscosity in interface jet dynamics. As for a bubble initiated at the interface, an extended Rayleigh–Plesset model is proposed that well predicts the asymmetric dynamics of the bubble, which accounts for a faster contraction of the bubble top and a downward liquid jet.

Using complementary experiments and direct numerical simulations, we study turbulent thermal convection of a liquid metal (Prandtl number $\textit {Pr}\approx 0.03$) in a box-shaped container, where two opposite square sidewalls are heated/cooled. The global response characteristics like the Nusselt number ${\textit {Nu}}$ and the Reynolds number $\textit {Re}$ collapse if the side height $L$ is used as the length scale rather than the distance $H$ between heated and cooled vertical plates. These results are obtained for various Rayleigh numbers $5\times 10^3\leq {\textit {Ra}}_H\leq 10^8$ (based on $H$) and the aspect ratios $L/H=1, 2, 3$ and $5$. Furthermore, we present a novel method to extract the wind-based Reynolds number, which works particularly well with the experimental Doppler-velocimetry measurements along vertical lines, regardless of their horizontal positions. The extraction method is based on the two-dimensional autocorrelation of the time–space data of the vertical velocity.

The subsonic, compressible, potential flow around a hypersphere can be derived using the Janzen–Rayleigh expansion (JRE) of the flow potential in even powers of the incident Mach number ${\mathcal {M}}_\infty$. JREs were carried out with terms polynomial in the inverse radius $r^{-1}$ to high orders in two dimensions, but were limited to order ${\mathcal {M}}_\infty ^{4}$ in three dimensions. We derive general JRE formulae for arbitrary order, adiabatic index and dimension. We find that powers of $\ln (r)$ can creep into the expansion, and are essential in the three-dimensional (3-D) sphere beyond order ${\mathcal {M}}_\infty ^{4}$. Such terms are apparently absent in the 2-D disk, as we verify up to order ${\mathcal {M}}_\infty ^{100}$, although they do appear in other dimensions (e.g. at order ${\mathcal {M}}_\infty ^{2}$ in four dimensions). An exploration of various 2-D and 3-D bodies suggests a topological connection, with logarithmic terms emerging when the flow is simply connected. Our results have additional physical implications. They are used to improve the hodograph-based approximation for the flow in front of a sphere. The symmetry-axis velocity profiles of axisymmetric flows around different prolate spheroids are approximately related to each other by a simple, Mach-independent scaling.

For an immiscible oil drop immersed in a stably stratified ethanol–water mixture, a downwards solutal Marangoni flow is generated on the surface of the drop, owing to the concentration gradient, and the resulting propulsion competes against the downwards gravitational acceleration of the heavy drop. In prior work of Li et al. (Phys. Rev. Lett., vol. 126, issue 12, 2021, 124502), we found that for drops of low viscosity, an oscillatory instability of the Marangoni flow is triggered once the Marangoni advection is too strong for diffusion to restore the stratified concentration field around the drop. Here we experimentally explore the parameter space of the concentration gradient and drop radius for high oil viscosities and find a different and new mechanism for triggering the oscillatory instability in which diffusion is no longer the limiting factor. For such drops of higher viscosities, the instability is triggered when the gravitational effect is too strong so that the viscous stress cannot maintain a stable Marangoni flow. This leads to a critical drop radius above which the equilibrium is always unstable. Subsequently, a unifying scaling theory that includes both the mechanisms for low and for high viscosities of the oil drops is developed. The transition between the two mechanisms is found to be controlled by two length scales: the drop radius $R$ and the boundary layer thickness $\delta$ of the Marangoni flow around the drop. The instability is dominated by diffusion for $\delta < R$ and by viscosity for $R<\delta$. The experimental results for various drops of different viscosities can well be described with this unifying scaling theory. Our theoretical description thus provides a unifying view of physicochemical hydrodynamic problems in which the Marangoni stress is competing with a stable stratification.

The Richardson-scaling law states that the mean square separation of a fluid particle pair grows according to t3 within the inertial range and at intermediate times. The theories predicting this scaling regime assume that the pair separation is within the inertial range and that the dispersion is local, which means that only eddies at the scale of the separation contribute. These assumptions ignore the structural organization of the turbulent flow into large-scale shear layers, where the intense small-scale motions are bounded by the large-scale energetic motions. Therefore, the large scales contribute to the velocity difference across the small-scale structures. It is shown that, indeed, the pair dispersion inside these layers is highly non-local and approaches Taylor dispersion in a way that is fundamentally different from the Richardson-scaling law. Also, the layer's contribution to the overall mean square separation remains significant as the Reynolds number increases. This calls into question the validity of the theoretical assumptions. Moreover, a literature survey reveals that, so far, t3 scaling is not observed for initial separations within the inertial range. We propose that the intermediate pair dispersion regime is a transition region that connects the initial Batchelor- with the final Taylor-dispersion regime. Such a simple interpretation is shown to be consistent with observations and is able to explain why t3 scaling is found only for one specific initial separation outside the inertial range. Moreover, the model incorporates the observed non-local contribution to the dispersion, because it requires only small-time-scale properties and large-scale properties.

The dependence of Reynolds shear stress on wave phase is investigated for initially isotropic turbulence distorted by a progressive surface wave through direct numerical simulations. A wave-following streamline coordinate frame is used to analyse the turbulence dynamics such that the information of the varying direction of wave orbital motions is embedded into the coordinate system, which helps capture the effect of flow curvature on the turbulence and quantify momentum exchange between the near-surface and deep regions in the wavy domain accurately. It is found that the Reynolds shear stress is enhanced under the backward slope of the wave and can be scaled by the wave steepness and the streamwise velocity fluctuations. Analyses of the budget of Reynolds shear stress indicate that such a variation with the wave phase is caused by the variation in the production of the Reynolds shear stress and the effect of pressure fluctuations. Further investigation shows that the production of the Reynolds shear stress is closely associated with the wave surface curvature. A model that includes a correction term for the curvature effect for the pressure–strain correlation term is examined and is found to agree reasonably well with the simulation result. The correction term is found to make an appreciable contribution to the model, further supporting our finding that the wave curvature plays an important role in the turbulence dynamics near the surface.

The energy spectrum is commonly used to describe the scale dependence of turbulent fluctuations in homogeneous isotropic turbulence. In contrast, one-point statistical quantities, such as the turbulent kinetic energy, are mainly employed for inhomogeneous turbulence models. Attempts have been made to describe the scale dependence of inhomogeneous turbulence using the second-order structure function and two-point velocity correlation. However, unlike the energy spectrum, expressions for the energy density in the scale space fail to satisfy the requirement of being non-negative. In this study, a new expression for the scale-space energy density based on filtered velocities is proposed to clarify the reasons behind the negative values of the energy density and to obtain a better understanding of inhomogeneous turbulence. The new expression consists of homogeneous and inhomogeneous parts; the former is always non-negative, while the latter can be negative because of the turbulence inhomogeneity. Direct numerical simulation data of homogeneous isotropic turbulence and a turbulent channel flow are used to evaluate the two parts of the energy density and turbulent energy. It was found that the inhomogeneous part of the turbulent energy shows non-zero values near the wall and at the centre of a channel flow. In particular, the inhomogeneous part of the energy density changes its sign depending on the scale. A concave profile of the filtered-velocity variance at the wall accounts for the negative value of the energy density in the region very close to the wall.

The scaling of different features of streamwise normal stress profiles $\langle uu\rangle ^+(y^+)$ in turbulent wall-bounded flows is the subject of a long-running debate. Particular points of contention are the scaling of the ‘inner’ and ‘outer’ peaks of $\langle uu\rangle ^+$ at $y^+\approxeq ~15$ and $y^+ ={O}(10^3)$, respectively, their infinite Reynolds number limit, and the rate of logarithmic decay in the outer part of the flow. Inspired by the thought-provoking paper of Chen & Sreenivasan (J. Fluid Mech., vol. 908, 2021, p. R3), two terms of an inner asymptotic expansion of $\langle uu\rangle ^+$ in the small parameter $Re_{\tau }^{-1/4}$ are constructed from a set of direct numerical simulations (DNS) of channel flow. This inner expansion is for the first time matched through an overlap layer to an outer expansion, which not only fits the same set of channel DNS within 1.5 % of the peak stress, but also provides a good match of laboratory data in pipes and the near-wall part of boundary layers, up to the highest $Re_{\tau }$ values of $10^5$. The salient features of the new composite expansion are first, an inner $\langle uu\rangle ^+$ peak, which saturates at 11.3 and decreases as $Re_{\tau }^{-1/4}$. This inner peak is followed by a short ‘wall log law’ with a slope that becomes positive for $Re_{\tau }$ beyond ${O}(10^4)$, leading up to an outer peak, followed by the logarithmic overlap layer with a negative slope going continuously to zero for $Re_{\tau }\to \infty$.

The wake of a notchback Ahmed body presenting a bi-stable nature is investigated by performing wind tunnel experiments and large-eddy simulations. Attention is confined to the Reynolds number ($Re$) influence on the wake state instability within $5\times 10^{4}\leq Re \leq 25\times 10^{4}$. Experimental observations suggest a wake bi-stability with low-frequency switches under low $Re$. The wake becomes ‘tri-stable’ with the increase of $Re$ with the introduction of a new symmetric state. The higher presence of the symmetric state can be considered as a symmetrization of the wake bi-stability with an increasing $Re$. The wake symmetry under high $Re$ attributed to the highly frequent switches of the wake is extremely sensitive to small yaw angles, showing the feature of bi-stable flows. The wake asymmetry is confirmed in numerical simulations with both low and high $Re$. The wake asymmetries are indicated by the wake separation, the reattachment and the wake dynamics identified by the proper orthogonal decomposition. However, the turbulence level is found to be significantly higher with a higher $Re$. This leads to a higher possibility to break the asymmetric state, resulting in highly frequent switches showing symmetry.

Fault zones have the potential to act as leakage pathways through low permeability structural seals in geological reservoirs. Faults may facilitate migration of groundwater contaminants and stored anthropogenic carbon dioxide (CO$_2$), where the waste fluids would otherwise remain securely trapped. We present an analytical model that describes the dynamics of leakage through a fault zone cutting multiple aquifers and seals. Current analytical models for a buoyant plume in a semi-infinite porous media are combined with models for a leaking gravity current and a new model motivated by experimental observation, to account for increased pressure gradients within the fault due to an increase in Darcy velocity directly above the fault. In contrast to previous analytical fault models, we verify our results using a series of analogous porous medium tank experiments, with good matching of observed leakage rates and fluid distribution. We demonstrate the utility of the model for the assessment of CO$_2$ storage security, by application to a naturally occurring CO$_2$ reservoir, showing the dependence of the leakage rates and fluid distribution on the fault/aquifer permeability contrast. The framework developed within this study can be used for quick assessment of fluid leakage through fault zones, given a set of input parameters relating to properties of the fault, aquifer and fluids, and can be incorporated into basin-scale models to improve computational efficiency. The results show the utility of using analytical methods and reduced-order modelling in complex geological systems, as well as the value of laboratory porous medium experiments to verify results.

Fully developed turbulent flow in channels with mild to strong longitudinal curvature is studied by direct numerical simulations. The Reynolds based on the bulk mean velocity and channel half-width $\delta$ is fixed at $20\,000$, resulting in a friction Reynolds number of approximately 1000. Four cases are considered with curvature varying from $\gamma = 2\delta /r_c = 0.033$ to 0.333, where $r_c$ is the curvature radius at the channel centre. Substantial differences between the mean wall shear stress on the convex and concave walls are already observed for $\gamma = 0.033$. A log-law region is absent and a region with nearly constant mean angular momentum develops in the channel centre for strong curvatures. Spanwise and wall-normal velocity fluctuations are strongly amplified by curvature in the outer region of the concave channel side. Only near the walls, where curvature effects are relatively weak, do the mean velocity and velocity fluctuation profiles approximately collapse when scaled by wall units based on the local friction velocity. Budgets of the streamwise and wall-normal Reynolds-stress equations are presented and turbulence structures are investigated through visualizations and spectra. In the case with strongest curvature, the flow relaminarizes locally near the convex wall. On the concave channel side, large elongated streamwise vortices reminiscent of Taylor–Görtler vortices develop for all curvatures considered. The maximum in the premultiplied two-dimensional wall-normal energy spectrum and co-spectrum shifts towards larger scales with increasing curvature. The large scales substantially contribute to the wall-normal velocity fluctuations and momentum transport on the concave channel side.