Using a combination of mean flow spatial linear stability and two-dimensional volume-of-fluid (VoF) simulations, the physics governing the instability of high-speed liquid sheets being injected into a quiescent gas environment is studied. It is found that the gas shear layer thickness $\delta _G$ plays an influential role, where for values $\delta _G/H\lesssim 1/8$, the growth of sinuous and varicose modes is nearly indistinguishable. Here, $H$ is the liquid sheet thickness. With larger values of $\delta _G/H$, a second peak develops in the lower wavenumber region of the dispersion relation, and becomes increasingly dominant. This second peak corresponds to a large-scale sinuous mode, and its critical wavelength $\lambda _{crit,sinuous}$ is found to scale as $\lambda _{crit,sinuous}/H = 14.26 (\delta _G/H)^{0.766}$. This scaling behaviour collapses onto a single curve for various combinations of the liquid-based Reynolds ($Re_L$) and Weber ($We_L$) numbers, provided that $\delta _G/H > O({10^{-1}})$. For the varicose modes, the shape of the dispersion relation does not change with variations in $\delta _G/H$, and the liquid shear layer thickness has an almost negligible influence on the growth of instabilities. Two-dimensional VoF simulations are employed to examine the validity of the linear stability assumptions. These simulations also show that the dominant sinuous mode remains active as the process transitions into the nonlinear regime, and that this mode is ultimately responsible for fragmenting the sheet. Based on an energy budget analysis, the most influential contributors to the growth of the sinuous mode are the gas Reynolds shear stress and the lateral working of pressure on the gas side.

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths, vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn., vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids, vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. (Phys. Fluids, vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

Supersonic flow over a hollow cylinder/flare with a free-stream Mach number of 2.25 is numerically investigated in this study. Axisymmetric computational fluid dynamics simulations and global stability analysis (GSA) are performed for a wide range of cylinder radii and flare deflection angles. The onset of incipient and secondary separation is delayed as the cylinder radius is decreased due to the axisymmetric effects. The GSA reveals that a decrease in cylinder radius also postpones the emergence of global instability. The GSA results agree well with the results of direct numerical simulations for a supercritical case in the linear stage. The saturated flow exhibits pairs of unsteady streamwise streaks downstream of reattachment. The criterion of the global stability boundary established for supersonic flow over a compression corner (Hao et al., J. Fluid Mech, vol. 919, 2021, A4) is extended to its axisymmetric counterpart.

Numerical simulations of multiphase flows are crucial in numerous engineering applications, but are often limited by the computationally demanding solution of the Navier–Stokes (NS) equations. The development of surrogate models relies on involved algebra and several assumptions. Here, we present a data-driven workflow where a handful of detailed NS simulation data are leveraged into a reduced-order model for a prototypical vertically falling liquid film. We develop a physics-agnostic model for the film thickness, achieving a far better agreement with the NS solutions than the asymptotic Kuramoto–Sivashinsky (KS) equation. We also develop two variants of physics-infused models providing a form of calibration of a low-fidelity model (i.e. the KS) against a few high-fidelity NS data. Finally, predictive models for missing data are developed, for either the amplitude, or the full-field velocity and even the flow parameter from partial information. This is achieved with the so-called ‘gappy diffusion maps’, which we compare favourably to its linear counterpart, gappy POD.

We present a model for the volume-averaged dissipation rate in linear unsteady flow through porous media. The model is derived by blending a new small-time asymptotic expression for the dissipation rate obtained from boundary layer theory with the known large-time asymptotic expression obtained from Darcy's law. The resulting model is a second-order Volterra functional of the volume-averaged acceleration. We validate the model with an analytical solution for transient flow through a porous medium composed of circular tubes and with numerical simulations of transient and oscillatory flow through a cylinder array and through a hexagonal sphere pack.

Two-dimensional compressible flows in radial equilibrium are investigated in the ideal dilute-gas regime and the non-ideal single-phase regime close to the liquid–vapour saturation curve and the critical point. Radial equilibrium flows along constant-curvature streamlines are considered. All properties are therefore independent of the tangential streamwise coordinate. A differential relation for the Mach number dependency on the radius is derived for both ideal and non-ideal conditions. For ideal flows, the differential relation is integrated analytically. Assuming a constant specific heat ratio $\gamma$, the Mach number is a monotonically decreasing function of the radius of curvature for ideal flows, with $\gamma$ being the only fluid-dependent parameter. In non-ideal conditions, the Mach number profile also depends on the total thermodynamic conditions of the fluid. For high molecular complexity fluids, such as toluene or hexamethyldisiloxane, a non-monotone Mach number profile is admissible in single-phase supersonic conditions. For Bethe–Zel'dovich–Thompson fluids, non-monotone behaviour is observed in subsonic conditions. Numerical simulations of subsonic and supersonic turning flows are carried out using the streamline curvature method and the computational fluid dynamics software SU2, respectively, both confirming the flow evolution from uniform flow conditions to the radial equilibrium profile predicted by the theory.

Shock-tube experiments on Richtmyer–Meshkov (RM) instability at a perturbed SF$_6$ layer surrounded by air, induced by a cylindrical divergent shock, are reported. To explore the effects of reverberating waves and interface coupling on instability growth, gas layers with various shapes are created: unperturbed inner interface and sinusoidal outer interface (case US); sinusoidal inner and outer interfaces that have identical phase (case IP); sinusoidal inner and outer interfaces that have opposite phase (case AP). For each case, three thicknesses are considered. Results show that reverberating waves inside the layer dominate the early-stage instability growth, while interface coupling dominates the late-stage growth. The influences of waves on divergent RM instability are more pronounced than the planar and convergent counterparts, which are estimated accurately based on gas dynamics theory. Both the wave influence and interface coupling depend heavily on the layer shape, leading to diverse growth rates: the quickest growth for case AP, medium growth for case US, the slowest growth for case IP. In particular, for the IP case, there exists a critical thickness below which the instability growth is suppressed by both the reverberating waves and interface coupling. This provides an efficient way to modulate the growth of divergent RM instability. It is found that divergent RM instability involves weak nonlinearity and strong interface coupling such that the linear theory of Mikaelian (Phys. Fluids, vol. 17, 2005, 094105) can well reproduce the instability growth at late stages for all cases. This constitutes the first experimental confirmation of the Mikaelian theory.

Theoretical investigation of the primary Mach reflection (MR) configuration on V-shaped blunt leading edges (VBLEs) forms the focus of this study. By ignoring the secondary interactions, a theoretical method based on a simplified form of the continuity relation is developed to predict the shock configurations, including the detached shock, the Mach stem, the transmitted shock and the triple point. The comparison of the theoretical results with both numerical and previous experimental results shows the reliability of the theoretical approach in predicting shock structures across a wide range of free stream and geometric parameters. The theoretical model provides a detailed comprehension of the occurrence mechanism of inverse MRs on VBLEs and the influence of the free stream and geometric parameters on primary MR configurations. Along with the primary MR configuration, the curved shock or compression waves generated by the crotch are solved and offer insight into the transition from the MR to the regular reflection from the same family (sRR). The increase of the ratio $R/r$ and the free stream Mach number $M_0$ appears to facilitate the transition, while the effect of the half-span angle $\beta$ is non-monotonic. The predicted shock positions allow for the identification of the transition boundary between the primary MR and sRR. It is found that $R/r$ below a threshold (for a set $M_0$ value) produces MR, irrespective of $\beta$. If this threshold is exceeded, the configuration can transition from the primary MR to sRR and then back to the primary MR as $\beta$ increases.

Using direct numerical simulations, we investigate the heat transport in bulk and boundary flows separately in rotating Rayleigh–Bénard convection in cylindrical cells. In the bulk we observe a steep scaling relationship between the Nusselt number ($Nu$) and the Rayleigh number ($Ra$), which is consistent with the results from simulations using periodic boundary conditions. For the boundary flow, we observe a power law $Nu_{BF}\sim (Ra/Ra_w)^1$ at the leading order, where $Nu_{BF}$ is the local Nusselt number of the boundary flow and $Ra_w$ is the onset Rayleigh number of the wall mode. We develop a model using the boundary layer marginal stability theory to explain this power law, and further show that a more precise description of the data can be obtained if a higher-order correction is introduced. A striking finding of our study is the observation of a sharp transition in flow state, manifested by a sudden drop in $Nu_{BF}$ with a corresponding collapse of the boundary flow coherency. After the transition, the boundary flow breaks into vortices, leading to a reduction in flow coherency and heat transport efficiency. As the physical properties of the vortices should not depend on the aspect ratio, $Nu_{BF}$ for all aspect ratios collapse together after the transition. Moreover, the centrifugal force helps trigger the breakdown of the coherent boundary flow state. For this reason, $Nu_{BF}$ for the cases with non-zero centrifugal force collapse together. We further develop a method that enables us to separate the contributions from the bulk and boundary flows in the global Nusselt number using only the global $Nu$ and it does not require the centrifugal force to be absent.

Inertial particles in wall-bounded turbulence are known to form streaks, but experimental evidence and predictive understanding of this phenomenon is lacking, especially in regimes relevant to atmospheric flows. We carry out wind tunnel measurements to investigate this process, characterizing the transport of microscopic particles suspended in turbulent boundary layers. The friction Reynolds number $Re_\tau = {O}(10^4)$ allows for significant scale separation and the emergence of large-scale motions, while the range of viscous Stokes number $St^+=18$–870 is relevant to the transport of dust and fine sand in the atmospheric surface layer. We perform simultaneous imaging of both carrier and dispersed phases along wall-parallel planes in the logarithmic layer, demonstrating that streamwise particle streaks largely overlap with large-scale low-speed flow regions. The fluid–particle slip velocity indicates that with increasing inertia, the particle streaks outlive the low-speed fluid streaks. Moreover, two-point statistics show that the width of the particle streaks increases linearly with Stokes number, bounded by the size of the coherent flow structures. Finally, the particle-sampled flow topology suggests that particle streaks reside between the legs of hairpin packets. From these observations, we infer a conceptual view of the formation of particle streaks in the frame of the attached eddy model. A scaling for the particle streaks’ width is derived as a function of $Re_\tau$ and $St^+$, which reproduces the measured trends and predicts widths ${O}(0.1)$ m in the atmospheric surface layer, comparable to aeolian streamers observed in the field.

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

We develop a high-fidelity multi-physics model and, by using it, study the energy conversion process of a piezohydroelastic flag. Instead of a regular flag with a clamped upstream leading edge, we use an inverted flag so as to make the best of fluid elastic instability for energy harvesting. Moreover, different from many previous studies where a resistor–capacitor (RC) circuit is usually used, we adopt a resistor–inductor–capacitor (RLC) circuit for electricity generation. The influences of several key parameters associated with fluid, structure and electric dynamics are studied. Significantly different response modes are identified, among which the symmetric- and asymmetric-flutter modes are most suitable for sustainable energy harvesting, both emerging with moderate bending stiffness. If only deploying the RC circuit, increasing the resistance makes the flag more stable. By adding an inductor to turn the RC circuit into an RLC one, we observe the occurrence of ‘lock-in’ between the flag frequency and the circuit frequency as first reported on a regular flag by Xia et al. (Phys. Rev. Appl., vol. 3, 2015, 014009). This phenomenon can significantly enhance the energy output, but it only happens when the circuit resistance is sufficiently large. From a derivation based on dynamic mode decomposition analysis, we further identify an optimal condition for maximizing the energy output, which can serve as a guideline to determine whether deploying an inductor can boost the performance, and, if yes, the required inductance. The findings from this study can better guide the design of flow-induced-vibration-based piezoelectric energy harvesters for microelectronic devices.

We present a method for accurately determining the stability characteristics of spatially modulated shear layers. The algorithm can handle arbitrary commensurate states, which are not accessible to classical direct-numerical-simulation-based approaches. It uses spectral discretization of the field equations to handle field modulations and the spectrally accurate immersed boundary conditions method to handle the geometry modulations. The algorithm can deal with pattern interaction effects driven by modulations of different physical origins. Various tests demonstrate that the algorithm delivers spectral accuracy for eigenvalues and eigenfunctions. The algorithm can be easily extended to analyse many sources and patterns of modulation with minimal commitment to the user's time.

Time-varying flow separation on an accelerating prolate spheroid has been studied at various angles of incidence. Instantaneous pressure and scanning stereoscopic particle image velocimetry were used to shed light on the evolution of cross-flow structures for the Reynolds number ($Re$) range of $1.0\times 10^6\leq Re \leq 1.5\times 10^6$. The movement of separation lines is examined for various model accelerations to investigate on the interplay between acceleration and flow separation. The results demonstrate that for axial accelerations, the streamwise pressure distribution in the rear part of the prolate spheroid switches from an adverse to a favourable pressure gradient. At the same time, the circumferential adverse pressure gradient present during steady motion vanishes during said accelerations. In contrast, both streamwise and circumferential adverse pressure gradients strengthen when the model is axially decelerated. These dynamic pressure distributions influence the location of the separation line, which in turn moves closer to the model meridian during accelerations while moving outwards during decelerations. The streamwise vorticity distribution and the streamwise circulation both show how the separation-line position impacts the vortex formation. A high-vorticity region near the model surface is established during acceleration. In contrast, a decelerating model leads to transport of high-vorticity fluid into the outer area of the cross-flow separation. We further assess the memory effects following the near-impulsive velocity changes. The cross-flow retains the memory of moving separation lines shortly after the acceleration. However, the separation recovers quickly to a steady state.