Deep-water surface wave breaking affects the transfer of mass, momentum, energy and heat between the air and sea. Understanding when and how the onset of wave breaking will occur remains a challenge. The mechanisms that form unforced steep waves, i.e. nonlinearity or dispersion, are thought to have a strong influence on the onset of wave breaking. In two dimensions and in deep water, spectral bandwidth is the main factor that affects the roles these mechanism play. Existing studies, in which the relationship between spectral bandwidth and wave breaking onset is investigated, present varied and sometimes conflicting results. We perform potential-flow simulations of two-dimensional focused wave groups on deep water to better understand this relationship, with the aim of reconciling existing studies. We show that the way in which steepness is defined may be the main source of confusion in the literature. Locally defined steepness at breaking onset reduces as a function of bandwidth, and globally defined (spectral) steepness increases. The relationship between global breaking onset steepness and spectral shape (using the parameters bandwidth and spectral skewness) is too complex to parameterise in a general way. However, we find that the local surface slope of maximally steep non-breaking waves, of all spectral bandwidths and shapes that we simulate, approaches a limit of $1/\tan ({\rm \pi} /3)\approx 0.5774$. This slope-based threshold is simple to measure and may be used as an alternative to existing kinematic breaking onset thresholds. There is a potential link between slope-based and kinematic breaking onset thresholds, which future work should seek to better understand.

The first long-lived turbulent structures observable in planar shear flows take the form of localized stripes, inclined with respect to the mean flow direction. The dynamics of these stripes is central to transition, and recent studies proposed an analogy to directed percolation where the stripes’ proliferation is ultimately responsible for the turbulence becoming sustained. In the present study we focus on the internal stripe dynamics as well as on the eventual stripe expansion, and we compare the underlying mechanisms in pressure- and shear-driven planar flows, respectively, plane-Poiseuille and plane-Couette flow. Despite the similarities of the overall laminar–turbulence patterns, the stripe proliferation processes in the two cases are fundamentally different. Starting from the growth and sustenance of individual stripes, we find that in plane-Couette flow new streaks are created stochastically throughout the stripe whereas in plane-Poiseuille flow streak creation is deterministic and occurs locally at the downstream tip. Because of the up/downstream symmetry, Couette stripes, in contrast to Poiseuille stripes, have two weak and two strong laminar turbulent interfaces. These differences in symmetry as well as in internal growth give rise to two fundamentally different stripe splitting mechanisms. In plane-Poiseuille flow splitting is connected to the elongational growth of the original stripe, and it results from a break-off/shedding of the stripe's tail. In plane-Couette flow splitting follows from a broadening of the original stripe and a division along the stripe into two slimmer stripes.

A microscale lubrication flow of a gas between eccentric circular cylinders is studied on the basis of kinetic theory. The dimensionless curvature, defined by the mean clearance divided by the radius of the inner cylinder, is small, and the rotation speed of the inner cylinder is also small. The Knudsen number, defined by the mean free path divided by the mean clearance, is arbitrary. The Boltzmann equation is studied analytically using the slowly varying approximation following the method proposed in the author's previous study (Doi, Phys. Rev. Fluids, vol. 7, 2022, 034201). A macroscopic lubrication equation, which is a microscale generalization of the Reynolds lubrication equation, is derived. To assess this, a direct numerical analysis of the Boltzmann equation in a bipolar coordinate system is conducted using the Bhatnagar–Gross–Krook–Welander kinetic equation. It is demonstrated that the solution of the derived lubrication equation approximates that of the Boltzmann equation over a wide range of the eccentricity and the whole range of the Knudsen number. It is also demonstrated that another lubrication equation derived by a formal application of the slowly varying approximation produces a non-negligible error of the order of the square root of the dimensionless curvature for large Knudsen numbers.

The third-order law links energy transfer rates in the inertial range of magneto- hydrodynamic (MHD) turbulence with third-order structure functions. Anisotropy, a typical property in the solar wind, challenges the applicability of the third-order law with the isotropic assumption. To shed light on the energy transfer process in the presence of anisotropy, we conducted direct numerical simulations of forced MHD turbulence with normal and hyper-viscosity under various strengths of the external magnetic field ($B_0$), and calculated three forms of third-order structure function with or without averaging of the azimuthal or polar angles with respect to $B_0$ direction. Correspondingly, three estimated energy transfer rates were obtained. The result shows that the peak of normalized third-order structure function occurs at larger scales closer to the $B_0$ direction, and the maximum of longitudinal transfer rates shifts away from the $B_0$ direction at larger $B_0$. Compared with normal viscous cases, hyper-viscous cases can attain better separated inertial range, thus facilitating the estimation of the energy cascade rates. We find that the widespread use of the isotropic form of the third-order law in estimating the energy transfer rates is questionable in some cases, especially when the anisotropy arising from the mean magnetic field is inevitable. In contrast, the direction-averaged third-order structure function properly accounts for the effect of anisotropy and predicts the energy transfer rates and inertial range accurately, even at very high $B_0$. With limited statistics, the third-order structure function shows a stronger dependence on averaging of azimuthal angles than the time, especially for high $B_0$ cases. These findings provide insights into the anisotropic effect on the estimation of energy transfer rates.

High-speed vehicles experience a highly challenging environment in which the freestream Mach number and surface temperature greatly influence aerodynamic drag and heat transfer. The interplay of these two parameters strongly affects the near-wall dynamics of high-speed turbulent boundary layers (TBLs) in a non-trivial way, breaking similarity arguments on velocity and temperature fields, typically derived for adiabatic cases. We present direct numerical simulations of flat-plate zero-pressure-gradient TBLs spanning three freestream Mach numbers $[2,4,6]$ and four wall temperature conditions (from adiabatic to very cold walls), emphasising the choice of the wall-cooling parameter to recover a similar flow organisation at different Mach numbers. We link qualitative observations on flow patterns to first- and second-order statistics to explain the decoupling of temperature–velocity fluctuations that occurs at reduced wall temperatures and high Mach numbers. For these cases, we discuss the formation of a secondary peak of thermal production in the viscous sublayer, which is in contrast with the monotonic behaviour of adiabatic profiles. We propose different physical mechanisms induced by wall-cooling and compressibility that result in apparently similar flow features, such as a higher peak in the streamwise velocity turbulence intensity, and distinct features, such as the separation of turbulent scales.

Turbulent flows above a solid surface are characterised by a hydrodynamic roughness that represents, for the far velocity field, the typical length scale at which momentum mixing occurs close to the surface. Here, we are theoretically interested in the hydrodynamic roughness induced by a two-dimensional modulated surface, the elevation profile of which is decomposed in Fourier modes. We describe the flow for a sinusoidal mode of given wavelength and amplitude with Reynolds-averaged Navier–Stokes equations closed by means of a mixing-length approach that takes into account a possible surface geometrical roughness as well as the presence of a viscous sublayer. It also incorporates spatial transient effects at the laminar–turbulent transition. Performing a weekly nonlinear expansion in the bedform aspect ratio, we predict the effective hydrodynamic roughness when the surface wavelength is varied and we show that it presents a non-monotonic behaviour at the laminar–turbulent transition when the surface is hydrodynamically smooth. Further, with a self-consistent looped calculation, we are able to recover the smooth–rough transition of a flat surface, for which the hydrodynamic roughness changes from a regime where it is dominated by the viscous length to another one where it scales with the surface corrugation. We finally apply the results to natural patterns resulting from hydrodynamic instabilities such as those associated with dissolution or sediment transport. We discuss in particular the aspect ratio selection of dissolution bedforms and roughness hierarchy in superimposed ripples and dunes.

A theoretical and experimental investigation of two-dimensional (2-D) liquid curtains (gravitationally thinning liquid sheets) is provided under conditions where the curtain issues from a thin slot whose centreline is inclined with respect to the vertical. This analysis is motivated in part by recent works where it has been proposed that oblique liquid curtains (those exiting a non-vertical slot) may bend upwards against gravity when the relevant Weber number at the slot is less than unity ($We <1$). By contrast, Weinstein et al. (J. Fluid Mech., vol. 876, 2019, R3) have proposed that such $We<1$ curtains must be vertical and downward falling regardless of the inclination of the slot. Under low-Reynolds-number ($Re$) conditions typical of liquid film coating operations, our experiments show that the curtain shape follows the classic ballistic (parabolic) trajectory in the supercritical regime ($We>1$). In subcritical conditions ($We<1$), experiments show that the downward-falling curtain is vertical except in a relatively small region near the slot, where the combined effects of viscosity and surface tension induce the so-called teapot effect. These experimental results are confirmed by 2-D numerical simulations, which predict the curtain behaviour ranging from highly viscous ($Re = O(1)$) to nearly inviscid conditions. The one-dimensional (1-D) inviscid model of Weinstein et al. is recast in a different form to facilitate comparisons with the 2-D model, and 1-D and 2-D results agree favourably for supercritical and subcritical conditions. Despite the large parameter range explored, we have found no evidence that upward-bending curtains exist in an oblique configuration.