Active suspensions encompass a wide range of complex fluids containing microscale energy-injecting particles, such as cells, bacteria or artificially powered active colloids. Because they are intrinsically non-equilibrium, active suspensions can display a number of fascinating phenomena, including turbulent-like large-scale coherent motion and enhanced diffusion. Here, using a recently developed active fast Stokesian dynamics method, we present a detailed numerical study of the hydrodynamic diffusion in apolar active suspensions of squirmers. Specifically, we simulate suspensions of active but non-self-propelling spherical squirmers (or ‘shakers’), of either puller type or pusher type, at volume fractions from 0.5 % to 55 %. Our results show little difference between pulling and pushing shakers in their instantaneous and long-time dynamics, where the translational dynamics varies non-monotonically with the volume fraction, with a peak diffusivity at around 10 % to 20 %, in stark contrast to suspensions of self-propelling particles. On the other hand, the rotational dynamics tends to increase with the volume fraction as is the case for self-propelling particles. To explain these dynamics, we provide detailed scaling and statistical analyses based on the activity-induced hydrodynamic interactions and the observed microstructural correlations, which display a weak local order. Overall, these results elucidate and highlight the different effects of particle activity versus motility on the collective dynamics and transport phenomena in active fluids.

The shear-induced diffusivity of non-Brownian spheres in monodisperse suspensions undergoing viscous flow was calculated using simulations that account for particle roughness and friction as independent parameters. The diffusivity increases significantly as the friction coefficient is increased, and the effect is largest on rougher particles. Roughness reduces the transverse diffusivities relative to smoother particles for sufficiently concentrated suspensions of frictionless and low-friction particles. However, the diffusivity of roughened particles is larger than smoother ones at high values of the friction coefficient. The increase of the diffusivity with friction is associated with a significant broadening of the variance of the rotational velocities. The most prevalent observation, when correlating the microstructure to changes in diffusivity for frictionless particles, is that less diffusive systems, with larger roughness, form layers along the flow direction. These results confirm previous experimental and simulation results that roughness can decrease diffusivity at large concentrations using a more detailed model. Also, comparisons of the simulation results with previously published experimental measurements indicate that friction improves the alignment of the results with experiments.

Settling velocity statistics for dilute, non-Brownian homogeneous suspensions of polydisperse spheres having a log-normal size distribution are generated from Stokesian dynamics simulations, as a function of the total volume fraction $\phi$ and normalised width $\alpha$ of the particle size distribution. Several hundred instantaneous configurations are averaged to obtain reliable statistics. The paper reports data for the average and fluctuating settling velocity of each particle class in a suspension that is widely polydisperse – previous work was limited to only two or three classes, and the average settling velocity of each particle class was in most cases not reported – and provides an assessment of the accuracy of the analytical models proposed by Batchelor, Richardson & Zaki, Davis & Gecol and Masliyah–Lockett–Bassoon in predicting the simulation data. A limited comparison with dynamic simulations in which the particle microstructure is allowed to evolve in time is also included.

The motion of a sphere freely rising or falling in a 5d (d is the diameter of the sphere) square tube was numerically studied for the sphere-to-fluid density ratio ranging from 0.1 to 2.3 (0.1 ≤ ρs/ρ ≤ 2.3, ρs is the density of spheres and ρ the fluid density) and Galileo number from 140 to 230 (140 ≤ Ga ≤ 230). We report that Hopf bifurcation occurs at Gacrit ≈ 157, where both the heavy and light spheres lose stability. The helical motion is widely seen for all spheres at Ga > 160 resulting from a double-threaded vortex interacting with the tube walls, which becomes irregular at Ga ≥ 190 where heavy spheres act differently from their counterparts; that is, heavy spheres change their helical directions alternately while light spheres exhibit helical trajectories with jaggedness in connection with the shedding of the double-threaded vortices. This is because of the difference in inertia between the heavy and light spheres. We also checked the oscillation periods for the helical motion of the spheres. They show opposite variations with ρs/ρ for the two types of spheres. Light spheres (ρs/ρ ≤ 0.7) reach a zigzagging regime at Ga ≥ 200 where a vortex loop (hairpin-like vortical structure) is formed which may develop into a vortex ring downstream at small ρs/ρ. This might be the first time a transition from the helical motion to the zigzagging motion for heavy spheres (ρs/ρ ≥ 1.8) has been reported. Finally, we examined the dependence of both the terminal Reynolds number and the drag coefficient of the spheres on the Galileo number.

Colloidal gels have strong industrial relevance as they can behave as liquids or solids. The latter allows them to support a buoyant weight against gravity. However, the system is intrinsically out of equilibrium, which means that the colloids must eventually settle out of the suspension. The process of settling has been captured theoretically, but the presence of a delay time during which the gel appears relatively unaffected by gravity has not. Here, we modify existing frameworks to capture this delay, by treating the gel as a continuum with viscous response that is based on the local bond density. We can solve our model numerically to obtain the evolution of the colloid density profile and recover qualitatively the accumulation of a dense layer on top of the settling gel, as is observed experimentally in depletion gels. This numerical study is complemented by a theoretical analysis that allows us to identify an emergent time and length scale that set the dynamics of the gel. Our model provides a solid foundation for future studies that incorporate hydrodynamic erosion and tackle industrially relevant geometries.

The featured article ‘Break-up of a falling drop containing dispersed particles’ (Nitsche and Batchelor, J. Fluid Mech., 1997, vol. 340, pp. 161–175) is G. K. Batchelor's last published paper with his former postdoctoral associate J. M. Nitsche. The objective of the study was to investigate the randomness of the velocities of interacting rigid particles falling under gravity through a viscous fluid at a small Reynolds number and its consequence for the breakup of a falling cloud of particles. The study focused on a quintessential problem of the collective dynamics of interacting particles and has been an inspiration for subsequent work.

Shear-induced self-diffusion is a fundamental mode of transport in granular flows. Yet its critical behaviour and dependence on the particle solid fraction are still unclear. Here, we rationalize these dependencies by performing two-dimensional pressure-imposed numerical simulations of dense non-Brownian frictional suspensions. Our results, combined with existing numerical data on inertial granular flows, show that the shear-induced diffusion coefficients of both systems can be captured by a single function of the distance to jamming. They further show that the grain diffusive behaviour is underpinned by a specific random walk process, having a constant elementary step length driven at a frequency that increases with the solid fraction. The proposed scaling laws pave the way for a better understanding of mixing processes in granular media.

The transient shear-induced particle migration of frictional non-Brownian suspensions is studied using particle-resolved simulations. The numerical method – the fictitious domain method – is well suited to heterogeneous flows thanks to a frame-invariant formulation of the subgrid (lubrication) corrections that does not involve the ambient flow (Orsi et al., J. Comput. Phys., vol. 474, 2023, 111823). The paper aims to give an accurate quantitative picture of the mass and momentum balance during the flow. The various assumptions and local constitutive laws that together form the suspension balance model (SBM) are thoroughly examined. To this purpose, the various quantities of interest are locally averaged in space and time, and their profile across the channel is extensively studied, with specific attention to the time evolution of the different contributions, either hydrodynamic in nature or from contact interactions, to the shear and normal stresses. The latter, together with the velocity gradient in the wall-normal direction and the volume fraction profile, yield the local constitutive laws, which are compared with their counterpart obtained in homogeneous shear flow. A fair agreement is observed except in a layering area at the boundaries and at the very centre of the channel. In addition, the main assumption of the SBM, i.e. the local relation between the hydrodynamic force on the particles and the particle flux, is meticulously investigated. The hydrodynamic force is found to be mainly a drag, except in the lower range of the probed volume fractions, where a non-drag contribution is observed.

The spreading of large viscous drops of density-matched suspensions of non-Brownian spheres on a smooth solid surface is investigated experimentally at the global drop scale. The focus is on dense suspensions with a solid volume fraction equal to or greater than $40\,\%$, and for drops larger than the capillary length, i.e. for which the spreading is governed by the balance of gravitational and viscous forces. Our findings indicate that all liquids exhibit a power-law behaviour typical of gravity-driven dynamics, albeit with an effective suspension viscosity that is smaller than the bulk value. When the height of the drop is of the order of the particle size, the power law breaks down as the particles freeze while the contact line continues to advance.

Many natural and industrial processes involve the flow of fluids made of solid particles suspended in non-Newtonian liquid matrices, which are challenging to control due to the fluid's nonlinear rheology. In the present work, a Taylor–Couette canonical system is used to investigate the flow of dilute to semi-dilute suspensions of neutrally buoyant spherical particles in highly elastic base polymer solutions. Friction measurement synchronized with direct flow visualization are combined to characterize the critical conditions for the onset of elasto-inertial instabilities, expected here as a direct transition to elasto-inertial turbulence (EIT). Adding a low particle volume fraction (${\leq }2\,\%$, dilute regime) does not affect the nature of the primary transition and reduces the critical Weissenberg number for the onset of EIT, despite a significant decrease in the apparent fluid elasticity. However, for particle volume fractions ${\geq }6\,\%$ (semi-dilute regime), EIT is no longer observed in the explored Reynolds range, suggesting an apparent relaminarization with yet not further decrease in fluid elasticity. Instead, a new regime, termed here elasto-inertial dissipative (EID), was uncovered. It originates from particle–particle interactions altering particle–polymer interactions and occurring under elasto-inertial conditions comparable to those of EIT. Increasing particle volume fraction in the semi-dilute regime and, in so, the particle contribution to the overall viscosity, delays the onset of EID similarly to what was observed previously for EIT in lower elasticity fluids. After this onset, a decrease in the pseudo-Nusselt number observed with increasing inertia and particle-to-polymer concentration ratio confirms a particle-induced alteration of energy transfer in the flow.

When a partially miscible binary mixture is quenched below its critical temperature, it transitions from its single-phase to a two-phase region, undergoing phase separation. The processes of formation and coalescence of droplets are driven by diffusive and convective phenomena, taking place isotropically in the system. The application of an external force field, which exerts a different contribution on the two species, breaks the symmetry of phase separation, leading to the segregation of two equilibrated phases separated by a single interface. This study investigates the dynamics of phase segregation under an external force. The effects of various force magnitudes, captured by the Bond number, in both high- and low-viscosity mixtures, distinguished by different fluidity numbers, are quantified via numerical simulations by using the phase field model. The intricate dynamics of formation, floating and coalescence of droplets towards complete segregation are described along with the quantification of the segregation time, revealing different patterns for high and low Bond numbers. Results show that in none of the cases, formation and floating can be regarded as strictly serial processes. A universal scaling between segregation time, Bond number, fluidity number and domain size is not possible, with a power-law dependence emerging only under the diffusion-dominated regime.

Laboratory experiments and particle-resolved simulations are employed to investigate the settling dynamics of a pair of rigidly connected spherical particles of unequal size. They yield a detailed picture of the transient evolution and the terminal values of the aggregate's orientation angle and its settling and drift velocities as functions of the aspect ratio and the Galileo number $Ga$, which denotes the ratio of buoyancy and viscous forces acting on the aggregate. At low to moderate values of $Ga$, the aggregate's orientation and velocity converge to their terminal values monotonically, whereas for higher $Ga$-values the aggregate tends to undergo a more complex motion. If the aggregate assumes an asymmetric terminal orientation, it displays a non-zero terminal drift velocity. For diameter ratios much larger than one and small $Ga$, the terminal orientation of the aggregate becomes approximately vertical, whereas when $Ga$ is sufficiently large for flow separation to occur, the aggregate orients itself such that the smaller sphere is located at the separation line. Empirical scaling laws are obtained for the terminal settling velocity and orientation angle as functions of the aspect ratio and $Ga$ for diameter ratios from 1 to 4 and particle-to-fluid density ratios from 1.3 to 5. An analysis of the accompanying flow field shows the formation of vortical structures exhibiting complex topologies in the aggregate's wake, and indicates the formation of a horizontal pressure gradient across the larger sphere, which represents the main reason for the emergence of the drift velocity.

The first experimental results on pattern transitions in the co-rotation regime (i.e. the rotation ratio $\varOmega = \omega _o/\omega _i > 0$, where $\omega _i$ and $\omega _o$ are the angular speeds of the inner and outer cylinders, respectively) of the Taylor–Couette flow (TCF) are reported for a neutrally buoyant suspension of non-colloidal particles, up to a particle volume fraction of $\phi = 0.3$. While the stationary Taylor vortex flow (TVF) is the primary bifurcating state in dilute suspensions ($\phi \leq ~0.05$), the non-axisymmetric oscillatory states, such as the spiral vortex flow (SVF) and the ribbon (RIB), appear as primary bifurcations with increasing particle loading, with an overall de-stabilization of the primary bifurcating states (TVF/SVF/RIB) being found with increasing $\phi$ for all $\varOmega \geq ~0$. At small co-rotations ($\varOmega \sim 0$), the particles play the dual role of stabilization ($\phi < 0.1$) and destabilization ($\phi \geq ~0.1$) on the secondary/tertiary oscillatory states. The distinctive features of the ‘particle-induced’ spiral vortices are identified and contrasted with those of the ‘fluid-induced’ spirals that operate in the counter-rotation regime.

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.

The clustering of debris floating on liquid interfaces such as water surfaces is a complex phenomenon that finds its applications in numerous examples from industrial processing and environmental systems. The recent paper by Shin & Coletti (J. Fluid Mech., vol. 984, 2024, R7) presents an experimental campaign investigating the three effects of turbulence, particle interactions and interfacial effects, to elucidate how the three force scales drag, capillary forces and lubrication give rise to three distinct regimes of clustering in dense suspensions. The study, hence, provides a useful systematic to categorize the clustering mechanisms. As an important finding, it is shown that, depending on volume fraction and non-dimensional turbulent shear, particles either tend to cluster into aggregate sizes larger than the Kolmogorov scale or can break into pieces that are as small as the primary particle size.

The orientational dynamics of a spherical magnetic particle in linear shear flow subjected to an oscillating magnetic field in the flow plane is analysed in the viscous limit. The shear is in the $X$–$Y$ plane, the magnetic field is in the $X$ direction and the vorticity is perpendicular to the flow in the $Z$ direction. The relevant dimensionless groups are $\omega ^\ast$, the ratio of the frequency of the magnetic field and the strain rate, and $\varSigma$, the ratio of the magnetic and hydrodynamic torques. As $\varSigma$ is decreased, there is a transition from in-plane rotation, where the rotation is in the flow ($X$–$Y$) plane, to out-of-plane rotation, where the orientation vector is not necessarily in the $X$–$Y$ plane and the dynamics depends on the initial orientation. The particle rotation is phase-locked for in-plane rotation with discrete odd rotation number (number of rotations in one period of magnetic field oscillation), while the orbits are quasi-periodic with non-integer rotation number for out-of-plane rotation. For $\varSigma \gg 1$, regions of odd rotation number $n_o$ are bound by the lines $8 (n_o-1) \varSigma \omega ^\ast = 1$ and $8 (n_o+1) \varSigma \omega ^\ast = 1$, and there are discontinuous changes in the rotation number and mean and root-mean-square torque at these lines. For $\varSigma \ll 1$, the domains of in-plane rotation of finite width in the $\omega ^\ast$–$\varSigma$ plane extend into downward cusps at $\omega ^\ast = {1}/{2 n_o}$. The orbits are quasi-periodic between these domains, where the rotation is out of plane.

Particulate flows at moderate particle Reynolds numbers are important in critical engineering and geological applications. This experimental study explores neutrally buoyant suspensions in an outer-rotating coaxial rheometer for solid fractions, $\phi$, from 0.1 to 0.5, and particle Reynolds number, $Re$, from 0.5 to 800, covering laminar, transitional and turbulent regimes; $Re$ is defined in terms of the square of the particle diameter and the shear rate. For $0.1 < \phi < 0.4$ and $0.5 < Re <10$, the direct torque measurements normalised by the laminar flow torque, $M/M_{lam}$, are independent of $Re$, but depend on $\phi$. For the same range of $\phi$ and for $10< Re<100$, the normalised torques depend on both $\phi$ and $Re$, and show an increasing dependence on $Re$. As $Re$ increases, the flow transitions to turbulence. Small particles delay the turbulent transition for $\phi \leqslant 0.3$, while large particles augment the transition. A modified Reynolds number, $Re^\prime$, that depends linearly on the particle diameter and the maximum velocity, $U_{o}$, is introduced for both laminar and turbulent flows and shows a better correlation of the results as compared with $Re$. For $\phi = 50\,\%$, the normalised torque minus the torque at zero rotational speed is nearly independent of $Re^\prime$. Rheological models based on $Re^\prime$ and the Krieger–Dougherty relative viscosity are proposed in the laminar regime for $10< Re^\prime <500$; in the turbulent regime, a correlation is proposed in terms of $Re^\prime$ and $\phi$ for $1000< Re^\prime < 6000$.

We develop a mean-field model to examine the stability of a ‘quasi-2-D suspension’ of elongated particles embedded within a viscous membrane. This geometry represents several biological and synthetic settings, and we reveal mechanisms by which the anisotropic mobility of particles interacts with long-ranged viscous membrane hydrodynamics. We first show that a system of slender rod-like particles driven by a constant force is unstable to perturbations in concentration – much like sedimentation in analogous 3-D suspensions – so long as membrane viscous stresses dominate. However, increasing the contribution of viscous stresses from the surrounding 3-D fluid(s) suppresses such an instability. We then tie this result to the hydrodynamic disturbances generated by each particle in the plane of the membrane and show that enhancing subphase viscous contributions generates extensional fields that orient neighbouring particles in a manner that draws them apart. The balance of flux of particles aggregating versus separating then leads to a wave number selection in the mean-field model.

We have performed numerical simulations to investigate the phenomenon of axial pattern formation exhibited by a non-neutrally buoyant concentrated suspension. Continuum modelling of the concentrated suspension is done using the suspension balance model to identify the underlying mechanism of the phenomenon. We demonstrate that axial concentration variations become amplified to axial bands owing to the influence of the second normal stress difference ($N_2$), and the first normal difference ($N_1$) accentuates the effect of $N_2$. We demonstrate that the end walls of the rotating cylinder are necessary to prevent the smearing out of axial bands but are not a direct cause of the phenomenon.