When a contaminated liquid evaporates from within a porous material, the impurities or dirt accumulate and deposit within the pore space. This occurs during the cleaning of filters and fouling of textiles, and is related to the ‘coffee-ring’ problem. To investigate how and where dirt is deposited in the pore space, we present a model for the motion of an evaporation front through a porous material, and the related accumulation, transport, and deposition of dirt, assuming that the liquid remains stationary. For physically relevant parameters, vapour transport out of the porous material is quasi-steady and we derive a single ordinary differential equation describing the motion of the evaporation front in time. Model solutions exhibit spatially non-uniform profiles of the deposited dirt-layer thickness through the porous material. The dirt accumulation and evaporation problems are coupled: deposited dirt hinders vapour transport through the porous material, slowing the evaporation. We identify two scenarios in which the porous material becomes clogged with dirt. Accumulation of suspended dirt at the evaporating interface along with slow dirt diffusion results in the deposited dirt layers clogging the pores at the evaporating interface, halting the drying and trapping liquid in the porous material. Alternatively, slow dirt deposition results in the suspended dirt being pushed far into the porous material by the evaporation, eventually leaving only dirt (with no liquid) in the pore space. We investigate the dynamics of both clogging scenarios, characterising the parameter regimes for which each occurs. Both clogging scenarios must be avoided in practice since they may be detrimental to future filter efficacy or textile breathability.

This paper presents a comprehensive theoretical analysis of the interaction between two quasi-steady burning droplets with differing temperatures, sizes and distances, building upon the mass-flux-potential model and flame-sheet assumption. In contrast to existing research, this study introduces a fresh perspective on droplet interactions by considering the different temperatures of the droplets. Utilizing the bispherical coordinate approach, theoretical solutions for the Stefan flow, scalar fields, droplet evaporation/burning rates, interaction coefficients and flame positions have been derived successfully. A comparison with extensive numerical simulations indicates a good agreement between the analytical and numerical results under a variety of conditions. It is revealed that proximity between the droplets causes non-uniform evaporation rates on their surfaces, and in some cases, leads to condensation on the cooler droplet. Notably, when the temperatures of the two droplets differ, this results in an uneven temperature distribution across the flame surface, and increasing the temperature of one droplet substantially elevates the temperature of the nearby flame. This study also establishes a criterion for the transition between different combustion modes, specifically between group and separated combustion. The findings of this study are crucial in deepening our understanding of evaporation and combustion processes, as well as the dynamics of flame spreading, local ignition, and extinction in systems involving multiple droplets.

The self-sustained interactions between a flexible film and periodic vortices epitomize the spirit of fish swimming and flag flapping in nature, involving intricate patterns of flow–structure coupling. Here, we comprehensively investigate the multiple coupling states of a film in the cylinder wake mainly with experiments, complemented by theoretical solutions and nonlinear dynamical analyses. Four regimes of film motion states are identified in the parameter space spanned by the reduced velocity and the length ratio. These regimes are (i) keeping stationary, (ii) deflection flutter, (iii) hybrid flutter and (iv) periodic large-amplitude flapping, each governed by a distinct coupling mechanism, involving regular and irregular Kármán vortices, local instability of the elongated shear layers and 2P mode vortex shedding. The film futtering in regimes (ii) and (iii) is substantiated to be chaotic and bears a resemblance to the ‘entraining state’ of fish behind an obstacle in the river. The periodic flapping in regime (iv) manifests itself in an amalgam of standing and travelling waves, and has intrinsic relations to the ‘Kármán gaiting’ of fish in periodic vortices. With the spatiotemporal reconstruction for the periodic flapping, we procure the energy distributions on the film, revealing the energy transfer processes between the film and the large-scale vortices. The findings unequivocally indicate that the flow–structure interaction during the energy-release stage of the film is more intense than that during the energy-extraction stage. Given the similarities, the mathematical and physical methods presented in this work are also applicable to the research on biological undulatory locomotion.

Considerable efforts have been devoted to the understanding of the small-scale characteristics in turbulent flows. While the universality of small-scale quantities has been established for incompressible flows, their extension to high-pressure transcritical flows remains an open area of research. To address this question, we investigate the real-fluid thermodynamic effects on small-scale velocity statistics of high-pressure transcritical wall-bounded turbulence. We show that in the locally isotropic region for transcritical flows, low-order moments of small-scale statistics collapse for all cases and Kolmogorov's (1941) theory holds. However, real-fluid thermodynamic effects introduce deviations in the tails of the probability density function for the velocity derivative and, consequently, high-order moments of velocity gradients and dissipation rate in transcritical flows cannot collapse in the locally isotropic region. Analysis of the intermittency shows that the low-order structure functions in transcritical flows follow extended self-similarity, while the dependence of the intermittency factor on real-fluid effects is observed for high-order structure functions. The real-fluid effects on intermittency are explained by turbulent structures related to rare events. Additionally, the dissipation rate moments for transcritical flows follow a universal scaling with Reynolds number, and the scaling exponents are different from those of incompressible flows. These results extend the small-scale universality in incompressible flows (Schumacher et al., Proc. Natl Acad. Sci. USA, vol. 111, 2014, pp. 10961–10965) to realistic flows with significant changes in thermodynamic properties, and provide a physical underpinning of the scaling laws of small-scale statistics at transcritical conditions.

The spontaneous formation of zonal jets is a distinctive feature of geostrophic turbulence with the phenomenon witnessed in numerous numerical studies. In such systems, strong rotation anisotropises the spectral evolution of the energy density such that zonal modes are favoured. In physical space, this manifests as eddies zonally elongating and forming into zonal jets. In the presence of large scale dissipation, the flow may reach statistical stationarity such that the zonal structure persists in the zonal and time mean, and is supported by a flux of eddy momentum. What is unclear is how the excitation of Rossby waves arranges the underlying eddy momentum stresses to support the mean flow structures. To study this, we examine a steady-state flow in the so-called ‘zonostrophic’ regime, in which characteristic scales of geostrophic turbulence are well separated and there are several alternating zonal jets that have formed spontaneously. We apply a geometric eddy ellipse formulation, in which momentum fluxes are cast as ellipses that encode information about the magnitude and direction of flux; the latter is described using the tilt angle. With the aid of a zonal filter, it is revealed that the scales responsible for providing the momentum fluxes associated with the jet structure are much smaller than the characteristic scales identified, and occupy a region of the energy spectrum that has been typically associated with isotropic dynamics.

The aeroelasticity of a panel in the presence of a shock is a fundamental issue of great significance in the development of hypersonic vehicles. In practical engineering, cavity pressure emerges as a crucial factor that influences the nonlinear dynamical characteristics of the panel. This study focuses on the aeroelastic bifurcation of a flexible panel subjected to both cavity pressure and oblique shock. To this end, a computational method is devised, coupling a high-fidelity reduced-order model for unsteady aerodynamic loads with nonlinear structural equations. The solution is meticulously tracked by continuous calculations. The obtained results indicate that cavity pressure plays a pivotal role in determining the bifurcation and stability characteristics of the system. First, the system exhibits hysteresis behaviour in response to the ascending and descending dynamic pressures. The evolution of hysteresis behaviour originates from the phenomenon of cusp catastrophe. Second, variations in cavity pressure induce three types of bifurcation phenomena, exhibiting characteristics akin to supercritical Hopf bifurcation, subcritical Hopf bifurcation and saddle-node bifurcation of cycles. The system's response at the critical points of these bifurcations manifests as long-period asymptotic flutter or explosive flutter. Lastly, the evolution of the dynamical system among these three types of bifurcations is an important factor contributing to the discrepancies observed in certain research results. This study enhances the understanding of the nonlinear dynamical behaviour of panel aeroelasticity in complex practical environments and provides new explanations for the discrepancies observed in certain research results.

A Marangoni surfer is an object embedded in a gas–liquid interface, propelled by gradients in surface tension. We derive an analytical theorem for the lower bound on the viscous dissipation by a Marangoni surfer in the limit of small Reynolds and capillary numbers. The minimum dissipation can be expressed with the reciprocal difference between drag coefficients of two passive bodies of the same shape as the Marangoni surfer, one in a force-free interface and the other in an interface with surface incompressibility. The distribution of surface tension that gives the optimal propulsion is given by the surface tension of the solution for the incompressible surface and the flow is a superposition of both solutions. For a surfer taking the form of a thin circular disk, the minimum dissipation is $16\mu a V^2$, giving a Lighthill efficiency of $1/3$. This places the Marangoni surfers among the hydrodynamically most efficient microswimmers.

The interaction of opposite-signed twist waves on vortex tubes can lead to vortex bursting, a process where the core expands into a double ring-like structure with strong swirling flows. Previous works have studied vortex bursting on rectilinear vortices by axially perturbing the initial core size to generate the twist waves, and observed largely axisymmetric bursting dynamics. In this work, we numerically study bursting on vortical structures with curved centrelines, analysing the interaction between the centreline dynamics, twist wave generation and propagation, and vortex bursting. We focus on axially perturbed helical vortex tubes with small radius-to-pitch ratios up to $0.0625$, as well as vortex rings with a large radius-to-core size ratio $10$, both at a circulation-based Reynolds number $5000$. The results show that though the initial twist wave propagation speeds are relatively unaffected by the curvature and torsion of the centreline, the bursting process is altered significantly compared with rectilinear vortices. The self-induced rotation of the centreline of the helical tube induces a non-axisymmetric distortion of the bursting structure, which rapidly breaks up the vortex core into small-scale helical structures. A similar destabilization of the bursting structure also occurs on vortex rings. The enstrophy increase and accelerated energy decay associated with bursting are predominantly determined by the twist wave strength, rather than the curvature and torsion of the centreline. Combined, our findings imply that bursting could play an important role in transferring and dissipating energy of vortical structures in wakes, and turbulent flows in general.

Shear-thinning viscosity is a non-Newtonian behaviour that active particles often encounter in biological fluids such as blood and mucus. The fundamental question of how this ubiquitous non-Newtonian rheology affects the propulsion of active particles has attracted substantial interest. In particular, spherical Janus particles driven by self-diffusiophoresis, a major physico-chemical propulsion mechanism of synthetic active particles, were shown to always swim slower in a shear-thinning fluid than in a Newtonian fluid. In this work, we move beyond the spherical limit to examine the effect of particle eccentricity on self-diffusiophoretic propulsion in a shear-thinning fluid. We use a combination of asymptotic analysis and numerical simulations to show that shear-thinning rheology can enhance self-diffusiophoretic propulsion of a spheroidal particle, in stark contrast to previous findings for the spherical case. A systematic characterization of the dependence of the propulsion speed on the particle's active surface coverage has also uncovered an intriguing feature associated with the propulsion speeds of a pair of complementarily coated particles not previously reported. Symmetry arguments are presented to elucidate how this new feature emerges as a combined effect of anisotropy of the spheroidal geometry and nonlinearity in fluid rheology.

We investigate the influence of vortices remote from the boundary on the near-wall flow dynamics in wall-bounded flows. A vortex ring with precisely controlled local twist is introduced into the outer layer of a channel flow at a moderate Reynolds number. We find that the minimum vorticity flux for triggering the transition to turbulence is significantly reduced from the initial disturbance of an untwisted vortex ring to that of a twisted ring. In particular, the latter disturbance can cause vortex bursting in the early transitional stage. The impact of vortex bursting on the transition process is characterised by the near-wall, wall-normal velocity with the rapid distortion theory. The wall-normal velocity grows during vortex bursting, and leads to streak formation and then the transition to turbulence. The notable wall-normal velocity is induced by the large di-vorticity generated in vortex bursting. We model the growing radial component of the di-vorticity in terms of the local twist, and demonstrate that its surge is due to the generation of highly twisted vortex lines in vortex bursting. Then, we derive that the generation of the di-vorticity in the outer layer enhances the wall-normal velocity in the inner layer via the Poisson equation with the image method and the multipole expansion. Thus, we elucidate that the vortex bursting can have an effect on the transition process.