Issue 
A&A
Volume 528, April 2011



Article Number  A78  
Number of page(s)  7  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201015372  
Published online  04 March 2011 
Stellar mixing
III. The case of a passive tracer^{⋆}
^{1}
NASA, Goddard Institute for Space Studies,
New York,
NY
10025,
USA
^{2}
Department of Applied Physics and Applied Mathematics, Columbia
University, New York,
NY
10027,
USA
email: vmc13@columbia.edu
Received:
9
July
2010
Accepted:
25
August
2010
In Papers I–II, we derived the expressions for the turbulent diffusivities of momentum, temperature T, and mean molecular weight μ. Since the scalar Tμ fields are active tracers (by influencing the density and thus the velocity field), whereas passive tracers such as ^{7}Li are carried along by the flow without influencing it, it would be unjustified to use the diffusivities of the Tμ fields to represent the diffusivity of passive tracers. In this paper, we present the first derivation of a passive tracer diffusivity. Some key results are: a) In the general 3D case, the passive tracer diffusivity is a tensor given in algebraic form; b) the diffusivity tensor depends on shear, vorticity, T, and μgradients, thus including double diffusion and differential rotation; c) in the 1D version of the model, the passive tracer diffusivity is a scalar denoted by K_{c}; d) in doubly stable regimes, ∇_{μ} > 0, ∇ − ∇_{ad} < 0, K_{c} is nearly the same as those of the Tμ fields; e) in semiconvection regimes, ∇_{μ} > 0, ∇ − ∇_{ad} > 0, K_{c} is larger than that of the μfield; f) in salt fingers regimes ∇_{μ} < 0,∇ − ∇_{ad} < 0, K_{c} is smaller than that of the μfield; and finally, g) in the only case we know of a direct measurement of a passive tracer diffusivity, the oceanographic North Atlantic Tracer Release Experiment (NATRE), the model reproduces the data quite closely.
Key words: diffusion / instabilities / turbulence / stars: interiors / convection / methods: analytical
© ESO, 2011
1. Introduction
Several people (e.g., Charbonnel & Talon 2005, 2007; Charbonnel & Zahn 2007) have discussed the importance of properly quantifying the mixing properties of tracers that behave as passive tracers, such as^{7}Li. Thus far, the lack of a reliable expression for the diffusivity of a passive tracer was remedied by the use of heuristic expressions which lack predictive power because they are adjusted to reproduce known data. In lieu of empirical relations, one might think of employing the diffusivities of either heat and/or mean molecular weight μ that were derived in Papers I and II. Such an assumption is, however, questionable, because a tracer such as^{7}Li is by definition passive since it does not affect the density of the fluid and in consequence, does not affect the velocity field. On the other hand, both T and μ affect the density and for that reason are called active tracers. There is therefore a need to a) construct a diffusivity for passive tracers; b) show how it relates to that of heat and molecular weight; and c) assess its validity before it is used. This is the task of this work.
While the RSM (Reynolds stress model) discussed in I and II provides a methodology to treat the first two items, provided care is taken to account for all the variables in the problem, item c) is more difficult since we are not aware of any direct measurement of a passive tracer diffusivity in astrophysical settings. Since we believe that an assessment is necessary and perhaps even indispensable for the credibility of the model, we show how our results compare with the only case we know of direct measurement of a passive tracer diffusivity. It is not an astrophysical setting but an oceanographic one, the socalled North Atlantic Tracer Release Experiment (NATRE, Ledwell et al. 1993, 1998) in which the tracer SF_{6} was released in the ocean. From the measurements in time of its dispersion, a value of the corresponding diffusivity was derived. We applied our model to NATRE and showed that the model predictions fare well with the measured data.
Some of the results that have emerged from this analysis can be summarized as follows:
a) the model is algebraic; b) in the general 3D case, the passive tracer diffusivity is a tensor that depends on stable, unstable stratification, shear (differential rotation), and double diffusion processes; c) in the 1D version of the model, in doubly stable regimes, ∇_{μ} > 0, ∇ − ∇_{ad} < 0, the passive tracer diffusivity is nearly the same as those of the Tμ fields; d) in semiconvection regimes, ∇_{μ} > 0, ∇ − ∇_{ad} > 0, the passive tracer diffusivity is larger than that of the μfield; and e) in salt fingers regimes,∇_{μ} < 0,∇ − ∇_{ad} < 0, the passive scalar diffusivity is smaller than that of the μfield.
2. Passive tracer flux
Calling c the passive tracer field, we begin with the general conservation equation (d/dt = ∂_{t} + u_{i}∂_{i}), (1a)where u_{i} is the total velocity field and the diffusion flux J_{i} is given by Chapman & Cowling (1970): (1b)with a_{,i} = ∂a / ∂x_{i} and where κ_{T},κ_{p} are dimensionless functions. Often, in the literature, J_{i} includes the density and has the opposite sign and χ_{c} is called D. We also note that different authors have used different notations, e.g., Eq. (2.1) of Aller & Chapman (1960), Eq. (3) of Schatzman (1969), Eq. (1) of Michaud (1970), and Eq. (1) of Vauclair & Vauclair (1982).

a)
mean cfield. In carrying out the mass averaging process (Favre 1969; Canuto 1997, 1999) defined (2a)we obtain (D_{t} = ∂_{t} + U_{i}∂_{i}) (2b)where U_{i},C are the mean velocity and mean concentration fields, and the first term on the rhs of (2b) is the turbulent c flux defined as (primes indicate turbulent variables): (2c)At first sight, it may look as though Eq. (2b) is unrelated to the Tμ fields or for that matter, to shear and vorticity but, as we show below, the turbulent cflux (2c) entails all those dependences.

b)
passive tracer fluxes. Multiplying (1a) by u_{i}, Eq. (6a) of Paper I by c and adding the two, we obtain (3a)where F_{i} was defined in Eq. (6a) of Paper I as (3b)where the last term on the rhs of (3b) is the kinematic viscosity tensor and the other symbols have their usual meaning. Next, we mass average (3a) with the following results (3c)where are the Reynolds stresses. Subtracting its massaveraged from the original equation, we obtain, after several steps, the equation governing the concentration flux: (4a)where 2S_{ij} = U_{i,j} + U_{j,i},2V_{ij} = U_{i,j} − U_{j,i}, are the shear and vorticity characterizing the mean velocity field. After some further algebraic steps regarding the last terms on the rhs of (4a), we obtain (4b)It is interesting to compare (4b) with the corresponding equations for the Tμ fluxes, Eqs. (9e,11b) of Paper I. The first two terms in each of these three equations are similar since they represent the sources of the fluxes given by the interaction of the Reynolds stresses R_{ij} with the gradients of the corresponding mean variables. The second terms in each of the three equations are also similar since they represent the interaction of the mean shear and vorticity with the fluxes under consideration. The third terms in each expression is where the differences occur. The Tμ fields have similar terms but with opposite signs since the Tμ fields contribute in opposite ways to the density, as Eq. (1b) of Paper I shows. Since the Tμ fields are active tracers, the square terms represent potential energies of the two fields while such a term does not exist in the case of a passive scalar, where we instead have a term that comes from the term g_{i}ρ in the momentum equation, see the second term in (3b); specifically, when we multiply g_{i}ρ by c′, we obtain the term (4c)where we have used Eq. (1b) of Paper I. Finally, since a passive tracer is characterized by a kinematic diffusivity χ_{c}, the latter gives rise to the last term in (4b). Thus, Eq. (4b) tells us that the cfield interacts not only with the largescale fields C,U_{i}, but also with the fluctuating parts of the Tμ fields, as Eq. (4c) clearly shows. To proceed, we need to determine the two fluxes entering the third term in (4b). In the local limit, the first one was given in Eq. (124) of C99, and reads as: (4d)where . Analogously, we derive (4e)Substituting these two relations into (4b) and taking the stationary limit, we obtain that the passive tracer flux is given by the following equation: (5a)where (5b)As for the dissipationrelaxation time scales (5c)we introduce the following notation and approximation: (5d)where τ = 2K / ε is the eddy turnover time scale and Kε are the eddy kinetic energy and its rate of dissipation. The second relation in (5d) rests on the assumption that a passive tracer has the same dissipation time scale as the active tracers T and μ. Equations (5b) then become (5e)The derivation of the equation for the cflux is now complete except for the two dimensionless variables α_{1,2} which we discuss in the next section. We note the following points. First, the passive tracer flux depends on shear, vorticity, T and μgradients. Second, using the HamiltonCayley theorem, one can rewrite (5a) in the following form: (5f)where the tensor diffusivity tensor (K_{c})_{ij} can be written down explicitly (see C99). Third, as for the functions α_{1,2}, relations (13i) of Paper I show that the forms of π_{1,4} defined in Eq. (10c) of Paper I conserve the Tμ symmetry. It is therefore tempting to identify α_{1} with either π_{1,4} because they are both potentially acceptable candidates. However, since we have no reason to chose one over the other, the construction of the function, (5g)is a non trivial problem that we discuss in the next section.
Fig. 1 Doubly Stable case, ∇_{μ} > 0, ∇ − ∇_{ad} < 0, R_{μ} < 0. The dimensionless mixing efficiencies, defined in Eq. (6d), for the case of a passive tracer, spice, and the ratios of the former with those of μ, heat, buoyancy, and spice. 
3. The functions α_{1,2}
To determine the functions α_{1,2}, we suggest the following procedure. Let us consider the 1D version of Eqs. (5). Using the notation , we derive that (6a)Using the relation in (Paper II, footnote 1), it is easy to show that (6b)The computation of e_{33} is a bit more involved but using relations (Paper II, footnote 1) and (Paper I, Appendix C), we obtain the following result (6c)where R_{μ} = ∇_{μ}(∇ − ∇_{ad})^{1} is the socalled density ratio and the functions A_{h,μ} are given in (Paper I, Appendix C). Since the form of K_{c} is also of the type (Paper II, Eqs. (6b, c)), namely (α = h,μ,c), (6d)the corresponding A_{c} is given by (6e)The dimensionless functions Γ in (6d) are known as mixing efficiencies, and they were discussed in Paper II, Sect. 3. To double check (6e), in Appendix A we rederive it using a 1D model from the very beginning. If we compare the forms of the functions A_{h,μ} with that of A_{c}, we notice that they are different, which is why we stated earlier that the identification of K_{c} with either K_{h} and/or K_{μ} is hard to justify. The point will be clearer when we plot such diffusivities and their ratios in Figs. 1–3.
Fig. 2 Same as in Fig. 1 but for the case of semiconvection ∇_{μ} > 0, ∇ − ∇_{ad} > 0, R_{μ} > 1. The values of R_{μ} are the reciprocal of those of the case in Fig. 3. 
Next, we consider the active tracer represented by the μfield and take the limit in which the α_{μ}contraction coefficient becomes vanishingly small (7a)Physically, (7a) implies that a change in μ does not affect the density, in which case the μfield behaves like a passive tracer. In this limit, we should therefore have: (7b)Since K_{c} depends on the functions α_{1,2} and K_{μ} depends on the π′s, the limit (7b) will allow us to establish a relation between the α_{1,2} and the π′s for which we have an explicit form (Paper I, Eq. (13i)). We begin with Eq. (6e), which now becomes (8a)where A_{h} is given in Appendix C of Paper I which, in this limit, reads as follows: (8b)We next consider the form of A_{μ} given in Appendix C of Paper I, (8c)It is easy to verify that the third relation (7b) is satisfied if we take (8d)By repeating the same procedure with the Tfield instead, we obtain (8e)which is identical to (8d), thus assuring the validity of the procedure. Next, since the shortest time scale dominates, we begin by writing (9)Introducing the π′s from their definition (10c) of Paper I and α_{1} from (5d), we choose A,B in (9) to satisfy the two requirements (8d,e). The result is (10)Admittedly, relation (10) is not the only possible functional dependence that satisfies the limits (8d,e), but we have chosen it for its simplicity. As for α_{2}, its identification with π_{2} seems unproblematic: (11)
Fig. 4 The zprofile of the predicted passive tracer diffusivity vs. the values measured at NATRE. 
4. Results
Using the previous results and the variables (τN)^{2} from solving Eqs. (9ce) of Paper II, Figs. 1–3 show results for different values of Ri and R_{μ}. In Fig. 1, corresponding to the DoublyStable case ∇_{μ} > 0, ∇ − ∇_{ad} < 0 and thus R_{μ} < 0, we present the dimensionless mixing efficiencies Γ; see Eq. (6d), for a passive tracer panel a) and in panel b) we present the case of spice (a term borrowed from oceanography) defined as σ ≡ α_{T}T + α_{μ}μ. Since σ is invariant under the exchange of the T and μ fields, it could be viewed as a possible candidate for a passive tracer. The corresponding mixing efficiency is (12)In the remaining panels c)–f) we exhibit the ratios of the mixing efficiencies of a passive tracer with those of μ, heat, buoyancy and spice. As one can observe from panels c)–f), for all practical purposes, there is little difference between the diffusivity of a passive tracer and that of μ, heat or spice. However, since a doublystable situation is uncommon, the case is rather academic. In Fig. 2 we present the results for the case of semiconvection, ∇_{μ} > 0, ∇ − ∇_{ad} > 0, and thus R_{μ} > 1. As one can observe from panel c), the passive tracer diffusivity is larger than that of μ by up to a factor of two and it is smaller than that of heat by the same amount. Finally, in Fig. 3 we present the results for the case of saltfingers, ∇_{μ} < 0,∇ − ∇_{ad} < 0 and thus R_{μ} < 1, in which case we notice that the passive tracer mixing efficiency is smaller than that of μ but larger than that of heat by roughly the same amount.
5. Assessment of the model prediction against the NATRE data
As an assessment of the model, we present a comparison of the model results with those of the only case we know in which a passive tracer diffusivity was measured directly. It is not an astrophysical case but an oceanographic one corresponding to the experiment NATRE (North Atlantic Tracer Release Experiment) in which the tracer SF_{6} was released in the middle of the Atlantic Ocean at the depth of approximately 400 m. It was monitored for several months, and a diffusivity was deduced from the dispersion in time. For details, the reader should consult the original papers by Ledwell et al. (1993, 1998). In the oceanographic case, the role of the μ field is played by salinity, and R_{μ} is denoted by R_{ρ} defined analogously as the ratio of the salinity to the temperature zgradients. In Fig. 4 we present the zprofile of the measured diffusivity vs. the one predicted by the present model in its 1D version presented in Appendix A. The model’s results were obtained using the zprofile of the density ratio R_{ρ} from the measurements themselves (Laurent & Schmitt 1999). The functions A_{h,}μ that enter the second of (6c) were constructed using relations (Paper II, Eq. (6d)) in which the heattosalt flux ratio “r” was computed from its definition, Paper II, Eq. (7b). The variable x in A_{h,}μ was obtained by solving Eqs. (Paper II, Eqs. (9c–e)), and the variable is given by Eq. (Paper II, Eq. (6h)). The model reproduces the measured data fairly well.
6. Summary and conclusions
As one may have guessed on physical grounds, the diffusivity of a passive tracer is not the same as that of active tracers such as temperature, buoyancy, and mean molecular weight. The derivation presented in this paper makes that clear, and the results shown in Figs. 1–3 quantify that relationship, specifically, the ratios of the diffusivities of passive vs. active tracers. The model results are algebraic, hence not difficult to implement in stellar structureevolution codes.
References
 Aller, L. H., & Chapman, S. 1960, ApJ, 132, 1 [Google Scholar]
 Canuto, V. M. 1997, ApJ, 482, 827 [Google Scholar]
 Canuto, V. M. 1999, ApJ, 524, 311 (cited as C99) [NASA ADS] [CrossRef] [Google Scholar]
 Canuto, V. M. 2011, A&A, 528, A76, A77 (Papers I and II) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chapman, S., & Cowling, T. G. 1970, The mathematical theory of non uniform gases (Cambridge Univ. Press) [Google Scholar]
 Charbonnel, C., & Talon, S. 2005, Science, 309, 2189 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Charbonnel, C., & Talon, S. 2007, Science, 318, 922 [CrossRef] [PubMed] [Google Scholar]
 Charbonnel, C., & Zahn, J. P. 2007, A&A, 467, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Favre, A. 1969, in Problems of Hydrodynamics and continuous mechanics (Philadelphia: SIAM), 231 [Google Scholar]
 Ledwell, J. R., Watson, A. J., & Law, C. S. 1993, Nature, 364, 701 [NASA ADS] [CrossRef] [Google Scholar]
 Ledwell, J. R., Wilson, A. J., & Law, C. S. 1998, J. Geophys. Res., 103, 21499 [NASA ADS] [CrossRef] [Google Scholar]
 Michaud, G. 1970, ApJ, 160, 641 [NASA ADS] [CrossRef] [Google Scholar]
 Schatzman, E. 1969, A&A, 3, 331 [NASA ADS] [Google Scholar]
 St. Laurent, L., & Schmitt, R. W. 1999, J. Phys. Oceanogr., 29, 1404 [NASA ADS] [CrossRef] [Google Scholar]
 Vauclair, S., & Vauclair, G. 1982, ARA&A, 30, 3 [Google Scholar]
Appendix A: 1D Passive tracer vertical flux
In order to derive the vertical flux of a passive tracer field denoted by c, we begin with the general equation (A.1)where the rhs represents molecular effects with a kinematic diffusivity χ_{c}. We split the c,u_{i} fields into an average and a fluctuating component, c = C + c′ and and proceed as follows. After introducing the above relations in (A.1), using the incompressibility condition, averaging the result, and subtracting it from the initial equation for the total cfield, algebraic steps lead to the following equation for the fluctuating component (D_{t} = ∂_{t} + U_{i}∂_{i}): (A.2)Using an analogous procedure applied to the velocity field described by the NavierStokes equations, we derive the equation for the field w′ which reads as (A.3)where ν is the kinematic viscosity. As expected, the averages of either (A.2–A.3) yield 0 = 0 since all fluctuating variables have a zero average. We also recall that , where T′, μ′are the fluctuating components of the T and μ fields. The next step consists of multiplying Eq. (A.2) by w′ and Eq. (A.3) by c′, averaging the results, and summing the two equations. The result is the dynamic equation for the vertical flux of the passive tracer : (A.4)where we have lumped the pressure correlation and the dissipation terms together into the dissipationrelaxation time scale τ_{pc}. Next, we must construct the two correlations . This is done by using the two equations for the fluctuating Tμ fields, which are obtained using a procedure analogous to the one used to arrive at (A.2). The results are (A.5)where the χ′s represent molecular diffusivities. If we multiply the first of (A.5) by c′ and the second by T′, summing the two relations after averaging them, we obtain the two equations \newpage \noindentwhere, as in Eq. (A.4), the two time scales τ_{cθ} and τ_{cs} represent the relaxationdissipation times scales of the corresponding correlations. For simplicity, we have neglected the thirdorder moments. In Eqs. (A.6, A.7), J_{h,}μ are the vertical fluxes of heat and μ whose form was derived in Eq. (C.1) of Paper I to be (A.8)where K_{h,}μ are the heat and μ diffusivities whose general form is given by relations (6b,c) of Paper II: (A.9)The explicit form of , which is not needed in this context, is given by Appendix C of Paper I. Taking the stationary limits in Eqs. (A.4, A.6, A.7) and using (A.8), we obtain the desired form of the vertical flux of a passive tracer (A.10)where the dimensionless function A_{c} is given by (A.11)which coincides with relation (6e) of the text, and A_{ρ} is defined in the second of relation (6c) of the text.
All Figures
Fig. 1 Doubly Stable case, ∇_{μ} > 0, ∇ − ∇_{ad} < 0, R_{μ} < 0. The dimensionless mixing efficiencies, defined in Eq. (6d), for the case of a passive tracer, spice, and the ratios of the former with those of μ, heat, buoyancy, and spice. 

In the text 
Fig. 2 Same as in Fig. 1 but for the case of semiconvection ∇_{μ} > 0, ∇ − ∇_{ad} > 0, R_{μ} > 1. The values of R_{μ} are the reciprocal of those of the case in Fig. 3. 

In the text 
Fig. 3 Same as in Fig. 1 but for the case of saltfingers, ∇_{μ} < 0,∇ − ∇_{ad} < 0, R_{μ} < 1. 

In the text 
Fig. 4 The zprofile of the predicted passive tracer diffusivity vs. the values measured at NATRE. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.