Droplet spreading is ubiquitous and plays a significant role in liquid-based energy systems, thermal management devices and microfluidics. While the spreading of non-volatile droplets is quantitatively understood, the spreading and flow transition in volatile droplets remains elusive due to the complexity added by interfacial phase change and non-equilibrium thermal transport. Here we show, using both mathematical modelling and experiments, that the wetting dynamics of volatile droplets can be scaled by the spatial–temporal interplay between capillary, evaporation and thermal Marangoni effects. We elucidate and quantify these complex interactions using phase diagrams based on systematic theoretical and experimental investigations. A spreading law of evaporative droplets is derived by extending Tanner's law (valid for non-volatile liquids) to a full range of liquids with saturation vapour pressure spanning from $10^1$ to $10^4$ Pa and on substrates with thermal conductivity from $10^{-1}$ to $10^3\ {\rm W}\ {\rm m}^{-1}\ {\rm K}^{-1}$. In addition to its importance in fluid-based industries, the conclusions also enable a unifying explanation to a series of individual works including the criterion of flow reversal and the state of dynamic wetting, making it possible to control liquid transport in diverse application scenarios.

Understanding the dynamics of staircases in salt fingering convection presents a long-standing theoretical challenge to fluid dynamicists. Although there has been significant progress, particularly through numerical simulations, there are a number of conflicting theoretical explanations as to the driving mechanism underlying staircase formation. The Phillips effect proposes that layering in stirred stratified flow is due to an antidiffusive process, and it has been suggested that this mechanism may also be responsible for salt fingering staircases. However, the details of this process, as well as mathematical models to predict the evolution and merger dynamics of staircases, have yet to be established. We generalise the theory of the Phillips effect to a three-component system (e.g. temperature, salinity, energy) and demonstrate a regularised nonlinear model of layering based on mixing length parametrisations. The model predicts both the inception of layering and its long-term evolution through mergers, while generalising, and remaining consistent with, previous results for double-diffusive layering based on flux ratios. Our model of salt fingering is formulated using spatial averaging processes, and closed by a mixing length parametrised in terms of the kinetic energy and the ratio of the temperature and salt gradients. The model predicts a layering instability for a bounded range of parameter values in the salt fingering regime. Nonlinear solutions show that an initially unstable linear buoyancy gradient develops into layers, which proceed to merge through a process of stronger interfaces growing at the expense of weaker ones. Our results indicate that these mergers are responsible for the characteristic increase of buoyancy flux through thermohaline staircases.

Theoretical analysis of attachment-line instabilities is performed for supersonic swept flows using the compressible Hiemenz approximation for the mean flow and the successive approximation procedures for disturbances. The theoretical model captures the dominant attachment-line modes in wide ranges of the sweep Mach number ${M_e}$ and the wall temperature ratio. It is shown that these modes behave similar to the first and second Mack modes in the boundary layer flow. This similarity allows us to extrapolate the knowledge gained for Mack modes to the attachment-line instabilities. In particular, we find that at sufficiently large ${M_e}$, the dominant attachment-line instability is associated with the synchronisation of slow and fast modes of acoustic nature. Point-by-point comparisons of the theoretical predictions with the experiments of Gaillard et al. (Exp. Fluids, vol. 26, 1999, pp. 169–176) demonstrate that at ${M_e} > 4$, the theory captures a significant drop of the transition onset Reynolds number, which is below the contamination criterion of Poll $({R_\mathrm{\ast }} = 250)$ at ${M_e} > 6$. This contradicts the generally accepted assumption that the attachment-line flow is stable for ${R_\mathrm{\ast }} \le 250$. The theoretical critical Reynolds numbers lie well below the experimental transition-onset Reynolds numbers. Stability computations using the Navier–Stokes mean flow and accounting for the leading-edge curvature effect do not eliminate this discrepancy. Most likely, in the experiments of Gaillard et al., we face with an unknown effect that does not fit to the concept of transition arising from linear instability.

Shear-induced migration of elongated micro-swimmers exhibiting anisotropic Brownian diffusion at a population scale is investigated analytically in this work. We analyse the steady motion of confined ellipsoidal micro-swimmers subject to coupled diffusion in a general setting within a continuum homogenisation framework, as an extension of existing studies on macro-transport processes, by allowing for the direct coupling of convection and diffusion in local and global spaces. The analytical solutions are validated successfully by comparison with numerical results from Monte Carlo simulations. Subsequently, we demonstrate from the probability perspective that symmetric actuation does not yield net vertical polarisation in a horizontal flow, unless non-spherical shapes, external fields or direct coupling effects are harnessed to generate steady locomotion. Coupled diffusivities modify remarkably the drift velocity and vertical migration of motile micro-swimmers exposed to fluid shear. The interplay between stochastic swimming and preferential alignment could explain the diverse concentration and orientation distributions, including rheological formations of depletion layers, centreline focusing and surface accumulation. Results of the analytical study shed light on unravelling peculiar self-propulsion strategies and dispersion dynamics in active-matter systems, with implications for various transport problems arising from the fluctuating shape, size and other external or inter-particle interactions of swimmers in confined environments.

A methodology to generate sparse Galerkin models of chaotic/unsteady fluid flows containing a minimal number of active triadic interactions is proposed. The key idea is to find an appropriate set of basis functions for the projection representing elementary flow structures that interact minimally one with the other, resulting in a triadic interaction coefficient tensor with sparse structure. Interpretable and computationally efficient Galerkin models are obtained, since a reduced number of triadic interactions are computed to evaluate the right-hand side of the model. To find the basis functions, a subspace rotation technique is used, whereby a set of proper orthogonal decomposition (POD) modes is rotated into a POD subspace of larger dimension using coordinates associated with low-energy dissipative scales to alter energy paths and the structure of the triadic interaction coefficient tensor. This rotation is obtained as the solution of a non-convex optimisation problem that maximises the energy captured by the new basis, promotes sparsity and ensures long-term temporal stability of the sparse Galerkin system. We demonstrate the approach on two-dimensional lid-driven cavity flow at $Re = 2 \times 10^4$ where the motion is chaotic and energy interactions are scattered in modal space. We show that the procedure generates Galerkin models with a reduced set of active triadic interactions, distributed in modal space according to established knowledge of scale interactions in two-dimensional flows. This property, however, is observed only if long-term temporal stability is included explicitly in the formulation, indicating that a dynamical constraint is necessary to obtain a physically consistent sparsification.

Bioprosthetic heart valves create turbulent flow during early systole which might be detrimental to their durability and performance. Complex mechanisms in the unsteady and heterogeneous flow field complicate the isolation of specific instability mechanisms. We use linear stability analysis and numerical simulations of the flow in a simplified model to study mechanisms initiating the laminar–turbulent transition. The analysis of a modified Orr–Sommerfeld equation, which includes a model for fluid–structure interaction (FSI), indicates Kelvin–Helmholtz and FSI instabilities for a physiological Reynolds number regime. Two-dimensional parametrized FSI simulations confirm the growth rates and phase speeds of these instabilities. The eigenmodes associated with the observed leaflet kinematics allow for decoupled leaflet oscillations. A detailed analysis of the temporal evolution of the flow field shows that the starting vortex interacts with the aortic wall leading to a secondary vortex which moves towards the shear layer in the wake of the leaflets. This appears to be connected to the onset of the shear layer instabilities that are followed by the onset of leaflet motion leading to large-scale vortex shedding and eventually to a nonlinear breakdown of the flow. Numerical results further indicate that a narrower aorta leads to an earlier onset of the shear layer instabilities. They also suggest that the growing perturbations of the shear layer instability propagate upstream and may initiate the FSI instabilities on the valve leaflets.

Steady shock reflection where the incident shock is free of interaction with other waves has been well studied. In this paper, we consider the less studied shock reflection problem where the incident shock interacts with the wedge trailing-edge expansion fan, which occurs when the wedge trailing-edge height surpasses a threshold. The influence of this interaction on the advance of transition from Mach reflection to regular reflection is quantified in terms of the wedge trailing-edge height ratio. The wave pattern, including primary and reflected Mach waves, for Mach reflection with interaction is clarified using computational fluid dynamics (CFD) and the method of characteristics. Those reflected Mach waves having an important effect on Mach reflection are identified. A simplified Mach stem model that accounts for the direct role of the interaction on the incident shock and its indirect role on the reflected shock and slipline is built up on a past model without interaction. Both theory and CFD show that the Mach stem height decreases nonlinearly with increasing trailing-edge height.

Shock waves are of great interest in many fields of science and engineering, but the mechanisms of their formation, maintenance and dissipation are still not well understood. While all transport processes existing in a shock wave contribute to its compression and irreversibility, they are not of equal importance. To figure out the roles of viscosity and heat conduction in shock transition, the existence of smooth shock solutions and the counter-intuitive entropy overshoot phenomenon (the specific entropy is not monotonically increasing and exhibits a peak inside the shock front) are theoretically and numerically investigated, with emphasis on the effects of viscosity and heat conduction. Instead of higher-order hydrodynamics, the Navier–Stokes formalism is employed for its stability and simplicity. Supplemented with nonlinear thermodynamically consistent constitutive relations, the Navier–Stokes equations are adequate to demonstrate the general nature of shock profiles. It is found that heat conduction cannot sustain strong shocks without the presence of viscosity, while viscosity can maintain smooth shock transition at all strengths, regardless of heat conduction. Hence, the critical role in shock compression is played by viscosity rather than heat conduction. Nevertheless, the dispensability of heat conduction would not compromise its essential role in the emergence of an entropy peak. It is the entropy flux resulting from heat conduction that neutralises the positive entropy production and thus prevents the decreasing entropy from violating the second law of thermodynamics. This mechanism of entropy overshoot has not been addressed previously in the literature and is revealed using the entropy balance equation.