Steep wave breaking on a vertical cylinder (a typical foundation supporting offshore wind turbines) will induce slam loads. Many questions on the important violent wave loading and the associated secondary load cycle remain unanswered. We use laboratory experiments with unidirectional waves to investigate the fluid loading on vertical cylinders. We use a novel three-phase decomposition approach that allows us to separate different types of nonlinearity. Our findings reveal the existence of an additional quasi-impulsive loading component that is associated with the secondary load cycle and occurs in the backwards direction against that of the incoming waves. This quasi-impulsive force occurs at the end of the secondary load cycle and close to the passage of the downward zero-crossing point of the undisturbed wave. Wavelet analysis showed that the impulsive force exhibits superficially similar behaviour to a typical wave-slamming event but in the reverse direction. To monitor the scattered wave field and extract run-up on the cylinder, we installed a four-camera synchronised video system and found a strong temporal correlation between the arrival time of the Type-II scattered wave onto the cylinder and the occurrence of this quasi-impulsive force. The temporal characteristics of this quasi-impulsive force can be approximated by the Goda wave impact model, taking the collision of the Type-II scattered waves at the rear stagnation point as the impact source.

The background potential energy (BPE) is the only reservoir that double diffusive instabilities can tap their energy from when developing from an unforced motionless state with no available potential energy (APE). Recently, Middleton and Taylor linked the extraction of BPE into APE to the sign of the diapycnal component of the buoyancy flux, but their criterion can predict only diffusive convection instability, not salt finger instability. Here, we show that the problem can be corrected if the sign of the APE dissipation rate is used instead, making it emerge as the most fundamental criterion for double diffusive instabilities. A theory for the APE dissipation rate for a two-component fluid relative to its single-component counterpart is developed as a function of three parameters: the diffusivity ratio, the density ratio, and a spiciness parameter. The theory correctly predicts the occurrence of both salt finger and diffusive convection instabilities in the laminar unforced regime, while more generally predicting that the APE dissipation rate for a two-component fluid can be enhanced, suppressed, or even have the opposite sign compared to that for a single-component fluid, with important implications for the study of ocean mixing. Because negative APE dissipation can also occur in stably stratified single-component and doubly stable two-component stratified fluids, we speculate that only the thermodynamic theory of exergy can explain its physics; however, this necessitates accepting that APE dissipation is a conversion between APE and the internal energy component of BPE, in contrast to prevailing assumptions.

Resistive tearing instabilities are common in fluids that are highly electrically conductive and carry strong currents. We determine the effect of stable stratification on the tearing instability under the Boussinesq approximation. Our results generalise previous work that considered only specific parameter regimes, and we show that the length scale of the fastest-growing mode depends non-monotonically on the stratification strength. We confirm our analytical results by solving the linearised equations numerically, and we discuss whether the instability could operate in the solar tachocline.

Whether colliding drops will merge with or bounce off each other is critical to numerous processes, and the physics involved is notoriously complex. In particular, experiments show that both sufficiently slow and fast head-on drop collisions lead to merging, but that there is often an intermediate regime in which bouncing is observed; these transitions in behaviour were recently discovered to be surprisingly sensitive to the radius of the drops and the ambient gas pressure. We show here that these transitions between bouncing and merging are governed by nanoscale phenomena; namely, gas-kinetic and disjoining pressure effects. To capture these crucial effects, a novel, open-source computational model is developed for the simulation of colliding drops. The model uses a hybrid approach, based on solving the Navier–Stokes equations in the drop with a lubrication approach for the unconventional physics of the gas film. Our simulations show remarkably good agreement with experiments of head-on collisions and also provide new experimentally verifiable predictions.

Linear unsteady aerofoil theory, while successfully used for the prediction of unsteady aerofoil lift for many decades, has yet to be proven adequate for predicting the propulsive performance of oscillating aerofoils. In this paper we test the hypothesis that the central shortcoming of linear small-amplitude models, such as the Garrick function, is the failure to account for the flow acceleration caused by aerofoil thrust. A new analytical model is developed by coupling the Garrick function to a cycle-averaged actuator disc model, in a manner analogous to the blade-element momentum theory for wind turbines and propellers. This amounts to assuming the Garrick function to be locally valid and, in combination with a global control volume analysis, enables the prediction of flow acceleration at the aerofoil. The new model is demonstrated to substantially improve the agreement with large-eddy simulations of an aerofoil in combined heave and pitch motion.

In this experimental study, we explore the dynamics of the thermal boundary layer in liquid metal Rayleigh–Bénard convection, covering the parameter ranges of $0.026 \leq$ Prandtl numbers $(Pr) \leq 0.033$ and Rayleigh numbers ($Ra$) up to $2.9\times 10^9$. Our research focuses on characterising the thermal boundary layer near the top plate of a cylindrical convection cell with an aspect ratio of 0.5, distinguishing between two distinct regions: the shear-dominated region around the centre of the top plate and a location near the side wall where the boundary layer is expected to be affected by the impact or ejection of thermal plumes. The dependencies of the boundary layer thickness on $Ra$ at these positions reveal deviating scaling exponents with the difference diminishing as $Ra$ increases. We find stronger fluctuations in the boundary layer and increasing deviation from the Prandtl–Blasius–Pohlhausen profile with increasing $Ra$, as well as in the measurements outside the centre region. Our data illustrate the complex interplay between flow dynamics and thermal transport in low-$Pr$ convection.

A high-resolution numerical simulation of an air–water turbulent upward bubbly flow in a pipe is performed to investigate the turbulence characteristics and bubble interaction with the wall. We consider three bubble equivalent diameters and three total bubble volume fractions. The bulk and bubble Reynolds numbers are $Re_{bulk}= u_{bulk} D/\nu _w = 5300$ and $Re_{bub}= (\langle u_{bub}\rangle - u_{bulk}) d_{eq}/\nu _w = 533\unicode{x2013}1000$, respectively, where $u_{bulk}$ is the water bulk velocity, $\langle u_{bub}\rangle$ is the overall bubble mean velocity, $D$ is the pipe diameter and $\nu _w$ is the water kinematic viscosity. The mean water velocity near the wall significantly increases due to bubble interaction with the wall, and the root-mean-square water velocity fluctuations are proportional to $\bar {\psi }(r)^{0.4}$, where $\bar {\psi } (r)$ is the mean bubble volume fraction. For the cases considered, the bubble-induced turbulence suppresses the shear-induced turbulence and becomes the dominant flow characteristic at all radial locations including near the wall. Rising bubbles near the wall mostly bounce against the wall rather than slide along the wall or hang around the wall without collision. Low-speed streaks observed in the near-wall region in the absence of bubbles nearly disappear due to the bouncing bubbles. These bouncing bubbles generate counter-rotating vortices in their wake, and increase the skin friction by sweeping high-speed water towards the wall. We also suggest an algebraic Reynolds-averaged Navier–Stokes model considering the interaction between shear-induced and bubble-induced turbulence. This model provides accurate predictions for a wide range of liquid bulk Reynolds numbers.

We consider an internally heated fluid between parallel plates with fixed thermal fluxes. For a large class of heat sources that vary in the direction of gravity, we prove that $\smash { \smash {{\langle {\delta T} \rangle _h}} } \geq \sigma R^{-1/3} - \mu$, where $\smash { \smash {{\langle {\delta T} \rangle _h}} }$ is the average temperature difference between the bottom and top plates, $R$ is a ‘flux’ Rayleigh number and the constants $\sigma,\mu >0$ depend on the geometric properties of the internal heating. This result implies that mean downward conduction (for which $\smash { \smash {{\langle {\delta T} \rangle _h}} }< 0$) is impossible for a range of Rayleigh numbers smaller than a critical value $R_0:=(\sigma /\mu )^{3}$. The bound demonstrates that $R_0$ depends on the heating distribution and can be made arbitrarily large by concentrating the heating near the bottom plate. However, for any given fixed heating profile of the class we consider, the corresponding value of $R_0$ is always finite. This points to a fundamental difference between internally heated convection and its limiting case of Rayleigh–Bénard convection with fixed-flux boundary conditions, for which $\smash {{\langle {\delta T} \rangle _h}}$ is known to be positive for all $R$.

Two-dimensional free-surface flow over localised topography is examined, with the emphasis on the stability of hydraulic-fall solutions. A Gaussian topography profile is assumed with a positive or negative amplitude modelling a bump or a dip, respectively. Steady hydraulic-fall solutions to the full incompressible, irrotational Euler equations are computed, and their linear and nonlinear stability is analysed by computing eigenspectra of the pertinent linearised operator and by solving an initial value problem. The computations are carried out numerically using a specially developed computational framework based on the finite-element method. The Hamiltonian structure of the problem is demonstrated, and stability is determined by computing eigenspectra of the pertinent linearised operator. It is found that a hydraulic-fall flow over a bump is spectrally stable. The corresponding flow over a dip is found to be linearly unstable. In the latter case, time-dependent simulations show that ultimately, the flow settles into a time-periodic motion that corresponds to an invariant solution in an appropriately defined phase space. Physically, the solution consists of a localised large-amplitude wave that pulsates above the dip while simultaneously emitting nonlinear cnoidal waves in the upstream direction and multi-harmonic linear waves in the downstream direction.

Adding polymers to turbulent Newtonian fluid flows can have dramatic effects. A well-known example is a significant drag reduction by flexible polymers in turbulent wall-bounded flows. In numerical studies of polymer drag reduction, polymers are often modelled as dumbbells of two beads connected by a finitely extensible nonlinear elastic (FENE) spring. There are natural queries whether this highly simplified coarse-grained model is adequate for describing a polymer macromolecule in turbulent flows. By carrying out Eulerian–Lagrangian simulations of polymers, described by different models, in a turbulent pipe flow, Serafini et al. (J. Fluid Mech., vol. 987, 2024, R1) have demonstrated that the FENE dumbbell model can accurately capture polymer extension statistics as compared with the realistic Kuhn chain model. Their work further reveals the surprising result that increasing the number of beads in a FENE chain worsens its accuracy in characterizing polymer spatial conformations at large Weissenberg numbers.

Many biological fluids are composed of suspended polymers immersed in a viscous fluid. A prime example is mucus, where the polymers are also known to form a network. While the presence of this microstructure is linked with an overall non-Newtonian response of the fluid, swimming cells and microorganisms similar in size to the network pores and polymer filaments instead experience the heterogeneous nature of the environment, interacting directly with the polymers as obstacles as they swim. To characterise and understand locomotion in these heterogeneous environments, we simulate the motion of an undulatory swimmer through three-dimensional suspensions and networks of elastic filaments, exploring the effects of filament and link compliance and filament concentration up to 20 % volume fraction. For compliant environments, the swimming speed increases with filament concentration to values approximately 10 % higher than in a viscous fluid. In stiffer environments, a non-monotonic dependence is observed, with an initial increase in speed to values 5 % greater than in a viscous fluid, followed by a dramatic reduction to speeds just a fraction of its value in a viscous fluid. Velocity fluctuations are also more pronounced in stiffer environments. We demonstrate that speed enhancements are linked to hydrodynamic interactions with the microstructure, while reductions are due to the filaments restricting the amplitude of the swimmer's propulsive wave. Unlike previous studies where interactions with obstacles allowed for significant enhancements in swimming speeds, the modest enhancements seen here are more comparable to those given by models where the environment is treated as a continuous viscoelastic fluid.

Bluff-body wakes generally become three-dimensional (3-D) and then turbulent when the Reynolds number exceeds a few hundred. Other than an alternate shedding of the spanwise vortices behind the body and a gradual decay and annihilation of the vortices with distance downstream, whether a secondary vortex street would develop in the relatively far wake has been a long-standing argument in the literature. This argument is addressed in the present study. Specifically, direct numerical simulations and transient growth analysis are performed to examine the two-dimensional and 3-D wakes of different bluff bodies, including circular cylinder, square cylinder, diamond cylinder and rectangular cylinders with different cross-sectional aspect ratios. We found that a secondary vortex street is absent for most 3-D and turbulent wakes. The root cause is the weakening of spanwise vortices by 3-D wake instability modes and streamwise circulation/vorticity. The weakened spanwise vortices induce reduced mean shear in the intermediate wake, which then induces much smaller perturbation energy growth that is below the threshold for the emergence of a secondary vortex street. This finding suggests that the 3-D and turbulence characteristics, and the momentum, mass and heat transport in the relatively far wake of bluff bodies, would not be influenced by extra anisotropy or inhomogeneity caused by a secondary vortex street.

We introduce a low-order dynamical system to describe thermal convection in an annular domain. The model derives systematically from a Fourier–Laurent truncation of the governing Navier–Stokes Boussinesq equations and accounts for spatial dependence of the flow and temperature fields. Comparison with fully resolved direct numerical simulations (DNS) shows that the model captures parameter bifurcations and reversals of the large-scale circulation (LSC), including states of (i) steady circulating flow, (ii) chaotic LSC reversals and (iii) periodic LSC reversals. Casting the system in terms of the fluid's angular momentum and centre of mass (CoM) reveals equivalence to a damped pendulum with forcing that raises the CoM above the fulcrum. This formulation offers a transparent mechanism for LSC reversals, namely the inertial overshoot of a forced pendulum, and it yields an explicit formula for the frequency $f^*$ of regular LSC reversals in the high-Rayleigh-number (Ra) limit. This formula is shown to be in excellent agreement with DNS and produces the scaling law $f^* \sim {Ra}^{0.5}$.

The influence of symmetry-breaking effects of ridge-type roughness on secondary currents in turbulent channel flow is investigated using direct numerical simulations. The ridges have triangular cross-section, which is systematically varied from isosceles to right-angled triangle, introducing an imbalance to the slopes of the ridges’ lateral surfaces while the streamwise homogeneity of the surfaces is maintained. In all cases, secondary current vortices are produced, but asymmetric ridge cross-sections break the symmetry of these vortices. As a result of the asymmetry-induced misalignment and imbalance in the secondary current vortices, net spanwise flow emerges. The magnitude of the spanwise flow increases with the slope ratio of the ridge lateral surfaces and significantly modifies the mean flow topology, leading to the merging of critical points in the case of the right-angled triangular ridge shape. Within the cavities, the net spanwise flow is accompanied by a non-zero mean spanwise pressure gradient, while from the perspective of the outer flow, the scalene ridge surfaces have a similar effect as a wall that is slowly moving in the spanwise direction. Overall, the present results suggest the existence of a special type of Prandtl's secondary currents of the second kind, namely those that result in net spanwise flow.

We present a systematic study on the effects of small aspect ratios $\varGamma$ on heat transport in liquid metal convection with a Prandtl number of $Pr=0.029$. The study covers $1/20\le \varGamma \le 1$ experimentally and $1/50\le \varGamma \le 1$ numerically, and a Rayleigh number $Ra$ range of $4\times 10^3 \le Ra \le 7\times 10^{9}$. It is found experimentally that the local effective heat transport scaling exponent $\gamma$ changes with both $Ra$ and $\varGamma$, attaining a $\varGamma$-dependent maximum value before transition-to-turbulence and approaches $\gamma =0.25$ in the turbulence state as $Ra$ increases. Just above the onset of convection, Shishkina (Phys. Rev. Fluids, vol 6, 2021, 090502) derived a length scale $\ell =H/(1+1.49\varGamma ^{-2})^{1/3}$. Our numerical study shows $Ra_{\ell }$, i.e. $Ra$ based on $\ell$, serves as a proper control parameter for heat transport above the onset with $Nu-1=0.018(1+0.34/\varGamma ^2)(Ra/Ra_{c,\varGamma }-1)$. Here $Ra_{c,\varGamma }$ represents the $\varGamma$-dependent critical $Ra$ for the onset of convection and $Nu$ is the Nusselt number. In the turbulent state, for a general scaling law of $Nu-1\sim Ra^\alpha$, we propose a length scale $\ell = H/(1+1.49\varGamma ^{-2})^{1/[3(1-\alpha )]}$. In the case of turbulent liquid metal convection with $\alpha =1/4$, our measurement shows that the heat transport will become weakly dependent on $\varGamma$ with $Ra_{\ell }\equiv Ra/(1+1.49\varGamma ^{-2})^{4/3} \ge 7\times 10^5$. Finally, once the flow becomes time-dependent, the growth rate of $Nu$ with $Ra$ declines compared with the linear growth rate in the convection state. A hysteresis is observed in a $\varGamma =1/3$ cell when the flow becomes time-dependent. Measurements of the large-scale circulation suggest the hysteresis is caused by the system switching from a single-roll-mode to a double-roll-mode in an oscillation state.

Liquid metal buoyant flow around two differentially heated horizontal cylinders in the presence of a uniform vertical magnetic field is investigated experimentally. While magneto-convection in pipes or ducts has been studied theoretically and experimentally in recent years, data for heat transfer at immersed obstacles are rare and, to our knowledge, detailed experimental investigations on this fundamental magnetohydrodynamic problem do not exist. In the present work, two horizontal cylinders inserted into an adiabatic rectangular cavity filled with gallium–indium–tin are kept at constant temperatures to establish a driving temperature gradient in the surrounding liquid metal. The buoyancy-driven flow, quantified by the Grashof number $Gr$, is varied in the range ${10^{6} \leq Gr \leq ~5\times 10^{7}}$. With increasing magnetic field, expressed via the Hartmann number $Ha$, different flow regimes are identified from measurements for $0 \leq Ha \leq ~3000$. The effect of the electromagnetic force primarily consists in suppressing turbulence and damping the convective flow. The heat transfer is quantified in terms of the non-dimensional Nusselt number $Nu$, and its dependence on $Gr/{Ha}^{2}$, which is identified as the important group governing the flow, is discussed.

Diffuse interface models are an important class of models used to describe multi-phase flows. In the case of incompressible viscous fluids there are a number of different diffuse interface models which have been known for many years. Nevertheless, a model based on mixture theory with a full set of conservation laws for the conservation of linear momentum of each constituent was not yet available. This gap was filled by ten Eikelder et al. (J. Fluid Mech., in press) recently and a first comparison with known models of Navier–Stokes/Cahn–Hilliard type is given. A detailed understanding of the relations between these models remains an important question.

The nature and behaviour of the drag coefficient $C_D$ of irregularly shaped grains within a wide range of Reynolds numbers $Re$ is discussed. The morphology of the grains is controlled by their fractal description, and they differ in shape. Using computational fluid dynamics tools, the characteristics of the boundary layer at high $Re$ has been determined by applying the Reynolds-averaged Navier–Stokes turbulence model. Both grid resolution and mesh size dependence are validated with well-reported previous experimental results applied in flow around isolated smooth spheres. The drag coefficient for irregularly shaped grains is shown to be higher than that for spherical shapes, also showing a strong drop in its value at high $Re$. This drag crisis is reported at lower $Re$ compared to the smooth sphere, but higher critical $C_D$, demonstrating that the morphology of the particle accelerates this crisis. Furthermore, the dependence of $C_D$ on $Re$ in this type of geometry can be represented qualitatively by four defined zones: subcritical, critical, supercritical and transcritical. The orientational dependence for both particles with respect to the fluid flow is analysed, where our findings show an interesting oscillatory behaviour of $C_D$ as a function of the angle of incidence, fitting the results to a sine-squared interpolation, predicted for particles within the Stokes laminar regime ($Re\ll 1$) and for elongated/flattened spheroids up to $Re=2000$. A statistical analysis shows that this system satisfies a Weibullian behaviour of the drag coefficient when random azimuthal and polar rotation angles are considered.

The mean flow in a turbulent boundary layer (TBL) deviates from the canonical law of the wall (LoW) when influenced by a pressure gradient. Consequently, LoW-based near-wall treatments are inadequate for such flows. Chen et al. (J. Fluid Mech., vol. 970, 2023, A3) derived a Navier–Stokes-based velocity transformation that accurately describes the mean flow in TBLs with arbitrary pressure gradients. However, this transformation requires information on total shear stress, which is not always readily available, limiting its predictive power. In this work, we invert the transformation and develop a predictive near-wall model. Our model includes an additional transport equation that tracks the Lagrangian integration of the total shear stress. Particularly noteworthy is that the model introduces no new parameters and requires no calibration. We validate the developed model against experimental and computational data in the literature, and the results are favourable. Furthermore, we compare our model with equilibrium models. These equilibrium models inevitably fail when there are strong pressure gradients, but they prove to be sufficient for boundary layers subjected to weak, moderate and even moderately high pressure gradients. These results compel us to conclude that history effects in mean flow, which negatively impact the validity of equilibrium models, can largely be accounted for by the material time derivative term and the pressure gradient term, both of which require no additional modelling.

The mobility of externally driven phoretic propulsion of particles is evaluated by simultaneously solving the solute conservation equation, interaction potential equation and the modified Stokes equation. While accurate, this approach is cumbersome, especially when the interaction potential decays slowly compared with the particle size. In contrast to external phoresis, the motion of self-phoretic particles is typically estimated by relating the translation and rotation velocities with the local slip velocity. While this approach is convenient and thus widely used, it is only valid when the interaction decay length is significantly smaller than the particle size. Here, by taking inspiration from Brady (J. Fluid Mech., vol. 922, 2021, A10), which combines the benefits of two approaches, we reproduce their unified mobility expressions with arbitrary interaction potentials and show that these expressions can conveniently recover the well-known mobility relationships of external electrophoresis and diffusiophoresis for arbitrary double-layer thickness. Additionally, we show that for a spherical microswimmer, the derived expressions relax to the slip velocity calculations in the limit of the thin interaction length scales. We also employ the derived mobility expressions to calculate the velocities of an autophoretic Janus particle. We find that there is significant dampening in the translation velocity even when the interaction length is an order of magnitude larger than the particle size. Finally, we study the motion of a catalytically self-propelled particle, while it also propels due to external concentration gradients, and demonstrate how the two propulsion modes compete with each other.