We investigate the effect of helicity on the scale-similar structures of homogeneous isotropic and non-mirror-symmetric turbulence based on the Lagrangian renormalised approximation (LRA), which is a self-consistent closure theory proposed by Kaneda (J. Fluid Mech., vol. 107, 1981, pp. 131–145). In this study, we focus on the time scale representing the scale-similar range. For the LRA, the Lagrangian two-time velocity correlation and response function determine the representative time scale. The LRA predicts that both the Lagrangian two-time velocity correlation and response function equation do not explicitly depend on helicity. We assume the extended scale-similar spectra and time scale by considering the helicity dissipation rate. Considering the small-scale structures, the requirements for the energy and helicity fluxes and response function equation to be scale similar, yield the conventional inertial-range power laws and provide the energy and helicity spectra $\propto k^{-5/3}$ and the time scale $\propto \varepsilon ^{-1/3} k^{-2/3}$, where $\varepsilon$ and $k$ denote the energy dissipation rate and wavenumber, respectively. Notably, energy flux can be scale similar only when $k^H /k \ll 1$, where $k^H = \varepsilon ^H/\varepsilon$ and $\varepsilon ^H$ denotes the helicity dissipation rate. This condition makes the energy cascade process in the scale-similar range completely independent of helicity. We also investigate the localness of the interscale interaction in the energy and helicity cascades for the LRA. We demonstrate that the helicity cascade is slightly non-local in scales compared with the energy cascade. This study provides a foundation on the modelling of non-mirror-symmetric turbulent flows.

We present theoretical results related to the experimental findings of Matsubara & Alfredsson (J. Fluid Mech., vol. 430, 2001, pp. 149–168) on the scaling of the energy spectra of the Klebanoff modes, i.e. streamwise-elongated vortical disturbances generated by free-stream turbulence in a flat-plate transitional boundary layer. The scaling is explained by a model that describes the streamwise evolution of the streamwise and spanwise energy spectra. The theoretical framework is based on the quasi-steady asymptotic solution of the boundary-region equations, on an axial-symmetric model of the free-stream spectrum, and on the spectral response of the boundary layer to the external perturbations.

Water particle orbits are key elements in the Lagrange wave formulation. The stochastic Miche implementation of the Lagrange model is a linear Gaussian two-dimensional or three-dimensional space–time model which exhibits typical nonlinear wave characteristics when transformed to Eulerian coordinates. This paper investigates the statistical relation between the degree of front–back asymmetry of individual waves and the orbit orientation for the particle located at the wave maximum at the point of observation. It is shown that, in the Lagrangian model with statistical front–back symmetry, for individual waves there is a clear connexion between the degree of individual wave asymmetry and the orientation of the randomly deformed elliptic orbit: a steep front correlates with upward tilt, a steep back is correlated with a downward tilt. This holds both for waves observed in time and in space, and the dependence is stronger for large amplitude waves than for smaller ones. The dependence is strongly dependent on the depth and on the significant steepness and spectral width. Inclusion of the average Stokes drift has a moderate effect on the dependence. For models with forced front–back asymmetry there is both a systematic dependence and a statistical correlation between asymmetry and tilt; for large amplitude waves the systematic relation dominates. The conclusions are based on Fourier simulations of Gauss–Lagrange waves of first and second order with a wind–sea Pierson–Moskowitz spectrum and a narrow swell JONSWAP spectrum.

We present direct numerical simulations of developing turbulent channel flows subjected to thermal expansion or contraction downstream of a heated or cooled wall. Using different constitutive relations for viscosity we analyse the response of variable property flows to streamwise acceleration/deceleration by separating the effect of streamwise acceleration/deceleration from the effect of wall-normal property variations. We demonstrate that, beyond a certain streamwise location, the flow can be considered in a state of ‘quasi-equilibrium’ regarding semilocally scaled variables. As such, we claim that the development of turbulent quantities due to streamwise acceleration/deceleration is localized to the region of impulsive heating/cooling, while changes in turbulence occurring farther downstream can be attributed solely to property variations. This finding allows us to study turbulence modulation in accelerating/decelerating flows using the semilocal scaling framework. By investigating the energy redistribution among the turbulent velocity fluctuations, we conclude that a change in bulk streamwise velocity has a non-local effect which originates from the change in mean shear and modifies the energy pathways through velocity-pressure-gradient correlations. On the other hand, the wall-normal property gradients have a local effect and act through the modification of the viscous dissipation. We show that it is possible to superimpose and compare the two different effects when using the semilocal scaling framework.

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.

We study concentration-driven natural convection boundary layers on horizontal surfaces, subjected to a weak, surface normal, uniform blowing velocity $V_i$ for three orders of range of the dimensionless blowing parameter $10^{-8}\le J=Re_x^3/Gr_x\le 10^{-5}$, where $Re_x$ and $Gr_x$ are the local Reynolds and Grashof numbers at the horizontal location $x$, based respectively on $V_i$ and ${\rm \Delta} C$, the concentration difference across the boundary layer. We formulate the integral boundary layer equations, with the assumption of no concentration drop within the species boundary layer, which is valid for weak blowing into the thin species boundary layers that occur at the high Schmidt number ($Sc \simeq 600$) of concentration-driven convection. The equations are then numerically solved to show that the species boundary layer thickness $\delta _d = 1.6\,x(Re_x/Gr_x)^{1/4}$, the velocity boundary layer thickness $\delta _v=\delta _d Sc^{1/5}$, the horizontal velocity $u = V_i(Gr_x/Re_x)^{1/4}f(\eta )$, where $\eta =y/\delta _v$, and the drag coefficient based on $V_i$, $C_D = 2.32/\sqrt {J}$. We find that the vertical profile of the horizontally averaged dimensionless concentration across the boundary layer becomes, surprisingly, independent of the blowing and the species diffusion effects to follow a $Gr_y^{2/3}$ scaling, where $Gr_y$ is the Grashof number based on the vertical location $y$ within the boundary layer. We then show that the above profile matches the experimentally observed mean concentration profile within the boundary layers that form on the top surface of a membrane, when a weak flow is forced gravitationally from below the horizontal membrane that has brine above it and water below it. A similar match between the theoretical scaling of the species boundary layer thickness and its experimentally observed variation is also shown to occur.

We present a local stability analysis of an idealized model of the stratified vortices that appear in geophysical settings. The base flow comprises an axisymmetric vortex with background rotation and an out-of-plane stable stratification, and a radial stratification in the thermal wind balance with the out-of-plane momentum gradient. Solving the local stability equations along fluid particle trajectories in the base flow, the dependence of short-wavelength instabilities on the Schmidt number $Sc$ (ratio between momentum and mass diffusivities) is studied, in the presence of curvature effects. In the diffusion-free limit, the well-known symmetric instability is recovered. In the viscous, double-diffusive regime, instability characteristics are shown to depend on three non-dimensional parameters (including $Sc$), and two different instabilities are identified: (i) a monotonic instability (same as symmetric instability at $Sc = 1$), and (ii) an oscillatory instability (absent at $Sc = 1$). Separating the base flow and perturbation characteristics, two each of base flow and perturbation parameters (apart from $Sc$) are identified, and the entire parameter space is explored for the aforementioned instabilities. In comparison with $Sc = 1$, monotonic and oscillatory instabilities are shown to significantly expand the instability region in the space of base flow parameters as $Sc$ moves away from unity. Neutral stability boundaries on the plane of $Sc$ and a modified gradient Richardson number are then identified for both these instabilities. In the absence of curvature effects, our results are shown to be consistent with previous studies based on normal mode analysis, thus establishing that the local stability approach is well suited to capturing symmetric and double-diffusive instabilities. The paper concludes with a discussion of curvature effects, and the likelihood of monotonic and oscillatory instabilities in typical oceanic settings.

The two-dimensional stability of vertically sheared inertial oscillations at ocean fronts is explored through a linear stability analysis and nonlinear simulations. Baroclinic effects reduce the minimum frequency of inertia-gravity waves to an extent determined by the balanced Richardson number ${{Ri}}$ of the front. Below a critical value of ${{Ri}}$, which depends on the strength of the inertial shear, the inertial oscillations become unstable to parametric subharmonic instability (PSI) resulting in growing perturbations that oscillate at half the inertial frequency $f$. Since the critical value is always greater than 1, PSI can occur at fronts stable to symmetric instability. Although modest in weak inertial shear, growth rates exceeding $f/2$ can be achieved for inertial shear greater than or equal to the thermal wind shear. Our formulation allows for non-hydrostatic perturbations and can be applied to initially unstratified geostrophic adjustment problems. We find that PSI will almost totally damp the transient oscillations that arise during geostrophic adjustment. The perturbations gain energy at the expense of the inertial oscillations through ageostrophic shear production. The perturbations then themselves become unstable to secondary Kelvin–Helmholtz instabilities creating a pathway by which the inertial oscillations can be dissipated rapidly. In contrast to symmetric and baroclinic instabilities that draw on a front's kinetic or potential energy, PSI acts to increase the energy stored in the balanced front as the convergence and divergence of the eddy-momentum fluxes set up a secondary circulation in the sense to stand up the front.

Energy transfer in turbulent fluids is non-Gaussian. We quantify non-Gaussian energy transfer between the atmosphere and bodies of water using a turbulent diffusion operator coupled with temporally self-affine velocity distributions and a recursive integration method that produce multifractal measures. The measures serve as input to a system of moment field equations (derived from Navier–Stokes) that generate and track high-frequency gravity waves that propagate through the water surface (as a result of the air–water interactions). The dimension of the support of the air–water turbulence produced by our methods falls within the range of theory and observation, and correspondingly, hindcast statistical measures of the water-wave surface such as significant water-wave height and wave period are well correlated to observational buoy data. Further, our recursive integration method can be used by spectral resolving phase-averaged models to interpolate temporal wind data to smaller scales to capture the non-Gaussian behaviour of the air–water interaction.

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.

Theoretical analysis of attachment-line instabilities is performed for supersonic swept flows using the compressible Hiemenz approximation for the mean flow and the successive approximation procedures for disturbances. The theoretical model captures the dominant attachment-line modes in wide ranges of the sweep Mach number ${M_e}$ and the wall temperature ratio. It is shown that these modes behave similar to the first and second Mack modes in the boundary layer flow. This similarity allows us to extrapolate the knowledge gained for Mack modes to the attachment-line instabilities. In particular, we find that at sufficiently large ${M_e}$, the dominant attachment-line instability is associated with the synchronisation of slow and fast modes of acoustic nature. Point-by-point comparisons of the theoretical predictions with the experiments of Gaillard et al. (Exp. Fluids, vol. 26, 1999, pp. 169–176) demonstrate that at ${M_e} > 4$, the theory captures a significant drop of the transition onset Reynolds number, which is below the contamination criterion of Poll $({R_\mathrm{\ast }} = 250)$ at ${M_e} > 6$. This contradicts the generally accepted assumption that the attachment-line flow is stable for ${R_\mathrm{\ast }} \le 250$. The theoretical critical Reynolds numbers lie well below the experimental transition-onset Reynolds numbers. Stability computations using the Navier–Stokes mean flow and accounting for the leading-edge curvature effect do not eliminate this discrepancy. Most likely, in the experiments of Gaillard et al., we face with an unknown effect that does not fit to the concept of transition arising from linear instability.

This study investigates the scenario of flow past a circular cylinder in proximity to a moving wall (or equally a body translating in still fluid parallel to a stationary wall). Fifty high-fidelity three-dimensional direct numerical simulations are performed over a parameter space of turbulent Reynolds numbers (Re) of 300 to 1000 combined with gap-to-diameter ratios (G/D) of 0.2 to 3. The flow, hydrodynamic and turbulence characteristics over the (Re, G/D) parameter space are examined in detail. Small-scale rib-like mode B structures and alternate vortex shedding are observed for all cases. The streamwise location for the vortex shedding (quantified by the wake recirculation length Lr) varies strongly with both Re and G/D. The variation of Lr with Re and G/D can be explained by the spanwise circulation $\varGamma_z$ fed into the wake, where the variation trends of Lr and $\varGamma_z$ are inversely correlated. The variations of the mean drag and fluctuating lift coefficients are also governed by the variations in Lr and $\varGamma_z$. The total kinetic energy in the wake region reduces drastically as G/D reduces below 0.8, which is contributed collectively by (i) reduction in the strength of the shed vortices, (ii) downstream movement of the location of vortex shedding and (iii) associated delayed generation of streamwise vortices. The present results on a moving wall also help to explain several flow and hydrodynamic characteristics reported in the literature for a stationary wall, because the moving wall eliminates the complex wall boundary layer and retains a ‘clean’ near-wall effect.

Saturated flow film boiling on a sphere has been numerically studied in this work for both vertical and horizontal flow configurations. The simulations were performed using a numerical methodology developed by the authors for boiling flows on three-dimensional unstructured meshes. For interface capturing, the coupled level set and volume of fluid method is used. The interface evolution, vapour wake dynamics and heat transfer have been thoroughly investigated by varying the saturated liquid flow velocity, sphere diameter and wall superheat. The relative importance of both the buoyancy and the inertial forces is described in terms of the Froude number $(Fr)$. The vapour bubble evolves periodically at low $Fr$ values, while a stable vapour column develops at high $Fr$ values. The interface evolution pattern obtained in the present work is in good agreement with the results of experimental studies available in the literature. For all the values of $Fr$, a stable vapour column develops for a large-diameter sphere and releases vapour bubbles of varying sizes. Furthermore, for a large-diameter sphere, surface capillary waves are observed at the interface, similar to the observations of some of the experimental studies available in the literature. The flow in the liquid and vapour wakes appears to be strongly coupled. The heat transfer in the present work is estimated using the spatially and temporally averaged Nusselt numbers. Finally, an fast Fourier transform analysis of the space-averaged Nusselt number reveals a strong interaction among the different forces.

The behaviour of turbulent, buoyant, planar plumes is fundamentally coupled to the environment within which they develop. The effect of a background stratification directly influences a plumes buoyancy and has been the subject of numerous studies. Conversely, the effect of an ambient co-flow, which directly influences the vertical momentum of a plume, has not previously been the subject of theoretical investigation. The governing conservation equations for the case of a uniform co-flow are derived and the local dynamical behaviour of the plume is shown to be characterised by the scaled source Richardson number and the relative magnitude of the co-flow and plume source velocities. For forced, pure and lazy plume release conditions the co-flow acts to narrow the plume and reduce both the dilution and the asymptotic Richardson number relative to the classic zero co-flow case. Analytical solutions are developed for pure plumes from line sources, and for highly forced and highly lazy releases from sources of finite width in a weak co-flow. Contrary to releases in quiescent surroundings, our solutions show that all classes of release can exhibit plume contraction and the associated necking. For entraining plumes, a dynamical invariance spatially only occurs for pure and forced releases and we derive the co-flow strengths that lead to this invariance.

We report a numerical investigation of a previously noticed but less explored flow state transition in two-dimensional turbulent Rayleigh–Bénard convection. The simulations are performed in a square domain over a Rayleigh number range of $10^7 \leq Ra \leq 2 \times 10^{11}$ and a Prandtl number range of $0.25 \leq Pr \leq 20$. The transition is characterized by the emergence of multiple satellite eddies with increasing $Ra$, which orbit around and interact with the main vortex roll in the system. Consequently, the main roll is squeezed to a smaller size compared with the domain and wanders around in the bulk region irregularly and extensively. This is in sharp contrast to the flow state before the transition, which is featured by a domain-sized circulatory roll with its vortex centre ‘condensed’ near the domain's centre. Detailed velocity field analysis reveals that there exists an abrupt increase in the energy fluctuations of the Fourier modes during the transition. Based on this phase-transition-like signal, the critical condition for the transition is found to follow a scaling relation as $Ra_t \sim Pr^{1.41}$ where $Ra_t$ is the critical Rayleigh number for the transition. This scaling relation is quantitatively explained by a phenomenological model grounded on the bistability behaviour (i.e. spontaneous and stochastic switching between the two flow states) observed at the edge of the transition. The model can also account for the effects of aspect ratio on the transition reported in the literature (van der Poel et al., Phys. Fluids, vol. 24, 2012).