We analyse the low-frequency dynamics of a high Reynolds number impinging shock-wave/turbulent boundary-layer interaction (SWBLI) with strong mean-flow separation. The flow configuration for our grid-converged large-eddy simulations (LES) reproduces recent experiments for the interaction of a Mach 3 turbulent boundary layer with an impinging shock that nominally deflects the incoming flow by $19.6^{\circ }$. The Reynolds number based on the incoming boundary-layer thickness of $Re_{\unicode[STIX]{x1D6FF}_{0}}\approx 203\times 10^{3}$ is considerably higher than in previous LES studies. The very long integration time of $3805\unicode[STIX]{x1D6FF}_{0}/U_{0}$ allows for an accurate analysis of low-frequency unsteady effects. Experimental wall-pressure measurements are in good agreement with the LES data. Both datasets exhibit the distinct plateau within the separated-flow region of a strong SWBLI. The filtered three-dimensional flow field shows clear evidence of counter-rotating streamwise vortices originating in the proximity of the bubble apex. Contrary to previous numerical results on compression ramp configurations, these Görtler-like vortices are not fixed at a specific spanwise position, but rather undergo a slow motion coupled to the separation-bubble dynamics. Consistent with experimental data, power spectral densities (PSD) of wall-pressure probes exhibit a broadband and very energetic low-frequency component associated with the separation-shock unsteadiness. Sparsity-promoting dynamic mode decompositions (SPDMD) for both spanwise-averaged data and wall-plane snapshots yield a classical and well-known low-frequency breathing mode of the separation bubble, as well as a medium-frequency shedding mode responsible for reflected and reattachment shock corrugation. SPDMD of the two-dimensional skin-friction coefficient further identifies streamwise streaks at low frequencies that cause large-scale flapping of the reattachment line. The PSD and SPDMD results of our impinging SWBLI support the theory that an intrinsic mechanism of the interaction zone is responsible for the low-frequency unsteadiness, in which Görtler-like vortices might be seen as a continuous (coherent) forcing for strong SWBLI.

A new amplitude equation is derived for high-frequency acoustic waves propagating through an incompressible vortical flow using multi-time-scale asymptotic analysis. The reduced model is derived without an explicit spatial-scale separation ansatz between the wave and vortical fields. As a consequence, the model is seen to capture very well the features of the wave field in the regime where the spatial scales of the wave and vortical fields are comparable, a regime for which an optimal reduced model does not seem to be available.

In this work, we study active particles with prescribed surface velocities in non-Newtonian fluids. We employ the reciprocal theorem to obtain the velocity of an active spherical particle with an arbitrary axisymmetric slip velocity in an otherwise quiescent second-order fluid. We then determine how the motion of a diffusiophoretic Janus particle is affected by complex fluid rheology, namely viscoelasticity and shear-thinning viscosity, compared to a Newtonian fluid, assuming a fixed slip velocity. We find that a Janus particle may go faster or slower in a viscoelastic fluid, but is always slower in a shear-thinning fluid as compared to a Newtonian fluid.

We investigate the spatial stability of compressible, viscous pipe flows with radius-dependent mean density profiles, subjected to solid body rotations. For a fixed Rossby number $\unicode[STIX]{x1D716}$ (inverse of the rotational speed), as the Reynolds number $Re$ is increased, the flow transitions from being stable to convectively unstable, usually leading to absolute instability. If flow compressibility is unimportant and $Re$ is held constant, there appears to be a maximum $Re$ below which the flow remains stable irrespective of any rotational speed, or a minimum azimuthal Reynolds number $Re_{\unicode[STIX]{x1D703}}$ $(=Re/\unicode[STIX]{x1D716})$ is required for any occurrence of absolute instabilities. Once compressible forces are significant, the effect of pressure–density coupling is found to be more severe below a critical $Re$, where as rotational speeds are raised, a stable flow almost directly transitions to an absolutely unstable state. This happens at a critical $Re_{\unicode[STIX]{x1D703}}$ which reduces with increased flow Mach number, pointing to compressibility aiding in the instability at these lower Reynolds numbers. However, at higher $Re$, above the critical value, the traditional stabilizing role of compressibility is recovered if mean density stratification exists, where the gradients of density play an equally important role, more so at the higher azimuthal modes. A total disturbance energy-based formulation is used to obtain mechanistic understanding at these stability states, where we find the entropic energy perturbations to dominate as the primary instability mechanism, in sharp contrast to the energy due to axial shear, known to play a leading role in incompressible swirling flows.

We report results from an experimental study on the formation of stable–streamlined and helical cavity wakes following the free-surface impact of Leidenfrost spheres. Similar to the observations of Mansoor et al. (J. Fluid Mech., vol. 743, 2014, pp. 295–326), we show that acoustic ripples form along the interface of elongated cavities entrained in the presence of wall effects as soon as the primary cavity pinch-off takes place. The crests of these ripples can act as favourable points for closure, producing multiple acoustic pinch-offs, which are found to occur in an acoustic pinch-off cascade. We show that these ripples pacify with time in the absence of physical contact between the sphere and the liquid, leading to extremely smooth cavity wake profiles. More importantly, the downward-facing jet at the apex of the cavity is continually suppressed due to a skin-friction drag effect at the colliding cavity-wall junction, which ultimately produces a stable–streamlined cavity wake. This streamlined configuration is found to experience drag coefficients an order of a magnitude lower than those acting on room-temperature spheres. A striking observation is the formation of helical cavities which occur for impact Reynolds numbers $Re_{0}\gtrsim 1.4\times 10^{5}$ and are characterized by multiple interfacial ridges, stemming from and rotating synchronously about an evident contact line around the sphere equator. The contact line is shown to result from the degeneration of Kelvin–Helmholtz billows into turbulence which are observed forming along the liquid–vapour interface around the bottom hemisphere of the sphere. Using sphere trajectory measurements, we show that this helical cavity wake configuration has 40 %–55 % smaller force coefficients than those obtained in the formation of stable cavity wakes.

In this study, the Bragg–Hawthorne equation (BHE) is extended in the context of a steady, inviscid and compressible fluid, thus leading to an assortment of partial differential equations that must be solved simultaneously. A solution is pursued by implementing a Rayleigh–Janzen expansion in the square of the reference Mach number. The corresponding formulation is subsequently used to derive a compressible approximation for the Trkalian model of the bidirectional vortex. The approximate solution is compared to a representative computational fluid dynamics simulation in order to validate the modelling assumptions under realistic conditions. The latter is found to exhibit an appreciable steepening of the axial velocity profile, which is accompanied by an axial dependence in the mantle location that is somewhat reminiscent of the radial shifting of mantles reported in some experimental trials and simulations. In this context, flows with a strong swirl intensity do not seem to be significantly affected by the introduction of compressibility. Rather, as the swirl intensity is reduced the effects of compressibility become more noticeable, especially in the axial and radial velocity components. It may also be realized that imparting a progressively larger swirl component stands to promote the axisymmetric distribution of flow field properties, and these include an implicit resistance to dilatational effects in the tangential direction. From a broader perspective, this study provides a viable approximation to the Trkalian motion associated with cyclonic flows, while serving as a limited proof of concept for the compressible Bragg–Hawthorne procedure applied to a steady, axisymmetric and inviscid fluid.

This paper describes a study of the generation of a plughole vortex and its consequences in a drainpipe during drainage of water from a stationary rectangular tank. The critical and minimum depths of water above the inlet of the drainpipe, where a surface dip starts to develop for drainpipes of various diameters, were examined parametrically. This study explored the following naturally occurring phenomena arising from a plughole vortex. (i) A plughole vortex initially causes a surface dip to develop towards the inlet of the drainpipe and as the surface dip approaches the inlet of the drainpipe it creates a droplet-shaped air bubble. (ii) A unique bubble transformation, i.e. from a droplet-shaped to a donut-shaped bubble ring, occurs just after the separation of the droplet-shaped air bubble from the surface dip. (iii) The donut-shaped bubble ring flows with the drain water and initially causes bubbly flow in the drainpipe. (iv) As the water head above the inlet of the drainpipe decreases, the droplet-shaped bubble size increases, and consequently, the bubble ring size increases and causes slug flow in the drainpipe. (v) As the slugs combine, the flow of the draining water eventually becomes annular flow in the drainpipe. Sounds, such as that of instantaneous fizz and bubble sink draining, were observed to be produced as a result of the bubble formation process. Temporal changes in the shape and size of the air bubbles were studied. Within the range of 0.45–0.6, the ratio of the bubble diameter to the bubble length was found to be linearly proportional to the ratio of the water depth to the diameter of the drainpipe. Several drainage cases were simulated numerically to observe the physics of these naturally occurring phenomena. The shapes and sizes of the vortices induced by plugholes have been visualised and analysed using the vortex core method.

We derive an alternative formulation of the turbulent kinetic energy equation for flows with strong near-wall density and viscosity gradients. The derivation is based on a scaling transformation of the Navier–Stokes equations using semi-local quantities. A budget analysis of the semi-locally scaled turbulent kinetic energy equation shows that, for several variable property low-Mach-number channel flows, the ‘leading-order effect’ of variable density and viscosity on turbulence in wall bounded flows can effectively be characterized by the semi-local Reynolds number. Moreover, if a turbulence model is solved in its semi-locally scaled form, we show that an excellent agreement with direct numerical simulations is obtained for both low- and high-Mach-number flows, where conventional modelling approaches fail.

An assessment of the turbulent boundary layer flow structure, which is coherent with the near-wall region, is carried out through a spectral coherence analysis. This spectral method is applied to datasets comprising synchronized two-point streamwise velocity signals at a near-wall reference position and a range of wall-normal positions spanning a Reynolds-number range $Re_{\unicode[STIX]{x1D70F}}\sim O(10^{3}){-}O(10^{6})$. Within each dataset, a self-similar structure is identified from the coherence between the turbulence in the logarithmic region and at the near-wall reference position. This self-similarity is described by a streamwise/wall-normal aspect ratio of $\unicode[STIX]{x1D706}_{x}/z\approx 14$, where $\unicode[STIX]{x1D706}_{x}$ and $z$ are the streamwise wavelength and wall-normal distance respectively.

Cavitation bubbles are nucleated at a high repetition rate in water by delivering a pulsed laser through a fibre optic. Continuous high-frequency cavitation drives a stream away from the fibre tip. Using high-speed photography and particle image velocimetry, the stream is characterised as a synthetic jet, generated by trains of vortices induced by non-spherical bubble collapse. At low laser power, the bubbles collapse before the arrival of a subsequent laser pulse. Yet, by increasing the laser power, a system of bubbles is formed which leads to complex bubble–bubble interactions. The synthetic jet is observed regardless of the bubble formation regime, demonstrating the stability of the phenomenon. Synthetic jet generation by repetitive bubble collapse extends the well-studied acoustic streaming from small-amplitude bubble oscillations.

Experimental measurements of the displacements of shock accelerated microparticles from shortly after shock interaction to the particle relaxation time show time-dependent drag coefficients ($C_{D}$) that are much higher than those predicted by quasi-steady and unsteady drag models. Nylon particles with mean diameter of $4~\unicode[STIX]{x03BC}\text{m}$, accelerated by one-dimensional normal shocks (Mach number $M_{s}=1.2$, 1.3 and 1.4), have measured $C_{D}$ values that follow a power-law behaviour. The drag is a function of the time-dependent Knudsen number, $Kn^{\ast }=M_{s}/Re_{p}$, where the particle Reynolds number ($Re_{p}$) is calculated using the time-dependent slip velocity. Some portion of the drag can be attributed to quasi-steady forces, but the total drag cannot be predicted by current unsteady force models that are based on the Basset–Boussinesq–Oseen equation and pressure drag. The largest contribution to the total drag is the unsteady component ($C_{D,us}$) until the particle attains $Kn^{\ast }\approx 0.5{-}1.0$, then the unsteady contribution decays. The quasi-steady component ($C_{D,qs}$) increases almost linearly with $Kn^{\ast }$, intersects the $C_{D,us}$ at $Kn^{\ast }\approx 2$ and becomes the primary contributor to the drag towards the end of the relaxation zone as $Re_{p}\rightarrow 0$. There are currently no analytical models that are able to predict the nonlinear behaviour of the shock accelerated particles during the relaxation phase of the flow.

Here, we examine the conditions for layer formation in a stratified shear flow. We reformulate the conditions for amplification of small perturbations to a uniform stratification first proposed by Phillips (Deep Sea Research and Oceanographic Abstracts, vol. 19, 1972, pp. 79–81, Elsevier) and Posmentier (J. Phys. Oceanogr., vol. 7 (2), 1977, pp. 298–300) using the sorted buoyancy coordinates introduced by Nakamura (J. Atmos. Sci., vol. 53 (11), 1996, pp. 1524–1537) and Winters & D’Asaro (J. Fluid Mech., vol. 317, 1996, pp. 179–193). We consider the possible dependence of the effective diffusivity on three non-dimensional parameters, the gradient Richardson number, the buoyancy Reynolds number and the Prandtl number, and obtain a simple expression for conditions favourable for layer formation. The new framework is applied to direct numerical simulations of stratified shear flow. We then apply a recent multi-parameter parameterization developed by Salehipour et al. (Geophys. Res. Lett., vol. 43 (7), 2016, pp. 3370–3379), which suggests that layer formation is favoured for large Prandtl numbers and moderate to large values of the gradient Richardson and buoyancy Reynolds numbers.