A&A 408, 409-414 (2003)
DOI: 10.1051/0004-6361:20031010

On non-linear hydrodynamic instability and enhanced transport in differentially rotating flows[*]

D. Richard 1-3

1 - NASA, Ames Research Center, MS 245-3, Moffett Field, California 94035, USA
2 - LUTH, Observatoire de Paris, 92190 Meudon Cedex, France
3 - GIT-SPEC, Commissariat à l'Énergie Atomique, l'Orme des Merisiers, 91190 Gif-sur-Yvette, France

Received 12 June 2002 / Accepted 17 June 2003

In this paper we argue that differential rotation can possibly sustain hydrodynamic turbulence in the absence of a magnetic field. We explain why the non-linearities of the hydrodynamic equations (i.e. turbulent diffusion) should not be neglected, either as a simplifying approximation or based on boundary counditions (Hawley et al. 1999). The consequences of lifting this hypothesis are studied for the flow stability and the enhanced turbulent transport. We develop a simple general model for the energetics of turbulent fluctuations in differentially rotating flows. By taking into account the non-linearities of the equations of motions, we give constraints on the mean flow properties for the possible development of shear instability. The results from recent laboratory experiments on rotating flows show - in agreement with the model - that the pertinent parameter for stability appears to be the Rossby number Ro. The laboratory experiments seem to be compatible with Ro<-1/2 and Ro > 1 in the inviscid or high rotation rates limit. Our results, taken in the inviscid limit, are coherent with the classical linear stability analysis, in the sense that the critical perturbation equals zero on the marginal linear stability curve. We also propose a prescription for turbulent viscosity which generalize the $\beta$-prescription derived in Richard & Zahn (1999).

Key words: hydrodynamics - instabilities - turbulence

1 Introduction

Differentially rotating flows are present in a wide variety of Astrophysical systems, including stellar interiors, accretion disks, or gaseous planets. Analytical studies teach us that differentially rotating flows are unstable, according to linear stability theory, whenever their angular momentum is decreasing outward (Rayleigh 1916). These flows have been extensively studied but they do not exist in astrophysical context. Nevertheless, linear stability does not mean stability. This is because turbulence is genuinely a non-linear phenomenon. Even when the primary instability is linear in nature, the self-sustenance of fully-developed turbulence requires for the non-linearities to enter the game at least - but not only - for the saturation of amplitudes. In that sense, linear stability analysis is limited to predict, for a restricted class of flows, their instability and only the very early development of the bifurcated solution. It also means that linear analysis predicts instabilities but does not predict unconditional stability. Numerical models fail to maintain hydrodynamical turbulence in Keplerian flows (Balbus et al. 1996). Based on past laboratory experiments from the 1930's (Taylor 1936; Wendt 1933), differentially rotating flows can become unstable for Reynolds numbers of order 10 000, a figure that to our knowledge has not been reached in published numerical simulations of Keplerian flows.

Richard & Zahn (1999) have shown, using laboratory experiments published in the 1930's, that differentially rotating flows exhibit finite-amplitude instabilities in the case where both angular velocity and angular momentum are increasing outward, a class of flows relevant for stellar interiors models. Recent laboratory studies (Richard et al. 2002), showed that the class of flows experiencing shear instabilities is much wider, as demonstrated by experimental results. We aim in this paper to find some hints about the physical mechanisms involved in these instabilities.

2 Background and motivations

In this section we present the basic equations for the evolution of turbulent fluctuations. We discuss the relevant approximations that can be done and their implications regarding the properties of the turbulence and the mean flow.

2.1 Turbulent transport equations

We consider the equation of motion for a viscous incompressible flow,

 \begin{displaymath}\partial_{\rm t} \vec{U} + (\vec{U} . \vec{\nabla}) \vec{U} = - {1
\over \rho} \vec{\nabla} P + \nu \vec{\nabla}^{2} \vec{U}
\end{displaymath} (1)

and we decompose the flow as its mean velocity and pressure fields and their fluctuations, namely

 \begin{displaymath}\vec{U}=\vec{u}+\vec{u{'}} \hspace{1cm} P = p +p{'}
\end{displaymath} (2)

where $\vec{U}$ and P are the total velocity and pressure fields, $\vec{u}$ and p are the time averaged velocity and pressure and $\vec{u'}$ and p' are their time dependent fluctuations. Multiplying the equations of motion by the fluctuation field and averaging over time, considering that

\begin{displaymath}\overline{\vec{u{'}}}= \vec{0} \hspace{1cm} \overline{p{'}} = 0,
\end{displaymath} (3)

(where $\overline{x}$ is the time average of x) we find the evolution equations for the turbulent fluctuations kinetic energy (i.e. the diagonal components of the Reynolds stress tensor) reported in Appendix A. Note that viscous diffusion has been neglected, but viscous dissipation has been kept. We consider circular motion $\overline{u_{r}} = \overline{u_{z}} = 0$), with axial and azimuthal symmetries ( $\partial_{z}= \partial_{\phi} = 0$) for the mean fields and define $\overline{u_{\phi} \over r } = \Omega $. The time averages of axial and azimuthal derivatives of fluctuating terms is neglected (i.e. net turbulent transport is only radial).
                           $\displaystyle {1\over 2} \partial_{\rm t} \overline{{u_{r}{'}}^{2}}$ = $\displaystyle 2 \Omega \; \overline{u_{r}{'}u_{\phi}{'}} - \overline{{u_r{'} \over r} \partial_r (r {u_r{'}}^2)} + \overline{u_r{'} {u_\phi{'}}^2 \over r}$  
    $\displaystyle + {1 \over \rho} \left( \overline{p{'}\partial_r u_r{'}} - \overline{\partial_r p{'} u_r{'}}\right)$  
    $\displaystyle - \nu \left( \overline{(\partial_r u{'}_{r})^2} + {1 \over 2} \ov...
...al_r u{'}_{\phi})^2} + {1 \over 2} \overline{(\partial_r u{'}_{z})^2 } \right),$ (4)

                            $\displaystyle {1\over 2} \partial_{\rm t}\overline{{u_{\phi}{'}}^{2}} =$ - $\displaystyle {\partial_{r} {\Omega r^{2}} \over r} \; {\overline{u_{r}{'}u_{\p...
...phi{'}\partial_r ( u_r{'} u_\phi{'})} - \overline{{u_r{'} u_\phi{'}}^2 \over r}$  
  - $\displaystyle \nu \left( {1 \over 2} \overline{(\partial_r u{'}_{\phi})^2} + (\overline{{u{'}_{r} / r}})^2
+{1 \over 2} (\overline{{u{'}_{\phi} / r}})^2 \right),$ (5)

 \begin{displaymath}{1\over 2} \partial_{\rm t}\overline{{u_{z}{'}}^{2}} = - \ove...
... u_z{'})} - {\nu \over 2} \overline{(\partial_r u{'}_{z})^2} .
\end{displaymath} (6)

This set of equations can be simplified by making the assumption that the turbulence is close to isotropy. This allows us to neglect the time averaged pressure fluctuations terms (see for example Kato & Yoshizawa 1997). We further assume that the characteristic spatial variation scale of the fluctuating components is smaller than the local radius (i.e terms of order  u'i u'j u'k / r can be neglected compared to the ones of order  $\partial_r u{'}_i u{'}_j u{'}_k$). The above equations then become,
$\displaystyle {1\over 2} \partial_{\rm t} \overline{{u_{r}{'}}^{2}}$ = $\displaystyle 2 \Omega \; \overline{u_{r}{'}u_{\phi}{'}} - \overline{u_r{'} \partial_r ( {u_r{'}}^2)}$$\displaystyle - \nu \left( \overline{(\partial_r u{'}_{r})^2} + {1 \over 2} \ov...
...al_r u{'}_{\phi})^2} + {1 \over 2} \overline{(\partial_r u{'}_{z})^2 } \right),$ (7)

 \begin{displaymath}{1\over 2} \partial_{\rm t}\overline{{u_{\phi}{'}}^{2}} = - {...
...i{'})} - {\nu \over 2} \overline{(\partial_r u{'}_{\phi})^2} ,
\end{displaymath} (8)

 \begin{displaymath}{1\over 2} \partial_{\rm t}\overline{{u_{z}{'}}^{2}} = -\over...
...} u_z{'})} - {\nu \over 2} \overline{(\partial_r u{'}_{z})^2}.
\end{displaymath} (9)

The total turbulent kinetic energy, defined as

\begin{displaymath}k = {1\over 2} \sum_{i=r,\phi,z} \overline{{u_{i}{'}}^{2}},
\end{displaymath} (10)

then obeys the following equation

 \begin{displaymath}\partial_{\rm t} k = - \overline{u_{r}{'}u_{\phi}{'}} r \part...
...al_r u_r{'} u_i{'}} - \nu \overline{(\partial_i u_i{'})^{2}} .
\end{displaymath} (11)

The first term of the rhs in Eqs. (7) and (8) represents the production of turbulence. We will also refer to them as "coupling" terms, as they couple the turbulent fluctuations to the mean flow. The second term of the rhs in the same equations and the first one of Eq. (9) describe the diffusion due to turbulence. The last term of Eqs. (7)-(9) is the viscous dissipation.

2.2 About linearized equations

In this section we will be considering the effect of neglecting the non-linear terms (i.e. turbulent diffusion) in Eqs. (7)-(9) and (11). There are several reasons why one would neglect the non-linearities of the hydrodynamic equations. The first one is to simplify the system by neglecting terms that are difficult to treat in an analytical stability analysis. Other reasons might include the belief that the physical system is such that the non-linearities actually have a negligible effect on its local or global evolution. Neglecting a priori such terms based on their local effect is rather difficult, considering that little is know on their actual amplitude which does depend closely on the intrinsic nature of the turbulence. Another option is to consider spatially averaged equations and make some hypothesis on the boundary conditions of the systems. For accretion disks, it has been argue that the average over the whole flow or part of the flow of the turbulent diffusion terms eventually vanishes. Making this hypothesis is equivalent to saying that there is no energy flux through the boundaries of the system. This is valid for an isolated system, but it seems unlikely for a disk, as it would not allow for energy transfer with its surroundings, in particular the central star. It does not seem realistic just by considering the mass flux between the two objects. It is likely that such energy transfer occurs through the boundary layer between the star and the disk. Imposing that condition also formally implies that the turbulent stress tensor vanishes on all boundaries, meaning that there is no turbulent viscosity at the edges of the disk. Furthermore, the radial velocity also has to vanish at the boundaries. We could also consider the properties of the boundaries beyond which the system cannot be considered as a continuous medium, hence cease to be described accurately by the equations of hydrodynamics. The description of the interactions between the fluid component of the disk and its outer non-continuous part is not a straightforward task. Our last argument on this matter is to point that these turbulent diffusion terms are the ones that are modeled by the anomalous turbulent viscosity. Neglecting them is equivalent to making the assumption that there is no turbulent viscosity, no enhanced transport. It is then an inevitable result that starting with such a truncated set of equations, one reaches the conclusion that there is no such transport.

Considering only the linear terms of the turbulent transport equations will give us insights - by definition - only on the linear stability of the flow. The linearized set of equations is:

                                 $\displaystyle {1\over 2} \partial_{\rm t} \overline{{u_{r}{'}}^{2}}$ = $\displaystyle 2 \Omega \; \overline{u_{r}{'}u_{\phi}{'}}$ $\displaystyle - \nu \left( \overline{(\partial_r u{'}_{r})^2} + {1 \over 2} \ov...
...al_r u{'}_{\phi})^2} + {1 \over 2} \overline{(\partial_r u{'}_{z})^2 } \right),$ (12)
$\displaystyle {1\over 2} \partial_{\rm t}\overline{{u_{\phi}{'}}^{2}}$ = $\displaystyle - {\partial_{r} {\Omega r^{2}} \over r} \; {\overline{u_{r}{'}u_{\phi}{'}} } - {\nu \over 2} \overline{(\partial_r u{'}_{\phi})^2},$ (13)
$\displaystyle {1\over 2} \partial_{\rm t}\overline{{u_{z}{'}}^{2}}$ = $\displaystyle - {\nu \over 2} \overline{(\partial_r u{'}_{z})^2},$ (14)
$\displaystyle \partial_{\rm t} k$ = $\displaystyle - \overline{u_{r}{'}u_{\phi}{'}} r \partial_r \Omega - \nu \overline{(\partial_i u_i{'})^{2}}.$ (15)

From Eq. (15) we conclude that a stationary turbulent state ( $\partial_{\rm t} k = 0$) can exist only if the following condition - depending on the sign of the angular velocity gradient - is satisfied
$\displaystyle \partial_{r} {\Omega } > 0 \;\;$ $\textstyle \Rightarrow$ $\displaystyle \;\; {\overline{u_{r}{'}u_{\phi}{'}} } < 0,$  
$\displaystyle \partial_{r} {\Omega } < 0 \;\;$ $\textstyle \Rightarrow$ $\displaystyle \;\; {\overline{u_{r}{'}u_{\phi}{'}} } > 0.$ (16)

For Rayleigh unstable flows (where $\partial_{r} {\Omega r^{2}} < 0$ and therefore $\partial \Omega <0$), Eqs. (12), (13) and (16) show that the linear coupling terms with the mean flow are always a source for turbulent fluctuations, a property that reflects the global linear instability of this class of flows. The non-linearities are needed there only for the saturation of amplitudes of the radial and azimuthal velocities along with the redistribution of energy towards axial motions of the new bifurcated flow. In the case of stable angular momentum stratification ( $\partial_{r} {\Omega r^{2}} > 0$), the first order coupling terms have opposite signs in Eqs. (7) and (8). This means that one of them is an energy sink for one of the component of the velocity fluctuations (Balbus et al. 1996). It then looks impossible to allow for the growth or to maintain turbulence as long as the angular momentum is increasing outward. This actually reflects the linear stability of such flows.

3 Considerations on the non-linear stability

In the following sections, will focus on the case of linearly stable flows (i.e. where $\partial_{r} {\Omega r^{2}} > 0$) and consider the energy equations including the turbulent diffusion.

3.1 Inviscid flow

We first consider the case where dissipation due to molecular viscosity can be neglected. We define u the characteristic turbulent velocity (e.g. the root-mean-square velocity) and $\lambda$ the length scale characteristic of the spatial variations of the velocity fluctuations, and we pose

 \begin{displaymath}{1\over 2} \partial_{\rm t} \overline{{u_{r}{'}}^{2}}
= 2 \Omega \; C^l_{r\phi} \; u^2 + C^{n}_r {u^3 \over \lambda} ,
\end{displaymath} (17)

 \begin{displaymath}{1\over 2} \partial_{\rm t}\overline{{u_{\phi}{'}}^{2}} = - {...
...over r} \; C^l_{r\phi} \; u^2+ C^{n}_\phi {u^3 \over \lambda},
\end{displaymath} (18)

 \begin{displaymath}{1\over 2} \partial_{\rm t}\overline{{u_{z}{'}}^{2}} = C^{n}_z {u^3 \over \lambda},
\end{displaymath} (19)

 \begin{displaymath}\partial_{\rm t} k = - r \partial_r \Omega \; C^l_{r\phi} u^2 + C^{n} {u^3 \over \lambda}\cdot
\end{displaymath} (20)

We have introduced the correlation coefficients Crnl, $C_\phi^{n}$, Czn, $C_{r\phi}^l$ and $C^{n} = C_r^{n}+ C_\phi^{n} +C_z^{n}$. By imposing simple constraints, we can deduce the signs of these parameters, which depend on the sign of the angular velocity gradient. In the case of a negative gradient for example, we know from Eq. (15) that the correlation product $\overline{u_{r}{'}u_{\phi}{'}}$ has to be positive, which translates into $C^l_{r\phi}>0$. It also implies that a stationary state can be reach only if Cnr is positive and $C^{n}_\phi$ is negative. Note that Cnz must be zero in that case because we neglected the viscous dissipation. These properties for all introduced coefficients are summarized in Table 1. From Eq. (20) we derive the amplitude of the fluctuations vorticity:

 \begin{displaymath}{u \over \lambda } = {C_{r\phi}^l \over C^{n}} r \partial_r \Omega.
\end{displaymath} (21)

This relation express that the vorticity extracted from the background flow by the fluctuations is proportional to the local shear.

\end{figure} Figure 1: Stability diagram of Couette-Taylor flow in the ( $\Omega _{\rm o}$, $\Omega _{\rm i}$) plane where $\Omega _{\rm o}$ and $\Omega _{\rm i}$ are the angular velocities (in Hz) of the outer and inner cylinder respectively; solid line: solid body rotation; dotted lines: viscous stability criteria; dashed lines: inviscid stability criteria (see text for detail). On the right (resp. left) of the solid body rotation curve, the angular velocity gradient is positive (resp. negative). The flows within the dark surface are Rayleigh unstable. The dots are the experimental instability boundaries within which the flows exhibit self-sustained turbulence; The lined surfaces stand for the predicted non-linearly instability. Experimental results are from Richard et al. (2002).


Table 1: Correlation coefficients properties (Inviscid case).
$\partial_r \Omega$ $C_{r\phi}^l$ Crn $C_\phi^{n}$ Czn Cn
>0 <0 >0 <0 0 <0
<0 >0 <0 >0 0 <0

We have seen that there always exists a sink term for one of the components of the energy fluctuations. A necessary condition for the existence of self-sustained turbulence is then that the non-linear terms overcome this negative production term (for the azimuthal component when $\partial_r \Omega <0$, for the radial one when $\partial_r \Omega >0$). Namely, from Eqs. (17) and (18)

                     $\displaystyle \partial_r \Omega > 0\!:$   $\displaystyle \;\; C^{n}_r {u^3 \over \lambda} \geqslant - 2 \Omega \; C^l_{r\phi} \; u^2,$ (22)
$\displaystyle \partial_r \Omega < 0\!:$   $\displaystyle \;\; C^{n}_\phi {u^3 \over \lambda} \geqslant {\partial_{r} {\Omega r^{2}} \over r} \; C^l_{r\phi} \; u^2,$ (23)

which reduces to,
              $\displaystyle \partial_r \Omega > 0\!:$   $\displaystyle \;\; {u \over \lambda} \geqslant \; - {C_{r\phi}^l \over C_r^{n}} \; 2 \Omega,$ (24)
$\displaystyle \partial_r \Omega < 0\!:$   $\displaystyle \;\; {u \over \lambda} \geqslant \; {C_{r\phi}^l \over C_\phi^{n}}\; {\partial_r \Omega r^2 \over r},$ (25)

Finally, using Eq. (21), we obtain the following relations,
                      $\displaystyle \partial_r \Omega > 0\!:$   $\displaystyle \;\; {C_{r\phi}^l \over C^{n}} r \partial_r \Omega \geqslant \; - {C_{r\phi}^l \over C_r^{n}} \; 2 \Omega \;,$ (26)
$\displaystyle \partial_r \Omega < 0\!:$   $\displaystyle \;\; {C_{r\phi}^l \over C^{n}} r \partial_r \Omega \geqslant \; {C_{r\phi}^l \over C_\phi^{n}}\; {\partial_r \Omega r^2 \over r}\cdot$ (27)

Introducing the dimensionless Rossby number, $Ro = r \partial_r \Omega / 2 \Omega$, and using the property $C^{n} = C_r^{n} + C_\phi^{n}$ we obtain the necessary conditions for instability:
                           $\displaystyle \partial_r \Omega > 0\!:$   $\displaystyle \;\; Ro \geqslant - {C^{n} \over C_r^{n}} = Ro_{\rm c} \;\;\;\;\;\;\; ( Ro_{\rm c} > 0 ) ,$ (28)
$\displaystyle \partial_r \Omega < 0\!:$   $\displaystyle \;\; Ro \leqslant - {C^{n} \over C_r^{n}}= Ro_{\rm c} \;\;\;\;\;\;\; ( Ro_{\rm c} < 0 ) .$ (29)

Here we have introduced the parameter $Ro_{\rm c}$. A similar result can be derived by simply posing $\partial_{\rm t}<{u_{i}{'}}^{2}> ~= 0$ and reducing Eqs. (17) and (19), which leads to Ro = - Cn / Crn. The quantity  $Ro_{\rm c}$ can be seen as constant by assuming that the ratio Crn/Cn is independent from the mean flow. This ratio quantify the redistribution of the energy extracted for the mean flow between the velocity components, by the non-linearities. Note that the coefficients Crn, $C_\phi^{n}$,Czn and Cn themself most likely depend on the mean flow (see Younis et al. 1999). Within this picture, $Ro_{\rm c}$ is a constant critical stability parameter, and we expect the flow to become turbulent for either $Ro > Ro_{\rm c}$ (when $\partial_r \Omega >0$) or $Ro < Ro_{\rm c}$ (when $\partial_r \Omega <0$). Townsend (1976) already noted that the Rossby number should approach a constant value when a linearly stable rotating shear flow becomes turbulent, by arguing that "it is conceivable and even likely that an asymptotic state can be reached with a constant ratio of transfer from the transverse motion to total energy transfer from the mean flow". The two actors in the balance are the gradient of angular velocity, shearing the velocity fluctuations, creating small scale vorticity, feeding the turbulent cascade, and the gradient of angular momentum, stabilizing the flow by damping the fluctuations through a "spring" mechanism. The energy extracted from the background flow into small scale vorticity must reach a critical value in order to overcome the constraint introduced by the stable gradient of angular momentum. The Rossby number, as a measure of the ratio between the shear and the angular momemtum gradient, appears a posteriori as a natural control parameter of the flow stability.

3.2 Viscous flow

Following the same path as in the inviscid case, we now add the viscous constraints to Eqs. (17)-(20):

 \begin{displaymath}{1\over 2} \partial_{\rm t} \overline{{u_{r}{'}}^{2}} = 2 \Om...
... {u^3 \over \lambda} - \nu C^{\nu}_{r} {u^2 \over \lambda^2} ,
\end{displaymath} (30)

 \begin{displaymath}{1\over 2} \partial_{\rm t}\overline{{u_{\phi}{'}}^{2}} = - {...
...u^3 \over \lambda} - \nu C^{\nu}_{\phi} {u^2 \over \lambda^2},
\end{displaymath} (31)

 \begin{displaymath}{1\over 2} \partial_{t}\overline{{u_{z}{'}}^{2}} = C^{n}_z {u^3 \over \lambda} - \nu C^{\nu}_{z} {u^2 \over \lambda^2} ,
\end{displaymath} (32)

 \begin{displaymath}\partial_{\rm t} k = - r \partial_r \Omega \; C^l_{r\phi} u^2 + C^{n} {u^3 \over \lambda}.- \nu C^{\nu} {u^2 \over \lambda^2}
\end{displaymath} (33)

where $C^{\nu} =C^{\nu}_{r}+C^{\nu}_{\phi}+C^{\nu}_{z}$. The properties of the correlation factors are reported in Table 2. The extraction of vorticity from the mean flow now becomes:

 \begin{displaymath}{u \over \lambda } = {C_{r\phi}^l \over C^{n}} r \partial_r \Omega + {\nu \over \lambda^2} {C^\nu \over C^n},
\end{displaymath} (34)

and the relations (24) and (25) translate into
                      $\displaystyle \partial_r \Omega > 0\!:$   $\displaystyle \;\; {u \over \lambda} \geqslant \; - {C_{r\phi}^l \over C_r^n} \; 2 \Omega + {\nu C^{\nu}_{r} \over C_r^n \lambda^2} \;,$ (35)
$\displaystyle \partial_r \Omega < 0\!:$   $\displaystyle \;\; {u \over \lambda} \geqslant \; {C_{r\phi}^l \over C_\phi^n}\...
...tial_r \Omega r^2 \over r} + {\nu C^{\nu}_{\phi} \over C_\phi^n \lambda^2}\cdot$ (36)

and after some rearrangement we finally obtain,
                 $\displaystyle \partial_r \Omega > 0\!:$   $\displaystyle \;\; Ro \geqslant Ro_{\rm c} \left( 1-{1 / Re^{+}} \right) ,$ (37)
$\displaystyle \partial_r \Omega < 0\!:$   $\displaystyle \;\; Ro \leqslant Ro_{\rm c} \left( 1+{1 / Re^{-}} \right) ,$ (38)

where we have introduced the following quantities,

\begin{displaymath}Re^+ = { 2 \Omega \lambda^2 \over \nu } (C_r^\nu /C_{r \phi}^l -{C^\nu C_r^n / C^n C_{r \phi}^l})^{-1} < 0 ,
\end{displaymath} (39)

\begin{displaymath}Re^- = { 2\Omega \lambda^2 \over \nu } ( C_\phi^\nu / C_{r \phi}^l-{C^\nu C_\phi^n / C^n C_{r \phi}^l})^{-1} > 0 .
\end{displaymath} (40)

The correction due to viscous dissipation lies in the 1 / Re+ and  1 / Re- terms. The necessary condition for instability reduces to the results found for the inviscid case, if $\nu$ is set to zero. For high rotation rates, the critical values of the Rossby number will eventually tend to the inviscid values.


Table 2: Correlation coefficients properties (Viscous case).
$\partial_r \Omega$ $C_{r\phi}^l$ Crn $C_\phi^{n}$ Czn Cn $C^\nu_i$ $C^\nu$
>0 <0 >0 <0 >0 <0 >0 >0
<0 >0 <0 >0 >0 <0 >0 >0

Considering all of the above, we expect that there exists a critical Rossby number for each class of flows - namely flows with $\partial_r \Omega >0$ and flows with $\partial_r \Omega <0$. This result can be verified by using experimental results on differentially rotating flows. This can be done by examining the value of the Rossby number at the onset of shear turbulence observed in laboratory Couette-Taylor experiment. The Couette-Taylor experiment consist of two coaxial cylinder between which the fluid is sheared, resulting in a differentially rotating azimuthal flow. Experimental results on stability (Richard et al. 2002) can actually be fitted with the values $Ro_{\rm c}=1$ for $\partial_r \Omega >0$, $Ro_{\rm c} = -1/2$, for $\partial_r \Omega <0$, $Re^+=-\Omega / 3$ and $Re^-=2 \Omega$. This is consistent with the picture of a constant critical Rossby number with a viscous correction vanishing for high rotation rates or low viscosity.

4 Turbulent transport and viscosity

Radial turbulent transport of angular momentum is quantified by the second-order correlation product $\overline{u_{r}{'}u_{\phi}{'}}$. The time evolution for this quantity after some rearrangement and spatial averaging over the azimuthal direction $\phi$, and the introduction of correlation coefficients, is given by:

 \begin{displaymath}\partial_{\rm t} \overline{u_{r}{'}u_{\phi}{'}}
=- C_{r\phi}...
...bda} - \nu C_{r \phi}^\nu {u^2 \over \lambda^2} + \Pi_{r\phi},
\end{displaymath} (41)

where $\Pi_{ij}$ is the pressure-strain correlation tensor. Following Kato & Yoshizawa (1997), in the case of isotropic turbulence, and with their notations, we have

\begin{displaymath}\Pi_{r\phi} = \frac{3}{2} C_2 {u{'}}^2 r \partial_r \Omega .
\end{displaymath} (42)

We can derive a formal expression for the turbulent diffusion of momentum, using the classical definition,

\begin{displaymath}\nu_{\rm t} = {\overline{u_{r}{'}u_{\phi}{'}} \over r \partial_r \Omega} = {C_{r\phi}^l u^2 \over r \partial_r \Omega},
\end{displaymath} (43)

For steady turbulence, in the limit where molecular viscosity in negligible compared to turbulent transport, from Eq. (41),

 \begin{displaymath}{\nu_{\rm t} \over \nu} = C_{r\phi}^l \left({{C_{r\phi}^l - (...
...ambda \over r}\right)^2 { r^3 \partial_r \Omega\over \nu}\cdot
\end{displaymath} (44)

Richard & Zahn (1999) derived an expression for the turbulent viscosity from Taylor (1936) and Wendt (1933) experimental data,

\begin{displaymath}{\nu_{\rm t} \over \nu} = \beta \; { \big\vert r^3 \partial_r \Omega \big\vert\over \nu} = \beta Re^*
\end{displaymath} (45)

where $\beta$ was found to be independent from the geometry and the background flow (see also Longaretti 2002 and Duschl et al. 2000 for other arguments in favor of this prescription). Identifying this expression with Eq. (44), one finds,

\begin{displaymath}\beta =\left\vert C_{r\phi}^l \left({{C_{r\phi}^l - (3/2) C_2...
...}}\right)^2 \left({\lambda \over r}\right)^2\right\vert \cdot
\end{displaymath} (46)

Richard & Zahn (1999) also identified the Re* parameter (named "gradient Reynolds number") as the control parameter for stability. This derivation was acheived by identifying the values of the Reynolds number for which the measured torque from laboratory experiments starts to differ from its theoretical value for the laminar Couette flow. For small gap Couette-Taylor apparatus, the classical Reynolds number defined as $Re = \Delta \Omega \Delta R R / \nu$ (where $\Delta R$ is the gap between the cylinders, $\Delta \Omega$ their differential angular velocity, and R the mean radius) was identified as the relevant control parameter. For gaps bigger than $\Delta R / R \gtrsim 1/20$ the critical Reynolds number behave like $(\Delta R / R )^2$, which means that Re* takes a constant value. To be more proper, this is a measure of the efficiency of the turbulent transport. The bifurcation identified from the torque data is then not directly linked to the stability of the flow, but is the threshold for which turbulent transport becomes appreciably more efficient than molecular viscosity. It should be read as the following condition

\begin{displaymath}{\nu_{\rm t} \over \nu} >> 1 \Rightarrow Re^* >> {1 \over \beta} \cdot
\end{displaymath} (47)

This is consistent with the numerical values from Richard & Zahn (1999) ( $Re^*_{\rm c} = 1.5 \times 10^{-5}$ and $\beta^{-1} = 1.5 \times 10^{-6}$).

5 Geometrical constraints

5.1 General case

In the case of a bounded flow, keeping up the assumption of isotropic turbulence, we have to assume that $\lambda$ will be of order of a macro-scale d. The stability criteria for a viscous flow is modified as $\lambda$ is constrained to smaller values than for a free shear flow (therefore being closer to the inviscid values of $Ro_{\rm c}$, according to Eqs. (37) and (38)) while for the inviscid case it remains unchanged. The radial transport will also be modified, as we have, according to Eq. (44)

 \begin{displaymath}{\nu_{\rm t} \over \nu} \propto d^2 { \vert r {\partial_r \Omega}\vert\over \nu}\cdot
\end{displaymath} (48)

5.2 Narrow gap Couette-Taylor experiment

This expression can be approximated in the case of narrow-gap Couette-Taylor experiment by taking for d the distance between the two cylinders $\Delta R$.

 \begin{displaymath}{\nu_{\rm t} \over \nu} \propto {\Delta \Omega \cdot \Delta R \cdot R \over \nu} = Re
\end{displaymath} (49)

The bifurcation threshold identified on torque measurement aforementioned will then scale as the classical Reynolds number, in agreement with the results from Richard & Zahn (1999).

6 Discussion

We have shown that a simple model of turbulent fluctuations energetics can explain laboratory experiment results on differentially rotating flow. The equations of motion show that shear instabilities can develop in rotating flows, as long as the non-linearities are taken into account. Making reasonable assumptions about the necessary properties of the velocity correlation products, we derived formal expressions for the control parameters for stability. Comparison with laboratory results confirms that the relevant parameter is the Rossby number. In the case of a viscous flow, the critical Rossby numbers exhibit a $\Omega^{-1}$ correction, therefore converging toward the same constant critical Rossby numbers as in the inviscid case, for large rotation rates. The turbulent viscosity inferred is consistent with the prescription previously proposed by Richard & Zahn (1999). The critical amplitude of perturbations are also consistent with linear stability theory results in the inviscid limit. We would like to stress the importance of the pressure fluctuations in the growth and self-sustaining of this kind of turbulence. Even though the time average of the related terms vanish from the velocity fluctuation equations, they can be seen as the main engine for isotropizing the turbulence, hence disappearing in time average when isotropy is acheived.

Complementary work is needed, in the laboratory or numerically, in order to seek whether the two particular values or the critical Rossby number matching the experimental results on the Couette-Taylor experiment are genuine to that system or are of more "universal" value.


The author would like to thank J.-P. Zahn, O. Dauchot, B. Dubrulle, F. Daviaud, S. Davis, J. Cuzzi, R. Young, A. Steinacker, F. Hersant and J.-M. Huré for stimulating discussions and support. This work was supported by the Programme National de Physique Stellaire (PNPS), the Commissariat à l'Énergie Atomique and by a Research Associateship from the National Research Council/National Academy of Sciences.


7 Online Material

Appendix A: Turbulent fluctuations equations

The turbulent transport equations are given here in cylindrical coordinates ($r,\phi,z$). Viscous diffusion is omitted

                                          $\displaystyle \partial_{\rm t}\overline{{u_{r}{'}}^{2}}
..._{r}{'}}^{2}}} - 2 \overline{u_{\phi} \over r}
\overline{u_{r}{'}u_{\phi}{'}} =$  
    $\displaystyle \quad-2 \left(
\overline{{u_{r}{'}}^{2}} {\partial_{r} \overline{...
...} \right) +
\overline{u_{r}{'}u_{z}{'}} {\partial_{z} \overline{u_{r}}} \right)$  
    $\displaystyle \quad- 2 \overline{ u_{r}{'} \left( {1 \over r}
\partial_{r}(r {u...
...phi}{'}) -{{u_{\phi}{'}}^{2}\over r}
+ \partial_{z}(u_{r}{'}u_{z}{'}) \right) }$  
    $\displaystyle \quad+2 {1\over \rho} \overline{p{'} \partial_{r}u_{r}{'}} - 2 {1 \over
\rho}\overline{ \partial_{r}(p{'} u_{r}{'}) }$  
    $\displaystyle \quad- \nu \Bigg( 2 \overline{\left( \partial_{r} u_{r} {'} \righ...
..._{\phi} {'} \right)^{2}} + \overline{\left( \partial_{r} u_{z} {'} \right)^{2}}$  
    $\displaystyle \quad +\overline{\left( {1 \over r}\partial_{\phi} u_{r} {'} - { ...
...ver r}\right)^{2}} +\overline{\left( \partial_{z} u_{r} {'} \right)^{2}} \Bigg)$ (A.1)

                                        $\displaystyle \partial_{\rm t}\overline{{u_{\phi}{'}}^{2}}
+ \overline{u_{r}}{\...
+ 2 \overline{u_{\phi} \over r} \overline{u_{r}{'}u_{\phi}{'}} =$  
    $\displaystyle \quad- 2 \left( \overline{u_{r}{'}u_{\phi}{'}}
{\partial_{r} \ove... + \overline{u_{\phi}{'}u_{z}{'}}
{\partial_{z} \overline{u_{\phi}}} \right)$  
    $\displaystyle \quad- 2 \overline{
u_{\phi}{'} \left(
{1 \over r} \partial_{\phi...
+ 2 {u_{r}{'}u_{\phi}{'} \over r} + \partial_{z}u_{\phi}{'}u_{z}{'}
    $\displaystyle \quad+2 {1 \over \rho} \overline{ p{'} \left( {1 \over r}\partial...
...r \rho} \overline{ \left( {1 \over r}\partial_{\phi}(p{'} u_{\phi}{'}) \right)}$  
    $\displaystyle \quad- \nu \Bigg( \overline{\left( \partial_{r} u{'}_{\phi}\right...
...hi} u{'}_{z}\right)^{2}} +\overline{\left( \partial_{z} u{'}_{\phi}\right)^{2}}$  
    $\displaystyle \quad +2\overline{\left( {1 \over r} \partial_{\phi} u{'}_{\phi} ...
...1 \over r} \partial_{\phi} u{'}_{r} + {u{'}_{\phi}\over r} \right)^{2}}} \Bigg)$ (A.2)

                                       $\displaystyle \partial_{\rm t}\overline{{u_{z}{'}}^{2}}
+ \overline{u_{r}}{\par...
+ \overline{u_{z}}{\partial_{z} \overline{{u_{z}{'}}^{2}}}
    $\displaystyle \quad-2
\overline{u_{r}{'} u_{z}{'}} {\partial_{r} \overli...
...}} \over r}
+ \overline{{u_{z}{'}}^{2}} {\partial_{z} \overline{u_{z}}}
    $\displaystyle \quad- 2 \overline{
u_{z}{'} \left(
{1 \over r} \partial_{r}(r u_...
...r r}\partial_{\phi}(u_{\phi}{'}u_{z}{'})
+ \partial_{z}{u_{z}{'}}^{2}
    $\displaystyle \quad+2 {1 \over \rho} \overline{p{'} \partial_{z}u_{z}{'}}
- 2 {1 \over \rho}\overline{ \partial_{z}(p{'} u_{z}{'}) }$  
    $\displaystyle \quad- \nu \Bigg( 2 \overline{\left( \partial_{z} u{'}_{z}\right)...
...{\phi} u{'}_{z}\right)^{2}} +\overline{\left( \partial_{r} u{'}_{z}\right)^{2}}$  
    $\displaystyle \quad+ \overline{\left( \partial_{z} u{'}_{r} \right)^{2}} +{\overline{\left( \partial_{z} u{'}_{\phi} \right)^{2}}} \Bigg)$ (A.3)

Copyright ESO 2003