The aeroacoustic characteristics of flying vehicles with pitch-fixed rotors differ from traditional helicopters with pitch-controlled rotor blades. Accurate predictions of rotor noise are still challenging because many uncertainty factors and unsteadinesses exist. This work investigates the aeroacoustic effects of rotational speed deviation, rotation speed fluctuation, blade vibration and blade geometric asymmetry. The analysis is based on the efficient computation of rotor noise under different working conditions. The mean aerodynamic variables are computed using the blade element moment theory, while small-amplitude fluctuations are introduced to account for the unsteadiness and uncertainty factors. It is shown that periodic rotation speed fluctuations and blade vibrations can produce significant extra tones. By contrast, if the fluctuations and vibrations are random, the noise level in a wide frequency range is increased. The intriguing result reminds us of the need to revisit the rotor broadband noise sources commonly attributed to turbulent flows. The influences are observer angle dependent, and the extra noise production is more significant in the upstream and downstream directions. The asymmetric blade geometry can cause extra tonal noise at the harmonics of the blade shaft frequency. The noise features of dual rotors are also investigated. Usually, the noise is sensitive to the initial phase difference and rotation directions due to the interference effect. However, the noise features are vastly altered if there are slight differences in the rotation speeds. Although the influences of some factors on rotor noise were already known, the present study provides a more comprehensive analysis of the problem. The results also highlight the need to consider these practical factors for accurate noise prediction of multi-rotor flying vehicles.

Immiscible and incompressible liquid–liquid flows are considered in a Taylor–Couette geometry and analysed by direct numerical simulations coupled with the volume-of-fluid method and a continuum surface force model. The system Reynolds number $Re \equiv r_i \omega _i d / \nu$ is fixed to $960$, where the single-phase flow is in the steady Taylor vortex regime, whereas the secondary-phase volume fraction $\varphi$ and the system Weber number $We \equiv \rho r_i^2 \omega _i^2 d / \sigma$ are varied to study the interactions between the interface and the Taylor vortices. We show that different Weber numbers lead to two distinctive flow regimes, namely an advection-dominated regime and an interface-dominated regime. When $We$ is high, the interface is easily deformed because of its low surface tension. The flow patterns are then similar to the single-phase flow, and the system is dominated mainly by advection (advection-dominated regime). However, when $We$ is low, the surface tension is so large that stable interfacial structures with sizes comparable to the cylinder gap can exist. The background velocity field is modulated largely by these persistent structures, thus the overall flow dynamics is governed by the interface (interface-dominated regime). The effect of the interface on the global system response is assessed by evaluating the Nusselt number $Nu_{\omega }$ based on the non-dimensional angular velocity transport. It shows non-monotonic trends as functions of the volume fraction $\varphi$ for both low and high $We$. We explain how these dependencies are closely linked to the velocity and interfacial structures.

The dynamics of an evaporating droplet in an unsteady flow is of practical interest in many industrial applications and natural processes. To investigate the transport and evaporation dynamics of such droplets, we present a numerical study of an isolated droplet in an oscillating gas-phase flow. The study uses a one-way coupled two-phase flow model to assess the effect of the amplitude and the frequency of a sinusoidal external flow field on the lifetime of a multicomponent droplet containing a non-volatile solute dissolved in a volatile solvent. The results show that the evaporation process becomes faster with an increase in the amplitude or the frequency of the gas-phase oscillation. The liquid-phase transport inside the droplet also is influenced by the unsteadiness of the external gas-phase flow. A scaling analysis based on the response of the droplet under the oscillating drag force is subsequently carried out to unify the observed evaporation dynamics in the simulations under various conditions. The analysis quantifies the enhancement in the droplet velocity and Reynolds number as a function of the gas-phase oscillation parameters and predicts the effects on the evaporation rate.

We experimentally investigate the evaporation of very volatile liquid droplets (Novec 7000 Engineered Fluid, chemical name hydrofluoroethers HFE-7000) in a turbulent spray. Droplets with diameters of the order of a few micrometres are produced by a spray nozzle and then injected into a purpose-built enclosed dodecahedral chamber filled with air containing various amounts of water vapour. The ambient temperature and relative humidity in the chamber are carefully controlled. We observe water condensation on the rapidly evaporating droplet, both for the spray and for a single acoustically levitated millimetric Novec 7000 droplet. We further examine the effect of humidity, and reveal that a more humid environment leads to faster evaporation of the volatile liquid, as well as more water condensation. This is explained by the much larger latent heat of water as compared with that of Novec 7000. We extend an analytical model based on Fick's law to quantitatively account for the data.

Well-resolved direct numerical simulations (DNS) have been performed of the flow in a smooth circular pipe of radius $R$ and axial length $10{\rm \pi} R$ at friction Reynolds numbers up to $Re_\tau =5200$ using the pseudo-spectral code OPENPIPEFLOW. Various turbulence statistics are documented and compared with other DNS and experimental data in pipes as well as channels. Small but distinct differences between various datasets are identified. The friction factor $\lambda$ overshoots by $2\,\%$ and undershoots by $0.6\,\%$ the Prandtl friction law at low and high $Re$ ranges, respectively. In addition, $\lambda$ in our results is slightly higher than in Pirozzoli et al. (J. Fluid Mech., vol. 926, 2021, A28), but matches well the experiments in Furuichi et al. (Phys. Fluids, vol. 27, issue 9, 2015, 095108). The log-law indicator function, which is nearly indistinguishable between pipe and channel up to $y^+=250$, has not yet developed a plateau farther away from the wall in the pipes even for the $Re_\tau =5200$ cases. The wall shear stress fluctuations and the inner peak of the axial turbulence intensity – which grow monotonically with $Re_\tau$ – are lower in the pipe than in the channel, but the difference decreases with increasing $Re_\tau$. While the wall value is slightly lower in the channel than in the pipe at the same $Re_\tau$, the inner peak of the pressure fluctuation shows negligible differences between them. The Reynolds number scaling of all these quantities agrees with both the logarithmic and defect-power laws if the coefficients are properly chosen. The one-dimensional spectrum of the axial velocity fluctuation exhibits a $k^{-1}$ dependence at an intermediate distance from the wall – also seen in the channel. In summary, these high-fidelity data enable us to provide better insights into the flow physics in the pipes as well as the similarity/difference among different types of wall turbulence.

A new time-dependent analysis of the global and local fluctuating velocity signals in grid turbulence is conducted to assess the scaling laws for non-equilibrium turbulence. Experimental datasets of static- and active-grid turbulence with different Rossby numbers $R_o({=}U/\varOmega M$: $U$ is the mean velocity, $\varOmega$ is the mean rotation rate and $M$ is the grid mesh size) are considered. Although the global (long-time-averaged) non-dimensional dissipation rate $C_\varepsilon$ is independent of the Reynolds number $Re_\lambda$ based on the global Taylor microscale, the local (short-time-averaged) non-dimensional dissipation rate $\left \langle C_\varepsilon (t_i) \right \rangle$ ($t_i$ is the local time) both in the static- and active-grid turbulence clearly show the non-equilibrium scaling $\left \langle C_\varepsilon (t_i)\right \rangle / \sqrt {Re_0} \propto \left \langle Re_\lambda (t_i) \right \rangle ^{-1}$ ($\left \langle Re_\lambda (t_i) \right \rangle$ and $Re_0$ are the Reynolds numbers based on the local Taylor microscale $\lambda (t_i)$ and the global integral length scale, respectively), which has only been confirmed for global statistics in the near field of grid turbulence. The local value of $\left \langle L(t_i) / \lambda (t_i) \right \rangle$ ($L(t_i)$ is the local integral length scale) shifts from the equilibrium to non-equilibrium scaling as $\left \langle Re_\lambda (t_i) \right \rangle$ increases, further confirming that the non-equilibrium scalings are recovered for local statistics both in the static- and active-grid turbulence. The local values of $\left \langle C_\varepsilon (t_i) \right \rangle$ and $\left \langle L(t_i) / \lambda (t_i) \right \rangle$ follow the theoretical predictions for global statistics (Bos & Rubinstein, Phys. Rev. Fluids, vol. 2, 2017, 022601).

Understanding and predicting the dynamics of dispersed micro-objects in microfluidics is crucial in numerous natural, industrial and technological situations. In this paper, we experimentally characterized the equilibrium velocity $V$ and lateral position $\varepsilon$ of various dispersed micro-objects, such as beads, bubbles and drops, in a cylindrical microchannel over an unprecedentedly wide range of parameters. By varying the dimensionless object size ($d \in [0.1; 1]$), the viscosity ratio ($\lambda \in [10^{-2}; \infty [$), the density ratio ($\varphi \in [10^{-3}; 2]$), the Reynolds number ($Re \in [10^{-2}; 10^2]$) and the capillary number ($Ca \in [10^{-3}; 0.3]$), we offer an exhaustive parametric study exploring various dynamics from the non-deformable viscous regime to the deformable inertial regime, thus enabling us to highlight the sole and combined roles of inertia and capillary effects on lateral migration. Experiments are compared and agree well with a steady three-dimensional Navier–Stokes model for incompressible two-phase fluids, including the effects of inertia and possible interfacial deformations. This model enables us to propose a correlation for the object velocity $V$ as functions of $d$, $\varepsilon$ and $\lambda$, obtained in the ${Re}={Ca}=0$ limit, but valid for a larger range of values of ${Re}$ and ${Ca}$ delimited by the validity of the linear regime. Next, we present stability maps for the centred position showing that non-deformable objects dominated by inertial effects are only stable if large enough, typically for $d \gtrsim 0.7$, whereas deformable objects dominated by capillary effects can be stable for much smaller sizes, provided the viscosity ratio is outside the range $0.7 \lesssim \lambda \lesssim 10$, in which deformability also plays a destabilizing effect, as for inertia.

When two drops collide, they may either exhibit complete coalescence or selectively generate secondary drops, depending on their relative sizes and physical properties, as dictated by a decisive interplay of the viscous, capillary, inertia and gravity effects. Electric field, however, is known to induce distinctive alterations in the topological evolution of the interfaces post-collision, by influencing a two-way nonlinear coupling between electro-mechanics and fluid flow as mediated by a topologically intriguing interfacial deformation. While prior studies primarily focused on the viscous-dominated regime of the resulting electro-coalescence dynamics, several non-intuitive features of the underlying morpho-dynamic evolution over the intertio-capillary regime have thus far remained unaddressed. In this study, we computationally investigate electrically modulated coalescence dynamics along with secondary drop formation mechanisms in the inertio-capillary regime, probing the interactions of two unequal-sized drops subjected to a uniform electric field. Our results bring out an explicit mapping between the observed topological evolution as a function of the respective initial sizes of the parent drops as well as their pertinent electro-physical property ratios. These findings establish electric-field-mediated exclusive controllability of the observed topological features, as well as the critical conditions leading to the transition from partial to complete coalescence phenomena. In a coalescence cascade, an electric field is further shown to orchestrate the numbers of successive stages of coalescence before complete collapse. However, an increase of the numbers of cascade stages with the electric field strength and parent droplet size ratio is non-perpetual, and the same is demonstrated to continue until only a threshold number of cascade stages is reached. These illustrations offer significant insights into leveraging the interplay of electrical, inertial and capillary-driven interactions for controllable drop manipulation via multi-drop interactions for a variety of applications ranging from chemical processing to emulsion technology.

We numerically study the melting process of a solid layer heated from below such that a liquid melt layer develops underneath. The objective is to quantitatively describe and understand the emerging topography of the structures (characterized by the amplitude and wavelength), which evolve out of an initially smooth surface. By performing both two-dimensional (achieving Rayleigh number up to $Ra=10^{11}$) and three-dimensional (achieving Rayleigh number up to $Ra=10^9$) direct numerical simulations with an advanced finite difference solver coupled to the phase-field method, we show how the interface roughness is spontaneously generated by thermal convection. With increasing height of the melt the convective flow intensifies and eventually even becomes turbulent. As a consequence, the interface becomes rougher but is still coupled to the flow topology. The emerging structure of the interface coincides with the regions of rising hot plumes and descending cold plumes. We find that the roughness amplitude scales with the mean height of the liquid layer. We derive this scaling relation from the Stefan boundary condition and relate it to the non-uniform distribution of heat flux at the interface. For two-dimensional cases, we further quantify the horizontal length scale of the morphology, based on the theoretical upper and lower bounds given for the size of convective cells known from Wang et al. (Phys. Rev. Lett., vol. 125, 2020, 074501). These bounds agree with our simulation results. Our findings highlight the key connection between the morphology of the melting solid and the convective flow structures in the melt beneath.

We experimentally and numerically study the dynamics of a liquid jet issued from a rotating orifice, whose breakup is regulated by a vibrating piezo element. The helical trajectory of the spiralling jet yields fictitious forces varying along the jet whose longitudinal projections stretch and thin the jet, affecting the growth of perturbations. We show that by quantifying these fictitious forces, one can estimate the jet intact length and size distribution of drops formed at jet breakup. The presence of the locally varying fictitious forces may render high-frequency perturbations, that would otherwise be stable in the abscence of stretching, unstable, as observed similarly in the case of straight jets stretching under gravity. The perturbation amplitude then dictates how strong the perturbation is coupled to the jet compared with random noise that is inherently present in any experimental set-up. In the present study we exploit the slenderness of the jet to separate the calculation of the base flow and the growth of perturbations. The fictitious forces calculated from the base flow trajectory are then used in a nonlinear slender-jet model, which treats the spiralling jet as a quasi-straight jet with locally varying body forces. We show both experimentally and numerically that jet breakup characteristics (e.g. intact length and drop size distribution) can be controlled by finite-amplitude perturbations created by mechanically induced pressure modulations. Finally, we revisit the integrated net gain approach developed for straight jets under gravity and we provide simple analogous relations for spiralling jets.

Despite significant advances in the field of man-made micro- and nanomotors, it remains a challenge to precisely control their motion in bounded environments. Here, we present a theoretical analysis of a thermally activated micromotor near a plane wall under the action of a background linear temperature field. The coupling between the autonomous and field-directed motions has been resolved using a combined analytical–numerical framework comprising general solutions in bispherical coordinates and the reciprocal theorem for creeping flows. Results reveal giant augmentation in swimming speeds, the controlling parameter zones for positive and negative thermotaxes and the flexibility of steering perpendicular to the field gradient for an isolated micromotor. Boundary-instigated thermo-fluidic modulations at different levels of confinements and preferential orientations cause directional switching of both the vertical translation and rotation parallel to the wall, thereby drastically altering the phase portraits of the swimmer dynamics. Contrasting trajectory characteristics, e.g. escape, attraction, are partitioned by unstable separatrices in the phase portraits, while competitive repulsion (attraction) after attraction (repulsion) characteristics emerge for different relative field strengths $\mathcal {S}$ and gradient orientations $\theta _T$. Below $\mathcal {S}=0.25$, highly counter-intuitive trajectories result when the micromotor is initially launched from an overlapping escape zone. Moreover, external-field-assisted microswimming can uniquely tune the directionality of wall-parallel translation, broadening the scope of dynamic regulation of self-propulsion. Thus, providing insights into a precisely controlled fuel-free actuation of micromotors near a physical obstacle, the present study stands as a step toward addressing the increasing demand for successful implementation of micromotors in futuristic clinical and environmental applications.

In this study, we revisit the spectral transfer model for the turbulent intensity in passive scalar transport (under large-scale anisotropic forcing), and a subsequent modification to the scaling of scalar variance cascade is presented. From the modified spectral transfer model, we obtain a revised scalar transport model using a fractional-order Laplacian operator that facilitates the robust inclusion of the non-local effects originating from large-scale anisotropy transferred across the multitude of scales in the turbulent cascade. We provide an a priori estimate for the non-local model based on the scaling analysis of the scalar spectrum, and later examine our developed model through direct numerical simulation. We present a detailed analysis on the evolution of the scalar variance, high-order statistics of the scalar gradient and important two-point statistical metrics of the turbulent transport to make a comprehensive comparison between the non-local model and its standard version. Finally, we present an analysis that seamlessly reconciles the similarities between the developed model with the fractional-order subgrid-scale scalar flux model for large-eddy simulation (Akhavan-Safaei et al., J. Comput. Phys., vol. 446, 2021, 110571) when the filter scale approaches the dissipative scales of turbulent transport. In order to perform this task, we employ a Gaussian process regression model to predict the model coefficient for the fractional-order subgrid model.

The heterogeneity of permeability in a porous medium considerably alters the behaviour of density-driven flows from what is observed in a medium of homogeneous permeability, and significantly enhances the mixing between the dense and light fluids during the flow. In this work, we present results from laboratory experiments performed in heterogeneous media consisting of horizontal layers of different permeabilities, investigating their effects on gravity current flows. We find that the mixing in our heterogeneous experimental set-ups can be ${O}(2)$ greater than that in a homogeneous medium of similar depth-averaged properties. The enhanced mixing in this setting is primarily because of transverse gravity-driven fingers and produced blunt front, which is the direct result of the layered structure. This enhanced rate of mixing dictates the gravity current height and length, making the current lose its long and thin shape much faster than a comparable current in a homogeneous medium. We discuss the experimental observations in detail and present relevant physical interpretations. Based on the experimental measurements and dimensionless modelling, we also derive semi-empirical formulas for predicting the gravity current length, height and mixing in a heterogeneous medium. Results from this work can be used in predicting the scale of mixing between different density fluids during contaminant transport in the subsurface environment.

The dynamics of sedimenting particles under gravity are surprisingly complex due to the presence of effective long-ranged forces. When the particles are polar with a well-defined symmetry axis and non-uniform density, recent theoretical predictions suggest that prolate objects will repel and oblate ones will weakly attract. We tested these predictions using mass polar prolate spheroids, which are composed of 2 mm spheres glued together. We probed different aspect ratios ($\kappa$) and centre of mass variations ($\chi$) by combining spheres of different densities. Experiments were done in both quasi-two-dimensional (2-D) and three-dimensional (3-D) chambers. By optically tracking the motion of single particles, we found that the dynamics were well described by a reduced mobility matrix model that could be solved analytically. Pairs of particles exhibited an effective repulsion, and their separation roughly scaled as $(\kappa - 1)/\chi ^{0.39}$, i.e. particles that were more prolate or had smaller mass asymmetry had stronger repulsion effects. In three dimensions, particles with $\chi >0$ were distributed more uniformly than $\chi =0$ particles, and the degree of uniformity increased with $\kappa$, indicating that the effective 2-body repulsion manifests for a large number of particles.

We investigate the Benjamin–Feir (or modulational) instability of Stokes waves, i.e. small-amplitude, one-dimensional periodic gravity waves of permanent form and constant velocity, in water of finite and infinite depth. We develop a perturbation method to describe to high-order accuracy the unstable spectral elements associated with this instability, obtained by linearizing Euler's equations about the small-amplitude Stokes waves. These unstable elements form a figure-eight curve centred at the origin of the complex spectral plane, which is parametrized by a Floquet exponent. Our asymptotic expansions of this figure-eight are in excellent agreement with numerical computations as well as recent rigorous results by Berti et al. (Full description of Benjamin–Feir instability of Stokes waves in deep water, 2021, arXiv:2109.11852) and Berti et al. (Benjamin–Feir instability of Stokes waves in finite depth, 2022, arXiv:2204.00809). From our expansions, we derive high-order estimates for the growth rates of the Benjamin–Feir instability and for the parametrization of the Benjamin–Feir figure-eight curve with respect to the Floquet exponent. We are also able to compare the Benjamin–Feir and high-frequency instability spectra analytically for the first time, revealing three different regimes of the Stokes waves, depending on the predominant instability.

The exhaust from combustion engines contains particulate matter (PM), which poses potential health risks to human lungs. Current emission laws place increasingly strict limitations on both PM and particle number, leading to the necessity of using wall-flow filters to separate out a significant amount of the introduced PM. As this leads to an increase in the filter's loading, it is regenerated continuously or periodically, leading to the rearrangement of individual particulate structures inside the filter channels. Such rearrangement events cause the formation of specific deposition patterns, which affect the filter's pressure drop, its loading capacity and the separation efficiency. In order to derive predictions on the formation of specific deposition patterns, the transient behaviour of individual particle structures needs to be examined. The present work investigates the detachment and transport of particle structures during filter regeneration with three-dimensional surface-resolved simulations using a lattice Boltzmann method. The goal of this work is the determination of relevant key quantities and their interpretation with respect to predictions regarding the resulting deposition patterns. In this context, it is shown that lift forces are not the predominant detachment forces for non-spherical particle structures, and that the stopping distance of such structures is too long to avoid back-end deposition.

Composite sweeping-enhanced resolvents, referred to as the ${\boldsymbol {R}}_s^2$ model, are proposed to predict the space–time statistics of large-scale structures in turbulent channel flows. This model incorporates two key mechanisms: (i) eddy damping is introduced to represent random sweeping decorrelation caused by nonlinear forcing, leading to a sweeping-enhanced resolvent ${{\boldsymbol {R}}_s}$; and (ii) the sweeping-enhanced resolvent ${{\boldsymbol {R}}_s}$ is composited into its iterations ${\boldsymbol {R}}_s^2$ to yield non-zero Taylor time microscales. The resulting ${\boldsymbol {R}}_s^2$ model can correctly predict the frequency spectra and two-point cross-spectra of large-scale structures. This model is compared numerically with eddy-viscosity-enhanced resolvent models. The latter are designed to represent energy transfer instead for time decorrelation, and thus underpredict the characteristic decay time scales. The ${\boldsymbol {R}}_s^2$ model correctly yields the characteristic decay time scales in turbulent channel flows.

Underwater stability of an air layer trapped in a micro-structure, plastron, is critical in drag reduction applications. Here, we investigate the wetting state and plastron stability of underwater superhydrophobic surfaces (SHS) under an intense acoustic drive. Flat surfaces and SHS are subjected to short acoustic pulses of different intensities. At low amplitude, the comparison between the results of various surfaces shows that plastron behaves like a water–air interface, whose presence can be detected from the phase of the reflected acoustic waves. At moderate intensity, a wetting transition towards a completely wetting state is observed and shown to be triggered by a sufficiently large acoustic radiation pressure. This wetting transition is well captured by a simplified model by balancing radiation pressure with the critical capillary pressure for the interface sliding. Cavitation clouds appear under strong excitation; their sizes and positions greatly depend on the surface acoustic boundary condition. For SHS in a Cassie–Baxter state (with an air layer), cavitation clouds appear at specific locations (from the solid surface) corresponding to the pressure anti-node of the transient standing wave generated by the reflection. This study unprecedentedly demonstrates the capability of acoustic waves to monitor and characterize plastron stability with low and moderate amplitudes, respectively.

Experiments were conducted to study the transition and flow development in a laminar separation bubble (LSB) formed on an aerofoil. The effects of a wide range of free-stream turbulence intensity ($0.15\,\%< Tu<6.26\,\%$) and streamwise integral length scale ($4.6\ {\rm mm}<\varLambda _{u}<17.2\ {\rm mm}$) are considered. The co-existence of modal instability due to the LSB and non-modal instability caused by streaks generated by free-stream turbulence is observed. The flow field is measured using hot-wire anemometry, which showed that the presence of streaks in the boundary layer modifies the mean-flow topology of the bubble. These changes in the mean flow field result in the modification of the convective disturbance growth, where an increase in turbulence intensity is found to dampen the growth of the modal instability. For a relatively fixed level of $Tu$, the variation of $\varLambda _{u}$ has modest effects. However, a slight advancement of the nonlinear growth of disturbances and eventual breakdown with the decrease in $\varLambda _{u}$ is observed. The data show that the streamwise growth of the disturbance energy is exponential for the lowest levels of free-stream turbulence and gradually becomes algebraic as the level of free-stream turbulence increases. Once a critical turbulence intensity is reached, there is enough energy in the boundary layer to suppress the laminar separation bubble, resulting in the non-modal instability taking over the transition process. Linear stability analysis is conducted in the fore position of the LSB. It accurately models incipient disturbance growth, unstable frequencies and eigenfunctions for configurations subjected to turbulence intensity levels up to 3 %, showing that the mean-flow modification due to the non-modal instability dampens the modal instability.

We consider the Nusselt–Rayleigh number problem of Rayleigh–Bénard convection and make the hypothesis that the velocity and thermal boundary layer widths, $\delta _u$ and $\delta _T$, in the absence of a strong mean flow are controlled by the dissipation scales of the turbulence outside the boundary layers and, therefore, are of the order of the Kolmogorov and Batchelor scales, respectively. Under this assumption, we derive $Nu \sim Ra^{1/3}$ in the high $Ra$ limit, independent of the Prandtl number, $\delta _T/L \sim Ra^{-1/3}$ and $\delta _u/L \sim Ra^{-1/3} Pr^{1/2}$, where $L$ is the height of the convection cell. The scaling relations are valid as long as the Prandtl number is not too far from unity. For $Pr \sim 1$, we make a more general ansatz, $\delta _u \sim \nu ^{\alpha }$, where $\nu$ is the kinematic viscosity and assume that the dissipation scales as $\sim u^3/L$, where $u$ is a characteristic turbulent velocity. Under these assumptions we show that $Nu \sim Ra^{\alpha /(3-\alpha )}$, implying that $Nu \sim Ra^{1/5}$ if $\delta _u$ were scaling as in a Blasius boundary layer and $Nu \sim Ra^{1/2}$ (with some logarithmic correction) if it were scaling as in a standard turbulent shear boundary layer. It is argued that the boundary layers will retain the intermediate scaling $\alpha = 3/4$ in the limit of high $Ra$.