On the stages of vortex decay in an impulsively stopped, rotating cylinder

The flow within an infinitely long cylinder exhibiting solid-body rotation (SBR) is impulsively stopped. The complete decay of the initial SBR is captured by means of direct numerical simulations for a wide range of Reynolds numbers ($Re$). Five distinct stages are identified during the decay process according to their flow structure and their underlying mechanisms of kinetic-energy dissipation. Initially, the laminar boundary layer undergoes a primary centrifugal instability, which causes the formation of coherent Taylor rolls. The flow then becomes turbulent, once the Taylor rolls are corrupted by secondary instabilities. Within the turbulent stage, two phases are distinguished. In the first turbulent phase, the SBR core is still intact and turbulence is sustained. The mean velocity profile is well described by the superposition of a near-wall region, a retracting SBR core and an intermediate region of constant angular momentum. In the latter region, the magnitude of angular momentum in viscous units $l^{+}(Re)$ is approximately constant in time. In the second turbulent phase, the SBR core breaks down, turbulence starts to decay exponentially and the kinetic energy of the mean flow decays logarithmically. Eventually, the flow relaminarises and the velocity profile of the analytical solution for purely laminar decay is recovered, albeit at an earlier temporal instant due to the net effect of transition and turbulence.


Introduction
Flows above concave walls have been studied for over a century due to the strong impact of curvature onto the properties of laminar, turbulent and transitional flows. For instance, curvature affects the stability of boundary layers (Rayleigh 1917), the mechanism of transition to turbulence through the occurrence of secondary instabilities (Swearingen & Blackwelder 1987), the statistical properties of turbulent boundary  I  II  III IV V FIGURE 1. Schematic describing the evolution of vorticity (colour coded) and azimuthal velocity profile (vectors) during the distinct stages of the spin-down process: initial condition (ic), laminar stage (I), instabilities and transition to turbulence (II), sustained turbulence with intact vortex core (III), corruption of vortex core and decay of turbulence (IV) and relaminarisation (V). Spirals illustrate the existence of turbulent fluctuations in stages III and IV.
layers (Meroney & Bradshaw 1975) and the decay of anisotropic turbulence in wall-bounded flows (Verschoof et al. 2016;Ostilla-Mónico et al. 2017). The present study describes a relatively simple flow which encompasses laminar, transitional, turbulent and decaying regimes under the influence of concave walls. An infinitely long cylinder of radius R has its axis aligned with the axial direction of a cylindrical coordinate system with radial, azimuthal and axial coordinates denoted respectively (r, ϕ, z). The cylinder is filled with an incompressible Newtonian fluid of kinematic viscosity ν and rotated with angular velocity Ω (see figure 1) until solid-body rotation (SBR) of constant axial vorticity ω = 2Ω is achieved. This flow is characterised by the following velocity field u ϕ (r) = Ωr, u r = u z = 0. (1.1a,b) The numerical experiment begins at the temporal instant t = 0 when, starting from the condition of SBR, the rotation of the cylinder wall is suddenly stopped. Following this event, a transient unsteady flow develops, referred to here as spin-down. Figure 1 presents the different stages of the spin-down process. Each stage is characterised by unique flow features, which are strongly influenced by the boundary conditions -specifically the concave walls of the cylinder. After the laminar boundary layer formation (stage I), centrifugal instabilities emerge as addressed experimentally by Euteneuer (1972) and Mathis & Neitzel (1985). These instabilities have also been studied analytically by Neitzel (1982) and Kim, Song & Choi (2008). However, only Euteneuer (1972) extended their work up to the nonlinear saturation of the primary instability (in stage II). Yet the subsequent stages of the spin-down process have not been investigated: the secondary instability (end of stage II), a stage of sustained turbulence (stage III), the decay of turbulence (stage IV) and the relaminarisation (stage V) itself.

Motivation and objectives
These later stages (stages II-V) are characterised by centrifugal instabilities and the onset of turbulence. The kinetic energy initially present in the SBR is not only dissipated through the viscous dissipation associated with a time-varying velocity profile of laminar spin-down, but also converted to turbulent fluctuations and eventually dissipated via turbulent viscous dissipation. When a large fraction of the total energy has been dissipated, turbulence does not self-sustain and a stage of viscous decay occurs yielding relaminarisation. As mentioned in studies on decaying Stages of vortex decay in an impulsively stopped cylinder 885 A6-3 Taylor-Couette (TC) flow by Verschoof et al. (2016) and Ostilla-Mónico et al. (2017), the rate at which energy is dissipated during the sustained and decaying turbulent regimes is not known a priori. In a similar manner, the details of production and dissipation of energy associated with turbulent fluctuations are also poorly understood. The statistical properties of the turbulent flow and the turbulence production and decay processes depend on the Reynolds number (Re) in at least two instances. On the one hand, there is the obvious impact of Re on the relative importance of viscous and inertial stresses, and thus on turbulent statistics. On the other hand, the value of Re determines the stability properties of the laminar boundary layer forming at the initial stages of spin-down, thereby determining the boundary-layer properties at the instant in which instability and transition to turbulence occur.
The objective of the present work is to describe the spin-down process throughout all of its phases from onset of centrifugal instabilities to the decay of turbulence. Particular focus lies on the flow stages that have not been discussed previously and on the analysis of the turbulent properties and the effect of Re. This study by no means strives to completely cover all aspects of the spin-down process. Rather the paper provides an initial overview of this complex transient flow and its phenomena. Each stage on its own has significant potential for further investigations and therefore makes the spin-down problem an interesting canonical flow to assess unsteady turbulence in the presence of concave walls.
Each stage of the flow inherits features of related canonical flows influenced by concave wall curvature. As such, the present study attempts to draw parallels to prior studies on concave boundary layers, which are therefore briefly reviewed in the following.

Review of concave-wall boundary layers
In an axisymmetric two-dimensional flow, the equilibrium between pressure and centrifugal forces is unstable if the magnitude of the circulation Γ (r) = 2π r 0 ωr dr (ω is the vorticity) decreases with increasing radius r. By identifying this inviscid centrifugal instability mechanism Rayleigh (1917) paved the way for research on the influence of curvature on wall-bounded flows. Subsequently, canonical flows with flat boundaries in the streamwise direction s were also assessed in their respective curved counterpart (figure 2). Examples include the spatially developing boundary layer (figure 2a, statistically steady, no streamwise pressure gradient ∂p/∂s = 0, spatially developing boundary-layer thickness ∂δ/∂s = 0), the fully developed channel flow (figure 2b, statistically steady, ∂p/∂s = 0, ∂δ/∂s = 0), the Couette flow (figure 2c, statistically steady, ∂p/∂s = 0, ∂δ/∂s = 0) and the temporally developing boundary layer (figure 2d, statistically unsteady, ∂p/∂s = 0, ∂δ/∂s = 0) after a sudden change of boundary conditions. In the following, we briefly review studies that modify these canonical flow scenarios to similar flows over concave walls. As it will become evident later, the present numerical experiment embodies aspects of all such flows.
Modifying the classic stability problem of a flat, spatially developing boundary layer to account for concave wall curvature (figure 2a) significantly changes its stability properties, as described by Floryan (1991) and Saric (1994). In flows above concave walls, pairs of streamwise vortices, i.e. Görtler vortices (Görtler 1941), are formed, which get corrupted further downstream by secondary instabilities. Experiments by Bippes (1972) and Swearingen & Blackwelder (1987) provide visualisations of two distinct secondary instability mechanisms: a sinuous mode, leading to spanwise meandering of the streamwise vortices and a varicose mode, resulting in hairpin-like structures. Linear stability theory was applied to assess the primary instability (Floryan & Saric 1982) and the secondary instabilities (Hall & Horseman 1991;Li & Malik 1995). Due to its high relevance in turbomachinery, recent work focuses on the receptivity of Görtler vortices towards roughness and free-stream turbulence (Schrader, Brandt & Zaki 2011;Wu, Zhao & Luo 2011), compressibility effects (Ren & Fu 2015) as well as the control of these instabilities (Sescu & Afsar 2018). The canonical, spatially developing boundary layer over flat plates has been extensively studied both in laboratory experiments and in recent numerical simulations (see Schlatter et al. 2009;Wu & Moin 2009;Sillero, Jiménez & Moser 2013). The sizeable computational cost limits the value of Re that can be achieved in numerical investigations. This limitation is shared also by studies of boundary layers over concave walls, which consist almost exclusively by laboratory experiments. Meroney & Bradshaw (1975), Hoffmann, Muck & Bradshaw (1985 and Barlow & Johnston (1988) allow transition to turbulence in a straight channel section and, before the flow becomes fully developed, a boundary layer of finite thickness enters a curved section of the channel. The studies revealed persistence of streamwise rolls with wavelengths similar to the boundary-layer thickness even in the turbulent stage, which result in enhanced Reynolds stresses.
Significantly lower computational effort and experimental complexity than in spatially developing boundary layers is required when turbulence is characterised in spatially confined, fully developed and pressure-driven flows such as channels and pipes (see figure 2b and amongst others Kim, Moin & Moser (1987)). Experimental (e.g. Hunt & Joubert 1979) and numerical (e.g. Moser & Moin 1987) investigations on fully developed curved channel flows also showed deviations in Reynolds stresses due to large-scale, streamwise and wall-parallel vortices with scales similar to the channel height (the so-called Dean instability - Dean (1928)).
Applying curvature to the classical Couette flow leads to a shear flow between two coaxial cylinders: the TC flow, first addressed by Taylor (1923) (figure 2c). As the system can be easily controlled and is satistically stationary, closed and symmetric, it allows accurate measurements with moderate complexity of the experimental apparatus as well as affordable numerical simulations even for large values of Re. Depending upon the relative and absolute rotational speeds, the radii of the two cylinders and fluid properties, a large variety of different flow structures can be produced. A broad body of literature deals with TC thanks to its simple set-up, the variety of competing physical phenomena occurring in the flow and the similarity with Rayleigh-Bénard convection. Recent reviews are given by Fardin, Perge & Taberlet (2014) and Grossmann, Lohse & Sun (2016).
The temporal evolution of a turbulent incompressible boundary layer after an impulsive acceleration of a flat wall -the so-called Stokes' first problem (figure 2d) -is the flat plate counterpart of the problem investigated here in the present study, at least during the first phases of the spin-down process. While the linear stability of the flow was already analysed by Luchini & Bottaro (2001) almost two decades ago, Kozul, Chung & Monty (2016) recently identified and closed a gap in the literature concerning the in-depth analysis of the turbulent state of such flows. Transferring Stokes' first problem to a flow with concave wall curvature, results in an azimuthally accelerated cylinder. According to the Rayleigh criterion the effects of centrifugal instabilities are only present in the case of a cylinder deceleration, which corresponds to the spin-down case investigated in the present study. As mentioned above, literature hereby is limited to the early stages of the flow, suggesting that there is merit in further investigations of the subsequent flow stages.

Outline
The paper is structured as follows. Section 2.1 contains a detailed description of the numerical method and the flow cases considered in the following. Particularly relevant are the description of the adopted Reynolds decomposition and the budget equations utilised to describe the temporal behaviour of the kinetic energy as described in § 2.2, the definition energy spectra (see § 2.4) and the details of the Lagrangian flow visualisation (see § 2.5). Starting with an overview over the temporal flow development in § 3, the different stages of the spin-down are discussed in detail for one of the simulated Re in § § 3.1-3.5. Finally, in § 4 the Re-influence is addressed by evaluating four simulations at different Re, ranging over almost one order of magnitude.

2.1.
Numerical procedure A newly created database of the turbulent spin-down process in cylinders is produced via direct numerical simulation (DNS). The code used in the present study is a mixed-discretisation parallel solver of the incompressible Navier-Stokes equations in cylindrical coordinates (Fabbiane 2011;Mascotelli 2016). Velocity and pressure fields are discretised via a Fourier-Galerkin approach along the two statistically homogeneous azimuthal (ϕ) and axial (z) directions, while second-order explicit compact finite-difference schemes (Lele 1992) based on a three-point computational stencil on an inhomogeneous grid are adopted in the radial direction (r). Spectral accuracy is therefore achieved for the discretisation of all differential operators in the statistically homogenous directions. The accuracies of the differential operators D 1 = ∂/∂r and D * = (∂/∂r)(r(∂/∂r)) operating in the radial direction are fourth and second order, respectively. The incompressibility constraint is enforced within machine accuracy by direct solution of the continuity equation, which is coupled with pressure through the radial component of the momentum equation for the collocation point in the fluid domain closest to the wall. The number of Fourier modes in the azimuthal direction decreases from the wall towards the cylinder axis as a linear function of r, so that the azimuthal resolution r∆ϕ is kept constant across the cylinder. The regularity boundary conditions (RBC, Lewis & Bellan 1990) are based upon the invariance of the solution with respect to the origin of the coordinate system. RBCs and N z are the number of maximum azimuthal and axial Fourier modes used to represent the flow field without accounting for the additional modes required to exactly remove the aliasing error. N r is the number of collocation points adopted in radial direction. The values of spatial and temporal resolutions are computed at the temporal instant of transition to turbulence, for which the friction velocity achieves its maximum value u max τ . Line colours are used in § 4.
are enforced at the cylinder axis, for all wavenumber pairs that exist throughout the cylinder cross-section, or at the radial position that represents the boundary for wavenumber pairs that only exist above certain values of r.
The governing equations are advanced in time starting from the initial condition of a fully established SBR. No pressure gradient is imposed in the axial direction. Temporal discretisation is achieved with an implicit second-order Crank-Nicholson scheme for the linear terms combined with an explicit third-order low-storage Runge-Kutta scheme for the nonlinear part of the governing equations. Random disturbances with constant energy of 10 −12 Ω 2 R 2 satisfying no-slip boundary conditions are superimposed to each wavenumber and velocity component of the initial velocity field in the whole cylinder volume. The first time step of the simulation forces the random disturbance to fulfil the continuity equation. In the resulting divergence-free field, the energy contained in each wavenumber and velocity component space is randomly distributed and is bound by 10 −12.55 Ω 2 R 2 and 10 −11.5 Ω 2 R 2 . Henceforth, governing equations and all variables are normalised via the cylinder radius R and the initial angular velocity Ω of the SBR. Four different numerical experiments are performed, characterised by different values of the Reynolds number Re = ΩR 2 /ν ∈ {3000, 6000, 12 000, 28 000}, where ν is the kinematic viscosity of the fluid. The discretisation parameters are summarised in table 1.
Spatial and temporal resolutions are set to fulfil the requirements for wall-bounded turbulence (see Kim et al. 1987) at all times of the temporal evolution of the flow. The resolutions are expressed in terms of viscous units, i.e. normalised via the kinematic viscosity ν of the fluid and the friction velocity u τ = √ τ w /ρ. Here, τ w is the spatially averaged wall shear stress and ρ is the fluid density. Normalisation in viscous units is indicated with the superscript +. The most stringent requirement for spatial resolution is achieved after onset of transition, when u τ reaches the maximum value u max τ . At this time instant, indicated by dashed black lines in figure 3(a), the azimuthal, axial and minimum radial resolutions are R + ∆ϕ ≈ 11, ∆z + ≈ 6 and ∆r + min ≈ 1.5, respectively. These values are computed without taking into account the additional modes used to exactly remove the aliasing error. It must be noted that such resolution is finer than the one required to correctly describe the onset of turbulent transition, as discussed in § 3. The resolution in viscous units improves for all other later time instants. The axial extent L cyl of the computational domain is a compromise between the need for accommodating several Taylor

Reynolds decomposition and energy budgets
The spin-down process is statistically unsteady, for which the Reynolds decomposition applied to the velocity field reads u i (r, ϕ, z, t) = u i ϕ,z (r, t) + u i (r, ϕ, z, t), (2.1) where index notation is used to indicate the ith velocity component in the respective direction of cylindrical coordinates (r, ϕ, z); · ϕ,z denotes averaging along the statistically homogeneous azimuthal and axial directions. It must be noted that the ensemble average of independent repetitions of the same experiment is applicable to unsteady problems. However, this has not been performed in the present study, for which spatial averaging resulted in sufficient statistical convergence. In equation (2.1), u i (r, ϕ, z, t) is the fluctuating velocity field about the average value u i ϕ,z (r, t). In the following, the notation u i = u i ϕ,z (r, t) and u i = u i (r, ϕ, z, t) is used for brevity. The temporal decay of the kinetic energy contained in the initial SBR is investigated in the present study. In the framework of the Reynolds decomposition, kinetic energy is split into mean kinetic energy K, associated with the averaged flow field u i , and turbulent kinetic energy k, associated with the fluctuating field u i . As u r = u z = 0 in the present flow, the mean kinetic energy is given by K = 1/2 u ϕ 2 . Its temporal evolution is governed by the following budget equation: where T m is the turbulent diffusion, V m the viscous diffusion and m the dissipation of K. The turbulence production term P couples the budget equation of K and k as it draws energy from the mean flow and transfers it to the fluctuating field. Following Mansour, Kim & Moin (1988) and Bilson & Bremhorst (2007) the budget equation for k = (1/2)( u z u z + u r u r + u ϕ u ϕ ) is given by Here, the viscous diffusion V t , the pressure diffusion Π d and the turbulent diffusion T t describe the transport of k, while t is its viscous dissipation. Beyond the averaging in the axial and azimuthal directions, the closed system allows averaging in the cylinder volume V, which is indicated in the following with the volume averaging operator [·]. The volume-averaged total kinetic energy [K] + [k] can be expressed as Due to the no-slip conditions at the decelerated cylinder walls, energy is constantly withdrawn from the system. The temporal change of kinetic energy [K] + [k] can be described by volume averaging and summation of the equations (2.2) and (2.3). All transport terms T t , V t , T m , V m and Π d contained in the two equations vanish by definition. Also the production P, which appears in both equations for [K] and [k] with an opposite sign, vanishes when the volume average of the kinetic energy [K] + [k] is evaluated. As a consequence, only the dissipation terms t and m remain (2.5) 2.3. Wall-based flow description To characterise turbulence in the developing boundary layer during spin-down, the classical friction Reynolds number is introduced (Jiménez et al. 2010), in which δ 99 is the boundary-layer thickness and u τ = √ τ w /ρ is the friction velocity based on the mean wall shear stress Note that in the investigated flow τ w , u τ , δ 99 and Re τ change with time t. Due to the lack of a constant outer velocity, the traditional definition of δ 99 is adapted by using the retracting vortex core in SBR as follows i.e. δ 99 is defined as the distance from the wall, at which the flow deviates 0.01ΩR from the initial SBR.
2.4. Energy spectra The appearance of Taylor-Görtler rolls will be assessed by the analysis of the power spectral density of k. The summands of k = 1/2 u i u i are split up into the contribution of different axial modes κ z of wavelength λ z = L cyl /κ z . The axial energy spectra are given by where R is the real part of a complex number,· indicates the Fourier coefficients and the superscript * denotes complex conjugation. By accumulating the energy of all modes except the base flow the summands of k are recovered, To gather information about the linear and subsequent nonlinear growth of the different modes κ z in the boundary layer, the energy spectrum Φ u ϕ u ϕ is spatially averaged over the fluid volume contained in the boundary layer V δ = π(2Rδ 99 − δ 2 99 )L cyl , resulting in the one-dimensional (1-D) spectrum for each time step; ξ ϕϕ thus describes the contribution to u ϕ u ϕ of each wavenumber κ z throughout the whole boundary layer.

Lagrangian flow visualisation: finite-time Lyapunov exponent
The identification of Lagrangian coherent structures (LCS) can provide deeper insight into the material transport of complex flows (Haller 2015). LCS separate flow regions of coherent movement. A classical method to calculate attracting LCS is the backward finite-time Lyapunov exponent (FTLE) σ b T int (t), which determines the local attraction rate of particle tracks over a finite time span. The first step to determine the backward FTLE is seeding massless tracers in the computational domain and tracking them backward in time over the timespan T int . In the present work, flow-map interpolation as suggested by Brunton & Rowley (2010) is applied to reduce the computational effort. The resulting flow map Ψ t−T int t is used to determine the right Cauchy-Green tensor (2.12) The terms of the deformation gradient ∇Ψ t−T int t are obtained by the numerical schemes described in § 2.1. By evaluating and normalising the maximum eigenvalue Λ max (C), 885 A6-10 F. Kaiser and others  I  II  III  IV  V  I  II  III  IV is determined. Shadden, Lekien & Marsden (2005) define thin ridges of the FTLE as LCS. For the present flow, the FTLE provides a clear visualisation of the complex three-dimensional interface between the boundary layer and the vortex core. The FTLE is particularly useful during the onset of three-dimensionality in the flow, as it captures the location and evolution of the secondary instability through an alternative visualisation of the flow field.

Flow stages
Dimensional analysis provides two dimensionless groups for the present flow: the Reynolds number Re = ΩR 2 /ν and a dimensionless time. Two distinct yet convertible representations of the dimensionless time are utilised in the following: the viscous time νt/R 2 and the outer time Ωt = Re · νt/R 2 , where Ωt = 2π represents a full revolution of the SBR. This section discusses the temporal development of the spindown flow at Re = 12 000, before the Re scaling is addressed in § 4.
Figures 4 and 5 provide a first overview of the flow evolution. Roman numbers I-V and background colouring indicate the different stages of the decay process. While each stage is characterised by unique features, the transition between stages occurs gradually and thus cannot be exactly localised in time.
During stage I, a stable, laminar boundary layer grows from the cylinder wall at a rate proportional to the viscous time scale √ νt (see figure 4a), while the friction velocity u τ decreases. Since the increase in boundary-layer thickness outweighs the decrease in u τ , Re τ increases mildly in time. Figure 4(b) shows the corresponding evolution of volume-averaged mean and turbulent kinetic energy. As expected, [K] decreases and the initial random disturbances contained in [k] also decay in the stable boundary layer. Stage II is characterised by the centrifugal instability associated with the emergence of Taylor vortices and subsequent breakdown to turbulence. The sudden increase of u τ and δ 99 (and thus Re τ ) is accompanied by the exponential growth of [k].
After transition to turbulence, stage III is entered. In this stage, [k] varies slowly in time. The flow consists of the superposition of a constantly retracting SBR core  I  II  III  IV  V   I  II  III  IV  at the centre of the cylinder and a turbulent boundary layer close to the cylinder wall. The latter is characterised by decreasing u τ and increasing δ 99 . The beginning of stage IV is marked by the breakdown of the SBR core, in which turbulent fluctuations become non-negligible and induce a reduction of u ϕ inside the vortex core. Two distinct phenomena occur during stage IV. First, the breakdown of the SBR core is accompanied by a sudden increase of δ 99 with no discernible effect on u τ or [k]. Then, after the SBR core has been completely eroded by turbulence, [k] and Re τ drop in time. Finally, the temporal evolution of the flow ends with the viscous decay of stage V.

Stage I: laminar boundary layer
The laminar stage of the flow, for which u i = 0 and thus u lam ϕ = u ϕ , can be described by the analytical solution of Neitzel (1982) where J 0 and J 1 are Bessel functions of the first kind and β i are roots of J 1 (β i ) = 0. For small t, while δ 99 /R is small and the local curvature is negligible, the present boundary layer above the concave wall is similar to Stokes' first problem (e.g. Schlichting 1979).
Thus, for small t, the growth rate of δ 99 coincides with the growth rate of the boundary 885 A6-12 F. Kaiser and others  (Euteneuer 1972) Euteneuer (1972) Euteneuer (1972) Euteneuer (1969) Kim & Choi (2006) Mathis & Neitzel (1985) Maxworthy ( Maxworthy (1971) and Mathis & Neitzel (1985) have been rescaled by a factor 0.5 due to their different definition of λ c , as noted by Kim & Choi (2006). (b) Temporal evolution of δ 99 and λ T (t) for the present numerical simulation compared against the measurement by Euteneuer (1972)  However, with increasing δ 99 /R the growth rate of δ 99 in the spin-down problem, computed exploiting the exact analytical solution (3.1), deviates from δ St 99 due to curvature effects. The deviation is such that δ 99 can be made directly proportional to √ νt, only if the proportionality coefficient a lam is a weak function of time.
For simplicity, in the following analyses the growth rate during stage I is characterised through a constant growth rate coefficient a lam , which is determined via a least square fit of the expression to the analytical solution for the spin-down process, which yields a lam ≈ 3.68. The approximation of neglecting the temporal variation of a lam in (3.3) introduces an error smaller than 0.01δ 99 for 0 < δ 99 < 0.28R. The good agreement of (3.3) (black line) with the numerical results is emphasised in figures 4(a) and 7(a).

3.2.
Stage II: emergence of Taylor rolls and laminar-to-turbulent transition If Re is large enough (the stability limit in the literature varies in the range of 128 < Re < 350), the boundary layer is linearly stable only until a critical time θ c = νt c /R 2 , after which the boundary layer undergoes linear primary centrifugal instability. The instability results in the emergence of radial plumes, which later evolve into Taylor rolls. The plumes occur at a characteristic time-dependent spacing λ T (t) in the axial direction, which at the critical time θ c is representative of the linearly most-amplified axial wavelength, the so-called critical wavelength λ c = λ T (θ c ). via linear stability analyses (i.e. Kim & Choi 2006) and laboratory experiments (Maxworthy 1971;Euteneuer 1969Euteneuer , 1972Mathis & Neitzel 1985). In the present work, λ c is measured as the most energetic wavelength of ξ ϕϕ (κ z , θ c ) (see § 2.4 for the definition of ξ ϕϕ ) during the early stage of the linear growth and shows excellent agreement with these existing studies. Figure 6(b) compares λ T (t) computed from the present numerical data against the experimental results of Euteneuer (1972) at Re = 12 000. Euteneuer (1972) estimated the critical wavelength by counting the number of toroidal Taylor rolls aligned along the axial dimension of his experimental set-up. It was found that the rate of change in the observed number of streamwise rolls changes abruptly. This change in growth rate was referred to as 'Knickstelle', the German word for 'kink', the origin of which will be discussed in the following. Figure 7 shows the initial evolution of δ 99 , u τ , [k], [K] and ξ ϕϕ . The flow at four selected time instances, marked with dashed black lines in figure 7, is visualised in figure 8. Figure 8 consists of a pre-multiplied spectrum κ z Φ u i u i , the 1-D spectrum ξ ϕϕ averaged over the boundary layer and a (r, z)-slice of u ϕ . The respective boundary-layer thickness δ 99 is depicted with an orange line. The complete    1, figure 8a) most modes decay and only a narrow band of axial wavenumbers κ z ≈ κ c is amplified. This is most clearly visible in figure 7(c), where the temporal evolution of ξ ϕϕ is reported for relevant wavenumbers. As the mean velocity profile changes due to the temporal growth of the boundary layer, the stability properties of the mean velocity profile also change, and so does the most amplified wavenumber. Already at Ωt = 2 wavenumbers κ z < κ c are amplified and carry more energy than κ c . This temporal change of the most amplified disturbances are represented in figure 7(d), in which ξ ϕϕ is reported as function of κ z at different temporal instances distinguishable through colour coding from the start of spin-down (bright colour) onwards.
Starting from Ωt = 2 the energy of almost all wavenumbers grows exponentially, leading to increasing [k] (figure 7b), until nonlinear saturation is reached and the growth stops at Ωt = 4.3. At this moment, the initially amplified wavenumbers achieve their maximum energy content. The fluctuations start to influence the mean profile u ϕ and thus δ 99 , u τ and [K]. The pre-multiplied spectrum Φ u z u z shows a distinct second peak at a radial position close to δ 99 , marked by an orange line in the left panels of figures 8(c) and 8(d). The second peak is related to the axial flow at the head of the radial jets, which starts the formation of plumes and streamwise vortices by redirecting the radial flow in axial direction. This process is accompanied by a simultaneous rapid growth of δ 99 , u τ and [k] (figures 7a and 7b). The friction velocity u τ reaches its maximum at Ωt = 5.3, right before the breakdown of this still coherent and quasi-axisymmetric flow through presence of a secondary instability. This effect will be discussed later in more detail. At this instance (figure 9e), the plumes are roughly twice as tall in the radial direction than they are wide in the axial direction, and begin to merge with adjacent plumes towards restoring an aspect ratio of unity. The sudden merging of consecutive plumes causes the dominant wavelength λ T to shift from the critical wavelength of early linear primary instability to λ T ≈ δ 99 (figure 7d). This, together with the rapidly growing δ 99 , explains the 'kink' observed by Euteneuer (1972). It can be seen in figure 6(b) that the numerical data reproduce the 'kink' around νt/R 2 ≈ 3.7 × 10 −4 (Ωt ≈ 4.4). However, during the later evolution of the flow, large deviations between the present data and the curve published by Euteneuer exist. These deviations are related to the onset of the secondary instabilities and breakdown to turbulence, as discussed in the following.
Eventually, the interaction of the slow fluid transported from the wall with the fast fluid in the core leads to secondary instabilities, which cause the breakdown of the quasi-axisymmetric flow patterns formed after the primary centrifugal instability. Two distinct competing mechanisms responsible for the onset of the secondary instability have been identified in prior studies on concave-wall flows (Saric 1994): the varicose and the sinuous mode, described in § 1. In order to observe which of the two mechanisms dominate the onset of secondary instabilities during spin-down, the LCS in the flow are visualised by means of the FTLE ( § 2.5). Figure 9 and supplementary material Movie2.mp4 show the temporal evolution of the three-dimensional flow during transition. Linear amplification of the random disturbances given by the initial condition (details in § 2.1) results in Taylor rolls that do not cover the whole circumference of the cylinder. As such, figure 9(a-c) present Taylor rolls with distinct start and ending points along the circumference. Under the action of shear, the two ends of the Taylor vortices roll up into hairpin-like vortices (figure 9d), which grow in time (figure 9e) and eventually lead to a turbulent flow populated by a richer range of vortical structures (figure 9f ). Note that the detected structures are similar to FTLE visualisations of an isolated hairpin vortex in a turbulent boundary layer (Green, Rowley & Haller 2007). Therefore, the varicose mode seems to be the dominant secondary instability mechanism of the spin-down process.
3.3. Stage III: turbulent regime with intact vortex core Stage III is entered once the plumes break down into a turbulent flow. The statistical flow quantities only change slowly and multiple characteristic features co-exist. Figure 10 shows two time instances at the beginning (figure 10a) and during the later evolution (figure 10b) of stage III. The radial plumes have merged, leading to an average axial plume spacing of approximately δ 99 . Moreover, velocity fluctuations start to reach outside the boundary layer into the SBR core (figure 10a). Figure 11 shows u ϕ and the mean axial vorticity ω z during stages I-III. The temporal development is colour coded and the black line depicts the end of the laminar-to-turbulent transition and the beginning of stage III. As soon as stage III is entered (Ωt > 9) a region of negligible mean axial vorticity ω z ≈ 0 is established at radial distances between the near-wall region and the SBR. A region where ω z = 0 implies a spatially constant angular momentum l = l(t) = u ϕ r and is thereby equivalent to the marginally stable case of the Rayleigh instability criterion. This region is well known in other flows on concave walls, such as the TC flow (e.g. Marcus 1984;Panton 1991;Ostilla-Mónico et al. 2016) and swirling pipe flow (Kitoh 1991). Recovering the region of spatially constant angular momentum in this unsteady problem implies that the time scale of self-organisation in the turbulent boundary layer is significantly smaller than the outer time scale of the temporally varying mean flow. This is in good agreement with recent laboratory experiments on unsteady and turbulent TC flow by Verschoof et al. (2018).
Examining the temporal development of the instantaneous value of Re τ during stage III (figure 4a) shows that the decline of u τ is compensated by the growth of δ 99 so that Re τ increases only slightly in time. Therefore, near-wall similarity is expected when turbulent statistics are normalised in time-varying viscous units, despite the statistical unsteadiness of the boundary layer. Figure 12(a) reports the mean velocity profile u + ϕ at different temporal instances. A clear collapse of u + ϕ against the wall-normal distance y + = (R − r)u τ /ν is visible as soon as stage III is entered. The approximate logarithmic behaviour is found for y + > 25, resembling what has been observed for the statistically steady TC flow, for which the constants M and N are related to the degree of curvature of the system (see, for example, Lathrop, Fineberg & Swinney (1991); Panton (1991)). This ratio, here described by δ 99 /R, varies in time. Interestingly, disregarding the significant increase of δ 99 /R during stage III, M and N are observed to be only a function of Re and do not depend on time.
An even clearer collapse within the boundary layer is obtained for the profiles of angular momentum in viscous units, as shown in figure 12(b). In addition to being constant throughout a range of wall-normal distances, l + is found to be approximately constant in time, indicating that l ∼ u τ during stage III. Profiles of both u + ϕ and l + at all values of Re are reported in the Appendix.  In order to understand the effects of turbulence on the decay of the SBR core, the terms of the budget equations for K and k ( § 2.2) are studied for an early (Ωt = 9) and a late (Ωt = 32) phase of stage III. Again, normalisation with the instantaneous wall units is adopted. Furthermore, the terms are pre-multiplied by the factor r/R to compensate for curvature effects and by y + to compensate for the use of a logarithmic scale (Hoyas & Jiménez 2008). Integrated over the cylinder volume, all transport terms are zero by definition. The aforementioned pre-multiplication allows one to capture this property visually in figure 13. As an additional benefit, this representation highlights the significant contribution of the outermost part of the boundary layer to the overall development of kinetic energy. Similarly to u + ϕ , the profiles only collapse in the near-wall region. Figure 13 from the outer buffer layer towards the wall. In the vorticity free region (y + > 100) all viscous terms except of t are negligible and the transport of K and k from the vortex core towards the buffer region is mainly driven by turbulence (T t and T m ). Starting from the buffer region (y + > 5), k is produced all through the boundary layer. In the logarithmic region, the sum Ξ = P + t of production and dissipation of k tends towards a local excess of dissipation, in contrast to fully developed flat boundary layers (Mansour et al. 1988). While + m is roughly unchanged between the two exemplary time instances, + t increases, implying that the relative contribution of [ t ] to the total dissipation [ tot ] increases during stage III. An overview of the temporal evolution of the statistics shown in figures 11-13 is given in supplementary material Movie3.mp4.

3.4.
Stage IV: vortex-core breakdown and decay of turbulence Regime III ends with the breakdown of the SBR core and the sudden increase of δ 99 (figure 4a). Figure 14 shows three time instances during stage IV. The region of constant angular momentum ( ω z ≈ 0) still exists during this flow stage. At the start of stage IV, u ϕ at the interface between SBR and boundary layer is higher than in the initial condition, as indicated by the dashed circle in figure 4. Eventually, the velocity profile within the whole SBR core substantially deviates from the initial condition, as indicated by the solid circle for the time instances Ωt = 90 and Ωt = 170. To understand how the breakdown of the SBR core occurs, the influence of the different K transport terms on the flow state at a time t is evaluated. The cumulative influence A of such terms is derived by integration of the respective term from t = 0 until the time t under consideration where T is a generic term in the K or k budget equation. Figure 14(b) shows A for the terms of the K budget equation at the same temporal instances shown previously in figure 14(a). Both the excess of mean kinetic energy at Ωt = 35 as well as the decaying core flow can be clearly related to the turbulent transport term A(T m ). At Ωt = 35 fast fluid from the outer areas of the cylinder has been transported inwards by turbulent fluctuations, inducing larger u ϕ in the SBR at the interface with the boundary layer. For the decaying core flow the mechanism is the opposite; A(T m ) transports K from the core and towards the wall. The azimuthal velocity profile suggests that the boundary layer does not yet reach the core region. Yet there are enough fluctuations in the core to allow a significant transport of K (see also figure 10b). To understand how fluctuations reach into the SBR core, A is also determined for the terms of the k budget equation (figure 15). The influence of turbulent transport A(T t ) on the propagation of k into the core flow is small. However, the impact of pressure diffusion A(Π d ) is significantly larger and not fully cancelled by A( t ).
Due to the closed geometry under consideration, the length scale of the boundary layer δ 99 is bounded by the cylinder radius R and tops off with the breakdown of the vortex core. Moreover, after the transition to turbulence, the wall shear stress and the friction velocity continuously decrease. As time passes this leads to small Re τ (figure 4a). While [P]/[ t ] is approximately constant and equals unity during stage III of the spin-down process, the ratio decreases during the later phases of stage IV. As a consequence also [k] starts to decline (figure 4b) and [ m ] becomes the main contributor to [ tot ] (see figure 5b).

Stage V: relaminarisation
Eventually, k becomes negligible and the flow relaminarises. The temporal evolution of u ϕ during transition from stage IV to V is shown in figure 16(a). The viscous nature of the flow is highlighted by means of the K budget equation for Ωt = 2400, shown in figure 16(b). Both P and T m are negligible compared to the other terms not involving turbulent quantities. With the decay of turbulence, the flow recovers the laminar state and u ϕ approaches the analytical solution for the laminar flow u lam ϕ , given by (3.1) and visualised as a dashed blue line in figure 16(a). However, the effect of turbulent dissipation results in a temporal shift ∆T when u ϕ is compared with the corresponding laminar solution u lam ϕ : u ϕ (t) = u lam ϕ (t + ∆T). As shown in figure 16(a), for Re = 12 000 the present solution at time Ωt = 2400 is compared to the laminar solution at Ωt = 4430, which results in effective temporal shift of Ω∆T = 2030.

Effect of Reynolds number
This section discusses the influence of varying Re in the range 3000 < Re < 28 000 on the temporal evolution of the flow. First, the Re-scaling of the transition onset is addressed and compared with classical Taylor-Couette flow in § 4.1. Second, the Re-dependent effects of transition on the evolution of δ 99 are evaluated in § 4.2 before the influence of Re on the evolution of Re τ is discussed in § 4.3. Furthermore, the scaling of energy decay (both [K] and [k]) is analysed in § 4.4. Eventually, a Re-independent flow feature is pointed out in § 4.5, that might facilitate future experimental investigations.

Critical Taylor number
In TC flow with a gap width d, a resting outer cylinder and a rotating inner cylinder of radius R i (angular velocity Ω i ) the onset of Taylor rolls typically occurs for a characteristic Taylor number of (4.1) see e.g. Fardin et al. (2014). In contrast, the stability of the spin-down flow is usually described by a critical time θ c and a critical wavelength λ c . The critical values can be obtained by stability analysis as performed by Neitzel (1982) and Kim & Choi (2006), which showed good agreement with present DNS results as reported in § 3.2. For the sake of comparison between TC and spin-down flow, the critical time θ c of the spin-down problem can also be translated into a Taylor number. Defining the SBR core as the inner cylinder and δ 99 as the gap width, the spin-down Taylor number results in Note that different flow events and their respective onset time can be analysed and related to a Taylor number. Examples include the visual onset of instabilities (Euteneuer 1972) or the time where the energy of the critical wavelength λ c first starts to increase. However, for the following evaluation, the critical time θ c is utilised, that  was obtained by the propagation theory of Kim & Choi (2006). Here, θ c is the point in time at which the energy growth rate of the most amplified disturbance is larger than the energy decay rate of the base flow. Its Re-dependence can be approximated by if Re is large enough. Equation (4.3) is inserted into the approximate expression for δ 99 (see (3.3)) to provide the critical boundary-layer thickness δ c ≈ 11.28R · Re −2/3 , which plugged into (4.2) leads to the critical Taylor number Ta c ≈ 37.9 R − δ c R = 37.9 1 − 11.28Re −2/3 . (4.4) Thus, similar to TC flow also for the spin-down problem a Taylor number Ta c can be formulated, that is related to the onset of Taylor rolls. Ta c tends to the constant value 37.9 for Re → ∞.

4.2.
Boundary-layer thickness δ 99 Figure 17(a) presents the evolution of δ 99 for different Re. The value of δ 99 collapses for all Re during the laminar stage if time is expressed in units of νt/R 2 . The thickness of the laminar boundary layer grows as described by (3.3) proportionally to √ νt with a proportionality constant a lam ≈ 3.68, until Taylor rolls start to emerge. The onset of Taylor rolls happens at an earlier time and thus at lower value of δ 99 with increasing Re. During stage II the fast increase of δ 99 is related to the growth of radial plumes as discussed in figure 8(d) in § 3.2. Subsequently, the growth of the plumes becomes nonlinear and eventually ceases. Secondary instabilities appear and initiate the breakdown to turbulence. Also in the turbulent flow of stage III the boundary layer grows approximately proportionally to √ νt (figure 17a), albeit with a different and Re-dependent proportionality constant a turb = a turb (Re). The ratio a turb /a lam is plotted in figure 17(b) for different values of Re. Utilising the present results, the ratio of these two quantities can be linked to Re by the empirical correlation a turb /a lam = 0.1Re 1/3 . (4.5) Stages of vortex decay in an impulsively stopped cylinder  Figure 18 shows the temporal history of Re τ for each value of Re considered in the present study. As expected, larger peak values of Re τ correspond to larger values of Re. For the same total energy of the initial perturbations, a larger Re leads to an earlier onset of the primary centrifugal instability, so that stage II is entered at an earlier time, which also yields a smaller value of δ 99 . For the two smallest values of Re considered in the present work, the simultaneous increase of δ 99 and decrease of u τ result in an almost constant value of Re τ during the boundary-layer growth with an intact vortex core in stage III. For larger values of Re, however, the growth of δ 99 increasingly outweighs the decrease of u τ , resulting in a mild increase of Re τ during phase III. The breakdown of the vortex core, which marks the beginning of stage IV and the related sudden increase of δ 99 , is associated with the peak values of Re τ for all considered values of Re. This is a consequence of the boundary-layer thickness being bounded by the cylinder radius (δ 99 /R = 1) and the monotonically decreasing u τ .

Decay of mean and turbulent kinetic energy
The decay of [K] and [k] at different Re is shown in figure 19. The data are normalised by the volume-averaged total kinetic energy of the initial condition [K 0 ] = [K](t = 0) = 1/4, which is identical for all considered Re. The temporal evolution of [K] and [k] are presented with a logarithmic (figures 19a and 19c) and a linear time scale (figures 19b and 19d) to emphasise the early and the later stages of the decay, respectively. Different non-dimensional time units are used, which are related by Ωt = Re · νt/R 2 .
The larger the value of Re, the slower the decay of [K] in units of Ωt, due to the expected less dominant role of viscosity; see figure 19(a). There is a clearly noticeable increase in the decay rate (i.e. the negative slope) of [K] when Taylor rolls emerge and the flow transitions to turbulence during stage II. Towards the end of stage III the decrease of [K] follows a logarithmic decay for all cases, which can be described by where h [K] is constant and g [K] (Re) depends on Re. The fact that the logarithmic decay of (4.6) is either approached from larger (low Re) or lower (high Re) values of [K] is best understood by considering the temporal evolution of [k] shown in figure 19(c).   figure 18). Roman numbers I-V and the background highlight the different stages of the decay.
Note that the initial decrease of [k] before the onset of primary instability, which was previously discussed in the context of figure 4( in stage V, if the latter is shifted by a Re-dependent ∆T = ∆T(Re) (see figure 19b), . (4.9) Note that the decay of [K lam ] and as such also of [K] is approximately exponential for large t (see (3.1)).

4.5.
Rankine-vortex analogy For all investigated Re, three distinct regions can be distinguished in the radial profile of u ϕ during stage III: the region of SBR, the region of spatially constant angular momentum ( ω z = 0, l = l(t)) and the shear layer at the wall. These are visualised in figure 20(a). The two innermost regions of the u ϕ profile are well approximated by a Rankine vortex u ran ϕ (l, r) ΩR = Ωr ∀ r √ l/Ω l/r ∀ r > √ l/Ω, (4.10) with angular momentum l in its irrotational portion chosen to match the ω x ≈ 0 portion of u ϕ . Rankine vortices at three different values of l = l i are presented in figure 20(a). At each value of l i a matching time instance is selected for each of the four investigated Re. As expected, the velocity profiles almost collapse. Only near the wall and at the interface between the region of SBR and ω z = 0 the profiles differ (see close ups in figure 20a). With increasing Re, the error made by the Rankine-vortex approximation decreases, as shown exemplarily for [K] and its approximation [K] ran estimated with 885 A6-26 F. Kaiser and others (4.10) in figure 20(b). This bears particular potential for experimental investigations, as now u ϕ and the connected quantities (i.e. angular momentum and K) can be estimated by determining the temporal velocity development u ϕ (r 1 , t) at a single radial position r 1 inside the region of constant angular momentum.

Conclusions
In the present work, we present the first DNS results for a complete spin-down process, which occurs when the rotation of an infinitely long cylinder containing fluid in SBR is suddenly stopped. The generated data, which cover almost a decade of Reynolds numbers up to Re = 28 000, allow for the observation of five different flow stages (I-V) during the decay of the initial SBR, each entailing its own characteristic phenomena. The following describes these five stages: (I) In the first stage of spin-down a laminar boundary layer, similar to Stokes' first problem, grows at the cylinder walls. (II) The formation and merging of Taylor rolls following linear and nonlinear growth of a primary centrifugal instability is the characteristic phenomenon of stage II. This causes a fast increase in boundary-layer thickness δ 99 , which is associated with a simultaneous shift of the dominant axial wavelength in the energy spectra. The subsequent transition to turbulence, initiated by secondary shear instabilities, is analysed by means of LCS. Hereby, the streamwise endpoints of the elongated streamwise vortices (Taylor rolls) are identified as the origin of evolving hairpin-like vortices and thus as the nucleus of transition. Larger values of the global Reynolds number Re are associated with earlier transition and smaller wavelengths of the primary instability. As such, the maximum of kinetic energy accumulated in the Taylor rolls decreases with increasing Re. (III) Once transition to turbulence is completed, the boundary-layer flow is turbulent while the retracting SBR vortex core is still intact. Therefore, the flow resembles turbulent TC flow between two concentric rotating cylinders, in which the inner cylinder in TC flow represents the SBR core in the present case. However, the unique feature of the present flow is that the radius of the SBR core varies in time, due to the growth of the boundary layer. In addition, turbulent fluctuations are allowed at the interface with the SBR core, whereas in TC flow a no-slip condition exists at the inner cylinder. Interestingly, in the investigated Reynolds number range, δ 99 still grows at the laminar growth rate √ νt, albeit with a different constant of proportionality. The simultaneous increase of δ 99 and decrease of u τ lead to a negligible temporal variation of Re τ and thus to a collapse of turbulence statistics in the near-wall region. In addition, the outer boundary layer organises as dictated by the marginal stability criterion (Rayleigh 1917). A slowly growing region of constant angular momentum l = u ϕ r forms, so that the mean velocity profile is well approximated by an appropriately defined Rankine vortex. It is found that the angular momentum l + = (r/R) u + ϕ in viscous units is constant in time. Similarly, the constants M and N describing the logarithmic behaviour of u + ϕ do not vary significantly during stage III, despite significant changes of the local curvature δ 99 /R in the boundary layer. During stage III and the early stage IV a logarithmic decay of [K] is observed. (IV) Eventually, the SBR core breaks down. Applying a temporal integration of k and K budget equations, the mechanism behind the core breakdown is identified: the pressure diffusion term Π d transports fluctuations into the core region.
These fluctuations in turn transport K from the core into the boundary layer (term T m ). Towards the end of stage IV turbulent dissipation t overtakes its production P and [k] starts to decrease exponentially. (V) As [k] diminishes, laminar flow is recovered and the velocity profile collapses with the analytical solution (Neitzel 1982). The net effect of transition and turbulence can be lumped into a single parameter, namely the temporal shift Ω∆T that needs to be added to the laminar solution in order to match the velocity profile of the actual flow.
While the present analysis is limited to a first overview of the observed flow regimes, many questions about the spin-down process remain that should be addressed in future works. For instance, the Re scaling of the primary instability should be further investigated. Very high Re values result in small δ 99 /R during the growth of the primary instability, and thus a diminishing influence of curvature. In TC flow the narrow-gap limit with negligible curvature converges towards the linearly stable Couette flow. However, for the present spin-down case, Stokes' first problem is approached when δ 99 /R is very small during transition. As such, for very large Re, Tollmien-Schlichting waves (Luchini & Bottaro 2001) could originate before or simultaneously with the centrifugal instability.
In § 3.2 the corruption of the primary centrifugal instability by secondary instabilities has been speculated to be related to non-axisymmetric Taylor rolls. This hypothesis could be further validated by enforcing the development of perfectly toroidal Taylor rolls, achievable via axisymmetric initial conditions. Furthermore, during the turbulent stage III, the similarity with TC flow suggests that models formerly applied to predict the torque scaling could also be extended to spin-down. Friction laws similar to the ones proposed by Panton (1991), Lathrop et al. (1991) could be derived by leveraging the fact that l + , M and N are temporally invariant and u ϕ is well modelled by a Rankine vortex. The analogy between TC flow and RB convection put forward by Eckhardt, Grossmann & Lohse (2007) could be extended to the present, statistically unsteady flow, in order to predict the dependence of u τ and δ 99 on Re and time.
As visualised via LCS in § 3.2, the turbulent non-turbulent interface (TNTI) in spin-down separates a growing turbulent boundary layer from a rotational retracting SBR core with negligible turbulent fluctuations. The presence of vorticity in the non-turbulent volume (see sketch in figure 2d) is a distinctive feature of TNTI in spin-down compared to other flows (da Silva et al. 2014), where the non-turbulent region is irrotational. Therefore, spin-down presents itself as an interesting canonical flow to assess the role of vorticity annihilation in the propagation of TNTI.   figure 12: Re ∈ {3000, 6000, 28 000}. Higher Re values result in larger values of u τ and thus in a wider spatial extend of both the logarithmic region as well as the region of constant angular momentum. While collapsing during stage III for the individual simulations, the constants M and N of the logarithmic region (3.4) as well as the normalised angular momentum l + (3.5) vary for different Re.