To get an insight into the dynamics of the oceanic surface boundary layer we develop an asymptotic model of the nonlinear dynamics of linearly decaying three-dimensional long-wave perturbations in weakly stratified boundary-layer flows. Although in nature the free-surface boundary layers in the ocean are often weakly stratified due to solar radiation and air entrainment caused by wave breaking, weak stratification has been invariably ignored. Here, we consider an idealized hydrodynamic model, where finite-amplitude three-dimensional perturbations propagate in a horizontally uniform unidirectional weakly stratified shear flow confined to a boundary layer adjacent to the water surface. Perturbations satisfy the no-stress boundary condition at the surface. They are assumed to be long compared with the boundary-layer thickness. Such perturbations have not been studied even in a linear setting. By exploiting the assumed smallness of nonlinearity, wavenumber, viscosity and the Richardson number, on applying triple-deck asymptotic scheme and multiple-scale expansion, we derive in the distinguished limit a novel essentially two-dimensional nonlinear evolution equation, which is the main result of the work. The equation represents a generalization of the two-dimensional Benjamin–Ono equation modified by the explicit account of viscous effects and new dispersion due to weak stratification. It describes perturbation dependence on horizontal coordinates and time, while its vertical structure, to leading order, is given by an explicit analytical solution of the linear boundary value problem. It shows the principal importance of weak stratification for three-dimensional perturbations.

A Lagrangian approach to both hydrostatic non-dispersive in the short-wave range and non-hydrostatic dispersive rotating shallow-water magnetohydrodynamics is developed, and used to analyse weakly and fully nonlinear waves described by the model. Hyperbolic structure in the non-dispersive case is displayed and Riemann invariants are constructed. Characteristic equations are used to establish criteria of breaking and formation of shocks by magneto-gravity waves, and conditions of the appearance of contact discontinuities in Alfvén waves. As in the case of non-magnetic rotating shallow water, rotation cannot prevent breaking. The Lagrangian equations of the model are reduced to a single partial differential ‘master’ equation, which is used to analyse the propagation of weakly nonlinear waves of both families, with or without weak rotation, and with or without weak short-wave dispersion. Corresponding modulation equations are constructed and their main properties sketched. The same master equation is used to obtain fully nonlinear finite-amplitude wave solutions in particular cases of no short-wave dispersion or no rotation.

We present an investigation of the fundamental physical processes involved in deep-water gravity wave breaking. Our motivation is to identify the underlying reason causing the deficiency of the eddy viscosity breaking model (EVBM) in predicting surface elevation for strongly nonlinear waves. Owing to the limitation of experimental methods in the provision of high-resolution flow information, we propose a numerical methodology by developing an EVBM enclosed standalone fully nonlinear quasi-potential (FNP) flow model and a coupled FNP plus Navier–Stokes flow model. The numerical models were firstly verified with a wave train subject to modulational instability, then used to simulate a series of broad-banded focusing wave trains under non-, moderate- and strong-breaking conditions. A systematic analysis was carried out to investigate the discrepancies of numerical solutions produced by the two models in surface elevation and other important physical properties. It is found that EVBM predicts accurately the energy dissipated by breaking and the amplitude spectrum of free waves in terms of magnitude, but fails to capture accurately breaking induced phase shifting. The shift of phase grows with breaking intensity and is especially strong for high-wavenumber components. This is identified as a cause of the upshift of the wave dispersion relation, which increases the frequencies of large-wavenumber components. Such a variation drives large-wavenumber components to propagate at nearly the same speed, which is significantly higher than the linear dispersion levels. This suppresses the instant dispersive spreading of harmonics after the focal point, prolonging the lifespan of focused waves and expanding their propagation space.

New experiments in highly turbulent, steady, subcritical and uniform water open-channel flows have been carried out to measure the mean turbulent kinetic energy (TKE) budget of sediment-laden boundary layer flows with two sizes (dp = 3 mm and 1 mm) of Plexiglas particles $({\rm relative\ density}\ = 1.192)$. The experiments covered energetic sediment transport conditions (Shields number of $0.35 < \theta < 1.2$) ranging from non-capacity to full-capacity flows in bedload-to-suspension-dominated transport modes (suspension number of $0.5 < {w_s}/{u_\ast } < 1.3$ where ${w_s}$ is the settling velocity and $u_*$ is the friction velocity) and for weakly to highly inertial, finite size turbulence-particle conditions (Stokes number of 0.1 < St < 3.5 and dp/η > 10 where $\eta$ is the Kolmogorov length scale). It was shown that the effects of sediments on the TKE budget are very pronounced in all large particle experiments for which a bedload layer of several grain diameter thickness is developed above the channel bed. When compared with the corresponding reference clear-water flows, the TKE shear-production rate for the 3 mm particle flows is strongly reduced in the wall region corresponding to the bedload layer. This turbulence damping is seen to increase with sediment load until full capacity for flows with constant Shields value, as well as with Shields number value. Inside this damped TKE shear-production zone, a distinct peak of maximal turbulence production appears to coincide with the upper edge of the bedload layer delimited by a sharp gradient in mean sediment concentration. This vertically upshifted peak of TKE production is accompanied by an enhanced net downward oriented TKE flux when compared with the reference clear-water flows. The downward diffused TKE is found to act in the bedload layer as a local energy source in reasonable balance with the sediment transport term. The mechanism behind this downward TKE transport was further analysed on the basis of coherent flow structure dynamics controlled by ejection- and sweep-type events. The agreement between the height of downward directed mean TKE flux and the height below which sweep-type events dominate the Reynolds shear-stress contribution over ejections, revealed the leading role played by sweeps in mean TKE transport. This agreement holds for all reference clear-water flows supporting the well-known wall-roughness-induced dominance of the sweep contribution in turbulent, rough clear-water boundary layer flows. Furthermore, for all 3 mm particle flows, the two referred to transition levels were significantly and similarly upshifted to the upper edge of the bedload layer. Only for these sediment-laden flows, the bedload layer thickness is seen to exceed the wall-roughness sublayer of the reference clear-water flows. This supports a strong analogy between wall-roughness effects in clear-water flows and bedload layer effects in sediment-laden flows, on the mean TKE budget induced by a similarly modified coherent flow structure dynamics. The bedload layer-controlled wall roughness is finally confirmed by the good prediction of the wall-roughness parameter ks of the logarithmic velocity distribution. An empirical formulation fitting the presented measurements is presented, valid over the range of Shields number values covered herein.

Direct numerical simulations of channel flow and temporal boundary layer at a Reynolds number $Re_{\tau } = 1500$ are used to assess the scale-by-scale mechanisms of wall turbulence. From the peak of turbulence production embedded at the small scales of the near-wall region, spatially ascending reverse cascades are generated that move through self-similar eddies growing in size with the wall distance. These fluxes are followed by spatially ascending forward cascades through detached eddies thus reaching sufficiently small scales where eventually scale energy is dissipated. This phenomenology is shared by both boundary layer and channel flow and is recognized as a robust physical feature characterizing wall turbulence in general. Specific features related to the flow configuration are indeed identified in the outer region. In particular, the central region of channels is characterized by a generalized Richardson energy cascade where large scales are in equilibrium with small scales at different wall distances through a combined forward cascade and spatial flux. On the contrary, the interface region of boundary layers is characterized by an almost two-dimensional physics where spatially ascending reverse cascades sustain long and wide interface structures with a forward cascade that survives only in the wall-normal scales. The overall scenario consists in a variety of scale motions that while protruding from the turbulent core towards the external region, squeeze at the interface thus sustaining vertical shear in a thin layer. The observed multidimensional physics sheds light on the complex interactions between outer entrainment and near-wall self-sustaining mechanisms with possible repercussions for theories.

The geometrical properties of streamlines, such as the curvatures, directions and positions, are studied in steady inviscid compressible flow fields via differential geometry theories and conservation laws. The influences of the streamline geometries on the flow speeds and pressures are also identified and discussed. By transforming the streamlines to fill the domain and satisfy the boundary conditions, a unified geometry-based solver, the streamline transformation method, is proposed for both subsonic and supersonic regions. The governing equations and boundary conditions along streamlines and shock waves are also derived. This method is verified by numerical results of three typical flow fields, including the subsonic channel flow, the supersonic downstream of attached shock waves and especially the subsonic/supersonic downstream of detached bow shock waves. Both two-dimensional planar and axisymmetric flow fields are considered. Compared with the results from computational fluid dynamics, good agreements are achieved by this method, while fewer computational resources, by an order of magnitude, are consumed. Features of these flow fields are also analysed from a geometrical perspective, such as flow speeds and pressures deviated by the wall curvatures, and three-dimensional effects in the after-shock flow fields. For a hyperbolic-shaped bow shock wave, the stand-off distances and the transitions from subsonic to supersonic regions are also discussed. As indicated by the accuracy, efficiency and applicability in a wide range of flow speeds, the streamline transformation method would be a potential candidate for the theoretical analysis and inverse design of high-speed flow fields, especially where the subsonic regions exist downstream of strong shock waves.