Coherent small-amplitude unsteadiness of the shock wave and the separation region over a canonical double cone flow, termed in literature as oscillation-type unsteadiness, is experimentally studied at Mach 6. The double cone model is defined by three non-dimensional geometric parameters: fore- and aft-cone angles ($\theta _1$ and $\theta _2$), and ratio of the conical slant lengths ($\varLambda$). Previous studies of oscillations have been qualitative in nature, and mostly restricted to a special case of the cone model with fixed $\theta _1 = 0^\circ$ and $\theta _2 = 90^\circ$ (referred to as the spike-cylinder model), where $\varLambda$ becomes the sole governing parameter. In the present effort we investigate the self-sustained flow oscillations in the $\theta _1$-$\varLambda$ parameter space for fixed $\theta _2 = 90^\circ$ using high-speed schlieren visualisation. The experiments reveal two distinct subtypes of oscillations, characterised by the motion (or lack thereof) of the separation point on the fore-cone surface. The global time scale associated with flow oscillation is extracted using spectral proper orthogonal decomposition. The non-dimensional frequency (Strouhal number) of oscillation is seen to exhibit distinct scaling for the two oscillation subtypes. The relationship observed between the local flow properties, instability of the shear layer, and geometric constraints on the flow suggests that an aeroacoustic feedback mechanism sustains the oscillations. Based on this understanding, a simple model with no empiricism is developed for the Strouhal number. The model predictions are found to match well with experimental measurements. The model provides helpful physical insight into the nature of the self-sustained flow oscillations over a double cone at high speeds.

Emptying–filling boxes have been studied in a wide range of configurations for decades, but the flow created in the box by two plumes rising from sources of arbitrary strength and elevation was previously unsolved. Guided by experiments and simplified analytical modelling, we reveal a rich array of two- and three-layer stratifications across seven possible flow regimes. The governing equations for these regimes show how the prevailing regime and stratification properties vary with three key parameters: the relative strength of the plumes, the height difference between their sources and a parameter characterising the resistance of the box to emptying. We observe and explain new behaviours not described in previous studies that are crucial to understanding emptying–filling boxes with multiple plumes. In particular, we demonstrate that the oft-assumed premise that $n$ plumes leads to a stratification with $n+1$ layers is not necessarily true, even in the absence of mixing. Two emptying–filling box models are developed: an analytical model addressing all combinations of the governing parameters and an extended model for three-layer stratifications that incorporates two mixing processes observed in the experiments. The predictions of these two models are in generally excellent agreement with measurements from the experimental campaign covering 69 combinations of the governing parameters. This study improves our understanding of emptying–filling boxes and could facilitate improvements to natural ventilation building design, as demonstrated by an example scenario in which occupants feel cooler upon the addition of a second source of heat.

The non-Oberbeck–Boussinesq effects on the stability of a vertical natural convection boundary layer are investigated using the linearised disturbance equations for air flows up to a temperature difference of $\Delta T=100\,{\rm K}$. Based on the linear stability results, the neutral curve is shown to be sensitive to the choice of reference temperature. When evaluated using the film temperature $T_f$, a lower film Grashof number is required to trigger the linear instability for larger $\Delta T$. The relative contributions of shear and buoyant production to the perturbation kinetic energy budget reveals that the marginally unstable modes are amplified based on different mechanisms: for lower wavenumbers at relatively small Grashof number, the instability is driven by buoyancy; whereas for higher wavenumbers and larger Grashof number, the flow becomes unstable due to a shear instability. The use of reference temperature is found to scale the shear- and buoyant-driven instabilities differently so that no single reference temperature definition would collapse the neutral curves. The linear stability result further demonstrates that at a given Grashof number a higher temperature difference would give a larger amplification rate of the perturbation, which then leads to an earlier onset of the nonlinearities when evaluated at $T_f$. Finally, by comparing the amplification rates obtained from direct numerical simulation and the linear stability results, the extent of the linear regime is determined for $\Delta T = 100\,{\rm K}$.

The evolution of the water-entry cavity affects the impact load and the motion of the body. This paper adopts the Eulerian finite element method for multiphase flow for simulations of the high-speed water-entry process. The accuracy and convergence of the numerical method are verified by comparing it with the experimental data and the results of the transient cavity dynamics theory. Based on the results, the representative characteristics of the cavity are discussed from the perspective of the cavity cross-section. It is found that the asymmetry of the cavity expansion and contraction durations is related to the motion of the free surface and the closure of the cavity. The uplift of the free surface suppresses cavity expansion, while the jet generated from free surface closure accelerates cavity contraction. The duration of the contraction of the cavity near the free surface is shorter than the expansion duration due to the change in the velocity distribution caused by the free surface motion. The necking phenomenon during deep closure leads to an increase in the internal pressure of the cavity, prolonging cavity contraction near the deep closure area. This work provides new insights into the cavity dynamics in high-speed water entry.

In this study we propose a novel data-driven reduced-order model for complex dynamics, including nonlinear, multi-attractor, multi-frequency and multiscale behaviours. The starting point is a fully automatable cluster-based network model (CNM) (Li et al., J. Fluid Mech., vol. 906, 2021, A21) that kinematically coarse grains the state with clusters and dynamically predicts the transitions in a network model. In the proposed dynamics-augmented CNM (dCNM) the prediction error is reduced with trajectory-based clustering using the same number of centroids. The dCNM is first exemplified for the Lorenz system and then demonstrated for the three-dimensional sphere wake featuring periodic, quasi-periodic and chaotic flow regimes. For both plants, the dCNM significantly outperforms the CNM in resolving the multi-frequency and multiscale dynamics. This increased prediction accuracy is obtained by stratification of the state space aligned with the direction of the trajectories. Thus, the dCNM has numerous potential applications to a large spectrum of shear flows, even for complex dynamics.

While the problem governing Stokes flow about a single particle that is subject to an external force is ill posed in two dimensions (the ‘Stokes paradox’), the related problem of two mutually repellent particles is well posed. Motivated by self-assembly phenomena in thin viscous membranes, we consider this problem in the limit of remote particles. Such limits are typically handled in the literature using reflection techniques, which provide successive approximations to the mutual hydrodynamic interactions. Since their starting point is a single particle in an unbounded fluid domain, these techniques are futile in the present two-dimensional problem. We show how this apparent contradiction is resolved via use of singular perturbations. We obtain a two-term approximation for the velocity acquired by circular disks, considering both rigid and free particle surfaces. We also illustrate our perturbation scheme for elliptic disks, deriving a renormalised single-particle velocity. The utility of our asymptotic scheme is illustrated in the general problem of hydrodynamic interaction between a cluster of remote disks.

The bubbly shock-driven partial cavitation in an axisymmetric venturi is studied with time-resolved two-dimensional X-ray densitometry. The bubbly shock waves are characterised using the vapour fraction and pressure changes across it, propagation velocity, and Mach number. The sharp changes in vapour fraction measured with X-ray densitometry, combined with high-frequency dynamic pressure measurements, reveal that the interaction of the pressure wave with the vapour cavity dictates the shedding dynamics. At the lowest cavitation number ($\sigma \sim 0.47$), the condensation shock front is the predominant shedding mechanism. However, as $\sigma$ increases ($\sigma \sim 0.78$), we observe an upstream travelling pressure discontinuity that changes into a condensation shock as it approaches the venturi throat. This coincides with the increasing strength of the bubbly shock wave as it propagates upstream, manifested by the increasing velocity of the shock front and the pressure rise across it. Consequently, the Mach number of the shock front increases and surpasses the critical value 1, favouring condensation shocks. Further, at higher $\sigma$ (${\sim }0.84\unicode{x2013}0.9$), both the re-entrant jet and pressure wave can cause cavity detachment. However, at such $\sigma$, the pressure wave likely remains subsonic. Hence cavity condensation is not favoured readily. This leads to the re-entrant jet causing the cavity detachment at higher $\sigma$. The shock front is accelerated as it propagates upstream through the variable cross-section of the venturi. This enhances its strength, favouring cavity condensation and eventual shedding. These observations explain the existence of shock fronts in an axisymmetric venturi for a large range of $\sigma$.

We explore predictions of two models of one-dimensional capillary rise in rigid and partially saturated porous media. One is an existing one from the literature and the second is a free-boundary model based on Richards’ equation with two moving boundaries of the evolving partially saturated region. Both models involve the specification of saturation-dependent functions for local capillary pressure and permeability and connect to classical models for saturated porous media. Existing capillary-rise experiments show two notable regimes: (i) an early-time regime typically well-described by classical capillary-rise theory in a fully saturated porous media, and (ii) a long-time regime that has anomalous dynamics in which the capillary-rise height may scale with a non-classical power law in time or have more complicated dynamics. We demonstrate that the predictions of both models compare well with experimental capillary-rise data over early- and long-time regimes gathered from three independent studies in the literature. The model predictions also shed light on recent scaling laws that relate the capillary pressure and permeability of the partially saturated media to the capillary-rise height. We use these models to probe computationally observed permeability relationships to capillary-rise height. We demonstrate that a recently proposed permeability scaling for the anomalous capillary-rise regime is indeed realized and is particularly apparent in the lower portion of the partially saturated media. For our free-boundary model we also compute capillary pressure measures and show that these reveal the linear relation between the capillary pressure and capillary-rise height expected for a capillarity–gravity balance in the upper portion of the partially saturated porous media.

We investigate two-dimensional (2-D) axisymmetric flow in toroidal geometry, with a focus on a transition between 2-D three-component flow and 2-D two-component flow. This latter flow state allows a self-organization of the system to a quiescent dynamics, characterized by long-living coherent structures. When these large-scale structures orient in the azimuthal direction, the radial transport is reduced. Such a transition, if it can be triggered in toroidally confined fusion plasmas, is beneficial for the generation of zonal flows and should consequently result in a flow field beneficial for confinement.

We present experiments on oscillating hydrofoils undergoing combined heaving and pitching motions, paying particular attention to connections between propulsive efficiency and coherent wake features extracted using modal analysis. Time-averaged forces and particle image velocimetry measurements of the flow field downstream of the foil are presented for a Reynolds number of $Re=11\times 10^3$ and Strouhal numbers in the range $St=0.16\unicode{x2013}0.35$. These conditions produce 2S and 2P wake patterns, as well as a near-momentumless wake structure. A triple decomposition using the optimized dynamic mode decomposition method is employed to identify dominant modal components (or coherent structures) in the wake. These structures can be connected to wake instabilities predicted using spatial stability analyses. Examining the modal components of the wake provides insightful explanations into the transition from drag to thrust production, and conditions that lead to peak propulsive efficiency. In particular, we find modes that correspond to the primary vortex development in the wakes. Other modal components capture elements of bluff body shedding at Strouhal numbers below the optimum for peak propulsive efficiency and characteristics of separation for Strouhal numbers higher than the optimum.

When a viscous fluid partially fills a Hele-Shaw channel, and is pushed by a pressure difference, the fluid interface is unstable due to the Saffman–Taylor instability. We consider the evolution of a fluid region of finite extent, bounded between two interfaces, in the limit that the interfaces are close, that is, when the fluid region is a thin liquid filament separating two gases of different pressure. In this limit, we derive a second-order ‘thin-filament’ model that describes the normal velocity of the filament centreline, and evolution of the filament thickness, as functions of the thickness, centreline curvature and their derivatives. We show that the second-order terms in this model, that include the effect of transverse flow along the filament, are necessary to regularise the instability. Numerical simulation of the thin-filament model is shown to be in accordance with level-set computations of the complete two-interface model. Solutions ultimately evolve to form a bubble of rapidly increasing radius and decreasing thickness.

A simple continuum model is proposed for discrete granular avalanches, as observed in slowly rotating drums. Subject to granular transfer across the basal interface, the model evolves the transient surface inclination and mid-slope velocity profile of the avalanches, from failure to arrest. This is done using new boundary conditions that allow entrainment or detrainment contingent on the basal shear rate. For entrainment to occur, the basal shear stress must overcome the erosion resistance of the jammed deposit, while for detrainment to take place the basal shear rate must vanish. The resulting avalanche dynamics is controlled by a single dimensionless number, the ratio of excess slope at failure to excess resistance to erosion, that can be determined from inclination history data. This number provides a measure of slope brittleness, from which the peak flow rate and arrest inclination can be predicted. Analytical techniques are used to derive model solutions, and clarify the resulting behaviour. Finally, the model is tested by comparing solutions with laboratory experiments and discrete particle simulations.

We explored the settling dynamics of vertically aligned particles in a quiescent, stratified two-layer fluid using particle tracking velocimetry. Glass spheres of $d=4\,{\rm mm}$ diameter were released at frequencies of 4, 6 and 8 Hz near the free surface, traversing through an upper ethanol layer ($H_1$), where H is height or layer thickess, varying from $10d$ to $40d$ and a lower oil layer. Results reveal pronounced lateral particle motion in the ethanol layer, attributed to a higher Galileo number ($Ga = 976$, ratio of buoyancy–gravity to viscous effects), compared with the less active behaviour in the oil layer ($Ga = 16$). The ensemble vertical velocity of particles exhibited a minimum just past the density interface, becoming more pronounced with increasing $H_1$, and suggesting that enhanced entrainment from ethanol to oil resulted in an additional buoyancy force. This produced distinct patterns of particle acceleration near the density interface, which were marked by significant deceleration, indicating substantial resistance to particle motion. An increased drag coefficient occurred for $H_1/d = 40$ compared with a single particle settling in oil; drag reduced as the particle-release frequency ($\,f_p$) increased, likely due to enhanced particle interactions at closer proximity. Particle pair dispersions, lateral ($R^2_L$) and vertical ($R^2_z$), were modulated by $H_1$, initial separation $r_0$ and $f_p$. The $R^2_L$ dispersion displayed ballistic scaling initially, Taylor scaling for $r_0 < H_1$ and Richardson scaling for $r_0 > H_1$. In contrast, $R^2_z$ followed a $R^2_z \sim t^{5.5}$ scaling under $r_0 < H_1$. Both $R^2_L$ and $R^2_z$ plateaued at a distance from the interface, depending on $H_1$ and $f_p$.

Droplet impact on oscillating substrates is important for both natural and industrial processes. Recognizing the importance of the dynamics that arises from the interplay between droplet transport and substrate motion, in this work, we present an experimental investigation of the spreading of a droplet impacting a sinusoidally oscillating hydrophobic substrate. We focus particularly on the maximum spread of droplets as a function of various parameters of substrate oscillation. We first quantify the maximum spreading diameter attained by the droplets as a function of frequency, amplitude of vibration, and phase at the impact for various impact velocities. We highlight that there can be two stages of spreading. Stage I, which is observed at all impact conditions, is controlled by the droplet inertia and affected by the substrate oscillation. For certain conditions, a Stage II spreading is also observed, which occurs during the retraction process of Stage I due to additional energies imparted by the substrate oscillation. Subsequently, we derive scaling analyses to predict the maximum spreading diameters and the time for this maximum spread for both Stage I and Stage II. Furthermore, we identify the necessary condition for Stage II spreading to be greater than Stage I spreading. The results will enable optimization of the parameters in applications where substrate oscillation is used to control the droplet spread, and thus heat and mass transfer between the droplet and the substrate.

Although the roll-streak (R-S) is fundamentally involved in the dynamics of wall turbulence, the physical mechanism responsible for its formation and maintenance remains controversial. In this work we investigate the dynamics maintaining the R-S in turbulent Poiseuille flow at $R=1650$. Spanwise collocation is used to remove spanwise displacement of the streaks and associated flow components, which isolates the streamwise-mean flow R-S component and the second-order statistics of the streamwise-varying fluctuations that are collocated with the R-S. This partition of the dynamics into streamwise-mean and fluctuation components facilitates exploiting insights gained from the analytic characterization of turbulence in the second-order statistical state dynamics (SSD), referred to as S3T, and its closely associated restricted nonlinear dynamics (RNL) approximation. Symmetry of the statistics about the streak centreline permits separation of the fluctuations into sinuous and varicose components. The Reynolds stress forcing induced by the sinuous and varicose fluctuations acting on the R-S is shown to reinforce low- and high-speed streaks, respectively. This targeted reinforcement of streaks by the Reynolds stresses occurs continuously as the fluctuation field is strained by the streamwise-mean streak and not intermittently as would be associated with streak-breakdown events. The Reynolds stresses maintaining the streamwise-mean roll arise primarily from the dominant proper orthogonal decomposition (POD) modes of the fluctuations, which can be identified with the time average structure of optimal perturbations growing on the streak. These results are consistent with a universal process of R-S growth and maintenance in turbulent shear flow arising from roll forcing generated by straining turbulent fluctuations, which was identified using the S3T SSD.

Pressure-driven flows of viscoelastic fluids in narrow non-uniform geometries are common in physiological flows and various industrial applications. For such flows, one of the main interests is understanding the relationship between the flow rate $q$ and the pressure drop $\Delta p$, which, to date, is studied primarily using numerical simulations. We analyse the flow of the Oldroyd-B fluid in slowly varying arbitrarily shaped, contracting channels and present a theoretical framework for calculating the $q-\Delta p$ relation. We apply lubrication theory and consider the ultra-dilute limit, in which the velocity profile remains parabolic and Newtonian, resulting in a one-way coupling between the velocity and polymer conformation tensor. This one-way coupling enables us to derive closed-form expressions for the conformation tensor and the flow rate–pressure drop relation for arbitrary values of the Deborah number ($De$). Furthermore, we provide analytical expressions for the conformation tensor and the $q-\Delta p$ relation in the high-Deborah-number limit, complementing our previous low-Deborah-number lubrication analysis. We reveal that the pressure drop in the contraction monotonically decreases with $De$, having linear scaling at high Deborah numbers, and identify the physical mechanisms governing the pressure drop reduction. We further elucidate the spatial relaxation of elastic stresses and pressure gradient in the exit channel following the contraction and show that the downstream distance required for such relaxation scales linearly with $De$.

We investigate by direct numerical simulations the fluid–solid interaction of non-dilute suspensions of spherical particles moving in triperiodic turbulence, at the relatively large Reynolds number of $Re_\lambda \approx 400$. The solid-to-fluid density ratio is varied between $1.3$ and $100$, the particle diameter $D$ is in the range $16 \le D/\eta \le 123$ ($\eta$ is the Kolmogorov scale) and the volume fraction of the suspension is $0.079$. Turbulence is sustained using the Arnold–Beltrami–Childress cellular-flow forcing. The influence of the solid phase on the largest and energetic scales of the flow changes with the size and density of the particles. Light and large particles modulate all scales in an isotropic way, while heavier and smaller particles modulate the largest scales of the flow towards an anisotropic state. Smaller scales are isotropic and homogeneous for all cases. The mechanism driving the energy transfer across scales changes with the size and the density of the particles. For large and light particles the energy transfer is only marginally influenced by the fluid–solid interaction. For small and heavy particles, instead, the classical energy cascade is subdominant at all scales, and the energy transfer is essentially driven by the fluid–solid coupling. The influence of the solid phase on the flow intermittency is also discussed. Besides, the collective motion of the particles and their preferential location in relation to properties of the carrier flow are analysed. The solid phase exhibits moderate clustering; for large particles the level of clustering decreases with their density, while for small particles it is maximum for intermediate values.

We consider two-dimensional flows above topography, revisiting the selective decay (or minimum enstrophy) hypothesis of Bretherton and Haidvogel. We derive a ‘condensed branch’ of solutions to the variational problem where a domain-scale condensate coexists with a flow at the (smaller) scale of the topography. The condensate arises through a supercritical bifurcation as the conserved energy of the initial condition exceeds a threshold value, a prediction that we quantitatively validate using direct numerical simulations. We then consider the forced–dissipative case, showing how weak forcing and dissipation select a single dissipative state out of the continuum of solutions to the energy-conserving system predicted by selective decay. As the forcing strength increases, the condensate arises through a supercritical bifurcation for topographic-scale forcing and through a subcritical bifurcation for domain-scale forcing, both predictions being quantitatively validated by direct numerical simulations. This method provides a way of determining the equilibrated state of forced–dissipative flows based on variational approaches to the associated energy-conserving system, such as the statistical mechanics of two-dimensional flows or selective decay.

Lubrication theory is adapted to incorporate the large normal stresses that occur for order-one Deborah numbers, $De$, the ratio of the relaxation time to the residence time. Comparing with the pressure drop for a Newtonian viscous fluid with a viscosity equal to that of an Oldroyd-B fluid in steady simple shear, we find numerically a reduced pressure drop through a contraction and an increased pressure drop through an expansion, both changing linearly with $De$ at high $De$. For a constriction, there is a smaller pressure drop that plateaus at high $De$. For a contraction, much of the change in pressure drop occurs in the stress relaxation in a long exit channel. An asymptotic analysis for high $De$, based on the idea that normal stresses are stretched by an accelerating flow in proportion to the square of the velocity, reveals that the large linear changes in pressure drop are due to higher normal stresses pulling the fluid through the narrowest gap. A secondary cause of the reduction is that the elastic shear stresses do not have time to build up to their steady-state equilibrium value while they accelerate through a contraction. We find for a contraction or expansion that the high $De$ analysis works well for $De>0.4$.

The Holton–Lindzen–Plumb (HLP) model of the quasi-biennial oscillation (QBO) is investigated in order to assess the impact of introducing intermittency in the wave forcing. Intermittency is introduced to HLP by allowing the amplitude of the waves which force the QBO to evolve according to a stationary random process, driven by a stochastic differential equation (SDE) with an associated time scale $\tau$. Provided that $\tau$ is much shorter than the QBO period, it is shown that the impact on the QBO of the intermittent forcing is captured by a single intermittency parameter $\lambda$, and the value of $\lambda$ is proportional to $\tau$ and otherwise depends upon the details of the SDE. Numerical simulations, using a family of mean-reverting Ornstein–Uhlenbeck processes as the choice of SDE, show that the effect of increasing the intermittency parameter is invariably to decrease the QBO amplitude and increase its period. Changes to the QBO amplitude and period are indeed found to collapse onto a single curve controlled by $\lambda$, as predicted by the theory, provided that $\tau$ is small enough for the approximations used to be valid. The extension to broadband forcing is discussed in the context of stochastic gravity wave parameterisation, with the eventual goal of developing a representation of source intermittency in the most general situation with close fidelity to the physics.