Table of Contents
ISRN Astronomy and Astrophysics
Volume 2012, Article ID 987275, 10 pages
Research Article

On Estimating Fluxes due to Small-Scale Turbulent Convection in a Rotating Star

Institute of Astronomy & Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK

Received 30 November 2011; Accepted 15 January 2012

Academic Editors: R. Avila and S. Bogovalov

Copyright © 2012 D. O. Gough. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The way in which turbulent fluxes are usually represented in computations of large-scale flow in the convection zones of the sun and other stars is briefly described. A model of an ensemble of eddies that is capable of generalization to circumstances more complicated than the usual essentially spherically symmetrical convection zone is outlined. Generalization usually requires the introduction of new postulates, and, in so doing, also lays bare some of the assumptions, often implicit, in the usual mixing-length formalisms.

1. Introduction

Computations of large-scale flow in stellar convection zones, be it a meridional circulation or simply rotation, require some representation of the fluxes of heat and momentum by turbulent motion on scales too small to be resolved. It is not uncommon in astrophysics to adopt the attitude that since the theory upon which any such representation is based (usually mixing-length theory) cannot be trusted, it is hardly necessary to exercise care in deriving the formulae for the fluxes. Readers can thus be faced with the prospect of learning of the results of complicated numerical computations, yet being unable to know precisely what those results mean.

The main point I wish to make is that the situation can be improved. The so-called theory used for calculating the turbulent fluxes can be refined, and the numerical machinery that has been built to model the large-scale flow can be used experimentally to test the refinements.

First I shall describe the broad principles behind many of the computations of turbulent fluxes in astrophysics. Then I shall outline a procedure whereby the representations of the small-scale motion might be developed in a coherent manner.

2. The Turbulent Fluxes

The first step in all computations of large-scale astrophysical flow, whether it is stated explicitly or not, is to imagine a separation of scales. Motion with characteristic length and time scales greater than some value, usually implicitly prescribed by the numerical technique, is treated by solving the governing equations directly, using some parametrized approximation to describe the motion with characteristic scales smaller than this value: ignoring the small-scale motion completely is a trivial form of such a parametrization. The large-scale fields of velocity and temperature, say, may be considered to derive from the total fields v and 𝑇 via the application of an averaging procedure, which I denote by an overbar. Thus𝐯=𝐯+u,𝑇=𝑇+𝜗.(1) I shall refer to fields such as 𝐯 and 𝑇 as mean quantities; the fields 𝐮 and 𝜗 are fluctuations. In this account, I have in mind particularly the convective envelopes of cool stars, in which turbulent motion is driven predominantly by buoyancy.

For simplicity of presentation I refer the fields to rectangular Cartesian coordinates 𝑥𝑖, with 𝑥3 vertical, and I shall ignore the curvature of the convection zone. This is a fair approximation if the scales to be parametrized are small enough. Equations of motion may then be written, in an obvious notation,𝜕𝜕𝑡𝜌𝑣𝑖+𝜕𝜕𝑥𝑗𝜌𝑣𝑖𝑣𝑗=𝜕𝑝𝜕𝑥𝑖𝑔𝜌𝛿3𝑖𝜕,(2)1𝜕𝑡𝜌𝜖+2𝜌𝑣𝑖𝑣𝑖+𝜕𝜕𝑥𝑖𝐹𝑖+𝜌𝑣𝑖+12𝜌𝑣𝑗𝑣𝑗𝑣𝑖=𝑔𝜌𝑣3,(3)𝜕𝜌+𝜕𝜕𝑡𝜕𝑥𝑖𝜌𝑣𝑖=0,(4) where 𝜌 and 𝑝 are density and pressure, 𝜖 and are specific internal energy and enthalpy (viz., internal energy and enthalpy per unit mass), 𝑔 is the magnitude of the gravitational acceleration 𝑔𝑖, 𝐹𝑖 is the radiative heat flux, 𝑡 is time, and 𝛿𝑖𝑗 is the Kronecker delta. Throughout my discussion I adopt the summation convention for repeated indices except where one of them is enclosed in parentheses. The total energy equation (3) is derived by combining the kinetic-energy equation, obtained by taking the (scalar) product of (2) with 𝑣𝑖 to relate the rate of change of the kinetic energy to the rate of working, with the thermal-energy equation𝜌𝑐𝑣𝜕𝑇𝜕𝑡+𝑣𝑖𝜕𝑇𝜕𝑥𝑖+𝑝𝜕𝑣𝑖𝜕𝑥𝑖+𝜕𝐹𝑖𝜕𝑥𝑖=0,(5) which is a statement of the first law of thermodynamics. In (5), 𝑐𝑣 is the specific heat at constant volume. These equations must, of course, be supplemented by an equation of state and an equation determining 𝐹𝑖. Viscous terms have been ignored; so too have perturbations to 𝑔.

To obtain equations for the large-scale flow the averaging procedure is applied to the full equations of motion. Usually some approximations are made in representing the fluctuations, the most common being the neglect of density fluctuations, except where they appear in the enthalpy flux. This constitutes part of the Boussinesq approximation, which has been justified by Spiegel and Veronis [1] for the case in which 𝑣𝑖=0, when the fluctuations are generated by buoyancy alone, and when the length scales of the fluctuations are all much smaller than the scale heights of density and pressure of the averaged state. Despite the fact that these conditions are believed not to be satisfied in most stars, I too shall nevertheless adopt the practice, obtaining𝜕𝜕𝑡𝜌𝑣𝑖+𝜕𝜕𝑥𝑗𝜌𝑣𝑖𝑣𝑗+𝜕𝑅𝑖𝑗𝜕𝑥𝑗𝜕=𝑝𝜕𝑥𝑖𝑔𝜌𝛿3𝑖,𝜕𝜕𝑡𝜌1𝜖+2𝜌𝑣𝑗𝑣𝑗+𝜕𝜕𝑥𝑖𝐹𝑖+𝜌𝑣𝑖+𝐹𝑐𝑖+12𝜌𝑣𝑗𝑣𝑗+𝑢𝑗𝑢𝑗𝑣𝑖+𝐾𝑡𝑖=𝑔𝜌𝑣3,𝜕𝜌+𝜕𝜕𝑡𝜕𝑥𝑖𝜌𝑣𝑖=0.(6) The only fluxes arising from the small-scale turbulent motion that need to be calculated for closing these equations are evidently the Reynolds stress 𝑅𝑖𝑗=𝜌𝑢𝑖𝑢𝑗, the turbulent heat flux 𝐹c𝑖=𝜌𝑢𝑖, where the prime denotes a turbulent fluctuation, and the turbulent kinetic-energy flux 𝐾𝑡𝑖=(1/2)𝜌𝑢𝑗𝑢𝑗𝑢𝑖.

In the absence of a large-scale mean flow the equations for the fluctuations are given in Boussinesq approximation by [1]𝜌𝜕𝑢𝑖𝜕𝑡+𝑢𝑗𝜕𝑢𝑖𝜕𝑥𝑗=𝜕𝑝𝜕𝑥𝑖+𝑔𝜌𝛿𝑇𝜗𝛿𝑖3,(7)𝜕𝜗𝜕𝑡+𝑢𝑖𝜕𝜗𝜕𝑥𝑖𝛽𝑖𝑢𝑖1=𝜌𝑐𝑝𝜕𝐹𝑖𝜕𝑥𝑖,(8)𝜕𝑢𝑖𝜕𝑥𝑖=0,(9) where 𝑐𝑝 is the specific heat at constant pressure and𝛽𝑖𝛿𝜌𝑐𝑝𝜕𝑝𝜕𝑥𝑖𝜕𝑇𝜕𝑥𝑖(10) is the superadiabatic lapse rate; 𝛿=(𝜕ln𝑝/𝜕ln𝜌)𝑇 is the inverse dimensionless isothermal compressibility. In what follows I shall assume that the perturbations are optically thick, so, ignoring, for simplicity, fluctuations in the thermal conductivity,𝐹𝑖=𝐾𝜕𝜗𝜕𝑥𝑖,(11) where 𝐾=4𝑎𝑐𝑇3/3𝜒𝜌, in which 𝑎 is the first radiation constant, 𝑐 is the speed of light, and 𝜒 is the Rosseland mean opacity. Various procedures have been suggested for treating optically thin perturbations (e.g., [25]); I do not discuss them here other than to suggest the procedure recommended by Gough [6], which is consistent with the procedure adopted here for dealing with the other fluctuation equations and which, unlike the other procedures, provides a smooth transition from optically thick to optically thin fluctuations.

3. Rough Estimates of the Turbulent Fluxes

In most cases some kind of diffusion approximation based on mixing-length theory is used to compute the turbulent fluxes. Let us consider first the convective heat flux. Since the Boussinesq approximation has already been used to derive (6), one might as well continue doing so and, in particular, neglect the contribution to the convective heat flux from pressure fluctuations. That flux can then be written𝐹c𝑖=𝜌𝑐𝑝𝜗𝑢𝑖.(12) (It is perhaps worth pointing out that for a perfect unionizing gas the pressure fluctuation makes precisely no contribution to when the latter is expressed in terms of the temperature fluctuation.) In this and all subsequent equations overbars are omitted from mean quantities, except where there is risk of confusion.

If there were no large-scale motion the star would be spherically symmetrical in the mean—I am ignoring the possibility of there being a substantial magnetic field—and only the vertical components 𝛽 and 𝐹𝑐 of the superadiabatic lapse rate and the convective flux would be nonzero. Moreover, the Reynolds stress would be diagonal; I denote its (3,3) component, the only component that enters hydrostatics when the spherical geometry is ignored, by 𝑝𝑡. In that case the simplest mixing-length estimates for the magnitudes of typical velocity and temperature fluctuations can be obtained by first ignoring the time derivatives and also ignoring the role of the pressure gradient in (7), so that the convective motion is considered to be purely vertical, with velocity 𝑤. Then one replaces 𝜕/𝜕𝑥𝑖 by 1, where , called the mixing-length, is some characteristic lengthscale of the largest turbulent eddies. For upwardly directed turbulent motion, this leads to𝑤2𝑔𝛿𝑇𝜗,1𝜗𝛽𝑤𝜅2𝜗,(13) where 𝜅 is the thermal diffusivity: 𝜅=𝐾/𝜌𝑐𝑝. In obtaining these estimates it was assumed that is much less than the scales of variation, 𝐻, of all the other variables appearing in (13) including 𝛽. That assumption is the basis of so-called local theories and allows fluctuations at any point to be expressed in terms of the mean state at that point. It is violated in any of the usual mixing-length models of stellar convective envelopes, for which the calibration of , usually to reproduce the solar radius, leads to a value of /𝐻 in excess of unity (e.g., [7, 8]).

To estimate the average 𝑤𝜗 one simply multiplies the estimates of 𝑤 and 𝜗 from (13), and introduces a scaling factor 𝛼𝐻 to account for the imperfect correlation between velocity and temperature fluctuations, giving𝐹𝑐18𝛼𝐻𝑆14𝑆+113𝐾𝛽,(14) where𝑆=𝑔𝛿𝑇1𝛽4𝜅2(15) is the product of the Prandtl number and a locally defined Rayleigh number based on the lengthscale . The factor 𝛼𝐻 depends on the way in which the mixing-length ideas have been employed to describe the dynamics; in a description of the turbulence more realistic than those assuming that all motion is vertical, it depends also on the supposed geometry of the turbulent eddies. The (3,3) component 𝑝𝑡 of the Reynolds stress tensor 𝑅𝑖𝑗 can be estimated similarly:𝑝𝑡𝛼𝑅𝜌𝑤214𝛼𝑅4𝑆+112𝜌𝜅2=𝛼𝑅11+4𝑆112𝑆1/22𝜌𝑔𝛿𝛽𝑇2,(16) in which 𝛼𝑅 plays a similar role to 𝛼𝐻 in (14). The turbulent kinetic-energy flux is usually ignored; it is proportional to 𝜌𝑤3, with a constant of proportionality that depends on the (uncertain) spatial correlations of the turbulent velocity field.

Deep in stellar convection zones 𝑆 is very large—in the solar convection zone it exceeds 1020. This implies that the coherent convective motion is essentially adiabatic. Microscopic diffusion of heat is negligible, and𝑤𝑔𝛿𝑇𝛽1/2,𝜗𝛽.(17) Then𝐹𝑐𝛼𝐻𝜌𝑐𝑝𝑔𝛿𝑇2𝛽3/2=𝛼𝐻𝑆1/2𝐾𝛽=𝐾𝑐𝑝𝛽,𝑡𝛼𝑅𝑔𝛿𝜌𝑇2𝛽,(18) which, of course, are independent of the thermal conductivity 𝐾.

The formula (18) for 𝐹𝑐 has been written as the product of a gradient 𝛽 and a transfer coefficient 𝐾𝑐. The reason for so doing is that it is this form that has been the most common basis for generalization to situations that are no longer spherically symmetrical in the mean. What is sometimes done (e.g., [9]) is to assume that (18) is the 3-component of a vector equation and that the transfer coefficient 𝐾𝑐 is a scalar, even though it appears to depend on the vectors 𝛽𝑖 and 𝑔𝑖. Thus one obtains𝐹𝑐𝑖=𝐾𝑐𝛽𝑖,(19) with the turbulent heat flux being parallel to the superadiabatic temperature gradient. One might more realistically envisage trying to generalize the formula to something like𝐹𝑐𝑖=𝐾𝑐𝑖𝑗𝛽𝑗,(20) since generally, as will become apparent later, one would not expect 𝐹𝑐𝑖 and 𝛽𝑖 to be parallel.

Accepting (19) exposes the difficulty of principle in deciding how the coefficient 𝐾𝑐 should be computed. Perhaps the simplest hypothesis would be to leave 𝑔 as it stands and consistently replace 𝛽 by the magnitude (𝛽𝑖𝛽𝑖)1/2 of 𝛽𝑖 in the formula (18) for 𝐾𝑐. But in practice even cruder approximations are often adopted, which at best are based on the value of 𝐾𝑐 computed from an initial trial model in which 𝑣𝑖 is constrained to be zero, even though the resulting available heat energy advected by the large-scale flow in the final model may turn out to be comparable with 𝐹𝑐𝑖.

The treatment of the Reynolds stress in the presence of shear is also motivated by the desire to reduce the formula to a linear Fick law. In most cases the usual scalar viscous law is adopted, with a turbulent shear viscosity derived in some way from the product of and a turbulent-velocity estimate such as that in (17).

It is convenient to consider the Reynolds stress tensor 𝑅𝑖𝑗 to be divided into two parts: a part 𝑅1𝑖𝑗 that vanishes when 𝑣𝑖=0 and the part 𝑅0𝑖𝑗 that does not. In the absence of a magnetic field or other perturbing agent, 𝑅0𝑖𝑗 is diagonal and acts as a (generally anisotropic) turbulent pressure 𝑝𝑡𝑖. It is usually ignored in computations of stellar structure. This may be dangerous, because the presence of any anisotropy in 𝑝𝑡𝑖 can modify large-scale flow (e.g., [10, 11]).

In modelling the remaining terms in the Reynolds stress, once again an assumption of localness is usually adopted so that it may be presumed that the mean flow may be replaced by a truncated Taylor expansion about any point of interest. Thus one can write𝜌𝑢𝑖𝑢𝑗=𝑝𝑡(𝑖)𝛿𝑖𝑗+𝑀𝑖𝑗𝑘𝑚𝜕𝑣𝑚𝜕𝑥𝑘+higher-orderterms.(21) There is no term depending on undifferentiated mean velocities 𝑣𝑖 because the formula must be invariant under Galilean transformation.

Gross simplifications of the physics have always been made in attempts to evaluate the viscosity tensor. For example, Wasiutyński [12] assumed fluid elements to travel with constant acceleration, relative to an inertial frame, from initial velocities having uncorrelated components. He obtained a formula which in Cartesian coordinates would have been, had spherical geometrical terms been ignored,𝑀𝑖𝑗𝑘𝑚=𝛿𝑖𝑛𝛿𝑗𝑚+𝛿𝑖𝑚𝛿𝑗𝑛𝑀𝑛𝑘,(22) where 𝑀nk is a symmetrical tensor with principal axes in the east-west, north-south, and vertical directions. Such a form, simplified further by assuming axisymmetry about the vertical, was adopted in some of the early discussions of the sun’s equatorial acceleration (e.g., [10, 12, 13]). The assumptions upon which the derivation is based imply, for example, that the convective velocities are not influenced significantly by the solar rotation. Some investigators assume complete isotropy of the turbulent motion, and no resistance to dilatation of the mean flow, which yields𝑀𝑖𝑗𝑘𝑚=𝛿𝑖𝑘𝛿𝑗𝑚+𝛿𝑖𝑚𝛿𝑗𝑘+𝛿𝑖𝑗𝛿𝑘𝑚𝜂,(23) where 𝜂 is a scalar. This is the ordinary shear viscosity law, which, unlike the anisotropic formulae, has a (no doubt erroneous) tendency to push the mean flow towards a state of rigid rotation.

4. The Eddy-Ensemble Approach

In order to refine the estimates of the turbulent fluxes one must consider more carefully the dynamics of the turbulent eddies. I adopt here an eddy-ensemble approach developed originally for pulsating stars by Gough [6, 14], who subsequently discussed a generalization to accommodate rotating flows [15] which this paper extends. The model representing the turbulent motion is much more carefully defined than in the previous section, and it is correspondingly more complicated to evaluate. Yet in the simple circumstance of convection in a fluid with no mean background flow it yields essentially the same results. The merit in the approach is that it lays bare the underlying assumptions, thereby making it clear how the theory can be extended to more complicated situations and what additional assumptions, if any, are needed to do so.

It is safer and simpler to compute the fluxes directly from the formulae 𝜌𝑢𝑖𝑢𝑗 and 𝜌𝑐𝑝𝜗𝑢𝑖 rather than to attempt an evaluation of transport coefficients such as 𝐾c𝑖𝑗 and 𝑀𝑖𝑗𝑘𝑚 or generalizations of them, as is more common (e.g., [16, 17]). The reason is that (20) and (21) depend on assumptions additional to those of the basic turbulence theory and must therefore predict fluxes that are less soundly based and more difficult to test. Moreover, the tensors 𝐾𝑐𝑖𝑗 and 𝑀𝑖𝑗𝑘𝑚 are of higher rank than the fluxes they determine, so one risks exacerbating the complexity of the calculation unnecessarily.

4.1. Axially Symmetric Convection

In order to establish the procedure, I consider first the simple case of a star that is spherically symmetrical. Then the turbulence is (statistically) axially symmetrical about the vertical, and 𝑅𝑖𝑗=𝑝𝑡(𝑖)𝛿𝑖𝑗 has only two independent components.

The approach I adopt here assumes that an eddy comes into being spontaneously (as a result of the convective instability acting on some random fluctuation in the medium) and subsequently grows in accordance with equations of motion linearized about the mean, background, state, and at all times there is some likelihood that it breaks up due to internal shear instability. The disruption is regarded as being instantaneous and occurs stochastically with probability proportional to the rms rate of strain in the eddy. It is assumed that the smaller-scale motion so generated contributes negligibly to the fluxes. I shall describe here only a local theory, using the Boussinesq approximation with 𝛽𝑖 constant within any eddy. The discussion is based on the review by Gough [18], where further elaborations and refinements, such as an approximate account of the smaller-scale motion, can be found.

To render the notation more transparent I rename the 𝑥𝑖 coordinates (𝑥,𝑦,𝑧), use vector notation where convenient, and denote the components 𝑢𝑖 of the fluctuating velocity field𝐮 by (𝑢,𝑣,𝑤). The convective flow is represented by a superposition of eigenfunctions of the linear problem, each of which is confined between two horizontal planes a distance ±(1/2) from its central level 𝑧0. I adopt stress-free isothermal boundary conditions, which yield the simplest formulae. Thus for 𝑧0(1/2)𝑧𝑧0+(1/2), a convectively unstable eddy created at time 𝑡0 with an initial vertical velocity amplitude 𝑊0 may be represented bŷ̂𝜗𝑊(𝑢,𝑣,𝑤,𝜗)=̂𝑢,𝑣,𝑤,0𝑒𝑞(𝑡𝑡0),(24) in whicĥ𝑢𝑥,𝑦,𝑧𝑧0=𝑘2𝑘𝑣𝑘sin𝑣𝑧𝑧0𝑓𝑥,̂𝑣𝑥,𝑦,𝑧𝑧0=𝑘2𝑘𝑣𝑘sin𝑣𝑧𝑧0𝑓𝑦,𝑤𝑥,𝑦,𝑧𝑧0𝑘=cos𝑣𝑧𝑧0̂𝜗𝑓,𝑥,𝑦,𝑧𝑧0=𝛽𝑞1𝑤𝑥,𝑦,𝑧𝑧0,(25) where the planform 𝑓(𝑥,𝑦) satisfies2𝑓+𝑘2𝑓=0,𝑓2=1,(26) and where 𝐤=(𝑘𝑥,𝑘𝑦,0) is the horizontal component, with magnitude 𝑘, of a wavenumber 𝐤 characterizing the eddy, having magnitude 𝑘; 𝑘𝑣=𝜋/ is its vertical component. Also 𝑓𝑥=𝜕𝑓/𝜕𝑥, and so forth, 2 is the horizontal Laplacian operator; the overbar denotes horizontal, or ensemble, average. Furthermore, the growth rate 𝑞 of the mode is given by1𝑞=2𝑘2𝜅+𝜇2+14𝑘4𝜅2=12𝑘4𝑆+112𝜅,(27) in which𝜇2=𝑔𝛿𝛽𝑘Φ𝑇,Φ=1+2𝑣𝑘2,(28) and now𝜇𝑆=2𝜅2𝑘4(29) (e.g., [19]), which differs from the definition (15) by just a constant geometrical factor. Where |𝑧𝑧0|>(1/2),𝑢=0and 𝜗=0 (although, of course, there are other eddies there). The formulae can easily be generalized to account for optically thin eddies [6]. Note that formally stable eddies also exist, with ̂𝜗=𝛽𝑞1𝑤, but their direct contribution to the fluxes is not significant; they formally contribute indirectly via a filling factor 𝑓 to be introduced below. The geometrical factor Φ and the planform 𝑓 define the degree of anisotropy of the turbulence.

An eddy is presumed to come into existence as a residue of the turbulent cascade from the breakup of previously existing eddies. The instant at which it is considered to have an existence of its own is not completely specified, although to be unstable it must have a temperature fluctuation that has a (randomly achieved) positive correlation with the vertical velocity. The actual criterion that is adopted to define the initial conditions is absorbed into the definition of the scaling factor Λ appearing in (35) and (37) for the turbulent fluxes, which in practice is subsequently calibrated by comparison with astronomical observation.

The turbulent fluxes can now be computed at height 𝑧 by first constructing the appropriate fluxes due to a single eddy and then averaging over all eddies that intersect the horizontal plane through 𝑧, weighting the average by the probability that the eddy has not been disrupted. The difficulties in performing this average arise from the computation of that probability, from assigning the distribution of initial amplitudes and the rate at which eddies are created, and for assigning a probability distribution to the eddy geometry characterized by the parameter Φ and, in the case of rotating flow discussed in the next subsection, the function 𝑓.

Let us consider first the computation of the disruption probability. Eddy disruption is thought to occur predominantly as a result of shear instability within the eddy; its probability of occurring is therefore proportional to the magnitude of that shear. The coefficient of proportionality depends on the geometry of the eddy and is difficult to estimate because shear instability is ill understood. Perhaps the most plausible first guess is to take as the measure of the internal shear the square root 𝑒=𝑒0exp[𝑞(𝑡𝑡0)] of the total squared rate of strain 𝑒2=𝑒𝑖𝑗𝑒𝑖𝑗, where 𝑒𝑖𝑗=(1/2)(𝜕𝑢𝑖/𝜕𝑥𝑗+𝜕𝑢𝑗/𝜕𝑥𝑖) and the angular brackets denote average over the volume occupied by the eddy. The ansatz is plausible, for at least it has the correct form in the limits of large and small 𝑘𝑣/𝑘 when the eddy motion approaches rectilinear shear. From the structure (25) of the velocity field, the disruption probability per unit time can then be evaluated to be1𝛾𝑒=2𝛾1Φ1𝑘𝑣𝑊0𝑞exp𝑡𝑡0,(30) where 𝛾 is some (constant) parameter to be calibrated later. The probability that an eddy created at time 𝑡0 still exists is therefore𝒫𝑡;𝑡0e=exp𝜆𝑞(𝑡𝑡0)1,(31) where 𝜆=𝛾𝑞1𝑒0. It is an important assumption of the theory that the eddy dynamics is dominated by its growth against the unstably stratified environment, which can be well approximated by linear theory. Then the duration of the process of eddy disruption can be regarded as being short compared with the mean eddy lifetime 𝜏, which permits disruption to be considered to be instantaneous. This implies that the mean eddy lifetime is given by𝜏=𝛾𝑒0𝑡0𝑡𝑡01𝒫d𝑡ln𝜆𝑞1,(32) the final, approximate, expression having been obtained as the leading term in an expansion for small 𝜆.

We come now to the distribution of 𝑊0 and the eddy creation rate. I shall adopt the common practice of assuming that the turbulent convective flow is dominated at any point by eddies with a unique size and shape. Since the flow in a spherically symmetrical (nonrotating) star is statistically axisymmetrical, all orientations of the planform 𝑓 are equally probable. Furthermore, since the theory is local and, for the present, the convection zone is considered to be time independent in the mean, the formula for the creation rate and the initial amplitude 𝑊0 depend on neither 𝑧0 nor 𝑡0 at any given value of 𝑧.

If 𝑛(𝑧0) is the rate of creation of eddies centred at 𝑧=𝑧0 per unit eddy-volume, then the proportion of the volume occupied by eddies at time 𝑡 is𝑓=𝑛𝑡𝒫d𝑡0=𝑛𝜏,(33) whence 𝑛=𝑓𝜏1; 𝑓 is essentially the filling factor of eddies with positive linear growth rates.

One is now in a position to evaluate the convective heat flux by integrating the contributions from all eddies. Formally𝐹𝑐=𝜌𝑐𝑝𝑛𝑡𝑤𝜗𝒫d𝑡0d𝑧0.(34) Only eddies that exist at height 𝑧 contribute to the integral over eddy locations 𝑧0, so the outer integral in this simple theory extends only from 𝑧/2 to 𝑧+/2. Moreover, consistent with the Boussinesq approximation, when evaluating that integral the spatial variation of 𝑛 can be ignored.

Note that if the horizontal and vertical wavenumbers are of the same order of magnitude, this expression for the probability is similar to that derived from the more literal interpretation of the mixing-length annihilation hypothesis suggested by Spiegel [20] in terms of vertically rising and falling fluid elements, which yields a breakup probability per unit time of 𝑤/, which is proportional to 𝑘v𝑊0. It is evident from (30) that the two expressions are rather different, especially for the highly elongated eddies having 𝑘𝑣𝑘 that are sometimes favoured by that literal interpretation [6], particularly if only those eddies with the greatest growth rates are considered, as Spiegel suggested. According to the formula proposed here, such eddies are prone to rapid disintegration by the shear in the highly elongated turbulent flow. Nevertheless, this approach is tantamount to no more than a sophisticated mixing-length formalism. It is evident that the manner in which the mixing-length hypothesis is interpreted affects the predicted anisotropy.

With these assumptions the averages can easily be performed. The analysis is essentially the same as that presented by Gough [6, 18]. The resulting fluxes are𝐹𝑐=ΛΦ5/2(Φ1)𝜌𝑐𝑝𝑇𝑞𝑔𝛿32=ΛΦ5/2(Φ1)1+4𝑆112𝑆13𝜌𝑐𝑝𝑔𝛿𝛽𝑇2𝛽,(35) which is essentially of the same form as (14). In this formula,4𝑓Λ=𝜋2𝛾2ln1/𝜆,(36) which demonstrates how the filling factor 𝑓, the amplitude at which the eddy is considered to come into existence (which determines the value of 𝑒0), and the scaling factor 𝛾 defining the eddy destruction rate combine together into a single (calibrateable) parameter. Note that if the assumption that is proportional to a scale height 𝐻 of the background state is made: =𝛼𝐻, the constant of proportionality 𝛼 combines similarly. Equation (35) also demonstrates that the fluxes depend only weakly on 𝜆, and therefore only weakly on when one considers the eddy to have come into existence. One can evaluate the mean-squared velocities for determining the Reynolds stress similarly: 𝑝𝑡=𝜌𝑤2=ΛΦ3(Φ1)𝜌2𝑞2=ΛΦ3(Φ1)1+4𝑆112𝑆12𝜌𝑔𝛿𝛽𝑇2,(37)𝑢2=𝑣2=12(Φ1)𝑤2.(38) Equations (35)–(37) with (29) are the same as (14)–(16), except for the constants multiplying the functional forms of 𝐹𝑐, 𝑝𝑡, and 𝑆, and they show the relation between 𝛼𝐻 and 𝛼𝑅, and Λ and Φ. Finally one can evaluate the flux of kinetic energy, 𝐾𝑡, likewise:𝐾𝑡=Λ𝐶Φ9/2(Φ1)3/2(Φ+3)1+4𝑆112𝑆13×𝜌𝑔𝛿𝛽𝑇3/23,(39) where 𝐶=82𝑓3/3𝜋2𝛾 is a velocity correlation constant. It requires yet another independent assumption to specify 𝐶. If, for example, it is assumed that 𝛾=1 and that every horizontal plane is tessellated with Christopherson hexagons (e.g., [21]), which can be represented by six wavenumbers 𝐤 of equal magnitude 𝑘 orientated uniformly in a circle, then 𝐶=4/33𝜋20.078 (cf. [22, 23]); if instead the plane is filled with closely packed cylinders with stagnant fluid in the interstices, as Böhm-Vitense [2] might prefer, then 𝐶0.034.

The analysis has not yet provided a method for computing the anisotropy parameter Φ. The simplified formalism presented here is incomplete, and it exposes important issues that must be resolved in order to close what remains a rather naive eddy-ensemble approach. However, before attempting a generalization, an observation of the structure of the formulae (35)–(37) is perhaps not out of place.

Although it has been assumed that turbulent eddies with a given value of Φ dominate the flow, it is appreciated that actually this Φ is merely meant to be representative of a spectrum of eddy shapes. What is meant by an eddy that dominates the flow is one that contributes maximally to a particular flux. Thus 𝑤2 is maximized at fixed when Φ=1.5, and 𝑢2 is maximized at Φ=3. If these values are used separately for the two fluxes, which is not an implausible procedure if the distribution of eddy shapes is rather flat, one finds𝑢2=12𝑤2.(40) This relation is also a property of a single eddy with Φ=2 (which is close to the geometric mean of the two stress-maximizing values, namely 3/2=2.1). So perhaps this intermediate value of Φ should be used to describe the representative eddy for computing the stress tensor. It is interesting, though perhaps merely coincidental, that the eddy shape that maximizes 𝐹𝑐 at fixed has Φ=5/3, which does not differ substantially from 2.

4.2. Convection in a Rotating Fluid

Generalization of the procedure outlined in the previous subsection in the presence of a general background mean flow is not straightforward, because it requires additional hypotheses. Nevertheless, some progress can be made.

The large-scale currents in stars envisaged here include both meridional circulation and zonal flow. I continue to assume that a meaningful separation of the scales of motion is possible, so that a local theory provides a fair first approximation to the turbulent fluxes. As usual only the leading terms in the Taylor expansion of 𝐯about a point will be considered.

Note first that the mean motion in the vicinity of a point can be decomposed into a translation, a rotation, a pure shearing motion, and an isotropic dilatation. The first does not affect the turbulent motion; the last modifies the turbulent pressure but does not add to it off-diagonal terms; it can be dealt with by generalizing the analysis by Gough [6] for sinusoidally pulsating stars. Rotation modifies the turbulence via the Coriolis forces, and pure shearing motion stretches the eddies; both of these do generate off-diagonal terms. Here I consider the special case where only the local rotation is appreciable; it is the simplest of the steady flows to deal with because there is no shear, and the modifications to the eigenfunctions that the mean flow imposes are relatively straightforward to calculate. Steps towards accounting for rectilinear shearing flow are being taken by Smolec et al. [24, 25].

In order to proceed I shall make two additional simplifying assumptions. Firstly, I assume that surfaces of constant potential temperature (or entropy) are horizontal. This may not be strictly the case, as Durney [9] has pointed out, although the absence of a large pole-equator temperature difference in the sun suggests that at least in that case the assumption may not be a bad first approximation, at least in the upper parts of the convection zone. Secondly, I assume that the vorticity of the mean flow is much less than the eddy growth rate, which (almost) circumvents the issue of how axisymmetry is broken by rotation. This assumption too is easily satisfied in the upper parts of the solar convection zone, but not so well for the slowest eddies that occupy much of the zone beneath. Thus, as with the (sometimes implicit) assumption that the mixing length is much less than the background scale heights, the formulae I derive are not strictly valid, and one should be aware of that when they are used in larger-scale calculations.

As before, the linear eigenfunctions of a plane Boussinesq layer are computed, but now the layer is presumed to rotate with angular velocity Ω=(1/2)curl𝐯. The Cartesian axes are orientated such that Ωlies in the 𝑥𝑧 plane. In view of the assumption Ω=|Ω|𝑞𝜇, there is no question of the rotation stabilizing the convection; the growth rate of the convective mode and its eigenfunction can be expanded in powers of 𝜀=Ω/𝜇, here up to second order:𝑞=𝑞01+𝜀̃𝑞1+𝜀2̃𝑞2,+̂𝑢=̂𝑢0+𝜀̂𝑢1+𝜀2̂𝑢2+,etc.(41) The leading-order solution is given by𝑤0𝑘=cos𝑣𝑧𝑧0𝑓=𝑞0𝛽1̂𝜗0,̂𝑢0=𝑘2𝑘𝑣𝑘sin𝑣𝑧𝑧0𝑓𝑥=𝑓𝑥𝑓𝑦1̂𝑣0,(42) which are equivalent to the formula (25) with 𝑞0 given by (27)–(29) for 𝑞. Thermal diffusion is important only in the immediately subphotospheric layers of the star. There the growth rate is high, 𝜀 is very small, and therefore so typically is the influence of the rotation. Deeper down where rotation may matter, the motion is close to being adiabatic. Therefore in calculating the corrections to the structure of the eddies it is expedient to assume adiabaticity, approximating 𝑞0 by 𝜇 and omitting some terms in the formulae for the eigenfunctions, in order to maintain a degree of lucidity. It would be quite straightforward to retain the nonadiabatic corrections were the necessity to arise, but I refrain from doing so here because the added complexity of the formulae would render them less easy to appreciate. The rotational modification to the convective mode can accordingly be represented bỹ𝑞1𝑤=0,1̂𝜗=0,1=0,̂𝑢1=2𝑘2𝑘cos𝜃sin𝑣𝑧𝑧0𝑓𝑦𝑘sin𝜃cos𝑣𝑧𝑧0𝑓𝑥𝑦,̂𝑣1=2𝑘2𝑘cos𝜃sinv𝑧𝑧0𝑓𝑥𝑘sin𝜃cos𝑣𝑧𝑧0𝑓𝑥𝑥,̃𝑞2=2𝑘2𝑓2𝑥sin𝜃+𝑘2𝑣cos2𝜃𝑤2=2sin2𝜃𝑧𝑧0𝑘cos𝑣𝑧𝑧0𝑓𝑥,̂𝜗2=𝜇1𝛽𝑘2sin2𝜃𝑓𝑥𝑥+𝑘2vcos2𝑤𝜃𝑓2,̂𝑢2=2𝑘2𝑘sin2𝜃cos𝑣𝑧𝑧0𝑘𝑣𝑧𝑧0𝑘sinv𝑧𝑧0𝑓𝑥𝑥=𝑓𝑥𝑥𝑓1𝑥𝑦̂𝑣2,(43) where 𝜃 is the angle of inclination of Ωto the vertical in the positive 𝑥 direction.

One can now go through the procedure described in Section 4.1 to evaluate the turbulent fluxes. That is quite straightforward except for one issue: the determination of the potential horizontal anisotropy of the turbulent velocity amplitudes resulting from the existence of a preferred direction introduced by the rotation. The assumption that the growth of the eddies satisfies linearized dynamics implies that anisotropy formally arises only from an anisotropy of initial conditions. Since the creation of the unstable convective eddy perturbations is not directly addressed by the formalism, an additional assumption must therefore be made. The assumption that I adopt here is simply that the filling factor 𝑓 is isotropic. When it was introduced in the discussion in Section 4.1 of axisymmetric turbulence I implicitly presumed it to be constant. But now that axisymmetry is lost, that presumption must be relaxed. Here, I adopt the principle of detailed balance, by applying (33) separately in every horizontal direction. That implies that the creation rate 𝑛 is proportional to 𝑞 and therefore shares with 𝑞 the same O(𝜀2) anisotropy. However, it transpires that the anisotropy it imparts to the velocity field is even smaller, reducing the O(𝜀2) terms by an O(ln1/𝜀) factor; it can therefore be neglected. After carrying out the averaging, there results𝐅𝑐=𝜀2sin2𝜃,𝜀sin𝜃,1𝜓𝐹𝑐,(44)𝑅𝑖𝑗=121(Φ1)+2𝜀2𝑅𝑥02𝜀201sin2𝜃2(1Φ1)+2𝜀2𝑅𝑦𝜀sin𝜃2𝜀2sin2𝜃𝜀sin𝜃1𝜓𝑝𝑡,(45) in which𝑅𝑥=Φ1,𝑅𝑦=(Φ1)cos23𝜃+4sin2𝜃,𝜓=12𝜀Φ1sin𝜃𝜀2Φ26(Φ1)cos2𝜃sin2𝜃,(46) and where 𝐹𝑐 and 𝑝𝑡 are given by (35) and (37). The general formula for 𝐊t is algebraically complicated, but in the case when the planform is composed of randomly orientated harmonic hexagons it simplifies to𝐊𝑡=1316𝜋𝜀122(Φ1)sin2𝜃,4(Φ1)𝜀sin𝜃,1𝜓𝐾𝑡,(47) where 𝐾𝑡 is given by (39).

The most obvious influence of the rotation on the convective fluxes is to rotate them principally out of the plane of 𝐠 and Ω, namely, in the 𝑦 direction in the coordinate system adopted here. The angle of rotation is proportional to 𝜀=Ω/𝜇. There are additional O(𝜀2) components generated in the plane of 𝐠 and Ω, and also the overall magnitudes of the fluxes are reduced somewhat, also by O(𝜀2). The rotation of the momentum fluxes generates off-diagonal terms in the Reynolds-stress tensor.

5. Discussion

Many assumptions have been necessary to arrive at a tractable procedure. Most notably, it has been assumed here that the turbulent medium can be represented as an ensemble of eddies that not only do not interact directly with each other—they do so implicitly via the background (mean) stratification, however—but that nonlinear advection processes and nonlinearities in the thermodynamics within an eddy (the latter being a consequence of the assumption that the theory is local, which is intimately linked with the Boussinesq assumption) can be ignored throughout almost the entire lifetime of each eddy. The only explicit nonlinearity is the eddy-disruption process, which is assumed to occur on a timescale short compared with the eddy lifetime, which is essential, although not sufficient, for the validity of the linearization assumption: each eddy is presumed to grow via the linear influence of the background state, although there is an implicit nonlinear backreaction on the mean state via the turbulent fluxes. That reaction has been presumed to preserve the alignment of the superadiabatic temperature gradient with gravity, although one can, with additional assumption about the preservation (or otherwise) of eddy shape, take a departure formally into account. Relating that to the fluxes, however, would require a subsequent global calculation.

The formalism is basically a mixing-length approach, which itself is also characterized by linear reasoning, either explicitly when describing the energy exchange with the background state or implicitly via the preservation of eddy size and shape, which is not predicted by the simplest forms of the theory; in any case, it always requires admitting the necessity of an additional assumption about eddy shape. Even though in linear theory the harmonic relation (26) is sufficient for determining the growth rate 𝑞 and the 𝑧 and 𝑡 dependence of the eigenfunctions, more detailed specification is necessary for relating the kinetic-energy flux to the heat flux and the Reynolds stress, in particular for determining the correlation constant 𝐶. Eddy size is basically the mixing length, which is not an integral part of the theory. Therefore nor is any procedure to define how it changes with changing vorticity (and, in more general circumstances, the changing shear) in the mean flow. Here I have assumed that the “filling factor” of eddies, a fundamental property of the eddy correlation, is unaffected, leaving the way clear to determine the effect on the eddy shape from the linearized dynamics. As always in such a theory, the value of the vertical extent of the eddy is left to other considerations. In particular, no attempt has been made to address how it might change with the properties—in the simple case considered in this paper, the value of the local rotation—of the background flow.

There remains much work to do before an even modestly reliable general theory emerges. Relaxing the restrictions that were imposed in the discussion above not only complicates the analysis but also requires one to face new problems whose resolution may depend on introducing new hypotheses. Accepting that rotation is not small compared with eddy growth rates, for example, raises the issue of the dependence of eddy creation on the orientation of the eddy. Imagining constant entropy surfaces to deviate from the horizontal introduces the possibility of baroclinic instability. How do these compete with the convective modes? It is easy to invent a set of superficially plausible hypotheses for tackling these problems, but some other no less implausible models may yield rather different results. The simple formulae (14)–(16), for example, result from a relatively wide variety of physical pictures of the turbulent motion, but the more subtle results such as (44)–(47) are more sensitive to the details of the assumptions. Can one even hope to develop a theory on the basis of linear eigenmodes? Or do nonlinearities play an essential role during the growth of an eddy? How inaccurate is the process of separating scales? It is certainly quite evident that assuming a unique eddy lengthscale at any location can be a poor recipe, especially near the edges on convection zones where local stability changes sense. One might improve the description of the small-scale turbulence by using a nonlocal theory, but how can one develop a procedure that is both tractable and realistic? Can the mixing-length formalism in any guise be expected to provide an adequate description of the turbulence?

It has been argued that because mixing-length theory is even less reliable at predicting fluxes in complicated situations than it is in the simplest of cases, the complexities of its consequences should be ignored and that it is not worthwhile attempting to generalize the theory. I disagree with this point of view. The mixing-length theory is still the only practical tool available for modelling stellar convection zones for evolutionary calculations—although work such as that by Trampedach and Augustson [26] may change that, but perhaps not before mixing-length theory has been calibrated more precisely [27]—and it should be utilized to fullest advantage until it is superseded. Why prefer a formula such as (23), which is surely incorrect, to one that attempts to include the effects of phenomena that are known to exist? Maybe the results obtained with a more intricate theory will at first be numerically no better at reproducing astronomical observations, but only by studying the consequences of rotation, shear, baroclinicity, and also of magnetic fields on the turbulence can some idea of the importance of these factors be acquired. Ogilvie [28], for example, has made some progress with the latter, using a rather different approach from the one adopted here. The introduction of each new factor tends to require new hypotheses, or to bring into prominence already existing hypotheses upon which the theory did not previously depend in an important way. With these come new parameters which need to be determined. In principle we have at our disposal laboratory experiments and numerical experiments performed with the solar modelling programmes to calibrate the formulae. Both should be utilized—indeed, some have already (e.g., [24, 25, 29]), but only to a limited extent—although eventually care must also be taken to ensure that sufficient checks remain to test the predictive power of the theory.


The author thanks Pierre Lesaffre for encouraging him to reformulate the theory and extend the expansion to second order in 𝜀 and then to publish it, and he thanks Günter Houdek for useful discussions; also Paula Younger for retyping the original paper. He is grateful to the Leverhulme Trust for an Emeritus Fellowship.


  1. E. A. Spiegel and G. Veronis, “On the Boussinesq approximation for a compressible fluid,” Astrophysical Journal, vol. 131, p. 442, 1960, correction: vol. 135, p. 665. View at Publisher · View at Google Scholar
  2. E. Böhm-Vitense, “Über die wasserstoffkonvektionszone in sternen verschiedener und effektivtemperaturen und leuchtkräfte,” Zeitschrift für Astrophysik, vol. 46, pp. 108–143, 1958. View at Google Scholar
  3. W. Unno, “Stellar radial pulsation coupled with the convection,” Publications of the Astronomical Society of Japan, vol. 19, p. 140, 1967. View at Google Scholar
  4. J. P. Cox, Principles of Stellar Structure, Gordon & Breach, London, UK, 1967.
  5. R. Kippenhahn and A. Weigert, Stellar Structure and Evolution, Springer, Berlin, Germany, 1990.
  6. D. O. Gough, “Mixing-length theory for pulsating stars,” Astrophysical Journal, vol. 214, no. 1, pp. 196–213, 1977. View at Google Scholar
  7. D. O. Gough and N. O. Weiss, “The calibration of stellar convection theories,” Monthly Notices of the Royal Astronomical Society, vol. 176, pp. 589–607, 1976. View at Google Scholar
  8. J. Christensen-Dalsgaard, W. Däppen, S. V. Ajukov et al., “The current state of solar modeling,” Science, vol. 272, no. 5266, pp. 1286–1292, 1996. View at Google Scholar
  9. B. R. Durney, “On the constancy along cylinders of the angular velocity in the solar convection zone,” Astrophysical Journal, vol. 204, no. 1, pp. 589–596, 1976. View at Publisher · View at Google Scholar
  10. R. Kippenhahn, “Differential rotation in stars with convective envelopes,” Astrophysical Journal, vol. 137, p. 664, 1963. View at Publisher · View at Google Scholar
  11. E. A. Spiegel and J.-P. Zahn, “The solar tachocline,” Astronomy and Astrophysics, vol. 265, no. 1, pp. 106–114, 1992. View at Google Scholar
  12. J. Wasiutyński, “Studies in hydrodynamics and structure of stars and planets,” Astrophisica Norvegica, vol. 4, 1946. View at Google Scholar
  13. W. J. Cocke, “On the solar differential rotation and meridional currents,” Astrophysical Journal, vol. 150, p. 1041, 1967. View at Publisher · View at Google Scholar
  14. D. O. Gough, “The pulsational stability of a convective atmosphere,” in Geophysical Fluid Dynamics II, pp. 49–84, Woods Hole Oceanographic Institution, 1965. View at Google Scholar
  15. D. O. Gough, “The estimation of fluxes due to small-scale turbulence,” in Proceedings of the Workshop on Solar Rotation, G. Belvedere and L. Paternò, Eds., pp. 337–350, Università di Catania, Catania, Italy, 1978.
  16. B. R. Durney and H. C. Spruit, “On the dynamics of stellar convection zones —the effect of rotation on the turbulent viscosity and conductivity,” Astrophysical Journal, vol. 234, pp. 1067–1078, 1979. View at Publisher · View at Google Scholar
  17. G. Rüdiger, P. Egorov, L. L. Kitchatinov, and M. Küker, “The eddy heat-flux in rotating turbulent convection,” Astronomy and Astrophysics, vol. 431, no. 1, pp. 345–352, 2005. View at Publisher · View at Google Scholar
  18. D. O. Gough, “The current state of stellar mixing-length theory,” in Proceedings of the International Astronomical Union Colloquium, E. A. Spiegel and J.-P. Zahn, Eds., pp. 15–56, Springer, Heidelberg, Germany, 1977.
  19. P. Ledoux, M. Schwarzschild, and E. A. Spiegel, “On the spectrum of trurbulent convection,” Astrophysical Journal, vol. 133, p. 184, 1961. View at Publisher · View at Google Scholar
  20. E. A. Spiegel, “A generalization of the mixing-length theory of turbulent convection,” Astrophysical Journal, vol. 138, p. 216, 1963. View at Publisher · View at Google Scholar
  21. S. Chandrasekhar, Hydrodynamic and Hydromagnetic Stability, Oxford University Press, Oxford, UK, 1961.
  22. P. H. Roberts, “On non-linear Bénard convection (with an Appendix by K. Stewartson),” in Non-Equilibrium Thermodynamics, Variational Techniques, and Stability, R. Donnelly, R. Herman, and I. Prigogine, Eds., p. 125, University of Chicago Press, 1966.
  23. D. O. Gough, E. A. Spiegel, and J. Toomre, “Modal equations for cellular convection,” Journal of Fluid Mechanics, vol. 68, no. 4, pp. 695–719, 1975. View at Google Scholar
  24. R. Smolec, G. Houdek, and D. O. Gough, “Modelling turbulent fluxes due to thermal convection in rectilinear shearing flow,” in Proceedings of the International Astronomical Union Symposium, N. H. Brummell, A. S. Brun, M. S. Miesch, and Y. Ponty, Eds., vol. 271, pp. 397–398, Cambridge University Press, 2011. View at Publisher · View at Google Scholar
  25. R. Smolec, G. Houdek, and D. O. Gough, “Nonlocal model for the turbulent fluxes due to thermal convection in rectilinear shearing flow,” in Proceedings of the 20th Stellar Pulsation Conference Series: Impact of New Instrumentation & New Insights in Stellar Pulsations, Granada, Spain, 2011.
  26. R. Trampedach and K. Augustson, “Using simulations of solar surface convection as boundary conditions on global simulations,” in Proceedings of the International Astronomical Union Symposium, N. H. Brummell, A. S. Brun, M. S. Miesch, and Y. Ponty, Eds., vol. 271, pp. 403–404, Cambridge University Press, 2011. View at Publisher · View at Google Scholar
  27. R. Trampedach and R. F. Stein, “The mass mixing length in convective stellar envelopes,” Astrophysical Journal, vol. 731, no. 2, article 78, 2011. View at Publisher · View at Google Scholar
  28. G. I. Ogilvie, “On the dynamics of magnetorotational turbulent stresses,” Monthly Notices of the Royal Astronomical Society, vol. 340, no. 3, pp. 969–982, 2003. View at Publisher · View at Google Scholar
  29. P. Tooth and D. O. Gough, “Towards an independent calibration of the mixing-length theory,” in Seismology of the Sun and the Sun-Like Stars, E. J. Rolfe, Ed., ESA, Noordwijk, The Netherlands, 1988. View at Google Scholar