We performed a linear stability analysis of a Newtonian ferrofluid cylinder surrounded by a Newtonian non-magnetic fluid in an azimuthal magnetic field. A wire is used at the centre of the ferrofluid cylinder to create this magnetic field. Isothermal conditions are considered and gravity is ignored. An axisymmetric perturbation is imposed at the interface between the two fluids and a dispersion relation is obtained allowing us to predict whether the flow is stable or unstable with respect to this perturbation. This relation is dependent on the Ohnesorge number of the ferrofluid, the dynamic viscosity ratio, the density ratio, the magnetic Bond number, the relative magnetic permeability and the dimensionless wire radius. Solutions to this dispersion relation are compared with experimental data from Arkhipenko et al. (Fluid Dyn., vol. 15, issue 4, 1981, pp. 477–481) and, more recently, Bourdin et al. (Phys. Rev. Lett., vol. 104, issue 9, 2010, 094502). A better agreement than the inviscid theory and the theory that only takes into account the viscosity of the ferrofluid is shown with the data of Arkhipenko et al. (Fluid Dyn., vol. 15, issue 4, 1981, pp. 477–481) and those of Bourdin et al. (Phys. Rev. Lett., vol. 104, issue 9, 2010, 094502) for small wavenumbers.

In this article, the objective is to characterize the influence of a continuous stratification on wake-induced vibrations of a circular cylinder. Experimental results are obtained by towing a horizontal cylinder in a horizontal direction perpendicular to its axes, at a constant speed in a linearly stratified fluid made of salty water. The cylinder can move vertically since it is fixed to free-to-rotate arms. The diameter of the cylinder, the length of the arms and the translating speed are varied. Two flow-induced vibrations modes are observed. The first one is analogous to a vortex-induced vibration (VIV) mode: it is associated with a resonance between the vortex shedding frequency and its frequency is proportional to the natural frequency of the mechanical system, as can be predicted by the VIV theory generalized to low mass ratios. Large amplitude oscillations of the cylinder are found to occur in a large range of velocities also in agreement with the low-mass-ratio VIV theory. The second mode is interpreted as a galloping mode. It has a low frequency and a large amplitude, and occurs for low Froude numbers and long arms. It corresponds to the regimes when the buoyancy forces are larger than the inertial and vertical drag forces. By computing the forces on the cylinder, it is shown that stratification is the source of a destabilizing lift when the cylinder departs from its horizontal motion. Only a weak effect of the Reynolds number on the stability characteristics has been observed in the considered range ($500 < Re < 15\,000$).

We consider a uniform ellipsoid of potential vorticity (PV), where we exploit analytical solutions derived for a balanced model at the second order in the Rossby number, the next order to quasi-geostrophic (QG) theory, the so-called QG+1 model. We consider this vortex in the presence of an external background shear flow, acting as a proxy for the effect of external vortices. For the QG model the system depends on four parameters, the height-to-width aspect ratio of the vortex, $h/r$, as well as three parameters characterising the background flow, the strain rate, $\gamma$, the ratio of the background rotation rate to the strain, $\beta$, and the angle from which the flow is applied, $\theta$. However, the QG+1 model also depends on the PV, as well as the Prandtl ratio, $f/N$ ($f$ and $N$ are the Coriolis and buoyancy frequencies, respectively). For QG and QG+1 we determine equilibria for different values of the background flow parameters for increasing values of the imposed strain rate up to the critical strain rate, $\gamma _c$, beyond which equilibria do not exist. We also compute the linear stability of this vortex to second-order modes, determining the marginal strain $\gamma _m$ at which ellipsoidal instability erupts. The results show that for QG+1 the most resilient cyclonic ellipsoids are slightly prolate, while anticyclonic ellipsoids tend to be more oblate. The highest values of $\gamma _m$ occur as $\beta \to 1$. For large values of $f/N$, changes in the marginal strain rates occur, stabilising anticyclonic ellipsoids while destabilising cyclonic ellipsoids.

Surfactant transport from foam film to foam film is an essential (yet poorly understood) aspect of the viscoplastic yielding behaviour of flowing foam. Recent experimental and modelling work by Bussonnière & Cantat (J. Fluid Mech., vol. 922, 2021, A25) has, however, helped to advance understanding of the relevant surfactant transport processes: the significance of that work is described herein.

The flow in and around a pair of plumes from sources that are vertically and horizontally offset is investigated. An analytical potential flow model is developed using an adapted version of the Milne–Thomson circle theorem to represent the flow due to adjacent circular and point sinks. This approximates the horizontal section through offset plumes, which have differing radii at the same vertical position. The predictions of this model are compared against Reynolds-averaged Navier–Stokes (RANS) simulations of the same system. These yield time-averaged results and so are particularly suitable for investigating the relatively weak entrainment field. Single and double plume results from RANS are also presented to compare with known results. Good agreement is found between the features in the analytical model and in the numerical solutions including the location of the stagnation point between the two plumes.

The maximum size and lifetime of an acoustically nucleated cavitation bubble is inversely proportional to the driving frequency and has achieved a limit of approximately 10 MHz. Smaller cavitation bubbles that are critical to microscopic applications require shorter lifetimes that correspond to higher oscillation frequencies. Here, we demonstrate that acoustic cavitation in the 100 MHz range and beyond can be achieved through wave propagation in a solid rather than in a liquid. The cavitation bubble is nucleated at a nano-sized fracture on a glass substrate, and its expansion is driven by a leaky Rayleigh wave, while the inertial collapse is induced by a trailing shock wave. As both waves travel at different velocities, the time interval between these two events is a function of the distance to the source. In this way, we demonstrate experimentally control of the lifetime of the bubbles in a range between 6 and 80 ns, corresponding to oscillation frequencies between 13 and 166 MHz. Our results agree with finite-volume fluid–structure interaction simulations.

The analysis of turbulence by way of higher-order spectral moments is uncommon, despite the relatively frequent use of such statistical analyses in other fields of physics and engineering. In this work, higher-order spectral moments are used to investigate the internal intermittency of the turbulent velocity and passive-scalar (temperature) fields. This study first introduces the theory behind higher-order spectral moments as they pertain to the field of turbulence. Then, a short-time Fourier-transform-based method is developed to estimate these higher-order spectral moments and provide a relative, scale-by-scale measure of intermittency. Experimental data are subsequently analysed and consist of measurements of homogeneous, isotropic, high-Reynolds-number, passive and active grid turbulence over the Reynolds-number range $35\leq R_{\lambda } \leq ~731$. Emphasis is placed on third- and fourth-order spectral moments using the definitions formalised by Antoni (Mech. Syst. Signal Pr., vol. 20 (2), 2006, pp. 282–307), as such statistics are sensitive to transients and provide insight into deviations from Gaussian behaviour in grid turbulence. The higher-order spectral moments are also used to investigate the Reynolds (Péclet) number dependence of the internal intermittency of velocity and passive-scalar fields. The results demonstrate that the evolution of higher-order spectral moments with Reynolds number is strongly dependent on wavenumber. Finally, the relative levels of internal intermittency of the velocity and passive-scalar fields are compared and a higher level of internal intermittency in the inertial subrange of the scalar field is consistently observed, whereas a similar level of internal intermittency is observed for the velocity and passive-scalar fields for the high-Reynolds-number cases as the Kolmogorov length scale is approached.

A rigorous bound on the energy dissipation in the flow around a sphere is derived analytically. The bound is allowed to depend on the viscosity of the fluid, the radius of the sphere, the free-stream velocity and the amplitude of the difference between the actual flow field and the irrotational flow around the sphere. The bound becomes asymptotically independent of viscosity as the Reynolds number tends to infinity. Provided that the velocity fluctuations are of the same order as the free-stream velocity, the scaling of the bound is consistent with a Reynolds-number-independent drag coefficient at high Reynolds numbers.

Granular flows are highly dissipative due to frictional resistance and inelasticity in collisions among grains. They are known to exhibit shock waves at velocities that are easily achieved in industrial and nature-driven flows such as avalanches and landslides. This experimental work investigates the formation of strong shock waves on triangular obstacles placed in a dry rapid granular stream in a confined two-dimensional set-up. Oblique attached shock waves are formed for mild turning angles and higher flow velocities, whereas strong bow shock waves are formed for higher turning angles and slower granular streams. A shadowgraph imaging technique elucidates interesting characteristics of the shock waves, especially in the vicinity of shock detachment. Velocity distributions in the form of scatter plots and probability distribution functions are calculated from the flow field data obtained by particle imaging velocimetry. The flow field around the granular shock wave region represents a bimodal distribution of velocities with two distinct peaks, one representing the supersonic flow within the free stream, and the other corresponding to the subsonic faction downstream of a shock wave. Connecting the two is a population that does not directly belong to either of the modes, constituting the non-equilibrium shock wave region. The effect of grain size and scaling, for fixed free-stream conditions and fixed channel width, on the shock detachment is presented. The mechanisms of the static heap formation and the shock detachment process in a confined environment are discussed.

We conducted a three-dimensional numerical simulation of bioconvection generated by oxygen-reactive chemotactic bacteria. This study investigated the bioconvection patterns, interference between plumes, and the wavelength of bioconvection patterns. In addition, we clarified the transport characteristics of cells and oxygen in the bioconvection. Multiple plumes occur in the suspension, and three-dimensional bioconvection is formed around the plumes by the cells with vortex rings arising around the plumes. Even if bioconvection at a high Rayleigh number is disturbed, the bioconvection is strongly stable with respect to disturbances, and the pattern does not change due to disturbances. Bioconvection changes depending on the physical properties of bacteria and oxygen, and in particular, the rate of oxygen consumption by bacteria significantly affects the strength of bioconvection. Bioconvection patterns with different plume arrangements and shapes are formed for different Rayleigh numbers or initial disturbances of the cell concentration. As a result, the wavelengths of the patterns also vary. As the Rayleigh number increases, interference between plumes is strengthened by the shortening of the pattern wavelength, so the velocities of both upward and downward flows increase. Many cells are located under the plumes, and a strong shear flow occurs in these regions. As the pattern wavelength decreases, the cells are affected by high shear stress. Then the convective transport of the entire suspension strengthens and the transport characteristics of cells and oxygen improve. When the chamber boundary is changed to side walls, bacteria adhere to the wall surface, and plumes are arranged regularly along the side walls.

Bounding walls or immersed surfaces are utilized in many industrial systems as the primary thermal source to heat a gas–solids mixture. Previous efforts to resolve the solids’ heat transfer near a boundary involve the extension of unbounded convection correlations into the near-wall region in conjunction with particle-scale theories for indirect conduction. Moreover, unbounded drag correlations are utilized in the near-wall region (without modification) to resolve the force exerted on a solid particle by the fluid. We rigorously test unbounded correlations and indirect conduction theory against outputs from direct numerical simulation of laminar flow past a hot plate and a static, cold particle. Here, local variables are utilized for consistency with unresolved computational fluid dynamics discrete element methods and lead to new unbounded correlations that are self-similar to those obtained with free-stream variables. The new drag correlation with local fluid velocity captures the drag force in both the unbounded system as well as the near-wall region while the classic, unbounded drag correlation with free-stream fluid velocity dramatically over-predicts the drag force in the near-wall region. Similarly, classic, unbounded convection correlations are found to under-predict the heat transfer occurring in the near-wall region. Inclusion of indirect conduction, in addition to unbounded convection, performs markedly better. To account for boundary effects, a new Nusselt correlation is developed for the heat transfer in excess of local, unbounded convection. The excess wall Nusselt number depends solely on the dimensionless particle–wall separation distance and asymptotically decays to zero for large particle–wall separation distances, seaming together the unbounded and near-wall regions.

The stability and transition in the bottom boundary layer under a solitary wave are analysed in the presence of finite-amplitude disturbances. First, the receptivity of the boundary layer is investigated using a linear input-output analysis, in which the environment noise is modelled as distributed body forces. The most ‘dangerous’ perturbations in a time frame until flow reversal are found to be arranged as counter-rotating streamwise-constant vortices. One of these vortex configurations is then selected and deployed to nonlinear equations, and streaks of various amplitudes are generated via the lift-up mechanism. By means of secondary stability analysis and direct numerical simulations, the dual role of streaks in the boundary-layer transition is shown. When the amplitude of streaks remains moderate, these elongated features remain stable until the adverse-pressure-gradient stage and have a dampening effect on the instabilities developing thereafter. In contrast, when the low-speed streaks reach high amplitudes exceeding 15 % of the free stream velocity at the respective phase, they become highly unstable to secondary sinuous modes in the outer shear layers. Consequently, a subcritical transition to turbulence, i.e. bypass transition, can be initiated already in the favourable-pressure-gradient region ahead of the wave crest.

The stability and dynamics of a falling liquid film over an anisotropic porous medium are studied using a one-domain approach. Our stability analysis shows a significant departure from the effective no-slip boundary condition in the isotropic case. Anisotropy does not affect the threshold of linear instability. However, a non-trivial dual effect of anisotropy on the film stability is observed depending on the permeability of the porous medium. This dual effect results from the net balance of the enhancement of viscous diffusion at the top Brinkman sublayer and the mitigation of viscous damping in the core Darcy sublayer. Three-equation models have been derived from the lubrication theory approximation in terms of the exact mass balance and averaged momentum balances in the porous and liquid layers. In the nonlinear regime, anisotropy has a dual effect by damping capillary waves at large permeabilities and enhancing them at low permeabilities. Anisotropy also affects wave speeds and shapes, modifies travelling-wave branches of solutions, affects the development of a time-periodic wavetrain by inlet forcing and alters the noise-driven dynamics of the flow. These effects result from the mitigation of mass exchange at the liquid–porous interface and the contribution of the cross-stream permeability in the Brinkman top sublayer to the viscous diffusion.

It is shown that thermal waves applied on the bounding surface of a horizontal slot generate a pumping effect. Reynolds stress developed by the change in the flow field due to the thermal inertia associated with the penetration of the wave into the fluid interior drives the net fluid movement in the horizontal direction. The induced flow rate increases with the wave speed, but excessive wave speed reduces it as convection becomes limited only in the near-wall area. The excessive increase of the wavelength and its excessive decrease reduce the flow rate. An increase in the wave amplitude increases the flow rate with a saturation limiting its growth. Judicious selection of the wave speed and wavelength of a thermal wave provides a means for optimizing the wave-induced pumping. The magnitude of the pumping increases with a reduction of Prandtl number.

This work reports experimental observation and theoretical explanation of the dynamics and morphology of a droplet passing through a soap film. During the process, the film undergoes four sequential responses: (1) film deformation upon droplet impact; (2) drop–film detachment; (3) coalescence of the film shell with the drop; (4) peel-off of the film shell. Physical models and the corresponding analytical expressions are developed to reveal the underlying physics for the observed four responses. It is identified that the film is an elongated catenoid under continuous stretch by the droplet, and that they separate at the fixed height of 5.8 times of the droplet radius while the detach point is located at the centre of the height. After separation, the droplet is wrapped with a film shell, which is then punctured by the ring tip of the converging surface wave at the impacting Weber number range of [45, 225]. The film shell then coalesces with the droplet, falls off with a fixed velocity and is eventually ejected as a bubble leaving the droplet with a transplanted surface of the soap solution.

We investigate the streamwise vorticity generation mechanism and distribution pattern in an unbounded steady inertial flow past a fixed Platonic polyhedron. Three angular positions are selected: an edge facing the flow (E), a face facing the flow (F) and a vertex facing the flow (V). We provide compelling evidence that the generation of the streamwise vorticity is primarily caused by the tilting of the transverse vorticity that originates from the particle front surface. Each inclined face on the front surface generates a pair of opposite-signed streamwise vortices. They are advected to the particle wake and form a chiral vorticity pattern which preserves the symmetry of the particle front surface. Two particles at dual angular positions exhibit highly similar vorticity patterns. Our study reveals a striking similarity between the vorticity patterns and the far-field optics diffraction pattern of a light beam past a polygonal aperture. We discover the deterministic vorticity generation mechanism to predict the streamwise vorticity patterns based on the distribution of edges and inclined faces on the particle front surface. Conversely, the vorticity patterns themselves can serve as a diagnostic tool to infer the geometry of the opaque particle front surface. Additionally, the vorticity patterns can be used to predict the stable angular position of a freely settling angular particle, which tends to be such that the number of streamwise vorticity pairs in the wake is maximized.

In this study, we performed a set of experiments of the thermal plume on an open cylinder heated from below and a simple scaling analysis. The flow structure of the plume was visualized using the shadowgraph technique, and the temperature at specific points was measured using a thermistor. Transient plumes in a developing stage, an equilibrating stage and a fully developed stage were described. The new scaling laws of the stem radius and the velocity of the rising plume were presented, which are different from those on a two-dimensional heated plate due to the cylindrical effect and agree with experimental results. In the fully developed stage, there is a transition route of the plume from a steady to chaotic state with an increase in the Rayleigh number, which involves a series of bifurcations. A reverse bifurcation from a periodic (Ra = 1.17 × 106) back to a steady (Ra = 1.30 × 106) state has been observed in the experiment, which is referred to as the period bubbling bifurcation. Thus, the present experimental results validate the previous numerical results presented by Zhang et al. (Phys. Fluids, vol. 33 (6), 2021, 064110). The bifurcation diagram, spectral analysis and attractor are used to characterize the transitional plume. In addition, the heat and mass transfer have also been quantified.

This study examines the effects of surface topography on the flow and performance of a self-propelled swimming (SPS) body. We consider a thin flat plate with an egg-carton roughness texture undergoing prescribed undulatory swimming kinematics at Strouhal number $0.3$ and tail amplitude to length ratio $0.1$; we use plate Reynolds numbers $Re=6$, 12 and $24\times 10^3$, and focus on $12\,000$. As the roughness wavelength is decreased, we find that the undulation wave speed must be increased to overcome the additional drag from the roughness and maintain SPS. Correspondingly, the extra wave speed raises the power required to maintain SPS, making the swimmer less efficient. To decouple the roughness and the kinematics, we compare the rough plates to equivalent smooth cases by matching the kinematic conditions. We find that all but the longest roughness wavelengths reduce the required swimming power and the unsteady amplitude of the forces when compared to a smooth plate undergoing identical kinematics. Additionally, roughness can enhance flow enstrophy by up to 116 % compared to the smooth cases without a corresponding spike in forces; this suggests that the increased mixing is not due to increased vorticity production at the wall. Instead, the enstrophy is found to peak strongly when the roughness wavelength is approximately twice the boundary layer thickness over the $Re$ range, indicating the roughness induces large-scale secondary flow structures that extend to the edge of the boundary layer. This study reveals the nonlinear interaction between roughness and kinematics beyond a simple increase or decrease in drag, illustrating that roughness studies on static shapes do not transfer directly to unsteady swimmers.