We derive from first principles analytic relations for the second- and third-order moments of $\boldsymbol{\mathsf{m}}$, the spatial gradient of fluid velocity $\boldsymbol{u}$, $\boldsymbol{\mathsf{m}} = \nabla \boldsymbol{u}$, in compressible turbulence, which generalize known relations in incompressible flows. These relations, although derived for homogeneous flows, hold approximately for a mixing layer. We also discuss how to apply these relations to determine all the second- and third-order moments of the velocity gradient experimentally for isotropic compressible turbulence.

A nonlinear ensemble-variational data assimilation is performed in order to estimate the unknown flow field over a slender cone at Mach 6, from isolated wall-pressure measurements. The cost functional accounts for discrepancies in wall-pressure spectra and total intensity between the experiment and the prediction using direct numerical simulations, as well as our relative confidence in the measurements and the estimated state. We demonstrate the robustness of the predicted flow by direct propagation of posterior statistics. The approach provides a unique first look at the flow beyond the sensor data, and rigorously accounts for the role of nonlinearity, unlike previous efforts that adopted ad hoc inflow syntheses. Away from the wall, two- and three-dimensional assimilated states both show rope-like structures, qualitatively similar to independent schlieren visualizations. Despite this resemblance, and even though the planar second modes are the most unstable upstream, three-dimensional waves must be included in the assimilation in order to accurately reproduce the wall-pressure measurements recorded in the AFRL Ludwieg Tube facility. The results highlight the importance of three-dimensionality of the field and of the base-state distortion on the instability waves in this experiment, and motivate future measurements that probe the three-dimensional nature of the flow field.

The dynamics of cloud cavitation about a three-dimensional hydrofoil are investigated experimentally in a cavitation tunnel with depleted, sparse and abundant free-stream nuclei populations. A rectangular planform, NACA 0015 hydrofoil was tested at a Reynolds number of $1.4\times 10^{6}$, an incidence of $6^{\circ }$ and a range of cavitation numbers from single-phase flow to supercavitation. High-speed photographs of cavitation shedding phenomena were acquired simultaneously with unsteady force measurement to enable identification of cavity shedding modes corresponding to force spectral peaks. The shedding modes were analysed through spectral decomposition of the high-speed movies, revealing different shedding instabilities according to the nuclei content. With no active nuclei, the fundamental shedding mode occurs at a Strouhal number of 0.28 and is defined by large-scale re-entrant jet formation during the growth phase, but shockwave propagation for the collapse phase of the cycle. Harmonic and subharmonic modes also occur due to local tip shedding. For the abundant case, the fundamental shedding is again large-scale but with a much slower growth phase resulting in a frequency of $St=0.15$. A harmonic mode forms in this case due to the propagation of two shockwaves; an initial slow propagating wave followed by a second faster wave. The passage of the first wave causes partial condensation leading to lower void fraction and consequent increase in the speed of the second wave along with larger-scale condensation. For a sparsely seeded flow, coherent fluctuations are reduced due to intermittent, disperse nuclei activation and cavity breakup resulting in an optimal condition with maximum reduction in unsteady lift.

A set of several exact coherent states in plane Couette flow is computed under spanwise wall oscillation control, with a range of wall oscillation amplitudes and periods $({A_w}, T)$. It is found that the wall oscillation generally stabilises the upper branch of the equilibrium solutions and achieves the corresponding drag reduction, while it influences modestly the lower branch. The stabilisation effect is found to increase with the oscillation amplitude with an optimal time period around ${T^{+}} \approx 100$. The exact coherent states reproduce some key dynamical behaviours of streaks observed in previous studies, while exhibiting the rich coherent structure dynamics that cannot be extracted from a phase average of turbulent states. Visualisation of state portraits shows that the size of the state space supporting turbulent solution is reduced by the spanwise wall oscillation, and the upper-branch equilibrium solutions become less repelling, with many of their unstable manifolds being stabilised. This change of the state space dynamics leads to a significant reduction in lifetime of turbulence. Finally, the main stabilisation mechanism of the exact coherent states is found to be the suppression of the lift-up effect of streaks, explaining why previous linear analyses have been so successful for turbulence stabilisation modelling and the resulting drag reduction.

We present direct numerical simulation results of a wave-current boundary layer in a current-dominated flow regime (wave driven to steady current ratio of 0.34) over bumpy walls for hydraulically smooth flow conditions (wave orbital excursion to roughness ratio of 10). The turbulent, wave-current channel flow has a friction Reynolds number of $350$ and a wave Reynolds number of $351$. At the lower boundary, a bumpy wall is introduced with a direct forcing immersed boundary method, while the top wall has a free-slip boundary condition. Despite the hydraulically smooth nature of the wave-driven flow, the phase variations of the turbulent statistics for the bumpy wall case were found to vary substantially when compared with the flat wall case. Results show that the addition of weak waves to a steady current over flat walls has a negligible effect on the turbulence or bottom drag. However, the addition of weak waves to a steady current over bumpy walls has a significant effect through enhancement of the Reynolds stress (RS) accompanied by a drag coefficient increase of $11\,\%$ relative to the steady current case. This enhancement occurs just below the top of the roughness elements during the acceleration portion of the wave cycle: Turbulent kinetic energy (TKE) is subsequently transported above the roughness elements to a maximum height of roughly twice the turbulent Stokes length. We analyse the TKE and RS budgets to understand the mechanisms behind the alterations in the turbulence properties due to the bumpy wall. The results provide a mechanistic picture of the differences between bumpy and flat walls in wave-current turbulent boundary layers and illustrate the importance of bumpy features even in weakly energetic wave conditions.

This paper investigates flow past a backward-facing step (BFS) in a duct at Reynolds number $Re = 5080$ based on step height, mean inflow velocity and fluid kinematic viscosity. The flow configuration matches a combustion experiment conducted by Pitz and Daily in 1983. High-resolution velocity fields are obtained in OpenFOAM by direct numerical simulation, and the flow field analysed by Lagrangian approaches. Trajectories of fluid particles in areas of interest are obtained by high-order numerical integration, and used to compute finite-time Lyapunov exponents (FTLEs) and polar rotation angles. Lagrangian coherent structures (LCSs) are extracted using geodesic theory, including hyperbolic LCSs and elliptic LCSs. We use complementary qualitative and quantitative LCS analyses to uncover the underlying flow structures. Notably, we find that a flow pathway in which fluid particles rarely diverge from adjacent particles is opened and closed by FTLE ridges determined by the periodic shedding of vortices from the BFS. Two dominant vortices with significant Lagrangian coherence, generated respectively by the separated boundary layer and shear layer, are self-sustaining and of comparable strength. Hyperbolic repelling LCSs act as transport barriers between the pathway and cores of the coherent vortices, thus playing a major part in the fluid entrainment process. Interactions between these different geometric regions partitioned by LCSs lead to intrinsic complexity in the BFS flow.

This paper describes a numerical investigation on the effects of small to large bluntness as well as the angle of incidence on the leading-edge separation in a laminar hypersonic flow. The results show that both bluntness and incidence strongly influence separation characteristics. A combined bluntness and incidence parameter is proposed and applied to plateau pressures in the separated region and is shown to demarcate effects of small and large bluntness at different angles of incidence. The investigation also showed profound changes in the internal structure of the separation bubble with changes in incidence and bluntness. The significance of the relationship between the two shear stress minima and the separation length is highlighted.

Based on experimental data acquired with particle image velocimetry, we examine turbulent boundary layers that are subjected to an abrupt change in wall roughness in the streamwise direction. Three different sandpapers (P24, P36 and P60) together with a smooth wall are used to form a number of different surface transition cases, including both R $\rightarrow$ S (where upstream surface is rough and second surface is either smooth or smoother compared with the upstream surface) and S $\rightarrow$ R (where upstream surface is smoother compared with the downstream surface; both surfaces are rough). This enables us to investigate the effect of the surface transition strength ($M = \ln [y_{02}/y_{01}]$, where $y_{01}$ and $y_{02}$ are the roughness lengths of the upstream and downstream surfaces, respectively) on the growth of the internal boundary layer (IBL) and the corresponding flow structure. The results show that for each surface transition group (i.e. R $\rightarrow$ S and S $\rightarrow$ R), the thickness of the IBLs is proportional to the strength of the surface transition, and that the IBLs are thicker for the S $\rightarrow$ R cases compared with their R $\rightarrow$ S counterparts for similar $|M|$, when normalised by the initial boundary layer thickness ($\delta _0$). The results also show that the growth rates of the IBLs could be represented by a power law, consistent with the previous studies. However, despite a wide range of scatter in the literature for the power-law exponent, an average value of $0.75$ (varies between $0.71$ and $0.8$ with no clear trend) is obtained in the present study considering all the surface transition cases. The pre-factor for the power-law fit, on the other hand, is found to be related to the strength of the surface transition. In addition to the variations in the velocity defect and diagnostic plots with the growth of the IBLs, sweep and ejection events appear to differ significantly (depending on the type of the step change). Two-point spatial correlations, moreover, show that the structure is more elongated in the wall-normal and streamwise directions, as the flow accelerates over the downstream surface (i.e. R $\rightarrow$ S cases). For the reverse transition cases (i.e. S $\rightarrow$ R, where the flow decelerates over the downstream rougher surface), however, the correlation coefficients shrink in size in both directions.

For one-way coupled (1WC) flows, Tom & Bragg (J. Fluid Mech., vol. 871, 2019, pp. 244–270) advanced the analysis of Maxey (J. Fluid Mech., vol. 174, 1987, pp. 441–465), which applied to weakly inertial particles, to particles of arbitrary inertia, and the new theoretical result revealed the role that different scales play in the preferential sweeping mechanism that leads to enhanced particle settling in turbulent flows. Monchaux & Dejoan (Phys. Rev. Fluids, vol. 2, 2017, 104302) showed using direct numerical simulations that, while for low particle loading the effect of two-way coupling (2WC) on the global flow statistics is weak, 2WC enables the particles to drag the fluid in their vicinity down with them, significantly enhancing their settling, and they argued that 2WC suppresses the preferential sweeping mechanism. We explore this further by considering the impact of 2WC on the contribution made by eddies of different sizes on the particle settling. In agreement with Monchaux & Dejoan, we show that even for low loading, 2WC strongly enhances particle settling, and we show how 2WC modifies the contribution from different flow scales. However, contrary to their study, we show that preferential sweeping remains important in 2WC flows. In particular, for both 1WC and 2WC flows, the settling enhancement due to turbulence is dominated by contributions from particles in straining regions of the flow, but for the 2WC case, the particles in these regions also drag the fluid down with them, leading to an enhancement of their settling compared with the 1WC case.

The impact of Strouhal number St (= 0.1–1.0), Reynolds number Re (= 50–2000) and dimensionless wavelength λ (= 0.5–2.0) on the hydrodynamic performance of a travelling wavy foil of a constant length is extensively investigated. The relationship of time-mean thrust with St, Re and λ is presented, suggesting that the propulsive force increases with increasing St, Re and λ. As such, the drag–thrust boundary advances as these parameters increase. A shorter λ makes the thrust steadier while a longer λ enhances the maximum instantaneous thrust. The latter is beneficial for prey to escape from a predator. The fluid added mass caused by the foil oscillation increases with St and λ but declines with Re (<500). Seven types of wake structures produced by the foil are identified, discussed and connected to thrust generation, showing how St, Re and λ affect the fluid dynamics, wake transition, vortex strength, wake jet, velocity, added mass, added damping, power input, efficiency and pressure profiles. The outcome of this work renders a physical basis for understanding the swimming of aquatic animals.

A mathematical model for the simulation of spherical bubble dynamics in binary alkane mixtures is presented. Detailed heat and mass transfer and phase transition are resolved, and air absorption and desorption are considered. As example mixtures, high-volatile heptane and low-volatile dodecane are investigated. The low-volatile component shows a convective counter-intuitive mass transport opposite to the diffusion flux. A staggered mass flux of heptane and dodecane is associated with a local segregation of mixture components in the surrounding liquid and a distinctive species distribution within the bubble. A comparative study with a pseudo-fluid where alkane species interdiffusion is absent and whose components cannot segregate reveals that for an oscillating bubble, the local segregation hardly affects bubble dynamics, while for a continuously growing bubble in a superheated liquid, the growth rate is considerably affected, particularly for a low heptane mixture percentage. This study demonstrates limitations of treating mixtures by a single-component surrogate fluid and may serve as a starting point for the development of multi-component cavitation models for computational fluid dynamics applications in real fluid mixtures as, e.g. fuels or hydraulic oils.

A fluid-filled cube rotating about an axis passing through the midpoints of opposite edges is subjected to small librations (i.e. modulation of the mean rotation). Low viscosity regimes, with Ekman number as small as $10^{-8}$ and equally small relative forcing amplitude, are explored numerically. The full inertial range of forcing frequencies, from 0 to twice the mean rotation rate are considered. The response flows are dominated by inertial wavebeams emitted from edges and/or vertices, depending on the forcing frequency. How these reflect on the cube's walls and focus onto edges and vertices lead to intricate patterns. Most of the results can be reconciled using linear inviscid ray-tracing theory with careful attention to wavebeam emissions and reflections. However, even at the low Ekman number and relative forcing amplitude considered, other effects are discernible which are not captured by ray tracing. These include a symmetry breaking due to viscous effects and a progressive wave, retrograde to the mean rotation and localized in the boundary layers of the cube, due to nonlinear effects and the librational forcing.

Thermo-responsive hydrogels are a promising material for creating controllable actuators for use in micro-scale devices, since they expand and contract significantly (absorbing or expelling fluid) in response to relatively small temperature changes. Understanding such systems can be difficult because of the spatially and temporally varying properties of the gel, and the complex relationships between the fluid dynamics, elastic deformation of the gel and chemical interaction between the polymer and fluid. We address this using a poro-elastic model, considering the dynamics of a thermo-responsive spherical hydrogel after a sudden change in the temperature that should result in substantial swelling or shrinking. We focus on two model examples, with equilibrium parameters extracted from data in the literature. We find a range of qualitatively different behaviours when swelling and shrinking, including cases where swelling and shrinking happen smoothly from the edge, and other situations that result in the formation of an inwards-travelling spherical front that separates a core and shell with markedly different degrees of swelling. We then characterise when each of these scenarios is expected to occur. An approximate analytical form for the front dynamics is developed, with two levels of constant porosity, that well-approximates the numerical solutions. This system can be evolved forward in time, and is simpler to solve than the full numerics, allowing for more efficient predictions to be made, such as when deciding dosing strategies for drug-laden hydrogels.

In this work, we investigate the problem of imbibition/drainage of a fluid in capillaries of arbitrary axisymmetric cross-sections filled initially with another immiscible one. The model predicts the location of the meniscus and its speed along the tube length with time. The two immiscible fluids may assume any density and viscosity contrasts. In addition, the axisymmetric profile of the tube maintains a relatively small angle of tangency to warrant that the axial velocity distribution assumes, approximately, a parabolic profile. The driving forces that may be encountered in this system include the capillary force, pressure force, gravitational force and an opposing viscous force. The orientation of the capillary force can be in the direction of the flow (e.g. during imbibition) or opposite to the flow (e.g. during drainage). Likewise, the gravitational force can be in the direction of the flow or opposite to it. In this work we account for all these possibilities. A differential equation is developed that defines the location of the meniscus with time. A fourth-order-accurate Runge–Kutta scheme has been developed to provide solutions for the different scenarios associated with this system. It is shown that the developed model reduces to those appropriate for straight tubes, which builds confidence in the modelling approach. The effects of changing the tangent along the profile of the tube, which influences the calculation of the radius of curvature of the meniscus, is also considered. Unlike the cases of straight capillary tubes, in tubes with arbitrary symmetric profiles, the friction force depends on the variations of the tube profile. Examples of converging/diverging capillary tubes that follow straight and power law profiles are investigated. In addition, the case of sinusoidal profiles has also been considered.

We investigate the dynamics of a fluid layer subject to a bottom heat flux and a top monotonically increasing temperature profile driving horizontal convection (HC). We use direct numerical simulations and consider a large range of flux-based Rayleigh numbers $10^6 \leq Ra_F \leq 10^9$ and imposed top horizontal to bottom vertical heat flux ratios $0 \leq \varLambda \leq 1$. The fluid domain is a closed two-dimensional box with aspect ratio $4\leq \varGamma \leq 16$ and we consider no-slip boundaries and adiabatic side walls. We demonstrate a regime transition from Rayleigh–Bénard (RB) convection to HC at $\varLambda \approx 10^{-2}$, which is independent of $Ra_F$ and $\varGamma$. At small $\varLambda$, the flow is organised in multiple overturning cells with approximately unit aspect ratio, whereas at large $\varLambda$ a single cell is obtained. The RB-relevant Nusselt number scaling with $Ra_F$ and the HC-relevant Nusselt number scaling with the horizontal Rayleigh number $Ra_L=Ra_F\varLambda \varGamma ^4$ are in good agreement with previous results from classical RB convection and HC studies in the limit $\varLambda \ll 10^{-2}$ and $\varLambda \gg 10^{-2}$, respectively. We demonstrate that the system is multi-stable near the transition $\varLambda \approx 10^{-2}$, i.e. the exact number of cells not only depends on $\varLambda$ but also on the system's history. Our results suggest that subglacial lakes, which motivated this study, are likely to be dominated by RB convection, unless the slope of the ice–water interface, which controls the horizontal temperature gradient via the pressure-dependence of the freezing point, is greater than unity.

We investigate the leading-edge separation, also called negligible boundary-layer thickness separation, induced by an impinging shock on a sharp flat plate. The canonical impinging shock wave/boundary-layer interaction configuration consisting of a wedge and a plate wall, placed in hypersonic free stream, is used for the investigations. We first construct a theoretical model for the leading-edge separated flow field (LESF) which predicts separation bubble geometry and surface pressure distribution as a function of three parameters: free-stream Mach number, wedge angle and the reattached flow turning angle. Markedly different predictions of the separated flow field are obtained for an oblique and a near normal reattachment on the plate surface. Experiments in a shock tunnel at a nominal Mach $6$ flow, with an impinging shock generated by a wedge of angle $26.6^{\circ }$, are used to validate the model. Schlieren flow visualization using a high-speed camera and surface pressure measurements using fast response piezoelectric sensors are the diagnostics employed. For a range of shock impingement locations, the LESFs are observed to be geometrically similar and in good agreement with the LESF model. When the impingement location gets closer to the leading edge, it is observed from experiments that the flow field is no longer geometrically similar, and the separation angle increases as the impingement gets closer to the leading edge beyond the range of similarity. The work thereby offers an elaborate description of the leading-edge separated flow when shock impingement occurs near the plate leading edge.

We report an experimental study of the viscous boundary layer (BL) properties of turbulent Rayleigh–Bénard convection in a cylindrical cell. The velocity profile with all three components was measured from the centre of the bottom plate by an integrated home-made particle image velocimetry system. The Rayleigh number $Ra$ varied in the range $1.82 \times 10^8 \le Ra \le 5.26 \times 10^9$ and the Prandtl number $Pr$ was fixed at $Pr = 4.34$. The probability density function of the wall-shear stress indicates that using the velocity component in the mean large-scale circulation (LSC) plane alone may not be sufficient to characterise the viscous BL. Based on a dynamic wall-shear frame, we propose a method to reconstruct the measured full velocity profile which eliminates the effects of complex dynamics of the LSC. Various BL properties including the eddy viscosity are then obtained and analysed. It is found that, in the dynamic wall-shear frame, the eddy viscosity profiles along the centre line of the convection cell at different $Ra$ all collapse on a single master curve described by $\nu _t^d / \nu = 0.81 (z / \delta _u^d) ^{3.10 \pm 0.05}$. The Rayleigh number dependencies of several BL quantities are also determined in the dynamic frame, including the BL thickness $\delta _u^d$ (${\sim } Ra^{-0.21}$), the Reynolds number $Re^d$ (${\sim }Ra^{-0.46}$) and the shear Reynolds number $Re_s^d$ (${\sim } Ra^{0.24}$). Within the experimental uncertainty, these scaling exponents are the same as those obtained in the static laboratory frame. Finally, with the measured full velocity profile, we obtain the energy dissipation rate at the centre of the bottom plate $\varepsilon _{w}$, which is found to follow $\langle \varepsilon _{w} \rangle _t \sim Ra^{1.25}$.

Adhesive contact between a thin elastic sheet and a substrate arises in a range of biological, physical and technological applications. By considering the dynamics of this process that naturally couples fluid flow, long-wavelength elastic deformations and microscopic adhesion, we analyse a sixth-order thin-film equation for the short-time dynamics of the onset of adhesion and the long-time dynamics of a steadily propagating adhesion front. Numerical solutions corroborate scaling laws and asymptotic analyses for the characteristic waiting time for adhesive contact and for the speed of the adhesion front. A similarity analysis of the governing partial differential equation further allows us to determine the shape of a fluid-filled blister ahead of the adhesion front. Finally, our analysis reveals a near-singular behaviour at the moving elastohydrodynamic contact line with an effective boundary condition that might be useful in other related problems.

This paper deals with the construction of the tip asymptotes for a hydraulic fracture deflating in a permeable elastic medium. Specifically, the paper describes the changing nature of the asymptotic fields during the arrest and recession phases following propagation of the fracture after fluid injection has ended. It shows that as the fracture deflates in the arrest phase, the region of dominance of the linear elastic fracture mechanics tip asymptote $w\sim x^{1/2}$ of the fracture aperture $w$ with distance $x$ from the front shrinks to the benefit of an intermediate asymptote $w\sim x^{3/4}$. Hence only the velocity-independent $3/4$ asymptote is left at the arrest–recession transition. Furthermore, with increasing receding velocity of the front, a linear asymptote $w\sim x$ develops progressively at the fracture tip, with $w\sim x^{3/4}$ again becoming an intermediate asymptote. These universal multiscale asymptotes for the arrest and recession phases are key to determining, in combination with a computational algorithm that can simulate the evolution of a finite fracture, the decaying stress intensity factor during arrest, the time at which the fracture transitions from arrest to recession, and the negative front velocity during recession.

The problem of propagation of steady periodic waves over a planar seabed is faced through the definition of a suitable semi-analytical iterative scheme. The latter is capable of describing highly nonlinear waves in deep, intermediate and shallow water conditions. Comparisons with the existing fifth-order theories show that the proposed model is accurate in all the regimes of motion and that it does not present any of the limitations affecting the Stokes and cnoidal wave solutions. Further, it also provides a reliable approximation of the dynamics of maximum-amplitude waves. The definition of the iterative scheme is preceded by a detailed study of the geometrical wave parameters. This latter part is aimed at the definition of a global scaling for water waves valid in all the regimes of motion.