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.

The analysis of Scott (J. Fluid Mech., vol. 741, 2014, pp. 316–349) is implemented numerically. Decaying turbulence is confined to a channel between two infinite, parallel, rotating walls. The Rossby and Ekman numbers are supposed small, the former condition making nonlinearity small, while the latter allows the turbulence to persist for the many rotational periods needed for the small nonlinearity to be effective. The flow is expressed as a combination of inertial waveguide modes, indexed by a two-dimensional wave vector $\boldsymbol{k}$ and an integer n. The $n = 0$ modes form a two-dimensional component of the flow, whereas the remainder is the wave component, on which attention is focused in this article. Assuming statistical axisymmetry and homogeneity in directions parallel to the walls, the second-order moments of the mode amplitudes yield a spectral matrix ${A_{nm}}(k,t)$ (where $k = |\boldsymbol{k} |$), of which the diagonal elements describe the distribution of energy over different modes. Wave-turbulence analysis provides an equation governing the time evolution of ${A_{nn}}$, $n \ne 0$, the wave spectra, which forms the basis for the present work. The initial distribution of energy is Gaussian and depends on a parameter $\varXi$, the initial spectral width. The problem has two other parameters, ${\beta _w}$ and ${\beta _v}$, which correspond to two distinct viscous dissipative mechanisms: wall damping due to boundary layers and volumetric damping by viscous effects throughout the flow. Results obtained by numerical solution include the time evolution of the total wave energy, E, and the detailed description of its distribution over k and n provided by ${A_{nn}}(k)$.

The reflection of a shock pulse at a liquid–gas interface occurs in many applications, from lithotripsy to underwater explosions and additive manufacturing. In linear theory, reflection and transmission at an interface depend only on the impedance difference, but this does not hold for a nonlinear pulse. This work develops an analytical framework for computing the reflection and transmission coefficients for an impulsive shock wave at a liquid–gas interface. The problem is treated analytically by considering idealised pulses and solving a series of consecutive Riemann problems. These correspond to the initial interaction with the interface and important subsequent wave interactions that enable a complete description of the process to be obtained. Comparisons with numerical and existing analytical approaches are made for the case of a water–air interface. In the acoustic limit, the method produces results identical to those of linear acoustic theory. As the pulse strength increases, the proposed method agrees well with numerical simulation results, whereas existing analytical methods that consider only the interface fail. We detail how a reflecting pulse can put water into tension without any incident negative pressure. It is further shown that the magnitude of the reflection coefficient decreases with increasing incident shock pressure, and the reflected pulse widens. Reflections of pulses with positive and negative pressures temporarily create negative pressure regions with greater magnitude than the incident pulse. Finally, we consider non-idealised waves. Comparisons with simulations show that the reflection characteristics can be explained qualitatively using the analytical method, and the reflection coefficients are predicted accurately.

Turbulent convection in the interiors of the Sun and the Earth occurs at high Rayleigh numbers $Ra$, low Prandtl numbers $Pr$, and different levels of rotation rates. To understand the combined effects better, we study rotating turbulent convection for $Pr = 0.021$ (for which some laboratory data corresponding to liquid metals are available), and varying Rossby numbers $Ro$, using direct numerical simulations in a slender cylinder of aspect ratio 0.1; this confinement allows us to attain high enough Rayleigh numbers. We are motivated by the earlier finding in the absence of rotation that heat transport at high enough $Ra$ is similar between confined and extended domains. We make comparisons with higher aspect ratio data where possible. We study the effects of rotation on the global transport of heat and momentum as well as flow structures (a) for increasing rotation at a few fixed values of $Ra$, and (b) for increasing $Ra$ (up to $10^{10}$) at the fixed, low Ekman number $1.45 \times 10^{-6}$. We compare the results with those from unity $Pr$ simulations for the same range of $Ra$ and $Ro$, and with the non-rotating case over the same range of $Ra$ and low $Pr$. We find that the effects of rotation diminish with increasing $Ra$. These results and comparison studies suggest that for high enough $Ra$, rotation alters convective flows in a similar manner for small and large aspect ratios, so useful insights on the effects of high thermal forcing on convection can be obtained by considering slender domains.

Oscillatory boundary layers over flat and rippled seabeds are well described in the literature. However, the presence of protruding vegetation stems has received no theoretical or experimental attention. The present work establishes an analytical constant viscosity model akin to the Stokes oscillatory boundary layer solution and a nonlinear varying-viscosity numerical model with a turbulence closure. The two models are used to describe the importance of vegetation and free stream velocity characteristics on spatially averaged oscillatory boundary layers: their friction factors, thickness and phase leads over the free-stream velocity. The models are periodic in time and resolve boundary and shear layers over the vertical, contrary to past efforts applying two-layer models. The models are extended to investigate the importance of finite wavelengths with steady streaming stresses and their associated mean velocity profile. Steady streaming is quantified both for the near-bottom streaming within the canopy and for the streaming in the shear layer above the canopy. Finally, akin to theoretical and experimental works on mean flows over unvegetated and flat seabeds due to oscillatory and nonlinear free-stream velocities, the numerical model investigates varying degree of nonlinearity for velocity- and acceleration-skewed velocity signals, and it is identified that the presence of vegetation stems gives rise to an additional contribution to the horizontal momentum balance which is not present for unvegetated conditions. Finally, it is discussed how the presence of a free surface, contrary to purely oscillatory conditions, alters the horizontal momentum balance within and above the canopy.

A permeable disk serves as a simplified model for the conversion of wind energy by a horizontal axis wind turbine. In this study, we investigate how inflow turbulence intensity (TI), $I_\infty$, and inflow turbulence integral length scale, $L_\infty$, influence the flow recovery in the wake, the capability of a permeable disk in extracting turbulence kinetic energy (TKE) of the incoming flow, and the statistics of wake-added turbulence using large-eddy simulation. The simulated inflows include various TIs (i.e. $I_\infty =2.5\,\%$–$25\,\%$) and integral length scales (i.e. $L_\infty / D =0.5$–$2.0$) for two thrust coefficients. Simulation results show that both inflow TI and integral length scale influence flow recovery via enhanced ejections and sweeps across the wake boundary, with the former strongly affecting the position where the wake starts to recover and the latter mainly on the recovery rate. Moreover, it is shown that increasing $I_\infty$ and $L_\infty$ increases the TKE extraction by the disk, occurring mainly at scales ($s$) greater than $0.5D$ and frequencies depending on the inflow integral length scale. As for the wake-added TKE, the inflow TI mainly affects its intensity, while the inflow integral length scale affects both its intensity and the sensitive frequencies, with the spectral distributions in scale space ($s$) being similar and the peak located around $s/D=1.0$ for the considered inflows.

Isolated-roughness-induced transitions controlled by local wall heating strips are studied via direct numerical simulation and BiGlobal linear stability analysis. The transition mechanisms are studied first with different wall temperatures. The separated shear layer–counter-rotating vortex system is found to be the main source for transitions. Symmetric and antisymmetric modes are observed in the wake, and the former is dominant. The local wall heating strip can delay the transition, and this effect is enhanced with higher heating temperature, wider strip and a combination of upstream and downstream control strips. The upstream strip lifts up the inlet flow and weakens the wake system in an indirect manner. The antisymmetric mode gradually vanishes, while the symmetric mode always exists but becomes weaker. The downstream strip exhibits a more effective transition delay by directly weakening the separated shear layer and vortex system in the wake. Vorticity transport analysis suggests that the downstream strip increases dissipation for streamwise vorticity and transfers it into wall-normal and spanwise vorticity. BiGlobal analyses indicate that the downstream strip shows less influence on the peak growth rate of the symmetric mode but significantly shrinks its unstable region. Analyses of the disturbance energy production indicate that the upstream strip wakens the wall-normal and spanwise shear at the same time, but the downstream strip mainly wakens the wall-normal one. More simulations are performed with different roughness heights, point-source disturbance and different roughness shapes. The results show that the current method remains effective enough in delaying transitions at a wide range of conditions.

High-pressure fluid transport in nanoporous media such as shale formations requires further understanding because conventional continuum approaches become inadequate due to their ultralow permeability and complexity of transport mechanisms. We propose a species-based approach for modelling two partially miscible, multicomponent fluids in nanoporous media – one that does not rely on conventional bulk fluid transport frameworks but on species movement. We develop a numerical model for species transport of partially miscible, non-ideal fluid mixtures using the chemical potential gradient as the driving force. The model considers the binary friction concept to include the friction between fluid molecules as well as between fluid molecules and pore walls, and incorporates the key multicomponent transport mechanisms – Knudsen, viscous and molecular diffusion. Under single-phase conditions, the system under consideration is quantified by introducing multicomponent Sherwood number (Sh), Péclet number (Pe) and fluid–solid friction modulus (φ). Despite the complexity of fluid transport in nanopores, the steady-state single-phase transport results reveal the contribution of diffusion in nanopores, where all parameters collapse on a set of master curves for the multicomponent Sh with a dependence on multicomponent Pe and φ. Unsteady state, two-phase transport modelling of the codiffusion process shows that light and intermediate alkanes are produced much higher than heavy alkanes when the vapour phase appears. We demonstrate that the pressure gradient is also crucial in promoting CO2 and alkane mixing during counterdiffusion processes. These results stress the need for a paradigm shift from classical bulk flow modelling to species-based transport modelling in nanoporous media.

Turbulent flows over rough surfaces can be encountered in a wide range of engineering applications. Despite the progress made after several decades of studies, the prediction of drag and roughness function from the surface geometrical parameters remains an open question. Several methods have shown encouraging results. However, they lack generality and present some scatter in the data. In this paper we propose a new parameter, the effective distribution ($ED$), which lays foundation on the effective slope with some changes to take into account the sheltering effect of large roughness elements and the drag induced by pinnacles higher than the average roughness elements. To develop this new correlation between geometrical features of the wall and the drag, we performed a set of simulations of the turbulent flow over a rough surface made of triangular elements varying their height and spatial distribution. The $ED$ correlates quite well both with the drag and the roughness function for a wide range of cases having different mean roughness height, skewness and kurtosis. To further validate the $ED$, and assessing how it can be generalized to real rough wall, an irregular wall made from the superposition of random sinusoidal function was considered. Results were consistent with the correlation here presented.

Flow control of a low-aspect-ratio flat-plate heaving wing at an average angle of attack of $10^{\circ }$ by a steady-blowing jet is numerically studied by using a feedback immersed boundary–lattice Boltzmann method. Blowing jets at the leading edge, mid-chord and trailing edge are considered. The wing enjoys the highest lift production with the trailing-edge downstream blowing jet, which improves the average lift by 50.0 % at $Re = 1000$ and 22.9 % at $Re = 5000$ through the enhancement of the tip vortex circulation caused by the increase in the mass flux of the shear layer at the wing tips. This increase in mass flux decreases as $Re$ increases from 1000 to 5000 due to its self-limiting mechanism. A mid-chord vertical blowing jet induces a middle vortex which enhances the lift production but the enhancement is smaller than that of trailing-edge downstream blowing jet. Other jet arrangements do not significantly increase the lift coefficient, but the mid-chord upstream blowing jet experiences a significant reduction in the drag coefficient, leading to an increase of 50.6 % in the average lift-to-drag ratio. The effectiveness of the flow control is not significantly affected by the aspect ratio.