Compound surfaces, consisting of periodic arrays of solid patches and free surfaces, exhibit hydrodynamic slipperiness which is quantified by their slip length. The limit of small solid fractions, where the slip length diverges, is of fundamental interest. This paper addresses longitudinal liquid flows over a periodic array of grooves which are partially invaded by the liquid. Assuming that the slats separating the grooves are infinitely thin, the solid fraction $\epsilon$ is set by the invasion depth. Inspired by the singular small-solid-fraction limit for non-invaded grooves (Schnitzer, Phys. Rev. Fluids, vol. 1, 2016, 052101R), we consider the idealised geometry of $90^{\circ }$ protrusion angles, where an integral force balance in the limit $\epsilon \to 0$ implies a slip length that scales as $\epsilon ^{-1}$. The problem exhibits a nested structure, where the liquid domain is conceptually decomposed into four distinct regions: a unit-cell region on the scale of the period, where the wetted portion of the slat appears as a point singularity; two regions on the scale of the wetted slat, where the flow essentially varies in one dimension; and a transition region about the tip of the slat. Analysing these regions using matched asymptotic expansions and conformal mappings yields the ratio of slip length to semi-period as $2\epsilon ^{-1} - (10/{\rm \pi} )\ln 2 + \cdots$.

Cross-waves are standing waves with crests perpendicular to a wave-maker; they are subharmonic waves excited by parametric instability. The modulational and chaotic behaviours of nonlinear cross-waves have been studied widely since the 1970s. Most of the previous work has focused on gravity waves where surface tension can be neglected. In this work we study cross-waves that are highly dependent on surface tension as well as gravity. By oscillating a planar wave-maker either vertically or horizontally with frequencies of 25 Hz through 40 Hz at one end of a rectangular basin, two-dimensional multi-component surface patterns are realized. Using the free-surface synthetic Schlieren technique to measure the surface elevations, multi-dimensional Fourier transforms are utilized to track the evolutionary spectrum of the water surface in both the temporal and spatial domains. Wavelet transforms are implemented to show the development of the various frequency components. Three-wave resonances with and without first subharmonics are observed for small nonlinearity. Three-dimensional oblique propagating cross-waves are generated at higher nonlinearity; unlike most previous cross-wave experiments, this staggered pattern propagates far downstream. Experimental evidence shows that two oblique propagating waves form a two-dimensional short-crested pattern, and that the lateral component of the waves develops into parametric sloshing modes corresponding to the width of the tank. Two regimes of nonlinear wave patterns, resonant triads and oblique propagating cross-waves, are delineated.

The role of the different helical components of the magnetic and velocity fields in the inverse spectral transfer of magnetic helicity is investigated through Fourier shell-to-shell transfer analysis. Magnetic helicity transfer analysis is performed on chosen data from direct numerical simulations of homogeneous isothermal compressible magnetohydrodynamic turbulence, subject to both a large-scale mechanical forcing and a small-scale helical electromotive driving. The root mean square Mach number of the hydrodynamic turbulent steady state taken as initial condition varies from 0.1 to about 11. Three physical phenomena can be distinguished in the general picture of the spectral transfer of magnetic helicity towards larger spatial scales: local inverse transfer (LIT), non-local inverse transfer (NLIT) and local direct transfer (LDT). A shell decomposition allows these three phenomena to be associated with clearly distinct velocity scales: the LDT is driven by large-scale velocity shear and associated with a direct magnetic energy cascade; the NLIT is mediated by small-scale velocity fluctuations which couple small- and large-scale magnetic structures; and the LIT by the intermediate spatial scales of the velocity field. The helical decomposition shows that like-signed helical interactions and interactions with the compressive velocity field are predominant. The latter has a high impact on the LDT and on the NLIT, but plays no role for the LIT. The locality and relative strength of the different helical contributions are mainly determined by the triad helical geometric factor, derived here in the compressible case.

Internal wave generation is fundamentally the conversion of barotropic to baroclinic energy that often occurs due to vertical acceleration of stratified flows over topographic features. Acceleration results in a phase lag between density (pressure) perturbations and the barotropic velocity. To estimate the conversion of barotropic to baroclinic energy, the density perturbation is often calculated using a time-invariant background density. Other phenomena, however, can also alter the phasing of density perturbations and vertical velocities, such as barotropic tidal heaving and internal wave interactions. Consequently, accurately accounting for these dynamics in energy budgets is important. Tidal averaging or modal decomposition are often used to isolate topographic energy conversion in the presence of these other phenomena. However, while effective, these methods do not provide insights into the dynamics of conversion either through time or over depth. Here, we present a new analytical approach to calculating barotropic to baroclinic conversion using a time-varying background density. Our method results in an additional term in the baroclinic energy budget that directly accounts for barotropic tidal heaving and internal wave interactions, depending on the formulation of the background density. The tidally averaged, domain-integrated conversion rate is consistent across methods. Isolation of topographic conversion demonstrates that conversion due to interactions between internal wave beams and barotropic tidal heaving lead to relatively small differences in the overall conversion. However, using a time-varying background density allows for full decomposition of barotropic to baroclinic conversion through time and the identification of regions where negative conversion related to mixing actually occurs.

Through boundary integral simulations and asymptotic analysis, we investigate the effect of a finite Navier slip length on the rheological proprieties of a dilute two-dimensional suspension of plate-like particles in the creeping flow limit. Specifically, we study the effects of Navier slip, particle thickness and Péclet number on the effective shear viscosity and average normal stress difference of an isolated two-dimensional plate-like particle in an unbounded shear flow field. We find that Navier slip induces a significant reduction in the effective viscosity and increases the average normal stress difference. The effect of slip becomes more enhanced as the thickness of the particle decreases and as the Péclet number increases. Remarkably, the analysis suggests that it is theoretically possible for a dilute suspension of slip plate-like particles at high Péclet numbers to have a shear viscosity smaller than that of the suspending fluid.

We report the results of a theoretical investigation of the stability of a hydrodynamic analogue of Landau levels, specifically circular orbits arising when a millimetric droplet self-propels along the surface of a vibrating, rotating liquid bath. Our study elucidates the form of the stability diagram characterising the critical memory at which circular orbits destabilise, and the form of instability. Particular attention is given to rationalising observations reported in prior experimental works, including the prevalence of resonant wobbling instabilities, in which the instability frequency is approximately twice the orbital frequency. We also explore the physical mechanism responsible for the onset of instability. Specifically, we compare the efficacy of different heuristic arguments proposed in prior studies, including propositions that the most unstable orbits arise when their radii correspond to the zeros of Bessel functions or when their associated wave intensity is extremised. We establish a new relation between orbital stability and the mean wave field, which supersedes existing heuristic arguments and suggests a rationale for the alternate wobbling and monotonic instabilities arising at onset as the orbital radius is increased progressively.

The spatial relationship between turbulent and particle concentration structures is investigated based on the turbulent velocity and particle concentration data obtained synchronously at the Qingtu Lake Observation Array site. In addition to the observation of particle concentration structures that contain not only large-scale coherence but also significant energy in the high-Reynolds-number atmospheric surface layer (ASL), the scale of turbulent motions that have the most significant coherence with particle concentration is found to follow a 1/2 power scaling law with the local height and ASL thickness. Moreover, large-scale turbulent velocity fluctuations have a significant amplitude modulation effect on particle concentration fluctuations, but the modulating influence is different for small dust particles and large saltating particles. Based on the interphase amplitude modulation, there exists a particle–turbulence structure phase difference that varies with height, which further makes the structure inclination angle of the particle concentration larger than that of the turbulence. In this scenario, a conceptual model reflecting the relationship between the two is proposed, and a quantitative formulation is further derived and found to be in good agreement with the experimental results. These findings and the proposed model contribute insights into particle–turbulence interactions, thereby providing theoretical support for a unified model of turbulence dynamics and particle kinematics.

The near-field characteristics of highly buoyant plumes, commonly referred to as lazy plumes, remain relatively poorly understood across a range of flow conditions, particularly compared with our understanding of far-field characteristics. Here, we perform fully resolved three-dimensional numerical simulations of round helium plumes to characterize the effects of different inlet Richardson, ${Ri}_0$, and Reynolds, ${Re}_0$, numbers on first- and second-order statistical moments as well as average vertical fluxes in the near field. For sufficiently high ${Re}_0$ at a particular ${Ri}_0$, heavy air can penetrate the core of the plume, reminiscent of spikes in the classical Rayleigh–Taylor instability. In the most turbulent simulation, this penetration becomes so strong that a recirculation zone forms along the centreline of the plume. Vertical fluxes are found to scale linearly with vertical distance from the plume inlet, consistent with experimental and numerical observations (Jiang & Luo, Flow Turbul. Combust., vol. 64, 2000, pp. 43–69; Kaye & Hunt, Intl J. Heat Fluid Flow, vol. 30, 2009, pp. 1099–1105). We analytically derive this linear scaling from the governing equations by making a radial entrainment hypothesis whereby ambient fluid is entrained, on average, only in the radial direction at a finite distance from the inlet. Through this derivation, we identify physical mechanisms that can cause these relationships to remain only approximately valid for the present simulations. Lastly, we identify near-field power-law scaling relations for the flux magnitudes based on ${Ri}_0$, and also examine vertical profiles of the non-dimensional Richardson number flux. Ultimately, insights from the present simulations are used to define near-, intermediate- and far-field regions in buoyant plumes.

Existing theoretical analyses of Faraday waves in Hele-Shaw cells rely on the Darcy approximation and assume a parabolic flow profile in the narrow direction. However, Darcy's model is known to be inaccurate when convective or unsteady inertial effects are important. In this work, we propose a gap-averaged Floquet theory accounting for inertial effects induced by the unsteady terms in the Navier–Stokes equations, a scenario that corresponds to a pulsatile flow where the fluid motion reduces to a two-dimensional oscillating Poiseuille flow, similarly to the Womersley flow in arteries. When gap-averaging the linearised Navier–Stokes equation, this results in a modified damping coefficient, which is a function of the ratio between the Stokes boundary layer thickness and the cell's gap, and whose complex value depends on the frequency of the wave response specific to each unstable parametric region. We first revisit the standard case of horizontally infinite rectangular Hele-Shaw cells by also accounting for a dynamic contact angle model. A comparison with existing experiments shows the predictive improvement brought by the present theory and points out how the standard gap-averaged model often underestimates the Faraday threshold. The analysis is then extended to the less conventional case of thin annuli. A series of dedicated experiments for this configuration highlights how Darcy's thin-gap approximation overlooks a frequency detuning that is essential to correctly predict the locations of the Faraday tongues in the frequency–amplitude parameter plane. These findings are well rationalised and captured by the present model.

Subgrid-scale (SGS) modelling is formulated using a local transport of spectral kinetic energy estimated by a wavelet multiresolution analysis. Using a spectrally and spatially local decomposition by wavelet, the unresolved inter-scale energy transfer and modelled SGS dissipation are evaluated to enforce explicitly and optimally their balance a priori over a range of large-eddy simulation (LES) filter widths. The formulation determines SGS model constants that optimally describe the spectral energy balance between the resolved and unresolved scales at a given cutoff scale. The formulation is tested for incompressible homogeneous isotropic turbulence (HIT). One-parameter Smagorinsky- and Vreman-type eddy-viscosity closures are optimised for their model constants. The algorithm discovers the theoretical prediction of Lilly (The representation of small-scale turbulence in numerical simulation experiments. In Proceedings of the IBM Scientific Computing Symposium on Environmental Sciences, pp. 195–210) at a filter cutoff scale in the inertial subrange, whereas the discovered constants deviate from the theoretical value at other cutoff scales so that the spectral optimum is achieved. The dynamic Smagorinsky model used a posteriori shows a suboptimal behaviour at filter scales larger than those in the inertial subrange. A two-parameter Clark-type closure model is optimised. The optimised constants provide evidence that the nonlinear gradient model of Clark et al. (J. Fluid Mech., vol. 91, issue 1, 1979, pp. 1–16) is prone to numerical instability due to its model form, and combining the pure gradient model with a dissipative model such as the classic Smagorinsky model enhances numerical stability but the standard mixed model is not optimal in terms of spectral energy transfer. A posteriori analysis shows that the optimised SGS models produce accurate LES results.

We perform direct numerical simulations (DNS) of the fully developed turbulent channel flow with a passive scalar subjected to constant time-averaged scalar fluxes at the wall. The first case (case I) considers a constant time-averaged scalar flux $\left \langle q_w \right \rangle$ along with the specific condition of a non-fluctuating scalar $\theta _w$ imposed at the wall implying zero wall scalar fluctuation $\theta '_w=0$ and zero variance (isoscalar boundary condition) whereas the second case (case II) accounts for a constant instantaneous scalar flux $q_w$ in time and space leading to zero gradient of scalar fluctuation along the normal to the wall $(\partial \theta '/ \partial x_n)_w = 0$ (isoflux boundary condition) implying a non-zero variance. The friction Reynolds number takes on the value $R_\tau =395$ and the molecular Prandtl numbers considered are $P_r=0.01$, 0.1, 1 and 10. The purpose is to investigate the effect of the wall scalar fluctuations on the scalar field. Emphasis is put on the mean passive scalar $\left \langle \theta \right \rangle$, the half-scalar variance $k_\theta$, the turbulent scalar fluxes $\tau _{i \theta }$, the correlation coefficients $R_{i \theta }$, the passive scalar to dynamic time-scale ratio $\mathcal {R}$, the turbulent Prandtl number $Pr_t$ and higher-order scalar statistics. Systematic comparisons between these two scalar fields are undertaken. As a result of interest, it is found that the mean scalar remains almost the same whatever the type of the boundary layer condition but not the scalar variance. The budgets of transport equations for the half-scalar variance and turbulent fluxes reveal that some scalar quantities such as the dissipation-rate and the molecular diffusion terms are highly modified in the near-wall region but not really the production, diffusion and the scalar-pressure gradient correlation terms in a first approximation. Visualization of the instantaneous scalar fields indicates that the topology of the structures is strongly modified as the Prandtl number increases from $P_r=0.01$ up to $10$. Finally, it is shown how to use this DNS database to devise and calibrate scalar flux equation models in the framework of both first and second moment closures. This study suggests that accounting for wall scalar fluctuations should be considered in the simulation of turbulent flows involving fluid and solid combinations at the interface and provides a useful high resolution DNS database.

Understanding particle motion in narrow channels can guide progress in numerous applications, from filtration to vascular transport. Thermal or active fluctuations of fluid-filled channel walls can slow down or increase the dispersion of tracer particles via entropic trapping in the wall bulges or hydrodynamic flows induced by wall fluctuations, respectively. Previous studies concentrated primarily on the case of a single Brownian tracer. Here, we address what happens when there is a large ensemble of interacting Brownian tracers – a common situation in applications. Introducing repulsive interactions between tracer particles, while ignoring the presence of a background fluid, leads to an effective flow field. This flow field enhances tracer dispersion, a phenomenon reminiscent of that seen for single tracers in incompressible background fluid. We characterise the dispersion by the long-time diffusion coefficient of tracers numerically and analytically with a mean-field density functional analysis. We find a surprising effect where an increased particle density enhances the diffusion coefficient, challenging the notion that crowding effects tend to reduce diffusion. Here, inter-particle interactions push particles closer to the fluctuating channel walls. Interactions between the fluctuating wall and the now-nearby particles then drive particle mixing. Our mechanism is sufficiently general that we expect it to apply to various systems. In addition, our perturbation theory quantifies dispersion in generic advection–diffusion systems.

In two-dimensional decaying homogeneous isotropic turbulence, kinetic energy and enstrophy are respectively transferred to larger and smaller scales. In such spatiotemporally complex dynamics, it is challenging to identify the important flow structures that govern this behaviour. We propose and employ numerically two flow-modification strategies that leverage the inviscid global conservation of energy and enstrophy to design external forcing inputs that change these quantities selectively and simultaneously, and drive the system towards steady-state or other late-stage behaviour. One strategy employs only local flow field information, while the other is global. We observe various flow structures excited by these inputs and compare them with recent literature. Energy modification is characterized by the excitation of smaller wavenumber structures in the flow than enstrophy modification.

Explosive dispersal of granular media widely occurs in nature across various length scales, also enabling engineering applications ranging from commercial or military explosive systems to the loss prevention industry. However, the complex particle–flow coupling makes the explosive dispersal behaviour of particles difficult to control or even characterize. Here, we study the central explosion-driven dispersal of dense particle layers using the coarse-grained computational fluid dynamics–discrete element method and present a comprehensive investigation of both macroscale dispersal behaviours and particle-scale pattern formation. Employing three independent dimensionless parameters that characterize the efficiency, homogeneity and completeness of explosive dispersal, we categorize the dispersal behaviours into ideal, partial, retarded and failed modes, and propose the corresponding thresholds. As the mass ratio of granular materials to central pressurized gases (M/C) spans four orders of magnitude, the dispersal mode transitions from ideal to partial, then to retarded and finally to failed mode. The transitions of dispersal modes correspond to the particle–flow coupling regime crossovers, which change from decoupling to weak, medium and finally to strong coupling as the dispersal mode undergoes corresponding transitions. We proceed to develop continuum models accounting for the shock compaction and the ensuing pulsation of the particle ring that are capable of identifying the ideal dispersal mode from various dispersal systems. We also provide insights into the origins of diverse particle-scale patterns that are strongly correlated with macroscale dispersal modes and critical for the accurate prediction of dispersal modes.

We present the analytical solution for the two-dimensional velocity and density fields within an approximation for laminar stratified inclined duct (SID) flows where diffusion dominates over inertia in the along-channel momentum equation but is negligible in the density transport equation. We refer to this approximation as the hydrostatic/gravitational/viscous in momentum and advective in density (HGV-A) approximation due to the leading balances in the governing equations. The analytical solution is valid for laminar flows in a two-layer configuration in the limit of long ducts. The non-dimensional volume flux within the HGV-A approximation is given by $Fr^* ={{Re}}_g/(AK)$, which is a control parameter with ${{Re}}_g$ the gravitational Reynolds number, $A$ the aspect ratio of the duct and $K$ a geometrical parameter that depends on the tilt of the duct and is obtained from the analytical solution. This analytical solution was validated against results from laboratory experiments, and allows us to gain new insight into the dynamics and properties of SID flows. Most importantly, constant values of $Fr^*$ describe, in both horizontal and inclined ducts, the transitions between increasingly turbulent flow regimes: from laminar flow, to interfacial waves, to intermittent turbulence and sustained turbulence.

This study explores the long-term impact of a weak initial departure from circular symmetry in coherent equivalent-barotropic vortices on their dynamics and evolution. An algorithm is developed which makes it possible to construct models of vortices that initially propagate with prescribed velocity. These solutions are used as the initial conditions for a series of numerical simulations. Simulations indicate that seemingly minor perturbations of dipolar form can control the propagation of vortices for extended periods, during which they translate over distances greatly exceeding their size. The numerical results are contrasted with the linear model, which assumes that the non-axisymmetric component of circulation in the vortex interior is relatively weak. The linear solutions reflect the self-propagation tendencies of coherent vortices to a much lesser degree, which underscores the role of fundamentally nonlinear mechanisms at play. The remarkable ability of quasi-monopolar vortices to retain the memory of weak initial perturbations helps to rationalize the wide range of the observed propagation velocities of coherent long-lived vortices in the ocean.

We present a theoretical and numerical stability analysis for a piston-driven planar shock against two-dimensional perturbations. The results agree with the well-established theory for isolated planar shocks: in the range of $h_c< h<1+2 {\mathcal {M}}_2$, where $h$ is the D'yakov–Kontorovich (DK) parameter related to the slope of the Rankine–Hugoniot curve, $h_c$ is its critical value corresponding to the onset of the spontaneous acoustic emission (SAE) and ${\mathcal {M}}_2$ is the downstream Mach number, non-decaying oscillations of shock-front ripples occur. The effect of the piston is manifested in the presence of additional frequencies occurring by the reflection of the sonic waves on the piston surface that can reach the shock. An unstable behaviour of the shock perturbation is found to be possible when there is an external excitation source affecting the shock, whose frequency coincides with the self-induced oscillation frequency in the SAE regime, thereby being limited to the range $h_c< h<1+2 {\mathcal {M}}_2$. An unstable evolution of the shock is also observed for planar shocks restricted to one-dimensional perturbations within the range $1< h<1+2 {\mathcal {M}}_2$. Both numerical integration of the Euler equations via the method of characteristics and theoretical analysis via Laplace transform are employed to cross-validate the results.

The adhesion of two surfaces relies on the propagation of an adhesion front. What is the dynamics of the front when both surfaces are coated with a thin layer of viscous liquid? Standard criteria from fingering instabilities would predict a stable front since viscous fluid pushes away air of low viscosity. Surprisingly, the front propagation may be unstable and generally leads to growing fingers. We demonstrate with model experiments where the two adhering surfaces are slightly tilted by an angle $\alpha$ that the origin of this interfacial instability relies on feeding the front from the surrounding thin film. We show experimentally that the typical wavelength of the instability is mainly dictated by the thickness of the oil layers $h$. In this wedge geometry, the propagation dynamics is found to follow a $t^{1/2}$ dependence and to saturate for an extension length of the order of $h/\alpha$.

The power exchange between fluid and structure plays a significant role in the force production and flight efficiency of flapping wings in insects and artificial flyers. This work numerically investigates the performance of flapping wings by using a high-fidelity fluid–structure interaction solver. Simulations are conducted by varying the hinge flexibility (measured by the Cauchy number, $Ch$, i.e. the ratio between aerodynamic and torsional elastic forces) and the wing shape (quantified by the radius of the first moment of area, $\bar {r}_1$). Results show that the lift production is optimal at $0.05 < Ch \leq 0.2$ and larger $\bar {r}_1$ where the minimum angle of attack is around $45^\circ$ at midstroke. The power economy is maximised for wings with lower $\bar {r}_1$ near $Ch=0.2$. Power analysis indicates that the optimal performance measured by the power economy is obtained for those cases where two important power synchronisations occur: anti-synchronisation of the pitching elastic power and the pitching aerodynamic and inertial powers and nearly in-phase synchronisation of the flapping aerodynamic power and the total input power of the system. While analysis of the kinematics for the different wing shapes and hinge stiffnesses reveals that the optimal performance occurs when the timing of pitch and stroke reversals are matched, thus supporting the effective transfer of input power from flapping to passive pitching and into the fluid. These results suggest that careful optimisation between wing shapes and hinge properties can allow insects and robots to exploit the passive dynamics to improve flight performance.