The elasto-inertial focusing and rotating characteristics of spheroids in a square channel flow of Oldroyd-B viscoelastic fluids are studied by the direct forcing/fictitious domain method. The rotational behaviours, changes in the equilibrium positions and travel distances are explored to analyse the mechanisms of spheroid migration in viscoelastic fluids. Within the present simulated parameters (1 ≤ Re ≤ 100, 0 ≤ Wi ≤ 2, 0.4 ≤ α ≤3), the results show that there are four kinds of equilibrium positions and six (five) kinds of rotational behaviours for the elasto-inertial migration of prolate (oblate) spheroids. We are the first to identify a new rotational mode for the migration of prolate spheroids. Only when the particles are initially located at a corner and wall bisector, some special initial orientations of the spheroids have an impact on the final equilibrium position and rotational mode. In other general initial positions, the initial orientation of the spheroid has a negligible effect. A higher Weissenberg number means the faster the particles migrate to the equilibrium position. The spheroid gradually changes from the corner (CO), channel centreline (CC), diagonal line (DL) and cross-section midline (CSM) equilibrium positions as the elastic number decreases, depending on the aspect ratio, initial orientation and rotational behaviour of the particles and the elastic number of the fluid. When the elastic number is less than the critical value, the types of rotational modes of the spheroids are reduced. By controlling the elastic number near the critical value, spheroids with different aspect ratios can be efficiently separated.

A Janus sphere with a stick–slip pattern can behave quite differently in its hydrodynamics compared with a no-slip or uniform-slip sphere. Here, using the Lorentz reciprocal theorem in conjunction with surface harmonic expansion, we rigorously derive the extended Faxén formula for the stresslet of a weakly stick–slip Janus sphere, capable of describing the anisotropic nature of the stresslet with an arbitrary axisymmetric stick–slip pattern in an arbitrary background flow. We find that slip anisotropy not only causes a variety of additional contributions to the stresslet, but also naturally renders a stresslet–rotation coupling that may turn a suspension of couple-free stick–slip Janus spheres into a dipolar one under the actions of an external couple. Moreover, to correctly account for the impacts of slip anisotropy on the stresslet, it is necessary to include at least the first four surface harmonic contributions. As a result, the anisotropies of both the stresslet and torque on the sphere in a linear flow field are purely reflected by a symmetric quadrupole and hexadecapole. These hydrodynamic quantities can be further mediated by an antisymmetric dipole and octupole due to the gradients of the imposed strain field. The average bulk stress and effective viscosity for a suspension of stick–slip spheres are also determined, showing characteristics quite distinct from those of a suspension of near spheres. If the spheres possess permanent dipole moments, in particular, additional stresslets and couplets can be generated by an applied external couple on each sphere and added into the bulk stress, accompanied by non-Jeffery orientational orbits of such dipolar stick–slip spheres. In addition to the above, the extended Faxén stresslet and torque relations found in this work will also provide the formulae needed for tackling problems involving hydrodynamically interacting stick–slip spheres on which small slip anisotropy may have profound impacts.

Numerical simulations of the interaction of internal solitary waves (ISWs) of opposite polarity are conducted by solving the incompressible Euler equations under the Boussinesq approximation. A double-pycnocline stratification is used. A method to determine when ISWs of both polarities exist is also presented. The simulations confirm previous work that the interaction of waves of the same polarity are soliton-like; however, here it is shown that when a fast ISW with the same polarity as a Korteweg–de Vries (KdV) solitary wave catches up and interacts with a slower ISW of opposite polarity, the interaction can be far from soliton-like. The energy in the fast KdV-polarity wave can increase by more than a factor of 5 while the energy in the slower negative-KdV-polarity wave can decrease by 50 %. Large trailing wave trains may be generated and in some cases multiple ISWs with KdV polarity may be formed by the interaction.

We present an analytical study of two-dimensional flow in a wedge driven by a time-dependent surface heat flux as a model problem to understand buoyancy-induced cross-shore flow. Besides the turbulent Prandtl number and the relevant Rayleigh number, both assumed to be of order unity, the solution is seen to depend on the geometry through a parameter $\beta$, representing the bottom slope. An analytic solution is sought in the asymptotic limit $\beta \ll 1$ for a water layer bounded by an adiabatic bottom surface subject to a harmonic heat flux on the upper surface. The analysis reveals that the motion at leading order can be expressed as the sum of a harmonic component and a steady component, the latter driven by nonlinear advection. This steady-streaming motion includes a near-shore vortex with associated recirculating motion that can affect cross-shore transport and dispersion in coastal environments. The analytical solution is compared with numerical solutions of the complete conservation equations for small values of $\beta$. Excellent quantitative agreement is found for values of the Rayleigh number below a critical value at which the periodic solution undergoes a period-doubling bifurcation, leading to the establishment of thermal-instability cells that dominate the offshore flow dynamics, while the near-shore dynamics remains well described by the analytical solution. The analysis illustrates that a periodic heat input that leads to a vertically inhomogeneous temperature distribution can result in residual motion, net heat fluxes and persistent temperature structure in the cross-shore direction.

The spontaneous emergence of structure is a ubiquitous process observed in fluid and plasma turbulence. These structures typically manifest as flows which remain coherent over a range of spatial and temporal scales, resisting statistically homogeneous description. This work conducts a computational and theoretical study of coherence in turbulent flows in the stochastically forced barotropic $\beta$-plane quasi-geostrophic equations. These equations serve as a prototypical two-dimensional model for turbulent flows in Jovian atmospheres, and can also be extended to study flows in magnetically confined fusion plasmas. First, analysis of direct numerical simulations demonstrates that a significant fraction of the flow energy is organized into coherent large-scale Rossby wave eigenmodes, comparable with the total energy in the zonal flows. A characterization is given for Rossby wave eigenmodes as nearly integrable perturbations to zonal flow Lagrangian trajectories, linking finite-dimensional deterministic Hamiltonian chaos in the plane to a laminar-to-turbulent flow transition. Poincaré section analysis reveals that Lagrangian flows induced by the zonal flows plus large-scale waves exhibit localized chaotic regions bounded by invariant tori, manifesting as Rossby wave breaking in the absence of critical layers. It is argued that the surviving invariant tori organize the large-scale flows into a single temporally and zonally varying laminar flow, suggesting a form of self-organization and wave stability that can account for the resilience of the observed large-amplitude Rossby waves.

Fluid-inertia torque remarkably affects the orientation of non-spherical particles in Newtonian flows whereas this torque induced by convective fluid inertia in particle-laden pseudo-plastic flows is still unknown. In the present study we numerically investigate the fluid-inertia torque on a neutrally buoyant spheroid in the Carreau-type pseudo-plastic fluid flows at finite Reynolds numbers with the immersed boundary method. The results show that compared with the fluid-inertia torque in Newtonian flows, the magnitude of the fluid-inertia torque on spheroids is remarkably attenuated by the shear-thinning rheology in pseudo-plastic fluid flows. The deviation of fluid-inertia torque between pseudo-plastic and Newtonian flows is more significant with decreasing Reynolds numbers, indicating the importance of the effect of shear-thinning rheology at small Reynolds numbers. Moreover, the spheroid rotation rate is reduced in pseudo-plastic fluids, and the equilibrium orientation of oblate spheroids changes non-monotonically with the shear-thinning effect in the linear shear flow of pseudo-plastic fluids. The present findings imply the importance of the effect of shear-thinning rheology on the torques of spheroids, which could be potentially applied for the control of particle orientations in pseudo-plastic fluids in the future.

We study the motion and deformation of a single bubble rising inside a cylindrical container filled with a viscoelastic material. We solve numerically the mass and momentum balances along with the constitutive equation for the viscoelastic stresses, without assuming axial symmetry to allow the growth of three-dimensional disturbances. Hence, we may predict the emergence of the notorious knife-edge shape of the bubble, which is a result of a purely elastic instability triggered in the locality of the trailing edge. Our results compare well with existing experiments. We visualize, for the first time to the best of our knowledge, the flow kinematics and dynamics that arise downstream of the bubble. We propose two quantities, one kinematical and one geometrical, for the determination of the onset of the instability. We demonstrate that extension-rate thinning in the constitutive law is necessary for the emergence of the instability. Moreover, our results indicate that increasing (a) the deformability of the bubble and (b) the initial extension rate hardening of the viscoelastic material, prior to thinning, triggers the instability earlier. These novel findings help us formulate and propose a mechanism that controls the onset of the instability and explain why the knife-edge shape is not encountered as frequently.

We study the application of Taylor's frozen hypothesis to the pressure fluctuations in turbulent channels by means of spatio-temporal data from direct numerical simulations with large computational domains up to the friction Reynolds number $R{e_\tau } = 2000$. The applicability of the hypothesis is quantitatively verified by comparing the wavenumber and Taylor (frequency) premultiplied spectra of the pressure fluctuations at each distance y from the wall. Using the local mean velocity $U(y)$ for the hypothesis leads to a large difference between both spectra with a value of $O(50\,{\%})$ for its maximum from the wall to $y/h \approx 0.2$, where h is the channel half-depth. Alternatively, the convection velocity of the pressure fluctuations ${C_p}(y)$, originally defined by Del Álamo & Jiménez (J. Fluid Mech., vol. 640, 2009, pp. 5–26) as a function of y, is investigated and adopted for the hypothesis. It is nearly equal to $U(y)$ from ${y^ + } = 20$ to the channel centre, where ${y^ + } = y{u_\tau }/\nu $, in which ${u_\tau }$ and $\nu $ represent the friction velocity and kinematic viscosity, respectively. For ${y^ + } \le 20$, ${C_p}(y)$ is almost constant with a value of around $12{u_\tau }$. Applying ${C_p}(y)$ for the hypothesis decreases the difference between both spectra down to a value of $O(10\,{\%})$ for its maximum from the wall to $y/h \approx 0.2$. Finally, the difference between the pressure wavenumber and frequency premultiplied spectra near the wall is reduced further via applying a wavenumber-dependent convection velocity. This wavenumber-dependent convection velocity is the linear combination of the convection velocities of the turbulent structures associated with the pressure field weighted by their relative contributions to the pressure variance.

At low Reynolds numbers, the flow through a pipe with a sudden expansion is characterized by the localized occurrence of flow instabilities, with the formation of a so-called turbulence puff. In the literature, physical experiments typically predict earlier occurrence of turbulence than computational fluid dynamics simulations. However, the behaviour of ‘natural’ transition to turbulence without perturbations, and the dependence of transition to turbulence on perturbations, are not yet fully understood, particularly for the simulations. The purpose of the present study is therefore to investigate this flow, including possible sources of perturbation in numerical simulations, and to evaluate their effect on transition to turbulence. Through the exploration of different flow rates, numerical settings and inlet perturbation amplitudes, and by evaluating coarse and refined simulations, insights into low-Reynolds-number transitional turbulent flows are obtained. The turbulence kinetic energy budget of the turbulence puff or slug characteristic of this flow is also evaluated. In conclusion, even when perturbations are not intentionally added, there can still be significant sources of numerical perturbation and error that trigger turbulence in simulations, but perturbations will need to be added in refined simulations in order to produce turbulence. Finally, the results agree with the notion that there may not be a scenario where the flow transitions naturally to turbulence without any perturbation.

Chiral fluids – such as fluids under rotation or a magnetic field as well as synthetic and biological active fluids – flow in a different way than ordinary ones. Due to symmetries broken at the microscopic level, chiral fluids may have asymmetric stress and viscosity tensors, for example giving rise to a hydrostatic torque or non-dissipative (odd) and parity-violating viscosities. In this article, we investigate the motion of rigid bodies in such an anisotropic fluid in the incompressible Stokes regime through the mobility matrix, which encodes the response of a solid body to forces and torques. We demonstrate how the form of the mobility matrix, which is usually determined by particle geometry, can be analogously controlled by the symmetries of the fluid. By computing the mobility matrix for simple shapes in a three-dimensional (3-D) anisotropic chiral fluid, we predict counterintuitive phenomena such as motion at an angle to the direction of applied forces and spinning under the force of gravity.

The properties of multiphase flows are challenging to measure, and yet effective properties are fundamental to modelling and predicting flow behaviour. The current study is motivated by rheometric measurements of a gas-fluidized bed using a coaxial rheometer in which the fluidization rate and the rotational speed can be varied independently. The measured torque displays a range of rheological states: quasistatic, dense granular flow behaviour at low fluidization rates and low-to-moderate shear rates; turbulent toroidal-vortex flow at high shear rates and moderate-to-high fluidization rates; and viscous-like behaviour with rate-dependent torque at high shear rates and low fluidization or at low shear rates and high fluidization. To understand the solid-like to fluid-like transitions, additional experiments were performed in the same rheometer using single-phase liquid and liquid–solid suspensions. The fluidized bed experiments are modelled as a Bingham plastic for low fluidization rates, and as a shear-thinning Carreau liquid at high fluidization rates. The suspensions are modelled using the Krieger–Dougherty effective viscosity. The results demonstrate that, by using the effective properties, the inverse Bingham number marks the transition from solid-like to viscous-flow behaviour; a modified gap Reynolds number based on the thickness of the shear layer specifies the transition from solid-like to turbulent vortical flow; and a gap Reynolds number distinguishes viscous behaviour from turbulent vortical flow. The results further demonstrate that these different multiphase flows undergo analogous flow transitions at similar Bingham or Reynolds numbers and the corresponding dimensionless torques show comparable scaling in response to annular shear.

Compound bubbles with a liquid coating in another continuous immiscible bulk phase are ubiquitous in a wide range of natural and industrial processes. Their formation, rise and ultimate bursting at the air–liquid interface play crucial roles in the transport and fate of natural organic matter and contaminants. However, the dynamics of compound bubbles has not received considerable attention until recently. Here, inspired by our previous work (Yang et al., Nat. Phys., vol. 19, 2023, pp. 884–890), we investigate the entrainment of daughter oil droplets in bulk water produced by a bursting oil-coated bubble. We document that the size of the entrained daughter oil droplet is affected by the oil coating fraction and the bulk liquid properties. We rationalize this observation by balancing the viscous force exerted by the extensional flow produced by bubble bursting with the capillary force resisting the deformation of the oil coating, and considering the subsequent end-pinching process which finally entrains the daughter oil droplets. We propose a scaling analysis for the daughter oil droplet size that well captures the experimental results for a wide range of oil coating fractions and Ohnesorge numbers of the bulk liquid. In addition, we discuss the non-monotonic variation of daughter droplet size with the Ohnesorge number, and show the eventual absence of daughter droplets because of the strong viscous effect in the high-Ohnesorge-number regime. Our findings may advance the fundamental understanding of compound bubble bursting and provide guidance and modelling constraints for bubble-mediated contaminant transport in liquids.

Electrolyte solutions play an important role in energy storage devices, whose performance relies heavily on the electrokinetic processes at sub-micron scales. Although fluctuations and stochastic features become more critical at small scales, the long-range Coulomb interactions pose a particular challenge for both theoretical analysis and simulation of fluid systems with fluctuating hydrodynamic and electrostatic interactions. Here, we present a theoretical framework based on the Landau–Lifshitz theory to derive closed-form expressions for fluctuation correlations in electrolyte solutions, indicating significantly different decorrelation processes of ionic concentration fluctuations from hydrodynamic fluctuations, which provides insights for understanding transport phenomena of coupled fluctuating hydrodynamics and electrokinetics. Furthermore, we simulate fluctuating electrokinetic systems using both molecular dynamics (MD) with explicit ions and mesoscopic charged dissipative particle dynamics (cDPD) with semi-implicit ions, from which we identify that the spatial probability density functions of local charge density follow a gamma distribution at sub-nanometre scale (i.e. $0.3\,{\rm nm}$) and converge to a Gaussian distribution above nanometre scales (i.e. $1.55\,{\rm nm}$), indicating the existence of a lower limit of length scale for mesoscale models using Gaussian fluctuations. The temporal correlation functions of both hydrodynamic and electrokinetic fluctuations are computed from all-atom MD and mesoscale cDPD simulations, showing good agreement with the theoretical predictions based on the linearized fluctuating hydrodynamics theory.

Viscous streaming is an efficient rectification mechanism to exploit flow inertia at small scales for fluid and particle manipulation. It typically entails a fluid vibrating around an immersed solid feature that, by concentrating stresses, modulates the emergence of steady flows of useful topology. Motivated by its relevance in biological and artificial settings characterized by soft materials, recent studies have theoretically elucidated, in two dimensions, the impact of body elasticity on streaming flows. Here, we generalize those findings to three dimensions, via the minimal case of an immersed soft sphere. We first improve existing solutions for the rigid-sphere limit, by considering previously unaccounted terms. We then enable body compliance, exposing a three-dimensional, elastic streaming process available even in Stokes flows. Such effect, consistent with two-dimensional analyses but analytically distinct, is validated against direct numerical simulations and shown to translate to bodies of complex geometry and topology, paving the way for advanced forms of flow control.

We have developed a parameter-free, two-phase, volume-averaged approach to predictively describe the spin-up flow of dilute, cluster-free ferrofluids excited by low-frequency rotating magnetic fields. Predictive validation of the model was performed through a thorough comparison with local velocity profile measurements, and it demonstrated its ability to capture the spin-up flow dynamics without the need for parameter tuning by carefully delineating the validity domain of the ferrofluid dilutedness conditions. To gain insight into the underlying flow mechanisms, we performed a systematic parametric analysis examining the effects of the induced magnetic field, the dipolar interactions between magnetic nanoparticles and the demagnetizing field. How these mechanisms shape the flow of dilute ferrofluids excited by low-frequency rotating fields in a standard spin-up flow geometry has been addressed using probabilistic nanoparticle orientational dynamics, combining Faxén's laws and the Smoluchowski equation to describe the transport of particle magnetic moments. Our findings revealed that the induced magnetic field is the primary driving force of ferrofluid spin-up flow. The dipole interactions and demagnetizing field, on the other hand, contribute only as secondary phenomena to the overall flow behaviour. Furthermore, we have discussed the potential extension of the two-phase approach, in particular with respect to the formation of chain-like aggregates under the influence of strong magnetic fields. Overall, our study provides valuable insights into the complex dynamics of ferrofluid flow and contributes to a comprehensive understanding of the key mechanisms governing the spin-up flow of dilute ferrofluids excited by low-frequency rotating magnetic fields.

The influence of surface roughness on transition to turbulence in a Mach 4.5 boundary layer is studied using direct numerical simulations. Transition is initiated by the nonlinearly most dangerous inflow disturbance, which causes the earliest possible breakdown on a flat plate for the prescribed inflow energy and Mach number. This disturbance primarily comprises two normal second-mode instability waves and an oblique first mode. When localized roughness is introduced, its shape and location relative to the synchronization points of the inflow waves are confirmed to have a clear impact on the amplification of the second-mode instabilities. The change in modal amplification coincides with the change in the height of the near-wall region where the instability wave speed is supersonic relative to the mean flow; the net effect of a protruding roughness is destabilizing when placed upstream of the synchronization point and stabilizing when placed downstream. Assessment of the effect of the roughness location is followed by an optimization of the roughness height, abruptness and width with the objective of achieving maximum transition delay. The optimization is performed using an ensemble-variational (EnVar) approach, while the location of the roughness is fixed upstream of the synchronization points of the two second-mode waves. The optimal roughness disrupts the phase of the near-wall pressure waves, suppresses the amplification of the primary instability waves and mitigates the nonlinear interactions that lead to breakdown to turbulence. The outcome is a sustained non-turbulent flow throughout the computational domain.