We report on the phenomenon of detonation synchronization and demonstrate the existence of the Arnold tongues and devil's staircases in the problem of gaseous detonation in a periodically inhomogeneous reactive medium. Universal properties of these dynamical structures – the Farey tree, fractal dimension and period-doubling bifurcations – are revealed.

Turbulent convection at low Prandtl numbers is in many aspects still terra incognita on the parameter map. One reason for this fact is that laboratory experiments on turbulent convection in this parameter regime are notoriously challenging as they require the use of opaque liquid metals. These working fluids prevent the application of typical optical imaging techniques such as particle image velocimetry. Recent experiments by Grannan et al. (J. Fluid Mech., vol. 939, 2022, R1) shed new light on the variety of regimes in liquid metal flows which include rotating convection, magnetoconvection and rotating magnetoconvection next to the classical Rayleigh–Bénard case. More importantly, the authors manage the seamless crossover from one regime into another. They were thus able to study low-Prandtl-number convection at different levels of complexity in a single experimental set-up. Their work provides new insights into the tight connections between characteristic large-scale flow behaviours and the resulting global heat transfer magnitudes. This has implications for convection in planetary cores and stellar convection zones and connected dynamo action.

The influence of nucleation on cavitation about a sphere from inception through to supercavitation at a transcritical Reynolds number of $1.5\times 10^{6}$ is investigated experimentally. Two extreme free-stream nuclei populations, deplete and abundant, were investigated. Unsteady surface pressures from two sensors on opposing sides of the sphere were acquired simultaneous with high-resolution high-speed photography at cavitation numbers between 1.0 and 0.3. High-resolution spectrograms derived from these measurements reveal principally bi-modal shedding in attached and detached regimes. Correlations between unsteady pressure measurements show the high modes to be axisymmetric and low modes asymmetric. Modal topology is also discerned from the high-speed imaging. The bi-modal shedding for lower cavitation numbers is driven by coupled re-entrant jet formation and upstream shockwave propagation. The attached regime is shown to have two sub-regimes. For the abundant case, the continuous supply of activated nuclei around the sphere periphery in the first bi-modal regime has the effect of driving the high symmetric mode preferentially over the asymmetric low mode compared with the deplete case. For the first bi-modal regime, frequencies were unaffected by nucleation changes although peak responses were centred at a cavitation number of about 0.8 for the deplete and 0.825 for the abundant. For the second attached regime, where cavity lengths are of the order of the sphere size, changes in nucleation altered frequencies and amplitudes of peak unsteady pressures. For the abundant case, the continuous nuclei supply significantly reduced coherence with modal peak amplitudes reduced by an order compared with the deplete case. Continuous nuclei activation increased the probability of the high mode over the low compared with the deplete case but to a lesser extent than the first regime. Nuclei activations also significantly reduced inter-cavity and cavity durations, but not growth and collapse phases, which increased modal frequencies compared with the deplete case. The second regime asymmetric low mode topology, for both nucleation cases, is shown to be alternate shedding of oblique vortices from diametrically opposing sides of the sphere similar to low Reynolds number shedding about spheres and other axisymmetric bluff bodies in single-phase flows.

We demonstrate dynamic stabilisation of axisymmetric Fourier modes susceptible to the classical Rayleigh–Plateau (RP) instability on a liquid cylinder by subjecting it to a radial oscillatory body force. Viscosity is found to play a crucial role in this stabilisation. Linear stability predictions are obtained via Floquet analysis demonstrating that RP unstable modes can be stabilised using radial forcing. We also solve the linearised, viscous initial-value problem for free-surface deformation obtaining an equation governing the amplitude of a three-dimensional Fourier mode. This equation generalizes the Mathieu equation governing Faraday waves on a cylinder derived earlier in Patankar et al. (J. Fluid Mech., vol. 857, 2018, pp. 80–110), is non-local in time and represents the cylindrical analogue of its Cartesian counterpart (Beyer & Friedrich, Phys. Rev. E, vol. 51, issue 2, 1995, p. 1162). The memory term in this equation is physically interpreted and it is shown that, for highly viscous fluids, its contribution can be sizeable. Predictions from the numerical solution to this equation demonstrate the predicted RP mode stabilisation and are in excellent agreement with simulations of the incompressible Navier–Stokes equations (up to the simulation time of several hundred forcing cycles).

Ion and water are transported by electroconvection near permselective membranes, resulting in complex phenomena associated with the flow–fines interaction. Sheltering the flow chaos by the shear flow is a common strategy in plasma fluids and has recently been successfully applied to control ionic fluids. The paper herein reveals the critical selection of shear velocity regarding the fluid from a chaotic to a steady state through numerical and theoretical analyses. For the shear sheltering, the dimensionless Debye length ${\lambda _D}$ with varying channel height is introduced to achieve a comprehensive discussion of the factors and laws affecting the shear vortex state. Based on an analysis of the vortex driving mechanism, the scaling of the slip velocity ${u_s}\sim {(\lambda _D^{ - 1}\Delta {\phi ^4})^{1/3}}$ is recommended as the critical selection factor for the steady and chaotic state under a fixed shear flow velocity, which involves the dimensionless Debye length ${\lambda _D}$ and voltage difference $\Delta \phi $. Furthermore, for ionic fluid control by shear flow, a critical shear velocity ${U_{HPC}}$ is proposed to distinguish the electroconvective flow from a chaotic state to a steady state. When the shear flow velocity ${U_{HP}} > {U_{HPC}}$, the shear flow shelters chaos, and the scaling law is also recommended for the regulation of the critical shear flow velocity ${U_{HPC}}$ jointly by ${\lambda _D}$ and $\Delta \phi $. The analysis is confirmed by direct numerical simulation and existing experimental data (J. Fluid Mech, vol. 813, 2017, pp. 799–823). This work provides a more comprehensive physical insight for shear sheltering and affects the design of electromembrane microfluidics.

We study the dynamics of atomization of a liquid column by a coaxial gas flow with varying gas pressures. Specifically, we analyse how the gas density increase associated with elevated gas pressures in the ambient and co-flowing gas jet influences the liquid destabilization and breakup process, as well as the resulting droplet formation and dispersion. We present new experimental results for a coaxial liquid–gas atomizer operating in a high-pressure environment, with gas–liquid momentum ratio in the range $M = 5\unicode{x2013}56$ and pressurized gas densities $\rho _g/\rho _0 = 1\unicode{x2013}5$, where $\rho _0$ is the ambient gas density at standard conditions. High-speed shadowgraphy images are used to quantify the spatially and temporally varying liquid–gas interface in the spray near-field. Liquid core lengths, spreading angles and other spray metrics are presented, and the influence of gas density is identified from the comparison with atomization at atmospheric conditions. In the spray mid-field, phase Doppler interferometry is used in conjunction with laser Doppler velocimetry to quantify the droplet size and velocities, as well as their radial variations across the spray. Results show an increase in droplet size at elevated ambient pressures, when keeping the gas–liquid momentum ratio constant. Finally, we show that these observations are in line with predictions from the Kelvin–Helmholtz and Rayleigh–Taylor instabilities, both of which are relevant to the gas–liquid atomization process.

The hydrodynamic mechanism of drag reduction by a flexible hairy coating was explored using the penalty immersed boundary method. A two-dimensional flexible hairy coating is constituted by multiple flexible filaments. A simulation of a cylinder without a hairy coating at a Reynolds number of 100 was also performed for comparison. The results of the simulations show good agreement with the experimental data by Niu & Hu (Phys. Fluids, vol. 23, 2011, 101701), where maximum drag reduction of 22% was attained at a particular length, bending rigidity, coating density and coating angle of the hairy coating. The hydrodynamic mechanism of drag reduction was characterized in terms of the wake pattern, shape deformation and kinetic energy of the hairy coating. The effect of a non-uniform bending rigidity of the hairy coating on drag reduction was explored. A stable streamline shape of the hairy coating was found to delay the vortex formation and stabilize the recirculation zone, resulting in decreased form drag. Active flapping of the hairy coating with enhanced vortex shedding is adverse to drag reduction. A hairy coating with a stiff base and flexible trailing edge is beneficial to maintaining a stable shape.

Budgets of turbulent kinetic energy (TKE) and turbulent potential energy (TPE) at different scales $\ell$ in sheared, stably stratified turbulence are analysed using a filtering approach. Competing effects in the flow are considered, along with the physical mechanisms governing the energy fluxes between scales, and the budgets are used to analyse data from direct numerical simulation at buoyancy Reynolds number $Re_b=O(100)$. The mean TKE exceeds the TPE by an order of magnitude at the large scales, with the difference reducing as $\ell$ is decreased. At larger scales, buoyancy is never observed to be positive, with buoyancy always converting TKE to TPE. As $\ell$ is decreased, the probability of locally convecting regions increases, though it remains small at scales down to the Ozmidov scale. The TKE and TPE fluxes between scales are both downscale on average, and their instantaneous values are correlated positively, but not strongly so, and this occurs due to the different physical mechanisms that govern these fluxes. Moreover, the contributions to these fluxes arising from the sub-grid fields are shown to be significant, in addition to the filtered scale contributions associated with the processes of strain self-amplification, vortex stretching and density gradient amplification. Probability density functions (PDFs) of the $Q,R$ invariants of the filtered velocity gradient are considered and show that as $\ell$ increases, the sheared-drop shape of the PDF becomes less pronounced and the PDF becomes more symmetric about $R=0$.

This paper theoretically studies the effect of eccentricity on the conditions of capillary emptying (determined by critical Bond number) in a horizontal annular tube in a downward gravity field. Experiments are conducted to compare with theoretical results. We find that non-horizontal eccentricity can lead to the occurrence of a re-entrant liquid-state transition (from liquid non-occlusion to liquid plug to liquid non-occlusion) with increasing Bond number, when the eccentricity (e) or inner-to-outer radius ratio (χ) is large enough, and the two liquid non-occlusion states correspond to different emptying mechanisms dominated by the gravity effect and the ‘wedge’ effect, respectively. Existence of the re-entrant transition is accompanied by occurrence of unconditional liquid non-occlusion at large enough or small enough contact angles regardless of Bond numbers. The critical Bond numbers at a contact angle γ for vertical upward eccentricity are equal to those at a contact angle 180° − γ for vertical downward eccentricity. In a parameter space (γ, e/(1 − χ)), the region with the re-entrant transition becomes larger with the eccentric angle varying from 0° (horizontal) to 90° (vertical). Optimization of geometrical parameters and inner and outer contact angles can lead to better effect of capillary emptying. This paper provides a very effective scheme for removing a liquid blockage from a capillary in optofluidics/microfluidics.

Wind gustiness in the marine atmospheric boundary layer affects significantly the dynamics of air–sea interaction. To understand the impacts of wind gust events, we perform large-eddy simulation of wind turbulence over a travelling wave to investigate the response of the wind field to an impulsive wind speed increase or decrease. It is found that the turbulence fluctuations and the terms in the turbulent kinetic energy budget equation have a delayed response to the change in the mean flow, while the response of the wave-coherent motions is quasi-stationary. The wave-coherent motions are investigated quantitatively through comparison with a viscous curvilinear model developed by Cao et al. (J. Fluid Mech., vol. 901, 2020, A27) and Cao & Shen (J. Fluid Mech., vol. 919, 2021, A38). We observe an asymmetric hysteresis between the growing wind and the decaying wind in the evolution of the form drag and the viscous drag. We find further that the variation of the wave growth rate during the wind gust is related closely to the contribution from the out-of-phase component of the vertical velocity. Our discoveries provide evidence for the necessity of improving non-equilibrium turbulence and wind input modelling to account for the wind gustiness effect in future studies.

This paper demonstrates experimentally that imposed periodic forcing can significantly alter the global flow characteristics of the flow over a double backward-facing step. The geometry consists of two equal height steps spaced up to eight step heights apart. A periodic zero-mass flux jet located at the first step's top corner was issued at frequencies ranging from below the step-mode instability frequency up to approximately five times the shear-layer instability frequency. Reattachment of the flow onto the first step was achieved for step separations as low as three single-step heights with imposed forcing; significantly shorter than the five single-step heights that occurred without forcing. A significant reduction in mean base pressure on the first step, and increase on the second step, occurred for low forcing frequencies. Even for large step separations, the effect of forcing on the flow persisted sufficiently far downstream to appreciably influence the development of the second recirculation zone. Importantly, previous forced single and unforced double backward-facing step flows provide reference cases to examine and discuss similarities and differences. This study offers insight into possibilities and potential outcomes of flow control for applications ranging from the drag reduction of ground vehicles such as pickup trucks, to enhanced mixing in industrial processes.

We consider as base flow the stationary vortex filament solution obtained by Castillo-Castellanos et al. (Phys. Rev. Fluids, vol. 6, 2021, 114701) in the far wake of a rotor with tip-splitting blades. The cases of a single blade and of two blades with a hub vortex are studied. In these solutions, each blade generates two closely spaced co-rotating tip vortices that form a braided helical pattern in the far wake. The long-wave stability of these solutions is analysed using the same vortex filament framework. Both the linear spectrum and the linear impulse response are considered. We demonstrate the existence of different types of instability modes. A first type corresponds to the local pairing of consecutive turns of the helical pattern, which is well described by the instability of a uniform helical vortex with a core size given by the mean separation distance of the vortices in the pair. A second type corresponds to the pairing of consecutive turns of the vortex pair and is observed only for densely braided patterns, which is well described by the instability of two interlaced helical vortices by straightening out the baseline helix. A third type of unstable modes modifies the separation distance between the vortices in each pair and amplifies specific (linear) wavelengths. These unstable modes also spread spatially with a weaker rate.

The turbulent kinetic energy and energy dissipation rate in the wake of a circular cylinder are examined at a Reynolds number of 1000. The turbulence characteristics are quantified using direct numerical simulation, which provides a comprehensive dataset that is almost impossible to acquire from physical experiments. The energy dissipation rate is decomposed into the components due to the mean flow, the coherent primary vortices and the remainder. It is found that the remainder component, which develops only in a three-dimensional turbulent wake and resides mainly in the regions of vortices, accounts for 95 % and 97 % of the total dissipation rate for 10 and 20 cylinder diameters downstream of the cylinder, respectively (while the remainder accounts for 62 % and 83 % of the total turbulent kinetic energy). Based on the remainder component, the validity of local isotropy, local axisymmetry, local homogeneity and homogeneity in the y–z plane for the turbulent dissipation in the wake is examined. The analysis reveals that the turbulent dissipation is largely locally homogeneous, but not locally isotropic or axisymmetric, even after the annihilation of the primary vortex street. In addition, the performances of the four corresponding surrogates to the true dissipation rate are evaluated. Owing to the general validity of local homogeneity, the surrogates of local homogeneity and homogeneity in the y–z plane perform well. Although local axisymmetry does not hold, the corresponding surrogate performs well, because errors from different terms largely cancel out. However, the surrogate of local isotropy generally underestimates the true dissipation rate.

The effect of a chordwise flexible aft-tail of a rigid heaving aerofoil on the dynamical transitions of the trailing-wake is studied here. The two-way coupled fluid–solid dynamics is simulated using an in-house fluid–structure interaction (FSI) platform, comprising a discrete forcing immersed boundary method based incompressible Navier–Stokes solver, weakly coupled with a finite difference method based structural solver. The FSI dynamics is studied in comparison to the corresponding rigid tail configuration. For the latter, mild jet-switching due to quasi-periodic movement of the wake vortices gives way to vigorous jet-switching as the dynamics transitions to a state of intermittency, where the quasi-periodic behaviour gets interspersed with chaotic windows. Introduction of a moderately flexible tail regularises this intermittent dynamics, eliminating jet-switching. The wake exhibits a deflected reverse Kármán pattern with fluctuating angles, governed by quasi-periodicity. With a highly flexible tail (very low rigidity), the wake shows almost a symmetric reverse Kármán street as periodicity is restored. Flexibility of the aft-tail is next controlled by changing its length, and flow is regularised and periodicity retained for moderate rigidity for increased length. Different dynamical states are established through robust nonlinear dynamical tools. The underlying flow-field behaviour, instrumental in suppressing the jet-switching phenomenon, is identified through a detailed investigation of the near-field vortex interactions dictated by the dynamics. A suite of measures has also been derived from the unsteady flow field to quantify the interactions of the key near-field vortices with a view to understanding the mechanism of switching and its subsequent suppression through flexibility.

For the Stokes number based on the Kolmogorov time scale ${St}_K$ up to $O(10^{2})$, the particle subgrid stress (particle stress) in the volume-average framework is studied by comparing the fluid residual stress, the particle Smagorinsky model and the particle scale-similarity model. To obtain the numerical database of the particle-laden turbulence, two-way coupling direct numerical simulation is carried out with isotropic and anisotropic forcing conditions. As the particle stress is related to the local flow structure, which is not reflected by ${St}_K$, a new Stokes number ${St}_R$ is introduced to extract the effect of the intensity of the fluid velocity fluctuation in the averaging volume. The degrees of agreement of the principal axes (eigenvectors) of the particle stress models to those of the fully resolved particle stress are regarded as functions of ${St}_R$ regardless of the averaging volume size. The fluid residual stress model shows the highest degree of agreement over a small ${St}_R$ range for both of the forcing cases, and similar predominance is also observed for the correlation coefficient reflecting the magnitude of the particle stress. The effects of ${St}_R$, ${St}_K$, the averaging volume size and the Reynolds number on the model coefficients are investigated based on the intensities of the deviatoric and isotropic parts of the fully resolved particle stress and its models. The Stokes number ${St}_R$ is found to be an essential factor to determine the model coefficients, as each effect is extracted reasonably by fixing ${St}_R$.

Pressure fluctuations play an essential role in the transport of turbulent kinetic energy and vibrational loading. This study focuses on examining the effect of wall cooling on pressure fluctuations in compressible turbulent boundary layers by high-fidelity direct numerical simulations. Pressure fluctuations result from the vorticity mode and the acoustic mode that are both closely dependent on compressibility. To demonstrate the effects of wall cooling at various compressibility intensities, three free-stream Mach numbers are investigated, i.e. $M_\infty =0.5$, 2.0 and 8.0, with real gas effects being absent for $M_\infty =8.0$ due to a low enthalpy inflow. Overall, opposite effects of wall cooling on pressure fluctuations are found between the subsonic/supersonic cases and the hypersonic case. Specifically, the pressure fluctuations normalized by wall shear stress $p^\prime _{rms}/\tau _w$ are suppressed in the subsonic and supersonic cases, while enhanced in the hypersonic case near the wall. Importantly, travelling-wave-like alternating positive and negative structures (APNS), which greatly contribute to pressure fluctuations, are identified within the viscous sublayer and buffer layer in the hypersonic cases. Furthermore, generating mechanisms of pressure fluctuations are explored by extending the decomposition based on the fluctuating pressure equation to compressible turbulent boundary layers. Pressure fluctuations are decomposed into five components, in which rapid pressure, slow pressure and compressible pressure are dominant. The suppression of pressure fluctuations in the subsonic and supersonic cases is due to both rapid pressure and slow pressure being suppressed by wall cooling. In contrast, wall cooling strengthens compressible pressure for all Mach numbers, especially in the hypersonic case, resulting in increased wall pressure fluctuations. Compressible pressure plays a leading role in the hypersonic case, mainly due to the APNS. Essentially, the main effects of wall cooling can be interpreted by the suppression of the vorticity mode and the enhancement of the acoustic mode.

We report analytical solutions of a problem involving a visco-elastic solid material layer sandwiched between two fluid layers, in turn confined by two long planar walls that undergo oscillatory motion. The resulting system dynamics is rationalized, based on fluid viscosity and solid elasticity, via wave and boundary layer theory. This allows for physical interpretation of elasto-hydrodynamic coupling, potentially connecting to a broad set of biophysical phenomena and applications, from synovial joint mechanics to elastometry. Further, obtained solutions are demonstrated to be rigorous benchmarks for testing coupled incompressible fluid–hyperelastic solid and multi-phase numerical solvers, towards which we highlight challenging parameter sets. Finally, we provide an interactive online sandbox to build physical intuition, and open-source our code-base.

The fluid–structure interaction problem composed of an elongated, finite-length, flexible cylinder falling in a fluid at rest is investigated experimentally. Tomographic reconstruction of the cylinder and three-dimensional particle tracking velocimetry of the surrounding fluid, based on the Shake-The-Box algorithm, are used jointly to capture both solid and fluid motions. Starting from the rectilinear vertical fall characterized by a steady wake, focus is placed on subsequent regimes involving, mainly in the horizontal direction, periodic rigid-body motions (RBM) of weak amplitude or periodic large-amplitude bending oscillations (BO). Two RBM regimes are explored: the TRA regime where the cylinder exhibits translational oscillations in a plane perpendicular to its axis, and the AZI regime in which the body displays an azimuthal oscillation around its centre. The associated unsteady wakes are composed of counter-rotating vortices bending near the body ends to connect with the adjacent vortex rows. Specific organizations of the vortical structures are uncovered, depending on the regime. In particular, in the AZI regime, they present an antisymmetrical distribution relative to the midspan point. For a sufficiently long cylinder, BO regimes emerge, resembling the structural modes of an unsupported beam. The associated wakes exhibit a cellular organization. Within each cell delimited by two deformation nodes, two counter-rotating vortex rows are shed per oscillation cycle. Flow velocity fluctuations are in phase opposition on each side of a deformation node. For both RBM and BO regimes, frequency and phase analyses of cylinder and wake behaviours, along the span, highlight the spatio-temporal synchronization of the unsteady flow and moving body.

A numerical investigation and reduced-order modelling of natural convection in a cavity with differentially heated sidewalls is discussed. The effect of conjugate boundary conditions on the cavity's heat transfer and natural flow circulation with varying aspect ratios and Rayleigh numbers is examined. Validation of the canonical differentially heated cavity reveals that a modification to the definition of the Rayleigh number (Ra) is required to reconcile the $Nu\sim Ra^{1/4}$ scaling (where Nu is the Nusselt number) and validity of previously proposed correlations for the heat transfer in vertical enclosures. Dynamic mode decomposition is used to uncover the underlying time-dependent flow structures and the results are compared with stability bifurcation studies in the literature. A flow mode previously associated with smaller aspect ratio cavities is identified as the unstable mode for a larger aspect ratio of 4. The effect of conjugate boundary conditions is scaled based on the ratio of the internal and external boundary layers, wherein higher external Reynolds numbers aid in the heat transfer as the Nusselt number approaches the isothermal limit. A reduced-order theoretical model is proposed to predict the Nusselt number for the conjugate boundary conditions. The performance of conjugate boundary conditions is connected to the flow stability, with the parallel-flow configuration acting in a destabilizing manner while the counterflow configuration has a stabilizing effect and results in the highest heat transfer. We also found that relatively large heat transfer can be achieved with substantially lower external actuation with intermediate aspect ratios, translating to less power in practice.

Richtmyer–Meshkov (RM) instability in convergent geometries (such as cylinders and spheres) plays a fundamental role in natural phenomena and engineering applications, e.g. supernova explosion and inertial confinement fusion. Convergent geometry refers to a system in which the interface converges and the fluids are compressed correspondingly. By applying a decomposition formula, the stretching or compression (S(C)) effect is separated from the perturbation growth as one of the main contributions, which is defined as the averaged velocity difference between two ends of the mixing zone. Starting from linear theories, the S(C) effect in planar, cylindrical and spherical geometries is derived as a function of geometrical convergence ratio, compression ratio and mixing width. Specifically, geometrical convergence stretches the mixing zone, while fluid compression compresses the mixing zone. Moreover, the contribution of geometrical convergence in the spherical geometry is more important than that in the cylindrical geometry. A series of cylindrical cases with high convergence ratio is simulated, and the growth of perturbations is compared with that of the corresponding planar cases. As a result, the theoretical results of the S(C) effect agree well with the numerical results. Furthermore, results show that the S(C) effect is a significant feature in convergent geometries. Therefore, the S(C) effect is an important part of the Bell–Plesset effect. The present work on the S(C) effect is important for further modelling of the mixing width of convergent RM instabilities.