This study investigates the correlation between the fluctuating wall heat flux, and the distribution and transport of Reynolds shear stress and turbulent heat flux in compressible boundary layers at Mach number 5.86 and friction Reynolds number 420, with a relatively weaker and a stronger wall cooling imposed. As illustrated from the probability density functions of the wall-heat-flux perturbations, with increasing wall cooling, the extreme wall heat flux is intensified and tends to be more negatively skewed. To examine the role of the extreme events in the transport of the momentum and heat, conditional analysis of the extreme positive and negative wall-heat-flux-perturbation events is conducted. In most regions of the boundary layer, the positive events are predominantly associated with an increase in Reynolds shear stress and a decrease in turbulent heat flux. Joint probability density functions of velocity and wall-heat-flux perturbations in the near-wall region indicate that the extreme positive events tend to be more correlated with ejections, which is particularly evident in the stronger wall-cooling case. To further shed light on the underlying mechanisms of the connections between wall heat flux and transport budgets, a transport equation for turbulent heat flux is derived, in a similar manner to that for Reynolds shear stress. The energy balance is inspected, with conditional analysis applied to budget terms and mean flow properties so as to quantify the correlation between wall-heat-flux fluctuations and energy evolution.

We study the stability and dewetting dynamics of a thin free-surface film composed of two miscible liquids placed on a solid substrate. Our study focuses on the development of a self-consistent model such that the mixture concentration influences both free-surface and wetting energies. By assuming a simple relation between these energies and the bulk and surface concentrations, we analyse their effect on the concentration distribution and dewetting down to the equilibrium film thickness determined by the fluid–solid interaction potential. The model, developed within the gradient dynamics formulation, includes the dependence of the free-surface energy on surface concentration leading to the Marangoni effect, while a composition-dependent Hamaker constant describes the wetting energy resulting from the fluid–solid interaction. We analyse the restrictions that must be fulfilled to ensure an equilibrium state for a flat film of a binary fluid. Then, we proceed by studying its linear stability. First, we consider the Marangoni effect while assuming that wetting energy depends only on the fluid thickness. Then, we include a dependence of wetting energy on concentration and study its effects. We find that the linear stability results compare very well with those of numerical simulations of the full nonlinear problem applied to the particular case of a binary melted metal alloy, even close to breakup times. Therefore, in practice, most of the evolution can be studied by using the linear theory, simplifying the problem considerably.

The boundary-layer stability on a section of a rotating wind turbine blade with an FFA-W3 series aerofoil at a chord Reynolds number of $3 \times 10^5$, with varying rotation and radii, is studied with direct numerical simulations and linear stability analyses. Low rotation does not significantly affect transition in the outboard blade region. The relative insensitivity to rotation is due to a laminar separation bubble near the leading edge, spanwise-deformed by a primary self-excited instability, promoting the secondary absolute instability of the Kelvin–Helmholtz (KH) vortices and rapid transition. Moderate increases in rotation, or moving inboard, stabilise the flow by accelerating the attached boundary layer and possibly inducing competition between cross-flow and KH modes. This delays separation and transition. Initially, for high rotation rates or radial locations close to the hub, transition is delayed. Nevertheless, strong stationary and travelling cross-flow modes are eventually triggered, spanwise modulating the KH rolls and shifting the transition line close to the leading edge. Cross-flow velocities as high as $56\,\%$ of the free stream velocity directed towards the blade tip are reached at the transition location. For radial locations farther from the hub, the effective angle of attack is decreased, and cross-flow transition occurs at lower rotation rates. The advance or delay of the transition line compared with a non-rotating configuration depends on the competing rotation effects of stabilising the attached boundary layer and triggering cross-flow modes in the separation flow region.

Stratified wakes past an isolated conical seamount are simulated at a Froude number of $Fr = 0.15$ and Rossby numbers of $Ro = 0.15$, 0.75 and $\infty$. The wakes exhibit a Kármán vortex street, unlike their unstratified, non-rotating counterpart. Vortex structures are studied in terms of large-scale global modes, as well as spatially localised vortex evolution, with a focus on rotation effects. The global modes are extracted by spectral proper orthogonal decomposition (SPOD). For all three studied $Ro$ ranging from mesoscale, submesoscale and non-rotating cases, the frequency of the SPOD modes at different heights remains coupled as a global constant. However, the shape of the SPOD modes changes from slanted ‘tongues’ at zero rotation ($Ro=\infty$) to tall hill-height columns at strong rotation ($Ro=0.15$). A novel method for vortex centre tracking shows that, in all three cases, the vortices at different heights advect uniformly at approximately $0.9U_{\infty }$ beyond the near wake, consistent with the lack of variability of the global modes. Under system rotation, cyclonic vortices and anticyclonic vortices (AVs) present considerable asymmetry, especially at $Ro = 0.75$. The vorticity distribution as well as the stability of AVs are tracked downstream using statistics conditioned to the identified vortex centres. At $Ro=0.75$, intense AVs with relative vorticity up to $\omega _z/f_{c}=-2.4$ (where $\omega_z$ is the vertical vorticity and $f_c$ is the Coriolis frequency) are seen with small regions of instability and they maintain large $\omega _z/f_{c}$ magnitude in the far wake. Recent stability analysis that accounts for stratification and viscosity is found to improve on earlier criteria and show that these intense AVs are stable.

A novel theoretical model for bubble dynamics is established that simultaneously accounts for the liquid compressibility, phase transition, oscillation, migration, ambient flow field, etc. The bubble dynamics equations are presented in a unified and concise mathematical form, with clear physical meanings and extensibility. The bubble oscillation equation can be simplified to the Keller–Miksis equation by neglecting the effects of phase transition and bubble migration. The present theoretical model effectively captures the experimental results for bubbles generated in free fields, near free surfaces, adjacent to rigid walls, and in the vicinity of other bubbles. Based on the present theory, we explore the effect of the bubble content by changing the vapour proportion inside the cavitation bubble for an initial high-pressure bubble. It is found that the energy loss of the bubble shows a consistent increase with increasing Mach number and initial vapour proportion. However, the radiated pressure peak by the bubble at the collapse stage increases with decreasing Mach number and increasing vapour proportion. The energy analyses of the bubble reveal that the presence of vapour inside the bubble not only directly contributes to the energy loss of the bubble through phase transition but also intensifies the bubble collapse, which leads to greater radiation of energy into the surrounding flow field due to the fluid compressibility.

The Taylor–Maccoll (T–M) equations are the governing equations for steady inviscid irrotational axisymmetric conical flow, and have been widely applied to the design of waveriders and intakes. However, only four classic solutions have been reported: external conical flow (ECF), Busemann flow and internal conical flow of types A and B (ICFA and ICFB). In this work, the analysis of the T–M equations clarifies all possible solutions and reveals their relations. The domain where elementary solutions exist is divided into four domains. The classic Busemann and ICFB solutions share the same elementary solution as the template in a domain called the pre-shock domain, while the classic ECF and ICFA solutions belong to a domain named the ECF domain. Two new solutions, the inner flow of ECF (IECF) and degenerate conical flow (DCF), are found in the domains named after the corresponding solutions, namely the IECF and DCF domains. The IECF behaves as the mass injection supporting the classic ECF on an imaginary cone surface, while the DCF behaves as the conical expansion of a uniform flow. Furthermore, possible combinations of pre-shock solutions and supersonic post-shock solutions are clarified. The classic solutions are special cases where the pre-/post-shock solutions are combined with uniform flows. In general, the Busemann and ICFB solutions can be combined with any post-shock solutions in accord with the shock relations, including the ECF, ICFA, IECF and DCF solutions. In addition, numerical analyses are conducted to verify the validity of the two new solutions, DCF, IECF and one combined solution Busemann–ECF.

Explosive dispersal of granular media widely occurs in nature across various length scales, enabling engineering applications ranging from commercial or military explosive systems to the loss prevention industry. However, the correlation between the explosive dispersal behaviour and the structure of dispersal system is far from completely understood, thereby compromising the prediction of the explosive dispersal outcome resulting from a specific dispersal system. Here, we investigate the dispersal behaviours of densely packed particle rings driven by the enclosed pressurized gases using coarse-grained computational fluid dynamics–discrete parcel method. Distinct dispersal modes emerge from the dispersal systems with vastly varying sets of the macro- and micro-scale structural parameters in terms of the dispersal completeness and the spatial uniformity of the dispersed mass. Further investigation reveals the variation in the dispersal modes arises from the collective effects of multiscale gas–particle coupling relationships. Specifically, the macroscale coupling dictates the cyclic momentum/energy transfer between gases and particle ring as an entirety. The mesoscale coupling relates to the inter-pore gas filtration through the thickness of the particle ring, leading to the mass/energy reduction of the explosive source. The microscale coupling involves the individual particle dynamics influenced by the local flow parameters. A persistent macroscale coupling results in an incomplete dispersal which takes the form of an aggregated annular band, whereas the meso- and micro-scale couplings alter the macroscale coupling to a different extent. By incorporating the effects of the variety of structural parameters on the multiscale gas–particle coupling relationships, a non-dimensional parameter referred to as the modified mass ratio is constructed, which shows an explicit correlation with the dispersal mode. We proceed to establish a dispersal ring model in the continuum frame which accounts for the macro and meso-scale coupling effects. This model proves to be capable of successfully predicting the ideal and validated failed dispersal modes.

Unlike in solids, heat transfer in fluids can be greatly enhanced due to the presence of convection. Under gravity, an unevenly distributed temperature field results in differences in buoyancy, driving fluid motion that is seen in Rayleigh–Bénard convection (RBC). In RBC, the overall heat flux is found to have a power-law dependence on the imposed temperature difference, with enhanced heat transfer much beyond thermal conduction. In a bounded domain of fluid such as a cube, how RBC responds to thermal perturbations from the vertical sidewall is not clear. Will sidewall heating or cooling modify flow circulation and heat transfer? We address these questions experimentally by adding heat to one side of the RBC. Through careful flow, temperature and heat flux measurements, the effects of adding side heating to RBC are examined and analysed, where a further enhancement of flow circulation and heat transfer is observed. Our results also point to a direct and simple control of the classical RBC system, allowing further manipulation and control of thermal convection through sidewall conditions.

To characterize fluid flow in the slip regime, the use of Navier–Stokes–Fourier (NSF) equations with slip boundary conditions is prevalent. This trend underscores the necessity of developing reliable and accurate slip boundary conditions. According to kinetic theory, slip behaviours are intrinsically linked to the gas scattering processes at the surface. The widely used Maxwell scattering model, which employs a single accommodation coefficient to describe gas scattering processes, reveals its limitations when the difference between accommodation coefficients in the tangential and normal directions becomes significant. In this work, we provide a derivation of velocity slip and temperature jump boundary conditions based on the Cercignani–Lampis–Lord scattering model, which applies two independent accommodation coefficients to describe the gas scattering process. A Knudsen layer correction term is introduced to account for the impact of the surface on the velocity distribution function, which is associated with the scattering model. The governing equation of the correction term is established based on the linearized Boltzmann equation. Additionally, two moments are derived to capture the collision effect in the Knudsen layer: a conserving moment of collision invariants, and an approximate higher-order conserving moment. These moments are then employed to determine the coefficients in the correction term. We demonstrate that the derived slip coefficients align closely with numerical results obtained by solving the Boltzmann equation in the Knudsen layer. Besides, we apply the derived slip boundary conditions within the framework of the NSF equations, yielding numerical results that exhibit excellent consistency with those obtained through molecular-level simulations.

The effects of reshock conditions, including the interface evolution state before reshock and the second shock intensity, on interface instability induced by two successive shocks propagating in the same direction are investigated via shock-tube experiments. It is observed that the reshock promotes the interface instability, and the post-reshock perturbation evolution relates to both the pre-reshock interface evolution state and second shock intensity. For the linear evolution of the twice-shocked interface, existing models perform poorly when either the pre-reshock interface shape effect or the secondary compression effect is pronounced, as current reduction factors fail to accurately describe these effects. Besides, the reshock-induced linear amplitude growth rate shows a non-monotonic dependence on the scaled pre-reshock amplitude, primarily due to the shape effect of the pre-reshock interface. For the post-reshock nonlinear evolution, the model proposed by Zhang & Guo (J. Fluid Mech., vol. 786, 2016, pp. 47–61) offers reasonable predictions when the second shock is weak. However, when the second shock is moderately strong, the model overestimates the bubble growth and underestimates the spike evolution under the influence of the significant secondary compression effect. Furthermore, empirical linear and nonlinear models capable of describing the dependence of the post-reshock evolution on reshock conditions are proposed based on the present experimental results and existing models.

Several transition scenarios are present in a hypersonic compression-ramp flow. In our previous work (Cao et al., J. Fluid Mech., vol. 941, 2022, p. A8), a complete transition process induced by the global instability of a compression-ramp flow was revealed. In a globally stable flow, however, the transition to turbulence can be promoted by convective instabilities, which is the focus of this work. The same flow conditions as in our previous work (Mach number 7.7, Reynolds number $8.6\times 10^5$ based on the flat-plate length) are considered here. Owing to a smaller ramp angle, a weakly separated flow forms on the compression ramp, which supports no global instability. Resolvent analysis identifies low-frequency streamwise streaks as the optimal response of base flow to upstream forcing. Local stability analysis reveals Mack's second mode in the boundary layer downstream of reattachment. By introducing random disturbances upstream of separation in direct numerical simulations, we observe breakdown to turbulence downstream of reattachment. Two transition scenarios are revealed, and they are highly dependent on the amplitude of upstream disturbances. For a large amplitude, strong streamwise streaks develop near the reattachment region, which break down to turbulence quickly. However, when the disturbance amplitude is reduced, the second-mode instability dominates the transition to turbulence.