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.

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.

Turbulent wakes are often characterized by dominant coherent structures over disparate scales. Dynamics of their behaviour can be attributed to interscale energy dynamics and triadic interactions. We develop a methodology to quantify the dynamics of kinetic energy of specific scales. Coherent motions are characterized by the triple decomposition and used to define mean, coherent and random velocity. Specific scales of coherent structures are identified through dynamic mode decomposition, whereby the total coherent velocity is separated into a set of velocities classified by frequency. The coherent kinetic energy of a specific scale is defined by a frequency triad of scale-specific velocities. Equations for the balance of scale-specific coherent kinetic energy are derived to interpret interscale dynamics. The methodology is demonstrated on three wake flows: (i) ${Re}=175$ flow over a cylinder; (ii) a direct numerical simulation of ${Re}=3900$ flow over a cylinder; and (iii) a large-eddy simulation of a utility-scale wind turbine. The cylinder wake cases show that energy transfer starts with vortex shedding and redistributes energy through resonance of higher harmonics. The scale-specific coherent kinetic energy balance quantifies the distribution of transport and transfer among coherent, mean and random scales. The coherent kinetic energy in the rotor scales and the hub vortex scale in the wind turbine interact to produce new scales. The analysis reveals that vortices at the blade root interact with the hub vortex formed behind the nacelle, which has implications for the proliferation of scales in the downwind near wake.

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.

Two-dimensional horizontally periodic Rayleigh–Bénard convection between stress-free boundaries displays two distinct types of states, depending on the initial conditions. Roll states are composed of pairs of counter-rotating convection rolls. Windy states are dominated by strong horizontal wind (also called zonal flow) that is vertically sheared, precludes convection rolls and suppresses heat transport. Windy states occur only when the Rayleigh number $Ra$ is sufficiently above the onset of convection. At intermediate $Ra$ values, windy states can be induced by suitable initial conditions, but they undergo a transition to roll states after finite lifetimes. At larger $Ra$ values, where windy states have been observed for the full duration of simulations, it is unknown whether they represent chaotic attractors or only metastable states that would eventually undergo a transition to roll states. We study this question using direct numerical simulations of a fluid with a Prandtl number of 10 in a layer whose horizontal period is eight times its height. At each of seven $Ra$ values between $9\times 10^6$ and $2.25\times 10^7$ we have carried out 200 or more simulations, all from initial conditions leading to windy convection with finite lifetimes. The lifetime statistics at each $Ra$ indicate a memoryless process with survival probability decreasing exponentially in time. The mean lifetimes grow with $Ra$ approximately as $Ra^4$. This analysis provides no $Ra$ value at which windy convection becomes stable; it might remain metastable at larger $Ra$ with extremely long lifetimes.

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.

We consider the role of gravity in solute transport when a thin droplet evaporates. Under the physically relevant assumptions that the contact line is pinned and the solutal Péclet number, ${Pe}$, is large, we identify two asymptotic regimes that depend on the size of the Bond number, ${Bo}$. When ${Bo} = O(1)$ as ${Pe}\rightarrow \infty$, the asymptotic structure of solute transport follows directly from the surface-tension-dominated regime, whereby advection drives solute towards the contact line, only to be countered by local diffusive effects, leading to the formation of the famous ‘coffee ring.’ In the distinguished limit in which ${Bo} = O({Pe}^{4/3})$ as ${Pe}\rightarrow \infty$, this interplay between advection and diffusion takes place alongside that between surface tension and gravity. In each regime, we perform a systematic asymptotic analysis of the solute transport and compare our predictions to numerical simulations. We identify the effect of gravity on the nascent coffee ring, providing quantitative predictions of the size, location and shape of the solute mass profile. In particular, for a fixed Péclet number, as the effect of gravity increases, the coffee ring is diminished in height and situated further from the contact line. Furthermore, for certain values of ${Bo}$, ${Pe}$ and the evaporation time, a secondary peak may exist inside the classical coffee ring. The onset of this secondary peak is linked to the change in type of the critical point in the solute mass profile at the droplet centre. Both the onset and the peak characteristics are shown to be independent of ${Pe}$.

Pattern-forming with externally imposed symmetry is ubiquitous in nature but little studied. We present experimental studies of pattern formation and selection by spatial periodic forcing in rapidly rotating convection. When periodic topographic structures are constructed on the heated boundary, they modulate the local temperature and velocity fields. Symmetric convection patterns in the form of regular vortex lattices are observed near the onset of convection, when the periodicity of the external forcing is set close to the intrinsic vortex spacing. We show that the new patterns arise as a dynamical process of imperfect bifurcation which is well described by a Ginzburg–Landau-like model. We explore the phase diagram of buoyancy strength and periodicity of external forcing to find the optimal experimental settings for which the vortex patterns best match that of the external forcing.

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.

A slip asymmetry can break the fore–aft symmetry of the local hydrodynamic force distribution on the surface of an otherwise no-slip or uniform-slip particle. Here, we use the Lorentz reciprocal theorem to demonstrate that such asymmetry, even in a fractional amount, can qualitatively alter the swimming characteristics of a self-propelled spherical squirmer, markedly different from those of no-slip or uniform-slip squirmers. Unlike the usual tangential squirming by the thrust-providing B1 mode and the type-determining B2 mode, we discover two unique features for a stick-slip squirmer. First, the squirmer can acquire a swimming velocity U without the B1 mode but simply by a symmetric extensile/contractile squirming from the B2 mode, which is able to reverse the swimming direction of the squirmer. Second, a stresslet $\boldsymbol{\mathsf{S}}$ can also be induced by a unidirectional squirming from the B1 mode, capable of inverting the squirmer's stresslet from extensile type to contractile type or vice versa to change the squirmer from puller to pusher or in a reverse manner. We further show that the two squirming modes can reinforce or compete with each other to enhance or diminish U and $\boldsymbol{\mathsf{S}}$ due to interplays between the asymmetric squirming forces on the stick and the slip faces. A phase diagram is also established to categorize a variety of newly emerging swimming states, such as an enhanced/degraded puller/pusher and a backward puller/pusher, depending on the relative strength of the squirming modes β = B2/B1, the direction of the stick-slip polarity and the degree of the slip disparity. As a result of such cooperative and competitive natures, a stick-slip squirmer can swim more or less efficiently than no-slip and uniform-slip ones. These distinctive features arising from stick-slip disparity can not only be made geometrically tuneable for steering the motion of a squirmer, but also provide new means for making efficient artificial microswimmers using amphiphilic Janus particles.

Different types of neural networks have been used to solve the flow sensing problem in turbulent flows, namely to estimate velocity in wall-parallel planes from wall measurements. Generative adversarial networks (GANs) are among the most promising methodologies, due to their more accurate estimations and better perceptual quality. This work tackles this flow sensing problem in the vicinity of the wall, addressing for the first time the reconstruction of the entire three-dimensional (3-D) field with a single network, i.e. a 3-D GAN. With this methodology, a single training and prediction process overcomes the limitation presented by the former approaches based on the independent estimation of wall-parallel planes. The network is capable of estimating the 3-D flow field with a level of error at each wall-normal distance comparable to that reported from wall-parallel plane estimations and at a lower training cost in terms of computational resources. The direct full 3-D reconstruction also unveils a direct interpretation in terms of coherent structures. It is shown that the accuracy of the network depends directly on the wall footprint of each individual turbulent structure. It is observed that wall-attached structures are predicted more accurately than wall-detached ones, especially at larger distances from the wall. Among wall-attached structures, smaller sweeps are reconstructed better than small ejections, while large ejections are reconstructed better than large sweeps as a consequence of their more intense footprint.

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.

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 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.