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

Astrophysical flows are often subject to both rotation and large-scale background magnetic fields. Individually, each is known to two-dimensionalize the flow in the perpendicular plane. In realistic settings, both of these effects are simultaneously present and, importantly, need not be aligned. In this work, we numerically investigate three-dimensional forced magnetohydrodynamic turbulence subject to the competing effects of global rotation and a perpendicular background magnetic field. We focus on the case of a strong background field and find that increasing the rotation rate from zero produces significant changes in the structure of the turbulent flow. Starting with a two-dimensional inverse energy cascade at zero rotation, the flow first transitions to a forward cascade of kinetic energy, then to a shear-layer dominated regime and finally to a second shear-layer regime where the kinetic energy flux is strongly suppressed and the energy transfer is mediated by the induced magnetic field. We show that the first two transitions occur at distinct values of the Rossby number, and the third occurs at a distinct value of the Lehnert number. The three-dimensional results are confirmed using an asymptotic two-dimensional, three-component model, which allows us to extend our results to the planetary-relevant case of an arbitrary angle between the rotation vector and guide field. More generally, our results demonstrate that, when considering the simultaneous limits of strong rotation and a strong guide field, the order in which those limits are taken matters in the misaligned case.

In an experimental context, the contamination of an air–liquid interface by ambient pollutants can strongly affect the dynamics and the stability of a given flow. In some configurations, the interfacial flow can even be blocked by surface tension effects. A cylindrical free-surface flow driven by a slow rotating disc is considered here as a generic example of such effects and is investigated both experimentally and numerically. We suggest here a simple numerical model, without any superficial transport of the pollutants, adaptable into any code for single-phase flows. For the stationary axisymmetric base flow, the radial velocity at the interface is set to zero whereas the usual stress-free boundary conditions are retained for the perturbations. The model does not feature any free parameter. For a geometrical aspect ratio of 1/4, known to display ambiguous behaviour regarding stability thresholds, the modal selection as well as a nonlinear stability island found in the experiments are well reproduced by the model, both qualitatively and quantitatively. The robustness of the model has also been validated by replacing the radial velocity profile by a more accurate experimental fit, with very little influence on the stability results.

The study of free-surface flows over vegetative structures presents a challenging setting for theoretical, computational and experimental analysis. In this work, we develop a multiple-scales asymptotic framework for the evolution of free-surface waves over rigid vegetation and a slowly varying substrate. The analysis quantifies the balance between the competing effects of vegetation and shoaling, and provides a prediction of the amplitude as the wave approaches a coastline. Our analysis unifies and extends existing theories that study these effects individually. The asymptotic predictions are shown to provide good agreement with full numerical simulations (varying depth) and published experimental results (constant depth).

The effect of the slip velocity in a bubbly jet ($Re = 3000$, initial void $\textrm {fraction} = 1.1\,\%$) is studied with various bubble slip ratios (0.07–2.1 times the inlet velocity, 0.01–0.31 m s$^{-1}$ in dimension). The governing equations, which are the conservation of mass, momentum and volume fraction, are solved by direct numerical simulations. A set of ordinary differential equations was derived by using a conventional one-dimensional integral framework on jets and by expressing the slip velocity using an exponential function. The one-dimensional analysis of bubbly jets successfully predicted the jet velocity radius, centreline velocity and the gas spread rate downstream of the bubbly jet. Second- and third-order statistics were also analysed to better understand the turbulent characteristics of the bubbly jet. A high slip ratio results in a rigid, narrower bubbly jet core region, where the turbulent kinetic energy is conserved along the centreline but does not diffuse towards the ambient region. The narrow, circular region around the core at high slip ratios decreased the turbulent kinetic energy. In the rigid bubbly jet core, turbulent diffusion of the gas phase is suppressed and has low correlations with the velocity fluctuations. Turbulence characteristics such as turbulence stress were compared with single-phase jet and plume results from existing literature. The magnitude of the turbulence characteristics at the lowest slip ratio is comparable with what is found in the literature, but there is a rapid transition from slip ratios of $0.14$ to $0.7$. Furthermore, the turbulent kinetic energy budget was analysed for each case. The production term had the largest contribution, and its magnitude was almost five times the buoyancy production downstream. High slip ratio cases show a positive peak of pressure strain and turbulent diffusion along the centreline.

Axisymmetric standing waves occur across a wide range of free surface flows. When these waves reach a critical height (steepness), wave breaking and jet formation occur. For travelling surface gravity waves, wave breaking is generally considered to limit wave height and reversible wave motion. In the ocean, the behaviour of directionally spread waves lies between the limits of purely travelling (two dimensions) and axisymmetric (three dimensions). Hence, understanding wave breaking and jet formation on axisymmetric surface gravity waves is an important step in understanding extreme and breaking waves in the ocean. We examine an example of axisymmetric wave breaking and jet formation colloquially known as the ‘spike wave’, created in the FloWave circular wave tank at the University of Edinburgh, UK. We generate this spike wave with maximum crest amplitudes of 0.15–6.0 m (0.024–0.98 when made non-dimensional by characteristic radius), with wave breaking occurring for crest amplitudes greater than 1.0 m (0.16 non-dimensionalised). Unlike two-dimensional travelling waves, wave breaking does not limit maximum crest amplitude, and our measurements approximately follow the jet height scaling proposed by Ghabache et al. (J. Fluid Mech., vol. 761, 2014, pp. 206–219) for cavity collapse. The spike wave is predominantly created by linear dispersive focusing. A trough forms, then collapses producing a jet, which is sensitive to the trough's shape. The evolution of the jets that form in our experiments is predicted well by the hyperbolic jet model proposed by Longuet–Higgins (J. Fluid Mech., vol. 127, 1983, pp. 103–121), previously applied to jets forming on bubbles.

Cross-flow, or vertical-axis, turbines are a promising technology for capturing kinetic energy in wind or flowing water and their inherently unsteady fluid mechanics present unique opportunities for control optimization of individual rotors or arrays. To explore the potential for beneficial interactions between turbines in an array, as well as to characterize important cycle-to-cycle variations, coherent structures in the wake of a single two-bladed cross-flow turbine are examined using planar stereo particle image velocimetry in a water channel experiment. There are three main objectives in the present work. First, the mean wake structure of this high chord-to-radius ratio rotor is described, compared with previous studies, and a simple explanation for observed wake deflection is presented. Second, the unsteady flow is then analysed via the triple decomposition, with the periodic component extracted using a combination of traditional techniques and a novel implementation of the optimized dynamic mode decomposition. The latter method is shown to outperform conditional averaging and Fourier methods, as well as uncover frequencies suggesting a transition to bluff-body shedding in the far wake. Third, vorticity and finite-time Lyapunov exponents are then employed to further analyse the oscillatory wake component. Vortex streets on both sides of the wake are identified, and their formation mechanisms and effects on the mean flow are discussed. Strong axial (vertical) flow is observed in vortical structures shed on the retreating side of the rotor where the blades travel downstream. Time-resolved tracking of these vortices is performed, which demonstrates that vortex trajectories have significant rotation-to-rotation variation within one diameter downstream. This variability suggests it would be challenging to harness or avoid such structures at greater downstream distances.

We report on the mobility and orientation of finite-size, neutrally buoyant, prolate ellipsoids (of aspect ratio $\varLambda =4$) in Taylor–Couette flow, using interface-resolved numerical simulations. The set-up consists of a particle-laden flow between a rotating inner and a stationary outer cylinder. The flow regimes explored are the well-known Taylor vortex, wavy vortex and turbulent Taylor vortex flow regimes. We simulate two particle sizes $\ell /d=0.1$ and $\ell /d=0.2$, $\ell$ denoting the particle major axis and $d$ the gap width between the cylinders. The volume fractions are $0.01\,\%$ and $0.07\,\%$, respectively. The particles, which are initially randomly positioned, ultimately display characteristic spatial distributions which can be categorised into four modes. Modes (i) to (iii) are observed in the Taylor vortex flow regime, while mode (iv) encompasses both the wavy vortex and turbulent Taylor vortex flow regimes. Mode (i) corresponds to stable orbits away from the vortex cores. Remarkably, in a narrow $\textit {Ta}$ range, particles get trapped in the Taylor vortex cores (mode (ii)). Mode (iii) is the transition when both modes (i) and (ii) are observed. For mode (iv), particles distribute throughout the domain due to flow instabilities. All four modes show characteristic orientational statistics. The focus of the present study is on mode (ii). We find the particle clustering for this mode to be size-dependent, with two main observations. Firstly, particle agglomeration at the core is much higher for $\ell /d=0.2$ compared with $\ell /d=0.1$. Secondly, the $\textit {Ta}$ range for which clustering is observed depends on the particle size. For this mode (ii) we observe particles to align strongly with the local cylinder tangent. The most pronounced particle alignment is observed for $\ell /d=0.2$ at around $\textit {Ta}=4.2\times 10^5$. This observation is found to closely correspond to a minimum of axial vorticity at the Taylor vortex core ($\textit {Ta}=6\times 10^5$) and we explain why.

Hypersonic laminar flow over a canonical 25°–55° double cone is studied using computational fluid dynamics and global stability analysis (GSA) with a free-stream Mach number of 11.5 and various unit Reynolds numbers. Axisymmetric simulations reveal that secondary separation occurs beneath the primary separation bubble beyond a critical Reynolds number. The numerical results agree well with existing experiments and the triple-deck theory with the axisymmetric effect on the incoming boundary layer treated by the Mangler transformation. The GSA identifies a three-dimensional global instability that is azimuthally periodic immediately prior to the emergence of secondary separation. The criterion of the onset of global instability in terms of a scaled deflection angle established for supersonic compression corner flows (Hao et al., J. Fluid Mech., vol. 919, 2021, A4) can be directly applied to double-cone flows. As the Reynolds number is further increased, the flow is strongly destabilized with the coexistence of multiple stationary and low-frequency oscillating unstable modes. Direct numerical simulations confirm that the supercritical double-cone flow is intrinsically three-dimensional, unsteady and exhibits strong azimuthal variations in the peak heating.

Rotating spheres in cylindrical channels roll or slide along the channel depending on the physical and geometric conditions. For a thorough investigation of the phenomenon, finite-element modelling is utilized to obtain the resistance coefficients for the motion of a sphere in a cylindrical channel, with an emphasis on near-wall motion. Extracted coefficients are compared with the data in the literature and utilized in exploring the conditions for rolling versus sliding along the channel. Sliding occurs due to the pressure build-up in the nip region between the wall and the rotating sphere in small confinement ratios, whereas rolling occurs when the shearing forces on the sphere are dominant in larger ratios. According to numerical results, a flow reversal takes place in the nip region and reduces the shear as well. Rolling versus sliding is demonstrated in experiments by using magnetic spherical particles, which are rotated by means of an external magnetic field inside cylindrical channels filled with viscous fluids. Faster axial velocities are observed in narrow channels while sliding than in wider channels while rolling for the same rotation rate of the sphere. Experiment observations are compared with the velocities evaluated from the resistance coefficients, showing that the distance between the sphere and the wall, which is dominated by roughness, plays an important role in the velocity of the sphere.

We propose a more conservative, physically-intuitive criterion, namely, the boundary enstrophy flux ($BEF$), to characterise leading-edge-type dynamic stall onset in incompressible flows. Our results are based on wall-resolved large-eddy simulations of pitching aerofoils, with fine spatial and temporal resolution around stall onset. We observe that $|BEF|$ reaches a maximum within the stall onset regime identified. By decomposing the contribution to $BEF$ from the flow field, we find that the dominant contribution arises from the laminar leading edge region, due to the combined effect of large clockwise vorticity and favourable pressure gradient. A relatively small contribution originates from the transitional/turbulent laminar separation bubble (LSB) region, due to LSB-induced counter-clockwise vorticity and adverse pressure gradient. This results in $BEF$ being nearly independent of the integration length as long as the region very close to the leading edge is included. This characteristic of $BEF$ yields a major advantage in that the effect of partial or complete inclusion of the noisy LSB region can be filtered out, without changing the $BEF$ peak location in time significantly. Next, we analytically relate $BEF$ to the net wall shear and show that its critical value ($=\max (|BEF|)$) corresponds to the instant of maximum net shear prevailing at the wall. Finally, we have also compared $BEF$ with the leading edge suction parameter ($LESP$) (Ramesh et al., J. Fluid Mech., vol. 751, 2014, pp. 500–538) and find that the former reaches its maximum value between $0.3^{\circ }$ and $0.8^{\circ }$ of rotation earlier.

We investigated the propagation of turbulent fronts in pipe flow at high Reynolds numbers by direct numerical simulation. We used a technique combining a moving frame of reference and an artificial damping to isolate the fronts in short periodic pipes, which enabled us to explore the bulk Reynolds number up to $Re=10^5$ with affordable computation power. We measured the propagation speed of the downstream front and observed that a fit of $1.971-(Re/1925)^{-0.825}$ (in unit of bulk speed) captures this speed above $Re\simeq 5000$ very well. The speed increases monotonically as $Re$ increases, in stark contrast to the decreasing trend above $Re\simeq 10\,000$ reported by Wygnanski & Champagne (J. Fluid Mech., vol. 59, 1973, pp. 281–335). The speed of the upstream front overall agrees with the former studies and $0.024+(Re/1936)^{-0.528}$ fits our data well, and those from the literature. Based on our analysis of the front dynamics, we proposed that both front speeds would keep their respective monotonic trends as the Reynolds number increases further. We show that, at high Reynolds numbers, the local transition at the upstream front tip is via high-azimuthal-wavenumber structures in the high-shear region near the pipe wall, whereas at the downstream front tip is via low-azimuthal-wavenumber structures in the low-shear region near the pipe centre. This difference is possibly responsible for the asymmetric speed scalings between the upstream and downstream fronts.

This study presents computational simulations of multicomponent and multiphase flows to reproduce the physical phenomena in the secondary atomization of a droplet induced by a hot temperature environment. The computational fluid dynamics model is based on the geometric volume of fluid method, with piecewise linear interface calculation reconstruction for accurate determination of the curvature and evaporation fluxes at the interface. The purpose of the model was to faithfully reproduce complex physical processes, such as internal gas cavity formation, liquid–vapour interface instability, cavity collapse and liquid jet ejection, and the pinch-off of a secondary droplet, leading to the microexplosion phenomenon that greatly enhances the evaporation rate of non-volatile liquid droplets. The solver was validated against the analytical solution in benchmark cases, and experimental data with bicomponent droplets reported in the literature. The developed model was used to predict the atomization of heavy fuel oil exposed at high temperatures under microgravity conditions. Different atomization regimes were identified, depending on the initial size of the internal bubbles. While small bubbles led to simple gas ejections, cavity collapse caused the larger bubbles to produce a jet formation. When the ratio between the bubble and droplet volumes was bigger than 0.7, microexplosions occurred. The results were found to be consistent with cases of bubble burst on flat surfaces, showing a strong dependence on the Ohnesorge number ($Oh$). Key observable quantities, particularly jet velocity and bubble cap drainage velocity, were found to agree with correlations reported in other studies. The similarities were also supported by studies extending over a wide range of simulations (4000 cases) at different $Oh$. An inversion in the dependence of the jet velocity on $Oh$ (above a critical value $Oh_c$) was observed.

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 consider miscible displacements in two-dimensional homogeneous porous media where the displacing fluid is less viscous and has a different density than the displaced fluid. We find that the dynamics evolve through nine possible regimes depending on the viscosity ratio, strength of density variations and the strength of the background flow, as characterized by the Péclet number. At early times the interface is dominated by longitudinal diffusion before undergoing a transition to a slumping regime where vertical flow is important. At intermediate times, vertical flow and diffusion can be neglected and there are three different limiting solutions: a fingering limit; an injection-driven gravity-current limit; and a density-driven gravity-current limit. Finally at late times, transverse diffusion becomes important and there is a transition from an apparent shutdown regime to a viscously enhanced Taylor-slumping regime. In each of the regimes, the dominant scalings are identified and reduced-order models for the evolution of the concentration field are developed. Lastly, three case studies are considered to illustrate the dominant physical balances in the geophysically relevant setting of geological $\textrm {CO}_2$ storage.

We study the dispersion characteristics of slightly buoyant droplets in a turbulent jet using large eddy simulations (LES). The droplet number density fields are represented using an Eulerian approach, with the dispersed phase modelled using the Fast-Eulerian method (Ferry & Balachandar, Intl J. Multiphase Flow, vol. 27, issue 7, pp. 1199–1226, 2001) that includes the droplet rise velocity. Radial concentration profiles and turbulent concentration fluxes for droplets of different sizes are analysed to quantify the ‘trajectory crossing effect’, when relative motions between particles and turbulent eddies tend to reduce turbulent diffusion. For finer LES grid resolutions, the model captures the differential, size-based dispersion characteristics of the droplets with the transverse dispersion of the larger droplet sizes suppressed, since trajectory crossing effects are explicitly resolved in LES. We examine a similarity solution model for the size-dependent radial concentration profiles based on a modified Schmidt number derived from the theory of turbulent diffusion of particles in the atmosphere proposed by Csanady (J. Atmos. Sci., vol. 20, pp. 201–208, 1963). The results are validated with the high resolution LES data and show good agreement. Then the size-dependent Schmidt number model is reformulated as a model for unresolved subgrid-scale trajectory crossing effects and used to calculate the subgrid concentration flux in a coarse LES of a turbulent jet, with slightly buoyant droplets injected at the centreline in the self-similar region of the jet. The results are compared to a simulation with higher grid resolution and a coarse simulation with a constant Schmidt number subgrid-scale model. We find that the subgrid model enhances the prediction accuracy of the concentration profiles and turbulent concentration flux for the coarse LES.

Atomization is often accompanied by phase change, which could significantly affect performance parameters such as the cooling efficiency and combustion efficiency of atomization. Nevertheless, the effect of phase change on jet atomization is rarely numerically studied due to the complexity of the coupling of the aerodynamics and the thermodynamics as well as the modelling difficulty caused by the cross-scale flow. In this study, comprehensive direct numerical simulations were carried out to evaluate the effects of phase change on the primary breakup and secondary atomization. Two methods dealing with phase-interface movement and mass transfer across the interface are built to meet the requirements of different modelling scales and Weber numbers. Simulation results indicate that phase change affects the flow behaviours and volume distribution of broken droplets in the primary breakup. In the secondary atomization, phase change leads to significantly different deforming morphologies of droplets with low Weber number and a more thorough breakup of droplets with high Weber number.

The dynamics and parameterization of mixing in temporally evolving turbulent open-channel flow is investigated through direct numerical simulations as the flow transitions from an initially neutral state to stable stratification. We observe three distinctly different mixing regimes separated by transitional values of turbulent Froude number $Fr$: a weakly stratified regime for $Fr >1$; an intermediate regime for $0.3< Fr<1$; and a saturated regime for $Fr<0.3$. The mixing coefficient $\varGamma =B/\epsilon _K$, (where $B$ is the buoyancy flux and $\epsilon _K$ is the dissipation rate of kinetic energy), is well predicted by the parameterization schemes of Maffioli et al. (J. Fluid Mech., vol. 794, 2016) and Garanaik & Venayagamoorthy (J. Fluid Mech., vol. 867, 2019, pp. 323–333) across all three regimes through instantaneous measurements of $Fr$ and the ratio $L_E/L_O$, where $L_E$ and $L_O$ are the Ellison and Ozmidov length scales, respectively. The flux Richardson number $R_f = B/(B+\epsilon _K)$ shows linear dependence on the gradient Richardson number $Ri_g$ up to a transitional value of $Ri_g =0.25$, past which it saturates again to a constant value independent of $Fr$ or $Ri_g$. By examining the flow as a balance of inertial, shear and buoyancy forces, we derive physically based scaling relationships to demonstrate that $Ri_g \sim Fr^{-2}$ and $Ri_g \sim Fr^{-1}$ in the weakly and moderately stratified regimes and that $Ri_g$ becomes independent of $Fr$ in the saturated regime. Our results suggest that the $L_E/L_O \sim Fr^{-1}$ scaling of Garanaik & Venayagamoorthy (J. Fluid Mech., vol. 867, 2019, pp. 323–333) in the intermediate regime manifests due to the influence of mean shear. The differences in the relationships between $Fr$ and $L_E/L_O$ for non-sheared flows within this regime are discussed.

Direct numerical simulations of particle-laden turbulent channel flows at friction Reynolds number $Re_\tau$ from $600$ to $2000$ have been performed to examine the near-wall particle streaks. Different from the well-observed small-scale particle streaks in near-wall turbulence of low $Re_\tau$, the present results show large-scale particle streaks through the computational domain formed for relatively high-inertia particles at high $Re_\tau$. Transferred by large-scale sweep and ejection events ($Q^-$), these high-inertia particles preferentially accumulate in near-wall regions beneath the large-scale low-speed flow streaks observed in the logarithmic region. The corresponding Stokes numbers are associated with the lifetime of large-scale $Q^-$ structures, which increases as the Reynolds number grows. The small-scale particle streaks with a typical Stokes number $St_\nu \approx 30$ are mainly driven by the $Q^-$ structures in the buffer layer, whose lifetime is approximately $30$ in viscous time unit. Therefore, we propose a new structure-based Stokes number normalized by the lifetime of $Q^-$ structures of different scales. The relevant flow scales that control the formation of the large-scale particle streaks are parameterized by the structure-based Stokes number. The small-scale (large-scale) particle streaks are most prominent when the buffer-layer (large-scale) structure-based Stokes number approaches unity. The present findings reveal that formation of near-wall particle streaks is governed by the $Q^-$ structures of different scales, and the particles with different inertia respond efficiently to the $Q^-$ structures of corresponding scales with respect to the particle translational motion.

We present an experimental study of the large-scale vortex (or large-scale circulation, LSC) in turbulent Rayleigh–Bénard convection in a $\varGamma =\text {diameter}/\text {height}=2$ cylindrical cell. The working fluid is deionized water with Prandtl number ($Pr$) around 5.7, and the Rayleigh number ($Ra$) ranges from $7.64\times 10^7$ to $6.06\times 10^8$. We measured the velocity field in various vertical cross-sectional planes by using the planar particle image velocimetry technique. The velocity measurement in the LSC central plane shows that the flow is in the single-roll form, and the centre of the single-roll (vortex) does not always stay at the centre of the cell; instead, it orbits periodically in the direction opposite to the flow direction of the LSC, with its trajectory in the shape of an ellipse. The velocity measurements in the three vertical planes in parallel to the LSC central plane indicate that the flow is in the vortex tube form horizontally filling almost the whole cell, and the centre line of the vortex tube is consistent with the so-called ‘jump rope’ form proposed by a previous study that combined numerical simulation and local velocity measurements in the low $Pr$ case (Vogt et al., Proc. Natl Acad. Sci. USA, vol. 115, 2018, pp. 12674–12679). In addition, we found that the oscillation of the local velocity in $\varGamma =2$ cells originates from the periodical orbiting of the vortex centre. Our velocity measurements further indicate that the vortex centre orbiting is absent in $\varGamma =1$ cells, at least in the $Ra$ range of our experiments.