This paper investigates the linear and nonlinear evolution of radiating modes under the influence of the spontaneously emitted Mach waves in a simple set-up of the supersonic boundary layers that develop in the entry region of a channel formed by two parallel semi-infinite flat plates. Two scenarios are considered. The first occurs in the boundary layers having identical wall conditions, where the Mach wave emitted by a radiating mode in one boundary layer influences the instability in the other. The second scenario takes place when the wall temperatures are different, in which case the spontaneously radiated Mach wave is reflected by the other boundary layer back to act on the radiating mode. Appropriate amplitude equations with the acoustic feedback effect being accounted for are derived. In each case, the effect of the spontaneously emitted sound contributes a linear term of delay type to the respective amplitude equation. For the first scenario, analytical and numerical studies of the amplitude equations show that due to the back action of the spontaneously radiated Mach wave, the amplitude exhibits rapid oscillations, and in the case of enhanced feedback effects, its envelope experiences near extinction followed by resurrection. The study of the coupled equations shows that the two modes with different initial amplitudes either undergo oscillations before attenuating, or terminate a finite-distance singularity at different locations. For the second scenario, the acoustic feedback produces similar effects in a broad range of wall temperature. The effects become pronounced, and the dependence on the wall temperature becomes more sensitive when the latter approaches the value corresponding to the resonance. Estimates suggest that such acoustic feedback is likely to be present in typical wind tunnel experiments and models for scramjet combustors.

Using laboratory experiments and numerical simulations, we examine the transfer of soluble material from small, spherical particles sinking in homogeneous turbulence at large Péclet number. A theoretical analysis predicts two distinct mechanisms of convective mass transfer: strain due to turbulence and slip due to gravitational settling. Their relative strength is parametrised by the sinking ratio, ${Sr} = w_0 \tau _\eta /a$, where $w_0$ is the quiescent settling velocity, $a$ is the particle radius and $\tau _\eta$ is the Kolmogorov time scale. This analysis predicts that the topology of the concentration wake changes from a symmetric topology at ${Sr} \ll 1$ to an asymmetric topology at ${Sr} \gg 1$ as the dominant mechanism of mass transfer changes. Particle tracking flow visualisations of small spheres releasing dye in turbulence confirm the existence of this change in mechanism at ${Sr} = O(1)$. We complement these experiments with numerical simulations of the mass transfer from sinking particles. The transfer rate predicted by the simulation is found to be in good agreement with literature data for mass transfer to turbulent suspensions of solid particles and is consistent with asymptotic expressions for mass transfer in uniform flow when ${Sr} \gg 1$. A decomposition of the convective fluxes confirms the transition in the transfer mechanism. At ${Sr} = O(1)$, both mechanisms provide comparable contributions to the transfer rate. Cross-correlation analysis reveals that particle-scale knowledge of both the recent strain and velocity history is required to predict the instantaneous transfer rate. Turbulence-induced particle rotation has a modest suppression effect upon convective transfer by sinking.

A surfactant-laden liquid film that lines the inside of an oscillating spherical cap is considered as a model of lung alveoli. Pulmonary surfactant solubility is described by Langmuir adsorption kinetics, modified by incorporating the intrinsic compressibility of the adsorbed monolayer. A novel boundary condition, supported by experimental data and scaling arguments, is applied at the rim. The condition enforces mass conservation of water and surfactant by matching the ‘large-scale’ dynamics of the alveolus to ‘small-scale’ equilibrium over mid-alveolar septa of small but finite thickness. Linear and weakly nonlinear analysis around the conditions in a non-oscillating cap indicates that the occurrence of shearing motion in the liquid is related to the non-zero film thickness over the rim, and shearing velocity at the interface is predicted an order-of-magnitude lower than the velocity of radial oscillation. Marangoni stresses dominate the interfacial dynamics, but capillary stresses affect significantly the interior flow field. In particular, they produce spatial modulations in flow rate, surface concentration of surfactant and wall shear stress, whose length scale varies with $Ca^{-1/3}$, i.e. is determined by a balance between capillary and viscous forces. Non-zero adsorption kinetics modifies at first order only the amplitude and phase of surface concentration, but affects all other variables at second order. In particular, it sets a steady drift of surfactant away from the alveolus and towards the rim. Finally, an attempt is made to relate the present predictions to physiological findings about air flow and particle deposition inside alveoli, and about shear stress-inflicted damage in diseased lungs.

This study presents a noise-robust closed-loop control strategy for wake flows employing model predictive control. The proposed control framework involves the autonomous offline selection of hyperparameters, eliminating the need for user interaction. To this purpose, Bayesian optimization maximizes the control performance, adapting to external disturbances, plant model inaccuracies and actuation constraints. The noise robustness of the control is achieved through sensor data smoothing based on local polynomial regression. The plant model can be identified through either theoretical formulation or using existing data-driven techniques. In this work we leverage the latter approach, which requires minimal user intervention. The self-tuned control strategy is applied to the control of the wake of the fluidic pinball, with the plant model based solely on aerodynamic force measurements. The closed-loop actuation results in two distinct control mechanisms: boat tailing for drag reduction and stagnation point control for lift stabilization. The control strategy proves to be highly effective even in realistic noise scenarios, despite relying on a plant model based on a reduced number of sensors.

Horizontal convection at large Rayleigh and Prandtl numbers is studied experimentally in a regime up to seven orders of magnitude larger in terms of Rayleigh numbers than previously achieved. To reach Rayleigh numbers up to $10^{17}$, the horizontal density gradient is generated using differential solutal convection by a differential input of salt and fresh water controlled by diffusion in a novel experiment in which the zero-net mass flux of water is ensured through permeable membranes. This set-up allows us to accurately measure the Nusselt number in solutal convection by carefully controlling the amount of salt water exchanged through the membranes. Combined measurements of density and velocity across more than five orders of magnitude in Rayleigh numbers show that the flow transitions from the Beardsley & Festa (J. Phys. Oceanogr., vol. 2, issue 4, 1972, pp. 444–455), Shishkina & Wagner (Phys. Rev. Lett., vol. 116, issue 2, 2016, 024302) regime to the Chiu-Webster et al. (J. Fluid Mech., vol. 611, 2008, pp. 395–426) regime and frames the present results within the scope of Shishkina et al. (Geophys. Res. Lett., vol. 43, issue 3, 2016, pp. 1219–1225), and the theory of Part 1 (Passaggia & Scotti, vol. 997, 2024, J. Fluid Mech., A5). In particular, we show that, even for large Prandtl numbers, the circulation eventually clusters underneath the forcing horizontal boundary, leaving a stratified core without motion. Finally, the previous regime diagrams (Hughes & Griffiths, Annu. Rev. Fluid Mech., vol. 40, 2008, pp. 185–208; Shishkina et al., Geophys. Res. Lett., vol. 43, issue 3, 2016, pp. 1219–1225) are extended by combining the present results at high Prandtl numbers, the results at low Prandtl numbers of Part 1, together with previous results from the literature. This work sets a new picture of the transition landscape of horizontal convection over six orders of magnitude in Prandtl number and sixteen orders of magnitude in Rayleigh number.

We show that rotating Rayleigh–Bénard convection, where a rotating fluid is heated from below, exhibits a non-Hermitian topological invariant. Recently, Favier & Knobloch (J. Fluid Mech., vol. 895, 2020, R1) hypothesized that the robust sidewall modes in rapidly rotating convection are topologically protected. By considering a Berry curvature defined in the complex wavenumber space, we reveal that the bulk states can be characterized by a non-zero integer Chern number, implying a potential topological origin of the edge modes based on the Atiyah–Patodi–Singer index theorem (Fukaya et al., Phys. Rev. D, vol. 96 2017, 125004; Yu et al., Nucl. Phys. B, vol. 916, 2017, pp. 550–566). The linearized eigenvalue problem is intrinsically non-Hermitian, therefore, the definition of Berry curvature generalizes that of the stably stratified problem. Moreover, the three-dimensional set-up naturally regularizes the eigenvector, avoiding the compactification problem in shallow water waves (Tauber et al., J. Fluid Mech., vol. 868, 2019, R2). Under the hydrostatic approximation, it recovers a two-dimensional analogue of the one which explains the topological origin of the equatorial Kelvin and Yanai waves (Delplace et al., Science, vol. 358, issue 6366, 2017, pp. 1075–1077). The non-zero Chern number relies only on rotation when the fluid is stratified, no matter whether it is stable or unstable. However, the neutrally stratified system does not support a topological invariant. In addition, we define a winding number to visualize the topological nature of the fluid. Our results represent a step forward for the topologically protected states in convection, but the bulk-boundary correspondence requires a further direct analysis for proof, and the robustness of the edge states under varying boundary conditions remains a question to be answered.

As most mathematically justifiable Lagrangian coherent structure detection methods rely on spatial derivatives, their applicability to sparse trajectory data has been limited. For experimental fluid dynamicists and natural scientists working with Lagrangian trajectory data via passive tracers in unsteady flows (e.g. Lagrangian particle tracking or ocean buoys), obtaining material measures of fluid rotation or stretching is an active topic of research. To facilitate frame-indifferent investigations in unsteady and sparsely sampled flows, we present a novel approach to quantify fluid stretching and rotation via relative Lagrangian velocities. This technique provides a formal objective extension of quasi-objective metrics to unsteady flows by accounting for mean flow behaviour. For extremely sparse experimental data, fluid structures may be significantly undersampled and the mean flow behaviour becomes difficult to quantify. We provide a means to maintain the accuracy of our novel sparse flow diagnostics in extremely sparse sampling scenarios, such as ocean buoy data and Lagrangian particle tracking. We use data from multiple numerical and experimental flows to show that our methods can identify structures beyond existing limits of sparse, frame-indifferent diagnostics and exhibit improved interpretability over common frame-dependent diagnostics.

The dynamics of the interaction of a system of two thin volatile liquid droplets resting on a soft viscoelastic solid substrate are investigated theoretically. The developed model fully considers the effect of evaporative cooling and the generated Marangoni stresses due to the induced thermal gradients, while also accounting for the effect of the gas phase composition and the diffusion of vapour in the atmosphere of the droplets. Using the framework of lubrication theory, we derive evolution equations for both the droplet profile and the displacement of the elastic solid, which are solved in combination with Laplace's equation for the vapour concentration in the gas phase. A disjoining-pressure/precursor-film approach is used to describe contact-line motion. The evolution equations are solved numerically, using the finite-element method, and we present a thorough parametric analysis to investigate the physical properties and mechanisms that affect the dynamics of droplet interactions. The results show that the droplets interact through both the soft substrate and the gas phase. In the absence of thermocapillary phenomena, the combined effect of non-uniform evaporation due to the increased vapour concentration between the two droplets and elastocapillary phenomena determines whether the drop–drop interaction is attractive or repulsive. The Marangoni stresses suppress droplet attraction at the early stages of the drying process and lead to longer droplet lifetimes. For substrates with intermediate stiffness, the emergence of spontaneous symmetry breaking at late stages of evaporation is found. The rich dynamics of this complex system is explored by constructing a detailed map of the dynamic regimes.

Simulating complex gas flows from turbulent to rarefied regimes is a long-standing challenge, since turbulence and rarefied flow represent contrasting extremes of computational aerodynamics. We propose a multiscale method to bridge this gap. Our method builds upon the general synthetic iterative scheme for the mesoscopic Boltzmann equation, and integrates the $k$–$\omega$ model in the macroscopic synthetic equation to address turbulent effects. Asymptotic analysis and numerical simulations show that the macroscopic–mesoscopic coupling adaptively selects the turbulence model and the laminar Boltzmann equation. The multiscale method is then applied to opposing jet problems in hypersonic flight surrounding by rarefied gas flows, showing that the turbulence could cause significant effects on the surface heat flux, which cannot be captured by the turbulent model nor the laminar Boltzmann solution alone. This study provides a viable framework for advancing understanding of the interaction between turbulent and rarefied gas flows.

Vortex rings are critical for thrust production underwater. In the ocean, self-propelled mesozooplankton generate vortices while swimming within a weakly stratified fluid. While large-scale biogenic transport has been observed during vertical migration in the wild and lab experiments, little focus has been given to the evolution of induced vortex rings as a function of their propagation direction relative to the density gradient. In this study, the evolution of an isolated vortex ring crossing the interface of a stable two-layer system is examined as a function of its translation direction with respect to gravity. The vortex ring size and position are visualized using planar laser-induced fluorescence (PLIF) and the induced vorticity field derived from particle image velocimetry (PIV) is examined. It is found that the production of baroclinic vorticity significantly affects the propagation of vortex rings crossing the density interface. As a result, any expected symmetry between vortex rings travelling from dense to light fluids and from light to dense fluids breaks down. In turn, the maximum penetration depth of the vortex ring occurs in the case in which the vortex propagates against the density gradient due to the misalignment of the pressure and density gradients. Our results have far-reaching implications for the characterization of local ecosystems in marine environments.

The linear stability of the incompressible flow in an infinitely extended cavity with rectangular cross-section is investigated numerically. The basic flow is driven by a lid which moves tangentially, but at yaw with respect to the edges of the cavity. As a result, the basic flow is a superposition of the classical recirculating two-dimensional lid-driven cavity flow orthogonal to a wall-bounded Couette flow. Critical Reynolds numbers computed by linear stability analysis are found to be significantly smaller than data previously reported in the literature. This finding is confirmed by independent nonlinear three-dimensional simulations. The critical Reynolds number as a function of the yaw angle is discussed for representative aspect ratios. Different instability modes are found. Independent of the yaw angle, the dominant instability mechanism is based on the local lift-up process, i.e. by the amplification of streamwise perturbations by advection of basic flow momentum perpendicular to the sheared basic flow. For small yaw angles, the instability is centrifugal, similar as for the classical lid-driven cavity. As the spanwise component of the lid velocity becomes dominant, the vortex structures of the critical mode become elongated in the direction of the bounded Couette flow with the lift-up process becoming even more important. In this case the instability is made possible by the residual recirculating part of the basic flow providing a feedback mechanism between the streamwise vortices and the streamwise velocity perturbations (streaks) they promote. In the limit when the basic flow approaches bounded Couette flow the critical Reynolds number increases very strongly.

We present the general analytical solution of the Riemann problem (decay of a jump discontinuity) for non-convex relativistic hydrodynamics. In convex dynamics, an elementary nonlinear wave, i.e. a rarefaction or a shock, originates at the discontinuity and travels towards one of the initial states. Between the left and right waves, an equilibrium state appears represented by a contact discontinuity. The exact solution to the Riemann problem in convex relativistic hydrodynamics was first addressed by Martí & Müller (J. Fluid Mech., vol. 258, 1994, pp. 317–333). In non-convex dynamics, two sequences of elementary nonlinear waves move towards the left and right initial states. Solving the Riemann problem involves determining the types of wave developing and the equilibrium state where they coincide. The procedure consists of constructing the wave curves associated with the nonlinear waves in the pressure–velocity phase space, where the intersection of the wave curves indicates the equilibrium state. We describe the relation between the wave curves, the explicit formulas for their calculation, and the outline of the process for a correct derivation and representation of the waves in the spatial domain. We present examples of the exact solution of a Riemann problem that illustrate the complex phenomena of non-convex dynamics by using the phenomenological non-convex equation of state proposed by Ibáñez et al. (Mon. Not. R. Astron. Soc., vol. 476, 2017, pp. 1100–1110).

We recently lost a valuable colleague and friend, Paul Steen, and the text below is a brief overview of his life.

The present study offers a twofold contribution on counter-gradient transport (CGT) of turbulent scalar flux. First, by examining turbulent scalar mixing through synchronized particle image velocimetry and planar laser-induced fluorescence on an inclined jet in cross-flow, we clarify the previously unexplained phenomenon of CGT, revealing key flow structures, their spatial distribution and modelling implications. Statistical analysis identifies two distinct CGT regions: local cross-gradient transport in the windward shear layer and non-local effects near the wall after injection. These behaviours are driven by specific flow structures, namely Kelvin–Helmholtz vortices (local) and wake vortices (non-local), suggesting that scalar flux can be decomposed into a gradient-type term for gradient diffusion and a term for large-eddy stirring. Second, we propose a new approach for reconstruction of turbulent mean flow and scalar fields using continuous adjoint data assimilation (DA). By rectifying model-form errors through anisotropic correction under observational constraints, our DA model minimizes discrepancies between experimental measurements and numerical predictions. As expected, the introduced forcing term effectively identifies regions where traditional models fall short, particularly in the jet centreline and near-wall regions, thereby enhancing the accuracy of the mean scalar field. These enhancements occur not only within the observation region but also in unseen regions, underscoring present DA approach's reliability and practicality for reproducing mean flow behaviours from limited data. These findings lay a solid foundation for adjoint-based model-consistent data-driven methods, offering promising potential for accurately predicting complex flow scenarios like film cooling.