Both Schwartz (J. Fluid Mech., vol. 62 (3), 1974, pp. 553–578) and Cokelet (Phil. Trans. R. Soc. Lond., vol. 286 (1335), 1977, pp. 183–230) failed to gain convergent results for limiting Stokes waves in extremely shallow water by means of perturbation methods, even with the aid of extrapolation techniques such as the Padé approximant. In particular, it is extremely difficult for traditional analytic/numerical approaches to present the wave profile of limiting waves with a sharp crest of $120^{\circ }$ included angle first mentioned by Stokes in the 1880s. Thus, traditionally, different wave models are used for waves in different water depths. In this paper, by means of the homotopy analysis method (HAM), an analytic approximation method for highly nonlinear equations, we successfully gain convergent results (and especially the wave profiles) of the limiting Stokes waves with this kind of sharp crest in arbitrary water depth, even including solitary waves of extreme form in extremely shallow water, without using any extrapolation techniques. Therefore, in the frame of the HAM, the Stokes wave can be used as a unified theory for all kinds of waves, including periodic waves in deep and intermediate depths, cnoidal waves in shallow water and solitary waves in extremely shallow water.

The impact of a deep-water plunging breaker on a finite height two-dimensional structure with a vertical front face is studied experimentally. The structure is located at a fixed horizontal position relative to a wave maker and the structure’s bottom surface is located at a range of vertical positions close to the undisturbed water surface. Measurements of the water surface profile history and the pressure distribution on the front surface of the structure are performed. As the vertical position, $z_{b}$ (the $z$ axis is positive up and $z=0$ is the mean water level), of the structure’s bottom surface is varied from one experimental run to another, the water surface evolution during impact can be categorized into three classes of behaviour. In class I, with $z_{b}$ in a range of values near $-0.1\unicode[STIX]{x1D706}_{0}$, where $\unicode[STIX]{x1D706}_{0}$ is the nominal wavelength of the breaker, the behaviour of the water surface is similar to the flip-through phenomena first described in studies with shallow water and a structure mounted on the sea bed. In the present work, it is found that the water surface between the front face of the structure and the wave crest is well fitted by arcs of circles with a decreasing radius and downward moving centre as the impact proceeds. A spatially and temporally localized high-pressure region was found on the impact surface of the structure and existing theory is used to explore the physics of this phenomenon. In class II, with $z_{b}$ in a range of values near the mean water level, the bottom of the structure exits and re-enters the water phase at least once during the impact process. These air–water transitions generate large-amplitude ripple packets that propagate to the wave crest and modify its behaviour significantly. At $z_{b}=0$, all sensors submerged during the impact record a nearly in-phase high-frequency pressure oscillation indicating possible air entrainment. In class III, with $z_{b}$ in a range of values near $0.03\unicode[STIX]{x1D706}_{0}$, the bottom of the structure remains in air before the main crest hits the bottom corner of the structure. The subsequent free surface behaviour is strongly influenced by the instantaneous momentum of the local flow just before impact and the highest wall pressures of all experimental conditions are found.

The wandering of a wing-tip vortex in free-stream turbulence was documented by analysis of multi-probe hot-wire measurements in a wind tunnel and flow visualisation and particle image velocimetry measurements in a water tunnel. An error-minimisation approach was applied to the hot-wire measurements to estimate the time history of the location of the vortex axis, whereas flow visualisation from two orthogonal views permitted the reconstruction of relatively long sections of the vortex axis. The amplitude of the wandering motion was found to scale with the turbulence intensity, the core radius and the vortex turnover time; this amplitude was insensitive to changes in the integral length and time scales of the turbulence. The period of the vortex wandering was distributed in the range between 1 and 10 vortex turnover times. The wavelength of wandering was distributed at a relatively long value, which scaled with the vortex turnover time. The velocity of vortex wandering depended on the vortex turnover time, but also contained an additional contribution that was consistent with motion induced by bending waves. The prevalence of the vortex turnover time as the scale for vortex wandering was interpreted as evidence that vortex-induced straining of the free-stream eddies bounds the interaction time between the two, thus limiting the time available for linear and angular momentum transfer.

It is well known that deterministic two-dimensional marine ice sheets can only be stable if the grounding line is positioned at a sufficiently steep, downward sloping bedrock. When bedrock conditions favour instabilities, multiple stable ice sheet profiles may occur. Here, we employ continuation techniques to examine the sensitivity of a two-dimensional marine ice sheet to stochastic noise representing short time scale variability, either in the accumulation rate or in the sea level height. We find that in unique regimes, the position of the grounding line is most sensitive to noise in the accumulation rate and can explain excursions observed in field measurements. In the multiple equilibrium regime, there is a strong asymmetry in transition probabilities between the different ice sheet states, with a strong preference to switch to the branch with a steeper bedrock slope.

Under suitable conditions, an immersed granular bed can be destabilized by local thermal forcing and the induced buoyant force. The destabilization is evident from the triggering and establishment of a dense fluid-like granular plume. Varying the initial granular layer average height $h$, a time series of the free layer surface is extracted, allowing us to dynamically compute the underlying volume of the granular layer. Different observed phenomena, namely the initial interface deformation, the lowering of the average granular interface (i.e. decrease of the granular layer volume) and the emission of a plume, are analysed. We show that the phenomenon is mainly driven by heat transfer, for large $h$ and also involves a variable height thermal boundary condition and Darcy flow triggering, for small $h$. Simple modelling with no adjustable parameters not only allows us to capture the observed scaling power laws but is also in quantitative agreement with the obtained experimental data.

We study granular chute flow using the classic Saint-Venant approach, with the shear stresses within the granular sheet being incorporated via a friction law due to Pouliquen & Forterre (J. Fluid Mech., vol. 453, 2002, pp. 113–151) and with the in-plane stresses (which are ignored in the traditional formulation for normal fluids) being represented by a viscous-like term recently derived by Gray & Edwards (J. Fluid Mech., vol. 755, 2014, pp. 503–534). On the basis of this model, we predict that the granular sheet is able to sustain monoclinal waves, i.e. travelling shock structures that monotonically connect a thick region of uniform flow to a thinner one. We examine the balance of forces that determine the shape of this particular waveform and give the precise window of system parameters for which monoclinal waves are expected to appear in experiments.

We investigate electrophoresis of a tightly fitting sphere of radius $R-h_{0}$ on the axis of a circular tube of radius $R$, using lubrication theory and ideas due to Schnitzer & Yariv (Phys. Fluids, vol. 26, 2014, 122002). The electrical charge clouds on both the cylindrical wall and the surface of the sphere are assumed thin compared to the gap between the sphere and cylinder, so that charge clouds do not overlap and ion exclusion effects are minimal. Nevertheless, non-uniform pumping of counter-ions within the charge clouds leads to a change in the ionic concentration outside the charge clouds in the narrow gap between sphere and cylinder. The electro-osmotic slip velocities at the two surfaces are modified, leading to a decrease in the electrophoretic velocity of the sphere at low Péclet numbers and an increase in the velocity at high Péclet numbers. When the field strength $E_{0}$ is low, it is known that the electrophoretic velocity $U_{0}$ is proportional to $E_{0}(\unicode[STIX]{x1D701}_{s}-\unicode[STIX]{x1D701}_{c})$ which is zero when the zeta potential $\unicode[STIX]{x1D701}_{s}$ on the sphere surface is equal to the zeta potential $\unicode[STIX]{x1D701}_{c}$ on the cylinder. The perturbation to the above low field strength electrophoretic velocity (at high Péclet number) is predicted to be proportional to $E_{0}^{3}(\unicode[STIX]{x1D701}_{s}+\unicode[STIX]{x1D701}_{c})^{2}(\unicode[STIX]{x1D70E}_{s}+\unicode[STIX]{x1D70E}_{c})$, where $\unicode[STIX]{x1D70E}_{s}$ and $\unicode[STIX]{x1D70E}_{c}$ are the surface charge densities on the sphere and cylinder. The choice of materials with similar or identical zeta potentials (and surface charge densities) for the cylinder and sphere should therefore facilitate the observation of velocities nonlinear in the field strength $E_{0}$, since the reference linear electrophoretic velocity will be small.

The motion of tracer particles is kinematically simulated around three elementary flow patterns; a Burgers vortex, a shear-layer structure with coincident vortices and a node-saddle topology. These patterns are representative for their broader class of coherent structures in turbulence. Therefore, examining the dispersion in these elementary structures can improve the general understanding of turbulent dispersion at short time scales. The shear-layer structure and the node-saddle topology exhibit similar pair dispersion statistics compared to the actual turbulent flow for times up to $3{-}10\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D702}}$, where, $\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D702}}$ is the Kolmogorov time scale. However, oscillations are observed for the pair dispersion in the Burgers vortex. Furthermore, all three structures exhibit Batchelor’s scaling. Richardson’s scaling was observed for initial particle pair separations $r_{0}\leqslant 4\unicode[STIX]{x1D702}$ for the shear-layer topology and the node-saddle topology and was related to the formation of the particle sheets. Moreover, the material line orientation statistics for the shear-layer and node-saddle topology are similar to the actual turbulent flow statistics, up to at least $4\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D702}}$. However, only the shear-layer structure can explain the non-perpendicular preferential alignment between the material lines and the direction of the most compressive strain, as observed in actual turbulence. This behaviour is due to shear-layer vorticity, which rotates the particle sheet generated by straining motions and causes the particles to spread in the direction of compressive strain also. The material line statistics in the Burgers vortex clearly differ, due to the presence of two compressive principal straining directions as opposed to two stretching directions in the shear-layer structure and the node-saddle topology. The tetrad dispersion statistics for the shear-layer structure qualitatively capture the behaviour of the shape parameters as observed in actual turbulence. In particular, it shows the initial development towards planar shapes followed by a return to more volumetric tetrads at approximately $10\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D702}}$, which is associated with the particles approaching the vortices inside the shear layer. However, a large deviation is observed in such behaviour in the node-saddle topology and the Burgers vortex. It is concluded that the results for the Burgers vortex deviated the most from actual turbulence and the node-saddle topology dispersion exhibits some similarities, but does not capture the geometrical features associated with material lines and tetrad dispersion. Finally, the dispersion around the shear-layer structure shows many quantitative (until 2–$4\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D702}}$) and qualitative (until $20\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D702}}$) similarities to the actual turbulence.

We investigate the onset of instability in a stably stratified two-component fluid in a vertical Hele-Shaw cell when the unstably stratified scalar has a settling velocity. This linear stability problem is analysed on the basis of Darcy’s law, for constant-gradient base states. The settling velocity is found to trigger a novel instability mode characterized by pairs of inclined waves. For unequal diffusivities, this new settling-driven mode competes with the traditional double-diffusive mode. Below a critical value of the settling velocity, the double-diffusive elevator mode dominates, while, above this threshold, the inclined waves associated with the settling-driven instability exhibit faster growth. The analysis yields neutral stability curves and allows for the discussion of various asymptotic limits.

The turbulent–non-turbulent interface (TNTI) of supersonic turbulent boundary layers is a fundamental but relatively unexplored physics problem. In this study, we present experimental results from fractal analysis on the TNTI of supersonic turbulent boundary layers, and test the applicability of the additive law for these flows. By applying the nanoparticle-tracer planar laser scattering (NPLS) technique in a supersonic wind tunnel, we obtain data covering nearly three decades in scale. The box-counting results indicate that the TNTI of supersonic turbulent boundary layers is a self-similar fractal with a fractal dimension of 2.31. By comparing data sets acquired from two orthogonal planes, we find that the scaling exponent does not depend on direction, consistent with the validity of the additive law for the TNTI of turbulent boundary layers in a scale range with the large-scale limit not exceeding approximately $0.05\unicode[STIX]{x1D6FF}$.