We have performed particle-resolved direct numerical simulations of many heavy non-spherical particles settling under gravity in the dilute regime. The particles are oblate spheroids of aspect ratio $1.5$ and density ratio $1.5$. Two Galileo numbers are considered, namely $111$ and $152$, for which a single oblate spheroid follows a steady vertical and a steady oblique path, respectively. In both cases, a strongly inhomogeneous spatial distribution of the disperse phase in the form of columnar clusters is observed, with a significantly enhanced average settling velocity as a consequence. Thus, in contrast to previous results for spheres, the qualitative difference in the single-particle regime does not result in a qualitatively different behaviour of the many-particle cases. In addition, we have carried out an analysis of pairwise interactions of particles in the well-known drafting–kissing–tumbling set-up, for oblate spheroids of aspect ratio $1.5$ and for spheres. We have varied systematically the relative initial position between the particle pair and we have considered free-to-rotate particles and rotationally locked ones. We have found that the region of attraction for both particle shapes, with and without rotation, is very similar. However, significant differences occur during the drafting and tumbling phases. In particular, free-to-rotate spheres present longer drafting phases and separate quickly after the collision. Spheroids remain close to each other for longer times after the collision, and free-to-rotate ones experience two or more collision events. Therefore, we have observed a shape-induced increase in the interaction time which might explain the increased tendency to cluster of the many-particle cases.

Through direct numerical simulations we investigate the effects of spatiotemporal intermittency as a result of stable stratification in surface heated stratified open channel flow. By adapting the density inversion criterion method of Portwood et al. [J. Fluid Mech., vol. 807, 2016, R2] for our flow, we demonstrate that the flow may be robustly separated into regions of active turbulence for which $Re_B \gtrsim {O}(10)$ and surrounding quiescent fluid, where $Re_B$ is the buoyancy Reynolds number. The intermittency in the flow spontaneously manifests as a deformed horizontal interface between the upper quiescent and lower turbulent flow, characterised by vigorous mixing from ‘overturning’ shear instabilities. The resulting vertical intermittency profile is accurately predicted by a local Monin–Obukhov length normalised by viscous wall units $\varLambda ^+$ such that the flow displays intermittency within the parameter range of $2.5 \lesssim \varLambda ^+ \lesssim 260$. By considering conditional averages of the ‘turbulent’ and ‘quiescent’ flow separately, we find the ‘turbulent’ flow within this region to be described by constant critical gradient Richardson and turbulent Froude numbers of $Ri_{g,c} \approx 0.2$ and $Fr_c \approx 0.3$. We find that the turbulent flow continues to display a $\varGamma \sim Fr^{-1}$ relationship when $Fr < Fr_c$, whereas the quiescent flow shows no correlation between $\varGamma$ and $Fr$, where $\varGamma$ is the flux coefficient. Hence, we demonstrate directly that for our flow, the emergence of an asymptotic ‘saturated’ $\varGamma$ regime in the limit of a low ‘global’ $Fr$ occurs due to intermittency and increasing contributions to measurements of $\varGamma$ from the quiescent flow.

Ambient air cushions the impact of drops on solid substrates, an effect usually revealed by the entrainment of a bubble, trapped as the air squeezed under the drop drains and liquid–solid contact occurs. The presence of air becomes evident for impacts on very smooth surfaces, where the gas film can be sustained, allowing drops to bounce without wetting the substrate. In such a non-wetting situation, Mandre & Brenner (J. Fluid Mech., vol. 690, 2012, p. 148) numerically and theoretically evidenced that two physical mechanisms can act to prevent contact: surface tension and nonlinear advection. However, the advection dominated regime has remained hidden in experiments as liquid–solid contact prevents rebounds being realised at sufficiently large impact velocities. By performing impacts on superheated surfaces, in the so-called dynamical Leidenfrost regime (Tran et al., Phys. Rev. Lett., vol. 108, issue 3, 2012, p. 036101), we enable drop rebound at higher impact velocities, allowing us to reveal this regime. Using high-speed total internal reflection, we measure the minimal gas film thickness under impacting drops, and provide evidence for the transition from the surface tension to the nonlinear inertia dominated regime. We rationalise our measurements through scaling relationships derived by coupling the liquid and gas dynamics, in the presence of evaporation.

This investigation characterises the time response and the transient turbulence dynamics undergone by rapidly decelerating turbulent pipe flows. A series of direct numerical simulations of decelerating flows between two steady Reynolds numbers were conducted for this purpose. The statistical analyses reveal that rapidly decelerating turbulent flows undergo four coherent, unambiguous transitional stages: inertial (stage I), a dramatic change of sign in the viscous force associated with the decay of the viscous shear stress at the wall together with a mild turbulence decay in the viscous sublayer; friction recovery (stage II), a recovery in viscous force and progressive decay in the turbulent inertia at the near-wall region; turbulence decay (stage III), a balanced decay in both turbulent inertia and viscous force at the near-wall and overlap regions; core relaxation (stage IV), slow turbulence decay at the core region. The FIK identity derived by Fukagata, Iwamoto and Kasagi (Phys. Fluids, vol. 14, 2002, L73–L76) was used to understand further how the flow dynamics influence the time response of the skin friction coefficient ($C_f$). The results show that although $C_f$ plateaus during the fourth stage, the turbulent contribution keeps decaying, undershoots and finally recovers to attain its final steady value. The time evolution of the azimuthal vorticity ($\omega _\theta$) flux reveals that as the flow is decelerated, a layer of negative $\omega _\theta$ is produced at the wall during the flow excursion. As time progresses, this negative vorticity propagates in the wall-normal direction, attenuating the pre-existing vorticity and producing a decay in the turbulence levels.

A comprehensive investigation, using experimental, computational and analytic methods, is reported on the motion of, and the forces on, spheres of different density ratios rolling freely down an incline in a fluid under gravity. The Reynolds number, based on sphere diameter and terminal velocity, ranged up to 1000 for the experiments, and up to 250 for the computer simulations. A modified Reynolds number, incorporating the density ratio, gravitational acceleration and angle of incline, was found to govern the saturated state of the flow. Transition from steady to unsteady flow was sensitive to mass ratio, with lighter spheres undergoing earlier transition. Indeed, positively buoyant spheres develop cross-slope oscillations prior to the onset of shedding. Also of interest, the transition to chaotic wake flow occurs at Reynolds numbers lower than for a sphere forced to roll at a constant speed. In addition to the average sphere motion, flow-induced vibrations were predicted and measured, with quasi-periodic lateral oscillations found to increase as the flow became more unstable, and to decrease with increased density ratio. The study confirms the time-averaged results of a previous experimental study, although closer inspection shows sensitivity to the relative surface roughness of the sphere and plane in experiments; this sensitivity is masked in typical log–log plots of drag against Reynolds number. Physical surface roughness appears to play a role analogous to the necessary imposed gap between the sphere and plane in computations, removing the singularity in drag that would prevent rolling for an incompressible fluid and perfectly smooth surfaces.

In Rayleigh Bénard convection, for a range of Prandtl numbers $4.69 \leqslant Pr \leqslant 5.88$ and Rayleigh numbers $5.52\times 10^5 \leqslant Ra \leqslant 1.21\times 10^9$, we study the effect of shear by the inherent large-scale flow (LSF) on the local boundary layers on the hot plate. The velocity distribution in a horizontal plane within the boundary layers at each $Ra$, at any instant, is (A) unimodal with a peak at approximately the natural convection boundary layer velocities $V_{bl}$; (B) bimodal with the first peak between $V_{bl}$ and $V_{L}$, the shear velocities created by the LSF close to the plate; or (C) unimodal with the peak at approximately $V_{L}$. Type A distributions occur more at lower $Ra$, while type C occur more at higher $Ra$, with type B occurring more at intermediate $Ra$. We show that the second peak of the bimodal type B distributions, and the peak of the unimodal type C distributions, scale as $V_{L}$ scales with $Ra$. We then show that the areas of such regions that have velocities of the order of $V_{L}$ increase exponentially with increase in $Ra$ and then saturate. The velocities in the remaining regions, which contribute to the first peak of the bimodal type B distributions and the single peak of type A distributions, are also affected by the shear. We show that the Reynolds number based on these velocities scale as $Re_{bs}$, the Reynolds number based on the boundary layer velocities forced externally by the shear due to the LSF, which we obtained as a perturbation solution of the scaling relations derived from integral boundary layer equations. For $Pr=1$ and aspect ratio $\varGamma =1$, $Re_{bs} \sim Ra^{0.375}$ for small shear, similar to the observed flux scaling in a possible ultimate regime. The velocity at the edge of the natural convection boundary layers was found to increase with $Ra$ as $Ra^{0.35}$; since $V_{bl}\sim Ra^{1/3}$, this suggests a possible shear domination of the boundary layers at high $Ra$. The effect of shear, however, decreases with increase in $Pr$ and with increase in $\varGamma$, and becomes negligible for $Pr\geqslant 100$ at $\varGamma =1$ or for $\varGamma \geqslant 20$ at $Pr=1$, causing $Re_{bs}\sim Ra_w^{1/3}$.

An artificial neural-network-based subgrid-scale (SGS) model, which is capable of predicting turbulent flows at untrained Reynolds numbers and on untrained grid resolution is developed. Providing the grid-scale strain-rate tensor alone as an input leads the model to predict a SGS stress tensor that aligns with the strain-rate tensor, and the model performs similarly to the dynamic Smagorinsky model. On the other hand, providing the resolved stress tensor as an input in addition to the strain-rate tensor is found to significantly improve the prediction of the SGS stress and dissipation, and thereby the accuracy and stability of the solution. In an attempt to apply the neural-network-based model trained for turbulent flows with a limited range of the Reynolds number and grid resolution to turbulent flows at untrained conditions on untrained grid resolution, special attention is given to the normalisation of the input and output tensors. It is found that the successful generalization of the model to turbulence for various untrained conditions and resolution is possible if distributions of the normalised inputs and outputs of the neural network remain unchanged as the Reynolds number and grid resolution vary. In a posteriori tests of the forced and the decaying homogeneous isotropic turbulence and turbulent channel flows, the developed neural-network model is found to predict turbulence statistics more accurately, maintain the numerical stability without ad hoc stabilisation such as clipping of the excessive backscatter, and to be computationally more efficient than the algebraic dynamic SGS models.

This work investigates the stability and transition to turbulence in a diverging channel subjected to a time-varying trapezoidal-shaped inflow boundary condition. Numerical simulations are performed for different deceleration rates and Reynolds numbers while maintaining a constant acceleration rate. The flow transition begins with two-dimensional primary instability with the formation of inflectional velocity profiles, followed by local separation and the emergence of an array of shear layer vortices. We divide simulation cases systematically into three categories based on the onset of secondary instability and the generation of streamwise vorticity. At low and medium Reynolds numbers (type I), the spanwise vortex rolls formed by inflectional instability remain two-dimensional and diffuse at the channel centre without exhibiting further instabilities. At high Reynolds numbers and deceleration rates (type II), the rolled shear layer exhibits secondary instability during the zero mean inflow phase, followed by local incipient turbulent structure formation. The streamwise vorticity that develops over the shear layer structures causes oscillations with a spanwise wavelength similar to those associated with the elliptic instability in a counter-rotating vortex pair. Using the Lamb–Oseen approximation of vortices in conjunction with the dynamic mode decomposition algorithm of the three-dimensional flow field, we captured successfully the characteristics of the secondary instability. The third category (type III) is characterized by periodic unsteady separation, secondary instability, and merging of shear layer vortices, which occurs when Reynolds numbers are high and deceleration rates are low.

Vertically neutral collapse of a pulsating bubble occurs when the boundaries above or below the bubble balance the buoyancy effect over a pulsation. In this study, the vertically neutral collapse of a bubble near a vertical rigid wall below the free surface is investigated. The boundary integral method (BIM) is employed to model the bubble dynamics with an open-domain free surface. Moreover, this method is validated against several buoyant bubble experiments. Bubble dynamics in such conditions are associated with three dimensionless parameters: the bubble-free surface distance $\gamma _{{f}}$, bubble–wall distance $\gamma _{{w}}$ and buoyancy parameter $\delta$. We derive the Kelvin impulse of a spherical bubble and the algebraic relationship for vertically neutral collapse, which proves to be accurate for predicting vertically neutral collapse when the bubble is relatively far from the boundaries. Four patterns of the vertically neutral collapse of the bubble for different $\gamma _{{w}}$ and $\gamma _{{f}}$ are identified: (i) formally downward jet; (ii) annular collapse; (iii) horizontal jet; and (iv) weak jet. Despite the downward jet shape, the ‘formally downward jet’ is in the vertically neutral collapse state in terms of the profile of toroidal bubbles and the orientation of local high-pressure zones around the bubble at jet impact. A bulge with a high curvature above the bubble in the ‘annular collapse’ pattern is formed during bubble collapse under two local high-pressure zones at the left and right extremities of the bubble. The ‘horizontal jet’ pattern has the greatest potential to attack the wall, and the power laws of the moment of the jet impact, jet velocity and bubble displacement with respect to the theoretical Kelvin impulse are discussed. In particular, we quantitatively illustrate the role of the free surface on bubble migration towards the wall through the variational power-law exponents of the bubble displacement with respect to $\gamma _{{w}}$.

The tidal flow through a channel connecting two basins with different tidal regimes can be optimally controlled by means of a turbine fence or array to maximise the extracted mechanical power. The paper gives the optimal control strategy as a function of the blockage ratio $\sigma$, i.e. the ratio of the turbine cross-section to the cross-section of the local passage of a turbine. The results presented are a physically consistent generalisation of the results of Garrett & Cummins (Proc. R. Soc. Lond. A, vol. 461, 2060, pp, 2563–2572), valid only for $\sigma =1$ and turbine efficiency of one, now for arbitrary blockage ratio $0 < \sigma \leqslant 1$. Published research over the past decade on the same topic has taken the momentum equation and the turbine drag force as a starting point. The new approach presented here, in contrast, takes the energy equation as the starting point and uses the relative volume flow as the control variable. As the work shows, this new approach has three advantages. First, starting with the energy equation allows us to derive an optimal flow control problem resulting in an Euler–Lagrange equation using the physically consistent and experimentally validated actuator disk model for the free surface flow of Pelz et al. (J. Fluid Mech., vol. 889, 2020) in a direct and formal way. The optimal control problem is solved (a) numerically and (b) analytically. In the latter case, the turbine characteristics are approximated by a rational function in the relevant design and operating range. The analytical solution (b) validated against the numerical solution (a) is surprisingly concise and easy to apply in practice, as shown by use cases. Second, instead of the induction factor, we use the volume flow that is the same for all turbines in a cascade, i.e. a row of turbines in the direction of flow, which significantly reduces the complexity of the optimal control task of turbine arrays. Third, we obtain a well-founded energy estimate, whereas previous methods overestimate the energy yield due to inconsistent turbine disc models (for the consistency and valid parameter ranges of different models, also in comparison with experiments, see Pelz et al., J. Fluid Mech., vol. 889, 2020). The results can be used for the conceptual design of turbine arrays, but also for a sound physically realistic and consistent resource assessment of tidal power for a system consisting of two basins, a channel and a turbine fence with $0<\sigma \leqslant 1$ and operated in a complete tidal cycle.

This paper proposes a new approach for evaluating numerically the forces and moments applied to a circular cylinder that is immersed in a fluid and which translates and rotates near a plane wall. Under the proposed approach, the flow is decomposed into inner and outer flows. The inner flow represents the flow in the thin interstice between the cylinder and the wall, and is obtained as an analytic expression using lubrication theory. The outer flow represents the flow far from the interstice, which does not depend on the magnitude of the gap between the cylinder and the wall, when the gap is small. The outer flow is obtained using numerical simulation as a function of both the Reynolds number and the slip coefficient. The force and moment coefficients are then obtained, as functions of the Reynolds number, slip coefficient and gap-to-diameter ratio, by combining the inner and outer solutions. Importantly, since the outer flow does not depend on the gap-to-diameter ratio, the parameter space to be explored by numerical simulations is greatly reduced compared to using finite gap ratio simulations. Moreover, the numerical difficulties associated with resolving the interstitial flow are avoided. The proposed approach can be extended to a wide range of rolling bodies, including spherical particles and wheels, and should significantly reduce the computational expense required to model the hydrodynamic forces and predict the subsequent motion of such bodies.

Surface buoyancy gradients over a quasigeostrophic fluid permit the existence of surface-trapped Rossby waves. The interplay of these Rossby waves with surface quasigeostrophic turbulence results in latitudinally inhomogeneous mixing that, under certain conditions, culminates in a surface buoyancy staircase: a meridional buoyancy profile consisting of mixed zones punctuated by sharp buoyancy gradients, with eastward jets centred at the sharp gradients and weaker westward flows in between. In this article, we investigate the emergence of this buoyancy staircase limit in surface quasigeostrophic turbulence and we examine the dependence of the resulting dynamics on the vertical stratification. Over decreasing stratification ($\mathrm {d}N(z)/\mathrm {d}z\leq 0$, where $N(z)$ is the buoyancy frequency), we obtain flows with a longer interaction range (than in uniform stratification) and highly dispersive Rossby waves. In the staircase limit, we find straight jets that are perturbed by eastward propagating along jet waves, similar to two-dimensional barotropic $\beta$-plane turbulence. In contrast, over increasing stratification ($\mathrm {d}N(z)/\mathrm {d}z\geq 0$), we obtain flows with shorter interaction range and weakly dispersive Rossby waves. In the staircase limit, we find sinuous jets with large latitudinal meanders whose shape evolves in time due to the westward propagation of weakly dispersive along jet waves. These along jet waves have larger amplitudes over increasing stratification than over decreasing stratification, and, as a result, the ratio of domain-averaged zonal to meridional speeds is two to three times smaller over increasing stratification than over decreasing stratification. Finally, we find that, for a given Rhines wavenumber, jets over increasing stratification are closer together than jets over decreasing stratification.

We investigate the Lagrangian properties of fingering convection in a parameter regime that is far beyond marginality, but that does not lead to staircase formation. We find that most fluid particles perform very short vertical excursions before switching direction. Only a small fraction of them travels vertically in the same direction for lengths substantially larger than the typical spatial scales associated with fingering instabilities. These long travels are associated with vertical velocities, buoyancy anomalies and advective fluxes constrained within a relatively narrow range of values. Despite their limited occurrence, they contribute the overwhelming majority of the up-gradient buoyancy flux. A description of fingering convection emerges that partitions the fluid into an incoherent background, where fluxes may be down-gradient, and a set of relatively long-lived, coherent, buoyancy-carrying structures that travel vertically, interacting with each other.

Reconstructing high-resolution flow fields from sparse measurements is a major challenge in fluid dynamics. Existing methods often vectorize the flow by stacking different spatial directions on top of each other, hence confounding the information encoded in different dimensions. Here, we introduce a tensor-based sensor placement and flow reconstruction method which retains and exploits the inherent multidimensionality of the flow. We derive estimates for the flow reconstruction error, storage requirements and computational cost of our method. We show, with examples, that our tensor-based method is significantly more accurate than similar vectorized methods. Furthermore, the variance of the error is smaller when using our tensor-based method. While the computational cost of our method is comparable to similar vectorized methods, it reduces the storage cost by several orders of magnitude. The reduced storage cost becomes even more pronounced as the dimension of the flow increases. We demonstrate the efficacy of our method on three examples: a chaotic Kolmogorov flow, in situ and satellite measurements of the global sea surface temperature and three-dimensional unsteady simulated flow around a marine research vessel.

By using direct numerical simulations (DNS) of bubbly flows with passive scalars, we show a transition in the scalar spectra from a $k^{-5/3}$ to a $k^{-3}$ scaling with the wavenumber $k$, in contrast with those of single-phase isotropic turbulence. For cases with a mean scalar gradient in the horizontal direction, the scalar spectrum decays faster than $k^{-3}$ at high wavenumbers. While the $k^{-3}$ scaling is well established in the bubbly flow velocity spectrum, the scalar spectrum behaviour is not fully understood. We find that the transition length scale of the scalar spectra is comparable to or below the bubble diameter and decreases with the molecular diffusivity of the scalar in the liquid. We use DNS to compute the scalar spectrum budget and show that the scalar fluctuations are produced by the mean scalar gradient at length scales above the bubble diameter, contrary to the velocity fluctuations. At length scales below the bubble diameter, the net scalar transfer scales as $k^{-1}$ inducing the $k^{-3}$ scaling of the scalar spectra. This finding is consistent with the hypothesis proposed by Dung et al. (J. Fluid Mech., vol. 958, 2023, p. A5) about the physical mechanism behind the $k^{-3}$ scaling. We also show dependencies of the bubble suspension's convective scalar diffusivity on the gas volume fraction and molecular diffusivity that differ based on the direction of the mean scalar gradient. For a mean scalar gradient in the vertical direction, we find and qualitatively explain a significant effect of the molecular diffusivity in the gas on the convective scalar diffusivity.

We report flow measurements in rotating Rayleigh–Bénard convection in the rotationally constrained geostrophic regime. We apply stereoscopic particle image velocimetry to measure the three components of velocity in a horizontal cross-section of a water-filled cylindrical convection vessel. At a constant, small Ekman number $Ek=5\times 10^{-8}$, we vary the Rayleigh number $Ra$ between $10^{11}$ and $4\times 10^{12}$ to cover various subregimes observed in geostrophic convection. We also include one non-rotating experiment. The scaling of the velocity fluctuations (expressed as the Reynolds number $Re$) is compared to theoretical relations expressing balances of viscous–Archimedean–Coriolis (VAC) and Coriolis–inertial–Archimedean (CIA) forces. Based on our results we cannot decide which balance is most applicable here; both scaling relations match equally well. A comparison of the current data with several other literature datasets indicates a convergence towards diffusion-free scaling of velocity as $Ek$ decreases. However, at lower $Ra$, the use of confined domains leads to prominent convection in the wall mode near the sidewall. Kinetic energy spectra point at an overall flow organisation into a quadrupolar vortex filling the cross-section. This quadrupolar vortex is a quasi-two-dimensional feature; it manifests only in energy spectra based on the horizontal velocity components. At larger $Ra$, the spectra reveal the development of a scaling range with exponent close to $-5/3$, the classical exponent for inertial range scaling in three-dimensional turbulence. The steeper $Re(Ra)$ scaling at low $Ek$ and development of a scaling range in the energy spectra are distinct indicators that a fully developed, diffusion-free turbulent bulk flow state is approached, sketching clear perspectives for further investigation.

Hydraulic fracturing for oil and gas production from shale formations, as well as natural geological phenomena, involve the propagation of thin viscous films within elastic media. For viscous fluids, stress diverges as the thickness of the film tends to zero, arresting the propagation of the film, and thus implying the contact line paradox. For free-surface films, this paradox is resolved by considering a precursor film, leading to Tanner's law. This approach was extended recently for viscous films between a thin elastic plate and a rigid solid, allowing calculation of the film propagation rate. In this work, we examine the effect of a pre-wetting layer on the rate of propagation of a viscous flow within an infinitely deep and long domain. We analyse the linear and nonlinear dynamic problems, and perform a self-similarity analysis. We find that peeling front propagation scales as time to the power of $1/9$ and $1/3$ for thin and thick pre-wetting layer limits, respectively. Our results contribute to the understanding of the contact line paradox in elastic media and the crucial role of the pre-wetting layer in resolving it.

With an immersed-boundary lattice-Boltzmann method, we consider the transit of a three-dimensional initially spherical capsule with a viscoelastic membrane through a cross-slot microchannel. The capsule is released with a small initial off-centre distance in the feeding channel, to mimic experiments where capsules or cells are not perfectly aligned with the centreline. Our main objective is to establish the phase diagram of the capsule's deformation modes as a function of the flow inertia and capsule membrane viscosity. We mainly find three deformation modes in the channel cross-slot. For a capsule with low membrane viscosity, a quasi-steady mode occurs at low Reynolds numbers ($Re$), in which the capsule can reach and maintain a steady ellipsoidal shape near the stagnation point, for a considerable time period. With $Re$ increasing to 20, an overshoot–retract mode is observed. The capsule deformation oscillates on an inertial–elastic time scale, suggesting that the dynamics is mainly driven by the balance of the inertial and membrane elastic forces. The membrane viscosity slows down the capsule deformation and suppresses the overshoot–retract mode. A capsule with high membrane viscosity undergoes a continuous-elongation mode, in which its deformation keeps increasing during most of its journey in the channel cross-slot. We summarise the results in phase diagrams, and propose a scaling model which can predict the deformation modes of a viscoelastic capsule in the inertial flow regime. We also discuss implications of the present findings for practical experiments for mechanical characterisation of capsules or cells.

To investigate the initial development of the Rayleigh–Bénard–Marangoni (RBM) instability in a relatively deep domain, direct numerical simulations for a large range of Marangoni and Rayleigh numbers were performed. In the simulations, the surface was assumed to be flat and surface cooling was modelled by a constant heat flux. The small-scale dynamics of the flow and temperature fields near the surface was fully resolved by using a non-uniform vertical grid distribution. A detailed investigation of the differences in physical mechanisms that drive the Rayleigh- and Marangoni-dominated instabilities is presented. To this end, various properties such as the maturation rate of convection cells, the fluctuating kinetic energy and the surface characteristic length scale were studied. It was confirmed that buoyancy forces and surface-temperature-gradient-driven Marangoni forces enhance one another in promoting the development of the RBM instability. When using a relevant measure of the effective thermal boundary layer thickness as length scale, both the critical Marangoni and Rayleigh numbers, obtained for the purely Marangoni- and purely Rayleigh-driven instabilities, were found to be in good agreement with the literature.

Turbulence in the restricted nonlinear (RNL) dynamics is analysed and compared with direct numerical simulations (DNS) of Poiseuille turbulence at Reynolds number $R=1650$. The structures are obtained by proper orthogonal decomposition (POD) analysis of the two components of the flow partition used in RNL dynamics: the streamwise mean flow and fluctuations. POD analysis of the streamwise mean flow indicates that the dominant POD modes, in both DNS and RNL dynamics, are roll-streaks harmonic in the spanwise direction. However, we conclude that these POD modes do not occur in isolation but rather are Fourier components of a coherent roll-streak structure. POD analysis of the fluctuations in DNS and RNL dynamics reveals similar complex structures consisting in part of oblique waves collocated with the streak. The origin of these structures is identified by their correspondence to POD modes predicted using a stochastic turbulence model (STM). These predicted POD modes are dominated by the optimally growing structures on the streak, which the STM predicts correctly to be of sinuous oblique wave structure. This close correspondence between the roll-streak structure and the associated fluctuations in DNS, RNL dynamics and the STM implies that the self-sustaining mechanism operating in DNS is essentially the same as that in RNL dynamics, which has been associated previously with optimal perturbation growth on the streak.