Search

It has been observed that flattening the mean velocity profile of pipe flow by body force can laminarise turbulence, a promising means to reduce frictional drag substantially. To explore whether there is a more efficient body force to eliminate turbulence, we consider time-independent active body forces with varying spatial dependencies. Results confirm that when using an active force, a flattened forced laminar profile is needed to eliminate turbulence, and that a purely streamwise body force is best for laminarisation. While these results required an expensive nonlinear optimisation, it was also observed that the optimal forced profile exhibits reduced linear transient growth (TG). To determine whether the reduction of linear TG alone is a sufficient target for the laminarisation of turbulence, a linear Lagrange Multiplier technique is used to minimise TG of perturbations to the forced laminar profile. The optimal velocity profiles reveal that TG for each azimuthal wavenumber strongly depends on the radial velocity gradient at some specific radial interval. The optimal velocity profile obtained by minimising TG of perturbations of azimuthal wavenumber m = 1 is shown to be able to eliminate turbulence at Re = 2400, but a more effective reduction of TG, including perturbations of higher m, is needed for laminarisation at higher Reynolds number Re = 3000. The streaks formed with the flattened forced laminar profile reveal the mechanism of laminarisation: it is found that the lift-up mechanism in the more flattened forced laminar profile creates a lower streak (i.e. closer to the wall) that is more stable, so that turbulence is harder to maintain through streak instability. Further numerical experiments by disturbing a periodic orbit validated that the breakdown of turbulence self-sustaining mechanisms is mainly caused by suppression of the formation of streaks.

It is well known that buoyancy suppresses, and can even laminarise turbulence in upward heated pipe flow. Heat transfer seriously deteriorates in this case. Through a new DNS model, we confirm that the deteriorated heat transfer within convective turbulence is related to a lack of near-wall rolls, which leads to a weak mixing between the flow near the wall and centre of pipe. Having surveyed the fundamental properties of the system, we perform a nonlinear nonmodal stability analysis. it is found that, the minimal seed becomes thinner and closer to the wall, with increase of buoyancy number C. Most importantly, we show that the critical initial energy required to trigger shear-driven turbulence keeps increasing, implying that attempts to artificially trigger it may not be an efficient means to improve heat transfer at larger C. The new minimal seed, found at C=6, is localised in streamwise direction and is active in the centre of pipe. To find this branch of optimal, we took advantage of a window of linear stability. While the nonlinear optimal causes transition to convective turbulence directly at this and larger C, transition via the linear instability passes via a travelling wave or periodic orbit solutions. Detailed analysis of the periodic solution reveals three stages: growth of the unstable eigenfunction, the formation of streaks, and the decay of streaks due to suppression of the instability. Flow visualization at C up to 10 also show similar features, suggesting that convective turbulence is sustained by these three typical processes.

Turbulence in the flow of fluid through a pipe can be suppressed by buoyancy forces. As the suppression of turbulence leads to severe heat transfer deterioration, this is an important and undesirable phenomenon in both heating and cooling applications. Vertical flow is often considered, as the axial buoyancy force can help drive the flow. With heating measured by the buoyancy parameter $C$, our DNS show that shear-driven turbulence may either be completely laminarised or transitions to a relatively quiescent convection-driven state. Buoyancy forces cause a flattening of the base flow profile, which in isothermal pipe flow has recently been linked to complete suppression of turbulence (K\"{u}hnen et al. Nat. Phys., 2018), and the flattened laminar base profile has enhanced nonlinear stability (Marensi et al. JFM, 2019). In agreement with these findings, the nonlinear lower-branch travelling-wave solution analysed here, which is believed to mediate transition to turbulence in isothermal pipe flow, is shown to be suppressed by buoyancy. A linear instability of the laminar base flow is responsible for the appearance of the relatively quiescent convection driven state for $C\gtrsim 4$ across the range of Reynolds numbers considered. In the suppression of turbulence, however, i.e. in the transition from turbulence, we find clearer association with the analysis of He et al. (JFM, 2016) than with the above dynamical systems approach, which describes better the transition to turbulence. The laminarisation criterion He et al. propose, based on an apparent Reynolds number of the flow as measured by its driving pressure gradient, is found to capture the critical $C=C_{cr}(Re)$ above which the flow will be laminarised or switch to the convection-driven type. Our analysis suggests that it is the weakened rolls, rather than the streaks, which appear to be critical for laminarisation., Comment: 30 pages, 21 figures

Motivated by the results of recent experiments (K\"uhnen et al., Flow Turb. Combust., vol. 100, 2018, pp. 919-943), we consider the problem of designing a baffle (an obstacle to the flow) to relaminarise turbulence in pipe flows. Modelling the baffle as a spatial distribution of linear drag $\mathbf{f}(\mathbf{x},t)=-\chi(\mathbf{x})\mathbf{u}_{tot}(\mathbf{x},t)$ within the flow ($\mathbf{u}_{tot}$ is the total velocity field and $\chi \ge 0$ a scalar field), two different optimisation problems are considered to design $\chi$ at a Reynolds number $Re=3000$. In the first, the smallest baffle defined in terms of a $L_1 $ norm of $\chi$ is sought which minimises the viscous dissipation rate of the flow. In the second, a baffle which minimises the total energy consumption of the flow is treated. Both problems indicate that the baffle should be axisymmetric and radially localised near the pipe wall, but struggle to predict the optimal streamwise extent. A manual search finds an optimal baffle one radius long which is then used to study how the amplitude for relaminarisation varies with $Re$ up to $15\,000$. Large stress reduction is found at the pipe wall, but at the expense of an increased pressure drop across the baffle. Estimates are then made of the break even point downstream of the baffle where the stress reduction at the wall due to the relaminarised flow compensates for the extra drag produced by the baffle., Comment: 32 pages, 24 figures

In the past two decades, our understanding of the transition to turbulence in shear flows with linearly stable laminar solutions has greatly improved. Regarding the susceptibility of the laminar flow, two concepts have been particularly useful: the edge states and the minimal seeds. In this nonlinear picture of the transition, the basin boundary of turbulence is set by the edge state's stable manifold and this manifold comes closest in energy to the laminar equilibrium at the minimal seed. We begin this paper by presenting numerical experiments in which three-dimensional perturbations are too energetic to trigger turbulence in pipe flow but they do lead to turbulence when their amplitude is reduced. We show that this seemingly counter-intuitive observation is in fact consistent with the fully nonlinear description of the transition mediated by the edge state. In order to understand the physical mechanisms behind this process, we measure the turbulent kinetic energy production and dissipation rates as a function of the radial coordinate. Our main observation is that the transition to turbulence relies on the energy amplification away from the wall, as opposed to the turbulence itself, whose energy is predominantly produced near the wall. This observation is further supported by the similar analyses on the minimal seeds and the edge states. Furthermore, we show that the time-evolution of production-over-dissipation curves provide a clear distinction between the different initial amplification stages of the transition to turbulence from the minimal seed., Comment: Under review for publication in Physical Review Fluids

In this short exposition, we describe equilibria and periodic orbits in terms of the flow map, {\Phi}, and discuss the essentials of the Jacobian-free Newton-Krylov (JFNK) method that can be used to find them. This method requires little more than calls to an existing time stepping code, which {\Phi} can be considered to represent. Fortran90 / MATLAB code is available to try it out for yourself, where, in the template/example the method is applied to the Lorenz system. This code is problem-independent and can be applied to large systems, having initially been developed to find periodic orbits in simulations of pipe flow., Comment: EPSRC Summer School on Modal decompositions in fluid mechanics. DAMTP, Cambridge, 5-8 August 2019. 16 pages. (v2: references added, hyperlinks corrected)

The transition from Tayor vortex flow to wavy-vortex flow is revisited. The Self-Sustaining Process (SSP) of Waleffe [Phys. Fluids 9, 883-900 (1997)] proposes that a key ingredient in transition to turbulence in wall-bounded shear flows is a three-step process involving rolls advecting streamwise velocity, leading to streaks which become unstable to a wavy perturbation whose nonlinear interaction with itself feeds the rolls. We investigate this process in Taylor-Couette flow. The instability of Taylor-vortex flow to wavy-vortex flow, a process which is the inspiration for the second phase of the SSP, is shown to be caused by the streaks, with the rolls playing a negligible role, as predicted by Jones [J. Fluid Mech. 157, 135-162 (1985)] and demonstrated by Martinand et al. [Phys. Fluids 26, 094102 (2014)]. In the third phase of the SSP, the nonlinear interaction of the waves with themselves reinforces the rolls. We show this both quantitatively and qualitatively, identifying physical regions in which this reinforcement is strongest, and also demonstrate that this nonlinear interaction depletes the streaks.

Steady states and traveling waves play a fundamental role in understanding hydrodynamic problems. Even when unstable, these states provide the bifurcation-theoretic explanation for the origin of the observed states. In turbulent wall-bounded shear flows, these states have been hypothesized to be saddle points organizing the trajectories within a chaotic attractor. These states must be computed with Newton's method or one of its generalizations, since time-integration cannot converge to unstable equilibria. The bottleneck is the solution of linear systems involving the Jacobian of the Navier-Stokes or Boussinesq equations. Originally such computations were carried out by constructing and directly inverting the Jacobian, but this is unfeasible for the matrices arising from three-dimensional hydrodynamic configurations in large domains. A popular method is to seek states that are invariant under numerical time integration. Surprisingly, equilibria may also be found by seeking flows that are invariant under a single very large Backwards-Euler Forwards-Euler timestep. We show that this method, called Stokes preconditioning, is 10 to 50 times faster at computing steady states in plane Couette flow and traveling waves in pipe flow. Moreover, it can be carried out using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes to these popular spectral codes. We explain the convergence rate as a function of the integration period and Reynolds number by computing the full spectra of the operators corresponding to the Jacobians of both methods., Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018)

We numerically compute the flow of an electrically conducting fluid in a Taylor-Couette geometry where the rotation rates of the inner and outer cylinders satisfy $\Omega_o/\Omega_i=(r_o/r_i)^{-3/2}$. In this quasi-Keplerian regime a non-magnetic system would be Rayleigh-stable for all Reynolds numbers $Re$, and the resulting purely azimuthal flow incapable of kinematic dynamo action for all magnetic Reynolds numbers $Rm$. For $Re=10^4$ and $Rm=10^5$ we demonstrate the existence of a finite-amplitude dynamo, whereby a suitable initial condition yields mutually sustaining turbulence and magnetic fields, even though neither could exist without the other. This dynamo solution results in significantly increased outward angular momentum transport, with the bulk of the transport being by Maxwell rather than Reynolds stresses.

The magnetorotational instability (MRI) is thought to be a powerful source of turbulence in Keplerian accretion disks. Motivated by recent laboratory experiments, we study the MRI driven by an azimuthal magnetic field in an electrically conducting fluid sheared between two concentric rotating cylinders. By adjusting the rotation rates of the cylinders, we approximate angular velocity profiles $\omega \propto r^{q}$. We perform direct numerical simulations of a steep profile close to the Rayleigh line $q \gtrsim -2 $ and a quasi-Keplerian profile $q \approx -3/2$ and cover wide ranges of Reynolds ($Re\le 4\cdot10^4$) and magnetic Prandtl numbers ($0\le Pm \le 1$). In the quasi-Keplerian case, the onset of instability depends on the magnetic Reynolds number, with $Rm_c \approx 50$, and angular momentum transport scales as $\sqrt{Pm} Re^2$ in the turbulent regime. The ratio of Maxwell to Reynolds stresses is set by $Rm$. At the onset of instability both stresses have similar magnitude, whereas the Reynolds stress vanishes or becomes even negative as $Rm$ increases. For the profile close to the Rayleigh line, the instability shares these properties as long as $Pm\gtrsim0.1$, but exhibits a markedly different character if $Pm\rightarrow 0$, where the onset of instability is governed by the Reynolds number, with $Re_c \approx 1250$, transport is via Reynolds stresses and scales as $Re^2$. At intermediate $Pm=0.01$ we observe a continuous transition from one regime to the other, with a crossover at $Rm=\mathcal{O}(100)$. Our results give a comprehensive picture of angular momentum transport of the MRI with an imposed azimuthal field.

Pipelines are used in a huge range of industrial processes involving fluids, and the ability to accurately predict properties of the flow through a pipe is of fundamental engineering importance. Armed with parallel MPI, Arnoldi and Newton--Krylov solvers, the Openpipeflow code can be used in a range of settings, from large-scale simulation of highly turbulent flow, to the detailed analysis of nonlinear invariant solutions (equilibria and periodic orbits) and their influence on the dynamics of the flow., Comment: 5 pages, 4 figures

Chaotic dynamics of low-dimensional systems, such as Lorenz or R\"ossler flows, is guided by the infinity of periodic orbits embedded in their strange attractors. Whether this also be the case for the infinite-dimensional dynamics of Navier--Stokes equations has long been speculated, and is a topic of ongoing study. Periodic and relative periodic solutions have been shown to be involved in transitions to turbulence. Their relevance to turbulent dynamics---specifically, whether periodic orbits play the same role in high-dimensional nonlinear systems like the Navier--Stokes equations as they do in lower-dimensional systems---is the focus of the present investigation. We perform here a detailed study of pipe flow relative periodic orbits with energies and mean dissipations close to turbulent values. We outline several approaches to reduction of the translational symmetry of the system. We study pipe flow in a minimal computational cell, and report a library of invariant solutions found with the aid of the method of slices. Detailed study of the unstable manifolds of a sample of these solutions is consistent with the picture that relative periodic orbits are embedded in the chaotic saddle and that they guide the turbulent dynamics., Comment: 26 pages, 11 figures, 1 table

Many transitional wall-bounded shear flows are characterised by the coexistence in state-space of laminar and turbulent regimes. Probing the edge boundary between the two attractors has led in the last decade to the numerical discovery of new (unstable) solutions to the incompressible Navier-Stokes equations. However, the iterative bisection method used to achieve this can become prohibitively costly for large systems. Here we suggest a simple feedback control strategy to stabilise edge states, hence accelerating their numerical identification by several orders of magnitude. The method is illustrated for several configurations of cylindrical pipe flow. Travelling waves solutions are identified as edge states, and can be isolated rapidly in only one short numerical run. A new branch of solutions is also identified. When the edge state is a periodic orbit or chaotic state, the feedback control does not converge precisely to solutions of the uncontrolled system, but nevertheless brings the dynamics very close to the original edge manifold in a single run. We discuss the opportunities offered by the speed and simplicity of this new method to probe the structure of both state space and parameter space., Comment: 12 pages, 5 figures

We report the results of three-dimensional direct numerical simulations for incompressible viscous fluid in a circular pipe flow with a sudden expansion. At the inlet, a parabolic velocity profile is applied together with a finite amplitude perturbation in the form of a vortex with its axis parallel to the axis of the pipe. At sufficiently high Reynolds numbers the recirculation region breaks into a turbulent patch that changes position axially depending on the strength of the perturbation. This vortex perturbation is believed to produce a less abrupt transition than in previous studies with a tilt perturbation, as the localized turbulence is observed via the formation of a wavy structure at a low order azimuthal mode, which resembles an optimally amplified perturbation. For higher amplitude, the localized turbulence remains at a constant axial position. It is further investigated using proper orthogonal decomposition, which indicates that the centre region close to the expansion is highly energetic., Comment: 12 pages, 7 figures, presented in BIFD, Paris 2015

Equilibrium solutions are believed to structure the pathways for ergodic trajectories in a dynamical system. However, equilibria are atypical for systems with continuous symmetries, i.e. for systems with homogeneous spatial dimensions, whereas relative equilibria (traveling waves) are generic. In order to visualize the unstable manifolds of such solutions, a practical symmetry reduction method is required that converts relative equilibria into equilibria, and relative periodic orbits into periodic orbits. In this article we extend the fixed Fourier mode slice approach, previously applied 1-dimensional PDEs, to a spatially 3-dimensional fluid flow, and show that is substantially more effective than our previous approach to slicing. Application of this method to a minimal flow unit pipe leads to the discovery of many relative periodic orbits that appear to fill out the turbulent regions of state space. We further demonstrate the value of this approach to symmetry reduction through projections (projections only possible in the symmetry-reduced space) that reveal the interrelations between these relative periodic orbits and the ways in which they shape the geometry of the turbulent attractor., Comment: 5 pages, 4 figures, 1 table

We report that many exact invariant solutions of the Navier-Stokes equations for both pipe and channel flows are well represented by just few modes of the model of McKeon & Sharma J. Fl. Mech. 658, 356 (2010). This model provides modes that act as a basis to decompose the velocity field, ordered by their amplitude of response to forcing arising from the interaction between scales. The model was originally derived from the Navier-Stokes equations to represent turbulent flows and has been used to explain coherent structure and to predict turbulent statistics. This establishes a surprising new link between the two distinct approaches to understanding turbulence., Comment: 21 figures

A new, fully-localised, energy growth optimal is found over large times and in long pipe domains at a given mass flow rate. This optimal emerges at a threshold disturbance energy below which a nonlinear version of the known (streamwise-independent) linear optimal (Schmid \& Henningson 1994) is selected, and appears to remain the optimal up until the critical energy at which transition is triggered. The form of this optimal is similar to that found in short pipes (Pringle et al.\ 2012) albeit now with full localisation in the streamwise direction. This fully-localised optimal perturbation represents the best approximation yet of the {\em minimal seed} (the smallest perturbation capable of triggering a turbulent episode) for `real' (laboratory) pipe flows., Comment: 7 pages, 6 figures

The aim in the dynamical systems approach to transitional turbulence is to construct a scaffold in phase space for the dynamics using simple invariant sets (exact solutions) and their stable and unstable manifolds. In large (realistic) domains where turbulence can co-exist with laminar flow, this requires identifying exact localized solutions. In wall-bounded shear flows the first of these has recently been found in pipe flow, but questions remain as to how they are connected to the many known streamwise-periodic solutions. Here we demonstrate the origin of the first localized solution in a modulational symmetry-breaking Hopf bifurcation from a known global travelling wave that has 2-fold rotational symmetry about the pipe axis. Similar behaviour is found for a global wave of 3-fold rotational symmetry, this time leading to two localized relative periodic orbits. The clear implication is that all global solutions should be expected to lead to more realistic localised counterparts through such bifurcations, which provides a constructive route for their generation., Comment: 5 pages, 6 figures

In stars and planets, magnetic fields are believed to originate from the motion of electrically conducting fluids in their interior, through a process known as the dynamo mechanism. In this Letter, an optimization procedure is used to simultaneously address two fundamental questions of dynamo theory: "Which velocity field leads to the most magnetic energy growth?" and "How large does the velocity need to be relative to magnetic diffusion?" In general, this requires optimization over the full space of continuous solenoidal velocity fields possible within the geometry. Here the case of a periodic box is considered. Measuring the strength of the flow with the root-mean-square amplitude, an optimal velocity field is shown to exist, but without limitation on the strain rate, optimization is prone to divergence. Measuring the flow in terms of its associated dissipation leads to the identification of a single optimal at the critical magnetic Reynolds number necessary for a dynamo. This magnetic Reynolds number is found to be only 15% higher than that necessary for transient growth of the magnetic field., Comment: Optimal velocity field given approximate analytic form. 4 pages, 4 figures

Symmetry reduction by the method of slices is applied to pipe flow in order to quotient the stream-wise translation and azimuthal rotation symmetries of turbulent flow states. Within the symmetry-reduced state space, all travelling wave solutions reduce to equilibria, and all relative periodic orbits reduce to periodic orbits. Projections of these solutions and their unstable manifolds from their $\infty$-dimensional symmetry-reduced state space onto suitably chosen 2- or 3-dimensional subspaces reveal their interrelations and the role they play in organising turbulence in wall-bounded shear flows. Visualisations of the flow within the slice and its linearisation at equilibria enable us to trace out the unstable manifolds, determine close recurrences, identify connections between different travelling wave solutions, and find, for the first time for pipe flows, relative periodic orbits that are embedded within the chaotic attractor, which capture turbulent dynamics at transitional Reynolds numbers., Comment: 24 pages, 12 figures

We propose a general strategy for determining the minimal finite amplitude isturbance to trigger transition to turbulence in shear flows. This involves constructing a variational problem that searches over all disturbances of fixed initial amplitude, which respect the boundary conditions, incompressibility and the Navier--Stokes equations, to maximise a chosen functional over an asymptotically long time period. The functional must be selected such that it identifies turbulent velocity fields by taking significantly enhanced values compared to those for laminar fields. We illustrate this approach using the ratio of the final to initial perturbation kinetic energies (energy growth) as the functional and the energy norm to measure amplitudes in the context of pipe flow. Our results indicate that the variational problem yields a smooth converged solution providing the amplitude is below the threshold amplitude for transition. This optimal is the nonlinear analogue of the well-studied (linear) transient growth optimal. At and above this threshold, the optimising search naturally seeks out disturbances that trigger turbulence by the end of the period, and convergence is then practically impossible. The first disturbance found to trigger turbulence as the amplitude is increased identifies the `minimal seed' for the given geometry and forcing (Reynolds number). We conjecture that it may be possible to select a functional such that the converged optimal below threshold smoothly converges to the minimal seed at threshold. This seems at least approximately true for our choice of energy growth functional and the pipe flow geometry chosen here., Comment: 27 pages, 19 figures, submitted to JFM

In most settings, from international pipelines to home water supplies, the drag caused by turbulence raises pumping costs many times higher than if the flow were laminar. Drag reduction has therefore long been an aim of high priority. In order to achieve this end, any drag reduction method must modify the turbulent mean flow. Motivated by minimization of the input energy this requires, linearly optimal forcing functions are examined. It is shown that the forcing mode leading to the greatest response of the flow is always of m=1 azimuthal symmetry. Little evidence is seen of the second peak at large m (wall modes) found in analogous optimal growth calculations, which may have implications for control strategies. The model's prediction of large response of the large length-scale modes is verified in full direct numerical simulation of turbulence ($Re=5300$, $Re_\tau\approx 180$). Further, drag reduction of over 12% is found for finite amplitude forcing of the largest scale mode, m=1. Significantly, the forcing energy required is very small, being less than 2% of that by the through pressure, resulting in a net energy saving of over 10%., Comment: 6 pages

Fully 3-dimensional computations of flow through a long pipe demand a huge number of degrees of freedom, making it very expensive to explore parameter space and difficult to isolate the structure of the underlying dynamics. We therefore introduce a `2+epsilon' dimensional model of pipe flow which is a minimal 3-dimensionalisation of the axisymmetric case: only sinusoidal variation in azimuth plus azimuthal shifts are retained, yet the same dynamics familiar from experiments are found. In particular the model retains the subcritical dynamics of fully resolved pipe flow, capturing realistic localised `puff'-like structures which can decay abruptly after long times, as well as global `slug' turbulence. Relaminarisation statistics of puffs reproduce the memoryless feature of pipe flow and indicate the existence of a Reynolds number about which lifetimes diverge rapidly, provided that the pipe is sufficiently long. Exponential divergence of the lifetime is prevalent in shorter periodic domains. In a short pipe, exact travelling-wave solutions are found nearby to flow trajectories on the boundary between laminar and turbulent flow. In a long pipe, the attracting state on the laminar-turbulent boundary is a localised structure which resembles a smoothened puff. This `edge' state remains localised even for Reynolds numbers where the turbulent state is global., Comment: 21 pages, 9 figures; as accepted, J. Fluid Mech

The laminar-turbulent boundary S is the set separating initial conditions which relaminarise uneventfully from those which become turbulent. Phase space trajectories on this hypersurface in cylindrical pipe flow look to be chaotic and show recurring evidence of coherent structures. A general numerical technique is developed for recognising approaches to these structures and then for identifying the exact coherent solutions themselves. Numerical evidence is presented which suggests that trajectories on S are organised around only a few travelling waves and their heteroclinic connections. If the flow is suitably constrained to a subspace with a discrete rotational symmetry, it is possible to find locally-attracting travelling waves embedded within S. Four new types of travelling waves were found using this approach., Comment: 24 pages, 14 figures. Accepted, Jou. Fluid Mech

This is a Reply to Comment arXiv:0707.2642 by Hof et al. on Letter arXiv:physics/0608292 which was subsequently published in Phys Rev Lett, 98, 014501 (2007). In our letter it was reported that in pipe flow the median time $\tau$ for relaminarisation of localised turbulent disturbances closely follows the scaling $\tau\sim 1/(Re_c-Re)$. This conclusion was based on data from collections of 40 to 60 independent simulations at each of six different Reynolds numbers, Re. In the Comment, Hof et al. estimate $\tau$ differently for the point at lowest Re. Although this point is the most uncertain, it forms the basis for their assertion that the data might then fit an exponential scaling $\tau\sim \exp(A Re)$, for some constant A, supporting Hof et al. (2006) Nature, 443, 59. The most certain point (at largest Re) does not fit their conclusion and is rejected. We clarify why their argument for rejecting this point is flawed. The median $\tau$ is estimated from the distribution of observations, and it is shown that the correct part of the distribution is used. The data is sufficiently well determined to show that the exponential scaling cannot be fit to the data over this range of Re, whereas the $\tau\sim 1/(Re_c-Re)$ fit is excellent, indicating critical behaviour and supporting experiments by Peixinho & Mullin 2006., Comment: 1 page, 1 figure