Flow between axially rotating concentric cylinders is well known to exhibit rich dynamics. Hence, Taylor instabilities have been studied, both experimentally and theoretically, for many years. Although usually studied in the abstract, such geometries arise in a range of practical situations including drilling, when a drilling fluid flow enters a well via a pipe that is the centre body and returns via the annulus between the pipe and the borehole wall. In drilling, the centre body rotates and the annular flow contains rock cuttings. Here, we report the development of an Eulerian-Eulerian solver, based on OpenFOAM, that solves for this cuttings transport problem in the presence of both gravity and Taylor vortices. To check the reliability of the solver, we conduct a set of experiments spanning a wide range of complex flow regimes. We show that the model successfully predicts, in all regimes, the observed complex redistribution of particulates. However, for suspension flows under viscously dominated conditions, high particle concentrations and in rectilinear flow, particle pressure and normal stress differences are sufficient to capture particle migration. Results show that in more complex flows exemplified by the Taylor–Couette flow studied here, more realistic predictions of non-Brownian particle migration require inclusion of forces arising through the relative velocity of the two phases including lift forces originating both from inertia and particle rotation.

Large eddy simulations (LES) are widely used to study the effects of surface morphology on turbulence statistics, exchange processes and turbulence topology in urban canopies. However, as LES are only approximations of reality, special attention is needed for the computational model set-up to ensure an accurate representation of the physical processes of interest. This paper shows that the choice of the numerical domain can significantly affect the accuracy of turbulent flow statistics, potentially causing a mismatch between numerical studies and experimental data. The study examines the influence of cross-stream aspect ratio (YAR), streamwise aspect ratio (XAR) and scale separation (SS) on first- and second-order flow statistics and turbulence topology. It is found that domains with a low YAR underestimate the velocity variance, while those with a low XAR overestimate the variance value. The study proposes a new approach based on the Buckingham Pi theorem to evaluate the effect of SS, as the existing method has major limitations for canopy flows. The results suggest that domains with small SS underpredict the variance value. To minimise the artificial impact of the numerical domain on turbulent flow statistics, the study recommends guidelines for future research, including a YAR of 3 or more, an XAR of 6 or more and an SS of 12 or more. Error tables are presented to allow researchers to select smaller domains than recommended, depending on their research interests in specific parts of the flow.

Large-eddy simulations of a flat-plate boundary layer, without a leading edge, subject to multiple levels of incoming free-stream turbulence are considered in the present work. Within an input–output model, where nonlinear terms of the incompressible Navier–Stokes equations are treated as an external forcing, we manage to separate inputs related to perturbations coming through the intake of the numerical domain, whose evolution represents a linear mechanism, and the volumetric nonlinear forcing due to triadic interactions. With these, we perform the full reconstruction of the statistics of the flow, as measured in the simulations, to quantify pairs of wavenumbers and frequencies more affected by either linear or nonlinear receptivity mechanisms. Inside the boundary layer, different wavenumbers at near-zero frequency reveal streaky structures. Those that are amplified predominantly via linear interactions with the incoming vorticity occur upstream and display transient growth, while those generated by the nonlinear forcing are the most energetic and appear in more downstream positions. The latter feature vortices growing proportionally to the laminar boundary layer thickness, along with a velocity profile that agrees with the optimal amplification obtained by linear transient growth theory. The numerical approach presented is general and could potentially be extended to any simulation for which receptivity to incoming perturbations needs to be assessed.

There are numerous examples of long-lived, surface-intensified anticyclones over submarine depressions and troughs in the ocean. These often co-exist with a large-scale cyclonic circulation. The latter is predicted by existing barotropic theory but the anticyclone is not. We extend one such theory, which minimizes enstrophy while conserving energy, to two fluid layers. This yields a bottom-intensified flow with cyclonic circulation over a depression. The solution is steady, an enstrophy minimum and stable. When the Lagrange multiplier, $\lambda$, is near zero, the total potential vorticity (PV) becomes homogenized, in both layers. For positive $\lambda$, the surface PV is anticyclonic and strongest at intermediate energies. In quasi-geostrophic numerical simulations with a random initial perturbation PV, the bottom-intensified cyclonic flow always emerges. Vortices evolve independently in the layers and vortex mergers are asymmetric over the depression; cyclones are preferentially strained out at depth while only anticyclones merge at the surface. Both asymmetries are linked to the topographic flow. The deep cyclones feed the bottom-intensified cyclonic circulation while the asymmetry at the surface is only apparent after that circulation has spun up. The result of the surface merger asymmetry is often a lone anticyclone above the depression. This occurs primarily at intermediate energies, when the surface PV predicted by the theory is strongest. Similar results obtain in a full complexity ocean model but with a more pronounced asymmetry in surface vortex mergers and, with bottom friction, significant bottom flow beneath the central anticyclone.

As presented in Annenkov & Shrira (Phys. Rev. Lett., vol. 102, 2009, 024502), when a surface gravity wave field is subjected to an abrupt perturbation of external forcing, its spectrum evolves on a ‘fast’ dynamic time scale of $O(\varepsilon ^{-2})$, with $\varepsilon$ a measure of wave steepness. This observation poses a challenge to wave turbulence theory that predicts an evolution with a kinetic time scale of $O(\varepsilon ^{-4})$. We revisit this unresolved problem by studying the same situation in the context of a one-dimensional Majda–McLaughlin–Tabak equation with gravity wave dispersion relation. Our results show that the kinetic and dynamic time scales can both be realised, with the former and latter occurring for weaker and stronger forcing perturbations, respectively. The transition between the two regimes corresponds to a critical forcing perturbation, with which the spectral evolution time scale drops to the same order as the linear wave period (of some representative mode). Such fast spectral evolution is mainly induced by a far-from-stationary state after a sufficiently strong forcing perturbation is applied. We further develop a set-based interaction analysis to show that the inertial-range modal evolution in the studied cases is dominated by their (mostly non-local) interactions with the low-wavenumber ‘condensate’ induced by the forcing perturbation. The results obtained in this work should be considered to provide significant insight into the original gravity wave problem.

The near wake of a small-scale wind turbine is investigated using particle image velocimetry experiments at different tip speed ratios ($\lambda$). The wind turbine model had a nacelle and a tower mimicking real-scale wind turbines. The near wake is found to be dominated by multiple coherent structures, including the tip vortices, distinct vortex sheddings from the nacelle and tower, and wake meandering. The merging of the tip vortices is found to be strongly dependent on $\lambda$. A convective length scale ($L_c$) related to the pitch of the tip vortices is defined that is shown to be a better length scale than turbine diameter ($D$) to demarcate the near wake from the far wake. The tower induced strong vertical asymmetry in the flow by destabilising the tip vortices and promoting mixing in the lower (below the nacelle) plane. The nacelle's shedding is found to be important in ‘seeding’ wake meandering, which, although not potent, exists close to the nacelle, and it becomes important only after a certain distance downstream ($x>3L_c$). A link between the ‘effective porosity’ of the turbine and $\lambda$ is established, and the strength and frequency of wake meandering are found to be dependent on $\lambda$. In fact, a decreasing trend of wake meandering frequency with $\lambda$ is observed, similar to vortex shedding from a porous plate at varying porosity. Such similarity upholds the notion of wake meandering being a global instability of the turbine, which can be considered as a ‘porous’ bluff body of diameter $D$.

The reflection of a rightward-moving oblique shock (RMOS) belonging to the first family, over an initially steady oblique shock wave (SOSW) produced by a wedge, is studied in this paper. To cover all possibilities, the problem is divided into a pre-shock reflection problem, for which the incident shock is assumed to reflect over the pre-interaction part of the SOSW, and a post-shock reflection problem, for which the incident shock is assumed to reflect over the post-interaction part. Such division, together with the definition of the equivalent problem defined on the reference frame co-moving with the nominal intersection point of the two shock waves, allows us to connect the reflection patterns with the six types of shock interference of Edney, which include type I–VI shock interferences depending on how an upstream oblique shock intersects a bow shock (types I and II are regular and Mach reflections of two shocks from the opposite sides; type III and type IV have two triple points or two Mach reflection configuration; type V and type VI are irregular and regular reflections of two shocks from the same side). We are thus able to identify all possible shock reflection types and find their transition conditions. Pre-shock reflection may yield IV, V and VI (of Edney's six types) shock interferences and post-shock reflection may yield I, II and III shock interferences. Pre- and post-shock reflections possibly occur at two different parts of the SOSW, and the complete reflection configuration may have one or both of them. Both transition condition study and numerical simulation are used to show how pre-shock reflection and post-shock reflection exist alone or coexist, leading to various types of combined pre-shock and post-shock reflections.

We revisit the problem of a solid sphere rising slowly in a rotating short container filled with a slightly viscous fluid, with emphasis on the drag force. The data of the classical experiments of Maxworthy (J. Fluid Mech., vol. 31, 1968, pp. 643–655) and recent experiments of Kozlov et al. (Fluids, vol. 8 (2), 2023, paper 49), and the available geostrophic and quasi-geostrophic theories, are subjected to a novel scrutiny by combined reprocessing and comparisons. The measured drag is, consistently, about 20 % lower than the geostrophic prediction (assuming that flow is dominated by the Ekman layers, while in the inviscid cores the Coriolis acceleration is supported by the pressure gradient). The major objective is the interpretation and improvement of the gap between data and predictions. We show that the data cover a small range of relevant parameters (in particular the Taylor number $T$ and the height ratio $H$ of cylinder to particle diameter) that precludes a thorough and reliable assessment of the theories. However, some useful insights and improvements can be derived. The hypothesis that the discrepancy between data and the geostrophic prediction is due to inertial effects (not sufficiently small Rossby number $Ro$ in the experiments) is dismissed. We show that the major reason for the discrepancy is the presence of relatively thick Stewartson layers about the cylinder (Taylor column) attached to the sphere. The $1/3$ layer displaces the boundary condition of the angular velocity ($\omega = 0$) outside the radius of the particle. This observation suggests a semi-empirical correction to the theoretical quasi-geostrophic predictions (which takes into account the Ekman layers and the $1/4$ Stewartson layers); the corrected drag is in fair agreement with the data. We demonstrate that the inertial terms are negligible for $Ro\,T^{1/2} <0.4$. We consider curve-fit approximations, and point out some persistent gaps of knowledge that require further experiments and simulations.

Three-dimensional simulations of the ordering of elastic capsule suspensions within planar Poiseuille flow channels are reported. The simulations utilize the immersed boundary method coupled with the lattice Boltzmann method to capture the complex flow-induced capsule deformations and hydrodynamic interactions within the suspensions. A parametric study is presented whereby the confinement ratio and the particle deformability are varied independently within a two-dimensional range relevant to this ordering phenomenon. The initial distribution of capsules is random, and the simulations evolve the system from a disordered state to an ordered one, while an order parameter that quantifies the fraction of capsules belonging to one-dimensional train assemblies is computed throughout time. A monotonic increase in ordering is observed with increasing deformability. However, an optimal confinement ratio is identified corresponding to a peak in the order parameter. This peak is attributed to the competition between increasing long-range capsule attractions and decreasing in-plane capsule density (with fixed volume fraction) as the confinement ratio increases. Simulations are also performed to understand how dispersity in capsule size and deformability impact the degree of ordering. It is shown that ordering is quite sensitive to dispersity in capsule size, and much less sensitive to dispersity in deformability. Overall, the results provide important insights for the design of microfluidic devices.

A computational study was performed of the transport of both non-adhesive and adhesive particles in a porous bed with a body-centred cubic (BCC) structure. Pore-scale simulation of the flow within the porous bed was achieved through combining the immersed boundary method and the lattice Boltzmann method. Particle transport is computed using an adhesive discrete-element method based on a multi-time-scale soft-sphere model. The fluid flow results are validated by comparison with experimental data for dimensionless permeability of flow in a porous bed of spheres. For computations with non-adhesive particles, the particles are observed to drift to the centre of ‘channels’ in the BCC array, within which most of the fluid flow occurs. The mechanism of this inward drift was found to be related to the phenomenon of oscillatory clustering, which is an inertial drift mechanism observed for particles in a corrugated channel. A measure for particle drift into these channels was developed, and the time rate of change of this measure was found to compare closely with an approximate theoretical prediction based on oscillatory clustering theory. The drift measure was observed to be limited at long time by hold-up of outlier particles caught in long-duration collisions with the fixed bed particles in regions of low fluid velocity magnitude. Simulations with adhesive particles exhibited marked increase in collision duration, as well as inhibition of the tendency to drift toward the flow channels due to adhesive hold-up.

This work studies two-dimensional fixed-flux Rayleigh–Bénard convection with periodic boundary conditions in both horizontal and vertical directions and analyses its dynamics using numerical continuation, secondary instability analysis and direct numerical simulation. The fixed-flux constraint leads to time-independent elevator modes with a well-defined amplitude. Secondary instability of these modes leads to tilted elevator modes accompanied by horizontal shear flow. For $Pr=1$, where $Pr$ is the Prandtl number, a subsequent subcritical Hopf bifurcation leads to hysteresis behaviour between this state and a time-dependent direction-reversing state, followed by a global bifurcation leading to modulated travelling waves without flow reversal. Single-mode equations reproduce this moderate Rayleigh number behaviour well. At high Rayleigh numbers, chaotic behaviour dominated by modulated travelling waves appears. These transitions are characteristic of high wavenumber elevator modes since the vertical wavenumber of the secondary instability is linearly proportional to the horizontal wavenumber of the elevator mode. At a low $Pr$, relaxation oscillations between the conduction state and the elevator mode appear, followed by quasi-periodic and chaotic behaviour as the Rayleigh number increases. In the high $Pr$ regime, the large-scale shear weakens, and the flow shows bursting behaviour that can lead to significantly increased heat transport or even intermittent stable stratification.

A computational analysis is performed to study the three-dimensional response of rectangular shear layers to plasma actuator-based control, in the context of sound mitigation of supersonic non-axisymmetric jets. A Mach $1.5$ rectangular jet with an aspect ratio $2:1$ is controlled using experimentally informed actuation patterns, referred to as M0, M1, M2, M3, M${\rm \pi}$ and M+/$-$1. While the first five progressively increase the phase difference between successive actuators thus enhancing three-dimensionality of the shear layer structures, the latter corresponds to the flapping mode of the jet. A preliminary linear analysis identifies that the frequency, $St\sim 1$, has a relatively high overall amplification within the baseline shear layer, and is hence utilized for control in the subsequent nonlinear simulations. Each actuation reveals unique near-field vortical and acoustic responses that have a profound impact on far-field noise levels. The M0 actuation induces circumferentially interconnected strong streamwise vortices, while M1 actuation enhances the circumferential variability in the coherent structures. The M2 actuation encompasses both these effects, and along with a very low tonal impact of forcing, produces the most desirable far-field noise mitigation (${\sim }2.6\,{\rm dB}$), contributed by a broadband reduction around the column-mode peak of the baseline jet. Beyond M2 actuation, effectiveness of control saturates, particularly along the direction of peak noise radiation. Through a near-field analysis of the acoustic component, the efficacy of M2 actuation is attributed to the attenuation of the radiative efficiency of the jet, including reduced energy in the supersonic phase speeds, and redistribution of energy into the higher helical modes. Further, it curtails the nonlinear difference interactions in the plume that energize column-mode frequencies, which often appear as strong intermittent sound-producing events. While the shear layer turbulent kinetic energy decreases with actuation, the controlled jets show minimal variations in mean flow properties, particularly under M2 actuation, suggesting this to be a promising small-perturbation-based noise control strategy.

Based on the multiple-scale expansion technique, a new set of extended nonlinear Schrödinger (ENLS) equations up to the third order is derived to account for the additional high-order bottom and dispersion effects as well as the nonlinear wave interaction on wave transformation over periodic sandbars of sinusoidal geometry. By employing the small-amplitude wave assumption, a closed-form analytical solution for Bragg scattering is obtained from the linearised ENLS equations, which demonstrates that a downshift of wave frequency of the maximum reflection is mainly due to the inclusion of the high-order bottom effect. The factors that affect the downshift of the resonant frequency are identified and a theoretical expression in parabolic form is derived to quantify the downshift magnitude. The fully ENLS equations are further analysed to reveal the additional wave nonlinear effects on Bragg scattering characteristics. Under the condition of infinitesimal sandbar amplitude, the ENLS equations render a theoretical expression of the critical value of $kh$ when the nonlinear wave self-modulation effect and the nonlinear wave cross-modulation effect are equal, whereas the former effect is responsible for wavenumber upshifting and the latter downshifting. When $kh$ is larger than the critical value, the increase of wave nonlinearity will enhance the downshift magnitude of the Bragg resonance, and vice versa. For finite amplitude of the bottom sandbar, the ENLS equations are solved numerically to examine the influence of both wave nonlinearity and sandbar amplitude on the characteristics of Bragg resonance. The results reveal that as the increase of sandbar amplitude, the critical $kh$ increases monotonically.

A computational model is developed to study the time-averaged mean dynamics of red blood cells (RBCs) driven by the time-averaged mean stress generated by two phase-shifted orthogonal ultrasonic standing waves in a viscous fluid. The cell is modelled as an ellipsoidal viscoelastic membrane enclosing the viscous fluid cytoplasm, the motion of which is described by the inclination angle of the ellipsoidal cell shape and the phase angle of the potential membrane cycle. Based on the acoustic perturbation method, the acoustic field and acoustic streaming field are solved to obtain the time-averaged mean stress, and then the temporal evolution equations of the inclination and phase angles of the cell are determined considering the torque balance and energy conservation. At a small acoustic pressure amplitude, this model reproduces the experimentally observed features of cell motion in orthogonal standing waves: the transition from steady stationary orientation to unsteady tumbling with the increase of the phase difference between the two standing waves. By turning up the acoustic pressure amplitude above a critical value, it is further predicted that the previously observed motions can be accompanied by the membrane tank-treading rotation. Observations of these motions, combined with the present computational model, can help to evaluate the mechanical properties of RBC membranes in an automated and high-throughput manner by acoustic methods.

Binary nanodroplet collisions have received increasing attention, whilst the identification of collision outcomes and the viscous dissipation mechanism have remained poorly understood. Using molecular dynamics simulations, this study investigates binary nanodroplet collisions over wide ranges of Weber number (We), Ohnesorge number (Oh) and off-centre distances. Coalescence, stretching separation and shattering are identified; however, bouncing, reflexive separation and rotational separation reported for millimetre-sized collisions are not observed, which is attributed to the enhanced viscous effect caused by the ‘natural’ high-viscosity characteristics of nanodroplets. Intriguingly, as an intermediate outcome, holes form in retracting films at relatively high We, arising from the vibration and thermal fluctuation of the films. Due to the combined effects of inertial, capillary and viscous forces, binary nanodroplet collisions fall into the cross-over regime, so estimating viscous dissipation becomes extremely important for distinguishing outcome boundaries. Based on the criterion that stretching separation is triggered only when the residual off-centre kinetic energy exceeds the surface energy required for separation, the boundary equation between coalescence and stretching separation is established. Here, viscous dissipation is calculated by the extracted flow feature from simulations, showing that the ratio of viscous dissipation to the initial kinetic energy depends only on Oh, not on We. Because of complex viscous dissipation mechanisms, the same boundary equation in the cross-over regime has also not been satisfactorily revealed for macroscale collisions. Therefore, the proposed equation is tested for wide data sources from both macroscale and nanoscale collisions, and satisfying agreement is achieved, demonstrating the universality of the equation.

A form of skin-friction drag decomposition is given based on the velocity–vorticity correlations, $\langle v\omega _z\rangle$ and $\langle -w\omega _y\rangle$, which represent the advective vorticity transport and vortex stretching, respectively. This identity provides a perspective to understand the mechanism of skin-friction drag generation from vortical motions and it has better physical interpretability compared with some previous studies. The skin-friction coefficients in incompressible turbulent channel flows at friction Reynolds numbers from 186 to 2003 are divided with this velocity–vorticity correlation-based identity. We mainly focus on the Reynolds number effects on the contributing terms, their scale-dependence and quadrant characteristics. Results show that the contributing terms and their proportions exhibit similarities and the same peak locations across the wall layer. For the first time, we find that the positive and negative regions in the spanwise pre-multiplied spectra of the turbulent inertia ($\langle v'\omega _z'\rangle +\langle -w'\omega _y'\rangle$) can be separated with a universal linear relationship of $\lambda _z^+=3.75y^+$. The linear relationship is adopted as the criterion to investigate the scale dependence of the velocity–vorticity coupling structures. It reveals that the negative and positive structures dominate the generation of friction drag associated with the advective vorticity transport and vortex stretching, respectively. Moreover, quadrant analyses of the velocity–vorticity correlations are performed to further examine the friction drag generation related to different quadrant motions.

By examining a systematic set of direct numerical simulations, we develop a model which captures the effect of roughness density on global and local heat transfer in forced convection. The surfaces considered are zero-skewed three-dimensional sinusoidal rough walls with solidities, $\varLambda$ (defined as the frontal area divided by the total plan area), ranging from low $\varLambda = 0.09$, medium $\varLambda = 0.18$ to high $\varLambda = 0.36$. For each solidity, we vary the roughness height characterised by the roughness Reynolds number, $k^+$, from transitionally rough to fully rough conditions. The findings indicate that, as the fully rough regime is approached, there is a pronounced breakdown in the analogy between heat and momentum transfer, whereby the velocity roughness function $\Delta U^+$ continues to increase and the temperature roughness function $\Delta \varTheta ^+$ attains a peak with increasing $k^+$. This breakdown occurs at higher sand-grain roughness Reynolds numbers ($k_s^+$) with increasing solidity. Locally, we find that the heat transfer can be meaningfully partitioned into two categories: exposed, high-shear regions experiencing higher heat transfer obeying a local Reynolds analogy and sheltered, reversed-flow regions experiencing lower and spatially uniform heat transfer. The relative contribution of these distinct mechanisms to the global heat transfer depends on the fraction of the total surface area covered by these regions, which ultimately depends on $\varLambda$. These insights enable us to develop a model for the rough-wall heat-transfer coefficient, ${C_{h,k}(k^+, \varLambda, Pr)}$, where $Pr$ is the molecular Prandtl number, that assumes different heat-transfer laws in exposed and sheltered regions. We show that the exposed–sheltered surface-area fractions can be modelled through simple ray tracing that is solely dependent on the surface topography and a prescribed sheltering angle. Model predictions compare well when applied to heat-transfer data of traverse ribs from the literature.