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

In this work we present experiments and simulations on the nucleation and successive dynamics of laser-induced bubbles inside liquid droplets in free-fall motion, i.e. a case where the bubbles are subjected to the influence of a free boundary in all directions. Within this spherical millimetric droplet, we have investigated the nucleation of secondary bubbles induced by the rarefaction wave that is produced when the shock wave emitted by the laser-induced plasma reflects at the drop surface. Interestingly, three-dimensional clusters of cavitation bubbles are observed. Their shape is compared with the negative pressure distribution computed with a computational fluid dynamics model and allows us to estimate a cavitation threshold value. In particular, we observed that the focusing of the waves in the vicinity of the free surface can give rise to explosive cavitation events that end up in fast liquid ejections. High-speed recordings of the drop/bubble dynamics are complemented by the velocity and pressure fields simulated for the same initial conditions. The effect of the proximity of a curved free surface on the jetting dynamics of the bubbles was qualitatively assessed by classifying the cavitation events using a non-dimensional stand-off parameter ${\Upsilon\hskip -1,05em -\,}$ that depends on the drop size, the bubble maximum radius and the relative position of the bubble inside the drop. Additionally, we studied the role of the drop's curvature by implementing a structural similarity algorithm to compare cases with bubbles produced near a flat surface to the bubbles inside the drop. Interestingly, this quantitative comparison method indicated the existence of equivalent stand-off distances at which bubbles influenced by different boundaries behave in a very similar way. The oscillation of the laser-induced bubbles promotes the onset of Rayleigh–Taylor and Rayleigh–Plateau instabilities, observed on the drop's surface. This phenomenon was studied by varying the ratio of the maximum radii of the bubble and the drop. The specific mechanisms leading to the destabilisation of the droplet surface were identified.

We study the gravity-induced collisions of charged spheres of dielectric materials dispersed in a gaseous medium. When the gap thickness between the surfaces of two spheres is shorter than the mean free path of the surrounding fluid medium, continuum assumptions for the hydrodynamics interactions are no longer valid, and the non-continuum lubrication interactions result in surface-to-surface contact in finite time. Two like-charged dielectric spheres attract each other at close separations for a wide range of size and charge ratio values. We use trajectory analysis to calculate the collision rate and, thus, explore the role of electrostatic interactions in the collision dynamics of a pair of like-charged dielectric spheres. We present the modifications of pair trajectories due to electrostatic forces and show how collision efficiencies vary with the non-dimensional parameter capturing the relative strength of the electrostatic force to gravity as well as the charge ratio and size ratio.

In direct numerical simulations (DNS) of turbulent Couette flow, the observation has been made that the long streamwise rolls increase in length with the Reynolds number (Lee & Moser, J. Fluid Mech., vol. 842, 2018, pp. 128–145). To understand this, we employ both linear stability theory and its extension to resolvent analysis. For this, we emphasise the high Reynolds number ($Re \rightarrow \infty$) and small streamwise wavenumbers ($\alpha \rightarrow 0$) limit, imposing the distinguished limit $Re_{\alpha }=Re \, \alpha = O(1)$. We find that in case of linear stability theory, $Re_{\alpha }$ acts as a global invariant in the resulting eigenvalue problem, while in case of resolvent analysis, $Re_{\alpha }$ acts as a local invariant in the behaviour of the energy of the system characterised through the first singular value $\sigma _1$ of the resolvent operator within the investigated asymptotic limit. In order to obtain constant streamwise structures for increasing Reynolds numbers, the respective streamwise wavenumber has to decrease, which verifies the observations from DNS studies of an increasing length of the streamwise structures with the Reynolds number. In linear stability theory, a parameter reduction is achieved for the above asymptotic limit, resulting in the modified Orr–Sommerfeld and Squire equations being dependent only on $Re_{\alpha }$. The behaviour of both the coherent structures obtained from linear stability theory and resolvent analysis are compared with each other and show similar behaviours over $Re_{\alpha }$.

It is known that the disintegration of vertical liquid curtains (sheets) is affected crucially by the amplification of free edge holes forming inside the curtain. This paper aims to investigate the influence of the hole expansion dynamics, driven by the so-called rim retraction, on the breakup of a liquid curtain, in both supercritical (Weber number $We > 1$) and subcritical ($We < 1$) conditions. The analysis is based on three-dimensional direct numerical simulations. For a selected supercritical configuration, the steady flow topology is first analysed. The investigation reveals the classic triangular shape regime of the steady curtain, due to the surface-tension-induced borders retraction towards its centre plane. The unsteady dynamics is then investigated as the curtain response to a hole perturbation introduced artificially in the steady flow configuration. The hole evolution determines a rim retraction phenomenon inside the curtain, which is influenced by both capillary and gravity forces. In supercritical conditions, the hole does not influence the curtain flow dynamics in the long-time limit. By reducing the Weber number slightly under the critical threshold ($We=1$), the initial amplification rate of the hole area increases, due to the stronger retraction effect of surface tension acting on the hole rims. The free hole expansion in fully subcritical conditions ($We < 1$) is investigated finally by simulating an edge-free curtain flow. As $We$ decreases progressively, the hole expands while it is convected downstream by gravity acceleration. In the range $0.4< We<1$, the subcritical curtain returns to the intact unperturbed configuration after the hole expulsion at the downstream outflow. For $We<0.4$, the surface tension force becomes strong enough to reverse the gravitational motion of the hole top point, which retracts upstream towards the sheet inlet section while expanding along the lateral directions. This last phenomenon causes finally the breakup of the curtain, which results in a columnar regime strictly resembling similar experimental findings of the literature.

The influence of surface roughness on transition to turbulence in a Mach 4.5 boundary layer is studied using direct numerical simulations. Transition is initiated by the nonlinearly most dangerous inflow disturbance, which causes the earliest possible breakdown on a flat plate for the prescribed inflow energy and Mach number. This disturbance primarily comprises two normal second-mode instability waves and an oblique first mode. When localized roughness is introduced, its shape and location relative to the synchronization points of the inflow waves are confirmed to have a clear impact on the amplification of the second-mode instabilities. The change in modal amplification coincides with the change in the height of the near-wall region where the instability wave speed is supersonic relative to the mean flow; the net effect of a protruding roughness is destabilizing when placed upstream of the synchronization point and stabilizing when placed downstream. Assessment of the effect of the roughness location is followed by an optimization of the roughness height, abruptness and width with the objective of achieving maximum transition delay. The optimization is performed using an ensemble-variational (EnVar) approach, while the location of the roughness is fixed upstream of the synchronization points of the two second-mode waves. The optimal roughness disrupts the phase of the near-wall pressure waves, suppresses the amplification of the primary instability waves and mitigates the nonlinear interactions that lead to breakdown to turbulence. The outcome is a sustained non-turbulent flow throughout the computational domain.

It is known that inertial lift forces can lead to particle focusing in channel flows; yet oscillatory straining effects have also been suggested as a mechanism for particle focusing in wavy channels. To explore the synergy between these two mechanisms, we analytically and experimentally investigate the focusing behaviour of rigid neutrally buoyant particles in a wavy channel. We decompose the particle-free channel flow into a primary Poiseuille flow and secondary eddies induced by the waviness. We calculate the perturbation of the particle on the particle-free flow and the resulting lateral lift force exerted on the particle using the method of matched asymptotic expansions. This yields a zeroth-order lift force arising from the Poiseuille flow and a first-order lift force due to the waviness of the channel. We further incorporate the inertial lift force into the Maxey–Riley equation and simulate particle trajectories in wavy channels. The interactions between the zeroth-order lift force and the particle-free flow largely determine the focusing locations. Experiments in wavy channels with varying amplitudes at channel Reynolds numbers ranging from 5 to 250 are consistent with the predictions of the focusing locations, which are mainly governed by the channel Reynolds number as well as the competition between the inertial lift and the oscillatory straining effects.

The formation process of the leading vortex ring in starting jets with uniform background co- and counter-flow has been studied numerically for $-0.5\leq R_v\leq 0.5$, where $R_v$ is the ratio of background velocity to jet velocity. For the cases with background counter-flow, the normal formation process of the leading vortex ring would be destroyed when $R_v<-0.4$, i.e. the trailing jet would overtake the leading vortex ring through the centre, a phenomenon reminiscent of vortex leapfrogging. As the velocity ratio $R_v$ increases, the formation number $F_{t^*}$ decreases from $9.6$ at $R_v=-0.4$ to $1.92$ at $R_v=0.5$. An analytical model based on the kinematic criterion has been developed so as to describe the relationship between the formation number $F_{t^*}$ and velocity ratio $R_v$. A linear relationship between the vortex core parameter and stroke ratio of starting jet ($\varepsilon \sim k_1L/D$) for the Norbury vortex ring has been established and used effectively to close the model. For co-flow with $0< R_v\leq 0.5$, the results from this model are consistent with the present numerical simulation and the experiments by Krueger et al. (J. Fluid Mech., vol. 556, 2006, pp. 147–166). For counter-flow, two different equations are proposed for $-0.4\leq R_v\leq -0.2$ and $-0.2< R_v<0$, respectively.

We investigate the dynamics of a low-density round jet, with a focus on the mechanisms governing the turbulent momentum and mass transfers as well as on the entrainment of ambient fluid. To that purpose, we combine a theoretical analysis, laboratory experiments and numerical simulations. The theoretical analysis relies on a general formulation of the entrainment decomposition for the case of large density differences, revealing the role of the processes contributing to the entrainment: turbulent kinetic energy production and variation in the shape of the mean velocity radial profiles. The spatial evolution of these terms has been evaluated by means of challenging experiments, providing a unique data set of combined velocity and density statistics of a low-density jet and an air jet. The same flows are investigated by means of large-eddy simulation (LES). Other than for providing complementary information on flow statistics, LES is here used to investigate the role of varying conditions imposed at the source, notably concerning the shape of the inlet velocity profile and the presence of a bottom wall surrounding the source. Experimental and numerical results provide clear insight on how a reduced density within the jet enhances the turbulent kinetic energy production (compared to an iso-density jet) and modifies the shape of the mean velocity profiles. Despite its clear influence on the flow statistics, the reduced density has overall little influence on the entrainment rate, which also shows little sensitivity to varying source conditions.

Two-dimensional numerical simulations with the particle tracking method were conducted to analyse the dispersion behind the detonation front and its mean structure. The mixtures were 2H$_2$–O$_2$–7Ar and 2H$_2$–O$_2$ of increased irregularity in ambient conditions. The detonation could be described as a two-scale phenomenon, especially for the unstable case. The first scale is related to the main heat release zone, and the second where some classical laws of turbulence remain relevant. The dispersion of the particles was promoted by the fluctuations of the leading shock and its curvature, the presence of the reaction front, and to a lesser extent transverse waves, jets and vortex motion. Indeed, the dispersion and the relative dispersion could be scaled using the reduced activation energy and the $\chi$ parameter, respectively, suggesting that the main mechanism driving the dispersion came from the one-dimensional leading shock fluctuations and heat release. The dispersion within the induction time scale was closely related to the cellular structure, particles accumulating along the trajectory of the triple points. Then, after a transient where the fading transverse waves and the vortical motions coming from jets and slip lines were present, the relative dispersion relaxed towards a Richardson–Obukhov regime, especially for the unstable case. Two new Lagrangian Favre average procedures for the gaseous detonation in the instantaneous shock frame were proposed and the mean profiles were compared with those from Eulerian procedure. The characteristic lengths for the detonation were similar, meaning that the Eulerian procedure gave the mean structure with a reasonable accuracy.

We experimentally study front propagation in a vortex lattice providing closed steady cellular flows and no mean flow. To this end, we trigger an autocatalytic reaction in a solution stirred by magnetohydrodynamic flows in a Hele-Shaw cell. We evidence a scale-invariant regime below some flow magnitude and a scale-dependent regime above, the scales referring here to the vortex scale and the front thickness. The transition between these regimes corresponds to a unitary Damköhler number $Da$: $Da=1$. The enhancement of the mean front velocity with the flow magnitude nicely agrees with the literature on numerical simulations and theoretical analyses in the scale-invariant regime $Da>1$, but displays noticeable discrepancies in the scale-dependent one $Da<1$. This shows that the transition between regimes is qualitatively sharp but quantitatively smooth.

Despite the widespread occurrence of pendant drops in nature, there is still a lack of combined studies on their dynamic and static stability. This study focuses on the dynamic and static stability of elongated drops with either a free or pinned contact line on a plane. We first examine static stability for both axisymmetric and non-axisymmetric perturbations subject to volume or pressure constraints. The stability limits for volume and pressure disturbances (axisymmetric) correspond to the maximum volume and pressure of the drops, respectively. Drops with free contact lines are marginally stable to non-axisymmetric perturbations because of their horizontal translational invariance, whereas pinned drops are stable. The linear dynamic stability is then investigated numerically through a boundary element model, restricted to volume disturbances. Results show that when the stability limit is reached, the first zonal mode has a zero frequency, suggesting that the thresholds for static and dynamic stability are essentially equivalent. Furthermore, natural frequencies experience sharp changes as the stability limit is approached. Another zero frequency mode associated with the horizontal motion of the centre of mass is also revealed by the numerical results, reflecting the horizontal translational invariance of drops with free contact lines. Finally, the frequency spectrum modified by gravity is explored, resulting in the identification of five gravity-induced frequency shift patterns. The frequency shifts break the spectral degeneracy for hemispherical drops with free contact lines, leading to various spectral orderings according to polar and azimuthal wavenumbers.

We investigate the dynamics and energy production capability of a flexible piezoelectric plate submerged close to the free surface and exposed to incident head gravity waves and current. A theoretical model is derived in which the flag and its wake are represented with a vortex line while the body of the fluid is considered to be inviscid. The model is employed to describe the hydrodynamic interactions between a flexible plate, its wake, gravity incident waves and the current. The model reveals two distinct vibration states of a piezoelectric device corresponding to almost similar optimal energy production levels. The first is associated with the cantilever fluttering mode of the plate, with limited dependency on the plate's flexibility across different Froude numbers and incoming wave frequencies. The other resembles the flow-induced flapping mode in more flexible plates, with the energy output showing a higher dependency on plate flexibility. The concurrent existence of these two energetic modes allows adjustment of the plate length to consistently achieve the maximum energy production level across different flow conditions. The role of the Froude number of the system's responses is explored and correlated to the appearance of gravity wave groups on the surface, each propagating with a different wavenumber. It is shown that a submergence depth of less than half of the body length is required to reach a high energetic condition in subcritical and critical flows. Finally, the optimal inductive and resistive values are related to proper matching between flow, mechanical and electrical time scales.