Next Article in Journal
Analysis of the Observation Results for Preferentially Oriented Particles in High-Level Clouds According to the EARLINET Lidar Network and MODIS Data
Next Article in Special Issue
Statistics of the Inertial Energy Transfer Range in d-Dimensional Turbulence (2 ≤ d ≤ 3) in a Lagrangian Renormalized Approximation
Previous Article in Journal
Bioaerosol Identification by Wide Particle Size Range Single Particle Mass Spectrometry
Previous Article in Special Issue
Non-Equilibrium Turbulent Transport in Convective Plumes Obtained from Closure Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

The Decay of Energy and Scalar Variance in Axisymmetric Turbulence

by
Peter A. Davidson
Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, UK
Atmosphere 2023, 14(6), 1019; https://doi.org/10.3390/atmos14061019
Submission received: 16 March 2023 / Revised: 31 May 2023 / Accepted: 5 June 2023 / Published: 13 June 2023

Abstract

:
We review the recent progress in our understanding of the large scales in homogeneous (but anisotropic) turbulence. We focus on turbulence which emerges from Saffman-like initial conditions, in which the vortices possess a finite linear impulse. Such turbulence supports long-range velocity correlations of the form u i u j = O ( r 3 ) , where u and u are separated by a distance r, and these long-range interactions dominate the dynamics of large eddies. We show that, for axisymmetric turbulence, the energy and integral scales evolve as u 2 ~ u / / 2 ~ t 6 / 5 and l ~ l / / ~ t 2 / 5 , where and // indicate directions that are perpendicular and parallel to the symmetry axis, respectively. These predictions are consistent with the evidence of direct numerical simulations. Similar results are obtained for the passive scalar variance, where we find that θ 2 ~ t 6 / 5 . The primary point of novelty in our discussion of passive scalar decay is that it is based in real (rather than spectral) space, making use of an integral invariant which is a generalization of the isotropic Corrsin integral.

1. Introduction

1.1. The Summer of 2002

In the summer of 2002, the author visited NCAR to discuss the behavior of the large scales in freely decaying, homogeneous turbulence with Jack Herring. It was to be the first of many memorable visits to NCAR. The primary focus at that time was on the behavior of Loitsyansky’s integral in isotropic turbulence, but it left the author with a longstanding interest in the dynamics of large scales in all types of homogeneous turbulence. Initially, that interest centered on large scales in so-called Batchelor turbulence, whose energy spectrum satisfies E ( k 0 ) ~ k 4 , and where Loitsyansky’s integral was shown to be conserved in fully developed turbulence (see Jackson et al. [1] and Ishida, Davidson, and Kaneda [2]). However, interest soon spread to include Saffman ( E ( k 0 ) ~ k 2 ) turbulence (Davidson, Okamoto, and Kaneda [3]; Krogstad and Davidson [4]), as well as the evolution of the large scales in rotating, stratified, and MHD turbulence (Davidson [5,6]). Finally, as part of an international network, the author began a study on the large scales in passive scalar mixing (Yoshimatsu and Kaneda [7]). Many of these studies have their roots in those discussions at NCAR in the summer of 2002.
This paper reviews some of the more recent developments that have been made in the large-scale behavior of freely decaying Saffman turbulence, with particular emphasis on energy decay and the kinematics of passive scalars. It particularly builds on the work of Davidson [5,6] and Davidson, Okamoto, and Kaneda [3], and on the fruitful discussions the author had with Yukio Kaneda and Katsunori Yoshimatsu in 2016 and 2017 as part of an international network.

1.2. Some Important Questions about the Decay of Passive Scalar Fluctuations

Let us start with some background information on the decay of passive scalar fluctuations in homogeneous turbulence, particularly the decay of scalar variance. Perhaps the first point to make is that, although there are points of contact with Richardson’s celebrated theory of relative dispersion (Richardson [8]), passive scalar decay and Richardson dispersion are distinct phenomena, with the former focusing on the integral scales and the latter focusing on the inertial sub-range (Monin and Yoglam [9]).
Perhaps the most important paper in passive scalar decay is that written by Corrsin [10]. Crucially, Corrsin noted that, under certain conditions, a passive scalar, θ , evolving in freely decaying isotropic turbulence, possesses an integral invariant of the two-point correlation function θ θ = θ ( x ) θ ( x ) . Specifically, he noted that
L θ =   V θ θ d r ,   r = x x ,
is conserved (when it exists), provided that the triple correlation u θ θ falls off faster than O ( r 2 ) at large separations, where r = r = x x . The possible existence of such an invariant immediately raises a host of questions, including the following.
(i)
Under what conditions is this integral convergent? In other words, when does L θ exist?
(ii)
If the integral exists, under what conditions would we expect to find u θ θ O ( r 3 ) , so that L θ is an invariant?
(iii)
If the invariant exists in isotropic turbulence, how does this generalize to anisotropic turbulence?
(iv)
Can we use the existence of the invariant to predict the decay rate of passive scalar variance?
We shall address all of these questions in this review. In particular, we generalize Corrsin’s integral to anisotropic and axisymmetric turbulence and use it to predict the decay rate of passive scalar variance.
Of course, the statistics of θ depend on those of the velocity field, and so the answers to the questions above must depend on the nature of the turbulence in which the scalar finds itself. The strongest long-range correlations are to be found in anisotropic Saffman turbulence, where the two-point velocity correlations fall off slowly as O ( r 3 ) , as discussed in Saffman [11,12] and Davidson, Okamoto, and Kaneda [3]. Consequently, the most stringent test for the existence of Corrsin’s integral invariant is in anisotropic Saffman turbulence, and so that is the case that we consider here. To focus our thoughts, we follow Davidson, Okamoto, and Kaneda [3] and restrict the discussion to statistically axisymmetric (but anisotropic) turbulence.
We shall see that, in axisymmetric Saffman turbulence, L θ is indeed convergent, though only conditionally so. It is also conserved. Moreover, the fact that L θ is conserved in Saffman turbulence strongly suggests that it is an invariant in all other forms of homogeneous turbulence. Finally, we combine the conservation of L θ with the self-similarity of large scales to show that θ 2 l / / θ l θ 2 is an invariant, where l / / θ and l θ are integral scales parallel and perpendicular to the axis of symmetry, respectively. This in turn demands that the passive scalar variance decays as θ 2 ~ t 6 / 5 , at least in Saffman turbulence.

1.3. The Structure of the Paper

The structure of the paper is as follows. We start, in Section 2, by reminding readers about the properties of certain integral invariants (Saffman invariants) associated with the two-point velocity correlations. Such invariants exist in so-called Saffman turbulence (turbulence in which the eddies possess a statistically significant amount of linear impulse) and are a statistical manifestation of linear momentum conservation. These invariants are important as their conservation constrains the way in which the integral scales can evolve, and this in turn has consequences for the rate of energy decay and the rate of growth of integral length scales.
Next, in Section 3 and Section 4, we generalize Corrsin’s passive scalar integral to anisotropic axisymmetric turbulence and note that its invariance reflects the conservation of the total amount of passive scalar within a large control volume embedded in the turbulence. We also note the analogy between our generalized Corrsin integral and the Saffman invariants.
The results of Section 2, Section 3 and Section 4 are then combined in Section 5, where we use these various integral invariants to determine the behavior of integral scales in axisymmetric turbulence. In particular, we deduce the rate of energy decay and passive scalar variance, as well as the evolution of various integral scales. We close, in Section 6 and Section 7, with the discussion and conclusions, respectively. As noted earlier, the primary point of novelty in the paper lies in the use of our generalized Corrsin integral to deduce the integral properties of a passive scalar.

2. Invariants of Freely Decaying, Axisymmetric Saffman Turbulence

As a prelude to our discussion of invariants in passive scalar mixing, it seems appropriate to review the large-scale dynamics of freely decaying Saffman turbulence. The review is based on the results of Saffman [11] for the most general anisotropic case, and on the findings of Davidson [5] and Davidson, Okamoto, and Kaneda [3], hereafter denoted as DOK12, for the case where the fields are statistically axisymmetric.

2.1. Anisotropic Saffman Turbulence

Saffman [11] considered homogeneous turbulence which emerges from an initial condition in which there are no long-range algebraic correlations in the vorticity field. Thus, he assumed that ω i ω j r and ω i ω j ω k r , as well as all higher-order long-range vorticity cumulants, are exponentially small at t = 0. (From now on, we shall abbreviate expressions such as ω i ( x ) ω j ( x ) r and ω i ω j ). This leads to a low-wavenumber energy spectrum of the form E ( k 0 ) ~ k 2 and so provides an alternative to the more traditional analysis of Batchelor and Proudman [13], who found E ( k 0 ) ~ k 4 . Interestingly, recent measurements taken by Krogstad and Davidson [4] support the idea that, at least for certain types of grid turbulence, the dynamics of freely decaying turbulence are indeed consistent with E ( k 0 ) ~ k 2 , rather than the more traditional assumption of E ( k 0 ) ~ k 4 .
The spectral tensor,
Φ i j ( k ) = 1 2 π 3 u i u j ( r ) e j k r d r ,
in Saffman turbulence takes the low-k form
Φ i j ( k 0 ) = δ i α k i k α k 2 δ j β k j k β k 2 M α β + O ( k ) ,
where M α β is symmetric and independent of k. Moreover, the coefficients M α β is related to the second moment of the two-point vorticity correlation, ω i ω j , by
( 2 π ) 3 M i j = 1 2 Ω i j 1 4 δ i j Ω k k ,   Ω i j = r 2 ω i ω j d r ,
(Saffman [11]). Evidently, Φ i j ( k ) is non-analytic at k = 0 as the value depends on the way in which we approach the origin. Expansion (2) yields
Φ i i ( k ) = δ α β k α k β k 2 M α β + O ( k ) ,   E ( k 0 ) = 4 3 π M α α k 2 ,
and the corresponding asymptotic form of u i u j is
u i u j = M α β π 2 δ i α 2 2 r i r α δ j β 2 2 r j r β r .
Note that u i u j = O ( r 3 ) .
The origin of the r 3 tails in u i u j reverts back to Biot–Savart law, in which an eddy with a finite linear impulse induces irrotational far-field velocity fluctuations of the form
u x = 1 4 π L eddy 1 / r + O ( r 4 ) ,
where L eddy is the linear impulse of the eddy, L eddy = 1 2   V eddy   x × ω d V . Thus, individual eddies cast a long shadow which falls off slowly as u = O ( r 3 ) , leading to the long-range algebraic correlations u i u j = O ( r 3 ) and u i u j u k = O ( r 3 ) . Since these far-field fluctuations are irrotational, they do not contribute to the leading-order terms in ω i ω j , which are of order ω i ω j O ( r 6 ) , rather than ω i ω j O ( r 5 ) , for t > 0 (Saffman [11]).
Since ω i ω j O ( r 6 ) , the integral Ω i j defined by (3) is absolutely convergent, and so M i j is well defined. On the other hand, the integrals
L i j = Lim V V u i u j d r
are only conditionally convergent. That is to say, they are convergent in the sense that they are independent of the size of V as V . They are also uniquely determined by M i j , since L i j may be written as a surface integral whose integrand is the far-field Expansion (5):
L i j = Lim V V u i u j d r = Lim S S u i u k r j d S k .
However, the numerical value of L i j depends on the shape of the surface S as it recedes to infinity. This conditional convergence is a direct consequence of the non-analytic form of Φ i j ( k ) at k = 0 . For the particular case of a large spherical volume of radius R and volume V R , Saffman [11] established that
L i j = Lim V R V R u i u j d r = ( 2 π ) 3 7 15 M i j + 1 15 δ i j M k k .
The integral L i j is particularly important because, as Saffman [11] suggested, L i j , and hence M i j , are invariants:
L i j = constant ,   M i j = constant .
This idea is as follows. Homogeneity tells us that
L i j = Lim V R 1 V R 1 / 2 V R u i d x 1 V R 1 / 2 V R u j d x
(see, for example, Davidson [6]), where V R is a large spherical control volume with surface S R and radius R. So, the dynamical behavior of L i j depends on that of the linear momentum P = u V R d V . However, P can only change as a result of pressure forces acting on the surface S R , or because momentum is convected across S R . Moreover, these are both random events spread over a large surface, and so the central limit theorem tentatively suggests that their cumulative effects are scaled as O S R 1 / 2 ~ R , i.e., the momentum change in any finite time is O V R 1 / 3 t . This in turn suggests that the pressure forces and momentum fluxes across S R are too weak to change P / V R 1 / 2 in the limit of V R , and that consequently L i j is an invariant in accordance with (11). This is essentially the argument of Saffman [11].
Unfortunately, it is rare that one satisfies the strict requirements for the central limit theorem to hold, especially in turbulence where the long-range interactions mean that remote events are rarely statistically independent. Consequently, we must regard the argument above as merely suggestive. However, formal proof of the invariance of L i j is given in the study conducted by Davidson [5,6].

2.2. Statistically Axisymmetric Saffman Turbulence

We now turn to the case of statistically axisymmetric Saffman turbulence, with or without reflectional symmetry. In particular, we need to summarize the findings of Davidson [5,6] and DOK12, as these directly lead to invariants of passive scalar mixing in axisymmetric turbulence. As usual, we use subscripts / / and to indicate quantities that are parallel and perpendicular to the symmetry axis, respectively, which we take to be the z axis. Here, it may be shown that M i j = 0   if i j , and so the only non-zero components of M i j are M / / and M x x = M y y = 1 2 M . This supports the idea that M / / and M are the only invariants of axisymmetric Saffman turbulence. Similarly, (3) tells us that Ω i j is diagonal, and that
( 2 π ) 3 M / / = 1 4 Ω / / 1 4 Ω ,   ( 2 π ) 3 M = 1 2 Ω / / .
Therefore, the most general form of Φ i j ( k 0 ) for Saffman turbulence is:
Φ / / = Φ z z = M / / k 4 k 4 + 1 2 M k z 2 k 2 k 4 + O ( k ) ,
Φ = Φ x x + Φ y y = M / / k z 2 k 2 k 4 + 1 2 M 1 + k z 4 k 4 + O ( k ) ,
(Davidson [5]), where k 2 = k 2 k z 2 , and we omit the off-diagonal terms. Thus, with regards to the spectral tensor, our two invariants can be written as follows:
M / / = Φ / / ( k z = 0 , k 0 ) = 1 ( 2 π ) 3 Lim k 0 e j k r u / / u / / d r = constant ,
1 2 M = Φ ( k z = 0 , k 0 ) = 1 ( 2 π ) 3 Lim k 0 e j k r u u d r = constant ,
M = Φ ( k = 0 , k z 0 ) = 1 ( 2 π ) 3 Lim k z 0 e j k z r z u u d r = constant ,
while
Φ / / ( k = 0 , k z 0 ) = 0 .
Similarly, the only non-zero components of L i j are L / / and L . As for the general anisotropic case, L / / and L are uniquely determined by M i j , although the relationship between L / / , L , M / / , and M depends on the shape of V in (7). For example, in the particular case of a large cylindrical control volume of radius R and height 2H, the following can be deduced:
L / / = ( 2 π ) 3 M / / 1 + 1 2 ( R / H ) 2 1 + ( R / H ) 2 3 / 2 + ( 2 π ) 3 M 1 4 ( R / H ) 2 1 + ( R / H ) 2 3 / 2 ,
L = ( 2 π ) 3 M 1 1 2 + 3 4 ( R / H ) 2 1 + ( R / H ) 2 3 / 2 + ( 2 π ) 3 M / / 1 2 ( R / H ) 2 1 + ( R / H ) 2 3 / 2 ,
(see DOK12), from which
L / / H / R = ( 2 π ) 3 M / / = constant , ( long cylinder ) ,
L H / R = ( 2 π ) 3 1 2 M = constant , ( long cylinder ) ,
L / / R / H = 0 , ( thin disc ) ,
L R / H = ( 2 π ) 3 M = constant , ( thin disc ) .
These are equivalent to the spectral results (15)→(17), where H / R corresponds to k z = 0       and       k 0 , while R / H is equivalent to k = 0       and       k z 0 . The constraint that L i j = constant is now reduced to
L / / = u / / u / / d r = constant ,
L = u u d r = constant ,
and this is true, irrespective of the shape of V. The invariance of L / / and L , or equivalently M / / and M , is confirmed in the direct numerical simulations of DOK12.

3. Passive Scalar Statistics in Axisymmetric Saffman Turbulence

Let us now turn to passive scalar mixing in Saffman turbulence. Let θ be the passive scalar which has a mean of zero, θ = 0 , and is governed by
D θ D t = α 2 θ .
We shall assume that, like the velocity field, its initial energy spectrum takes the form E θ ( k 0 ) ~ k 2 . (As we shall see, such an initial condition ensures that the spectrum E θ ( k 0 ) ~ k 2 then persists for all time.) In such a case, the spectral function
Φ θ ( k ) = 1 2 π 3 θ θ ( r ) e j k r d r ,
takes the low-k form
Φ θ ( k 0 ) = δ α β k α k β k 2 N α β + O ( k ) ,
and
E θ ( k 0 ) = 4 3 π N α α k 2 ,
where N α β is symmetric and independent of k. This has the same structural form as Φ i i ( k 0 ) in (4), provided that we substitute N α β for M α β , and so we conclude that θ θ = O ( r 3 ) .
In statistically axisymmetric turbulence, N i j = 0   if i j , and so the only non-zero components of N i j are N / / and N x x = N y y = 1 2 N . The low-k form of the spectral function Φ θ ( k ) is then simplified to
Φ θ ( k 0 ) = N / / + 1 2 N N / / 1 2 N k / / 2 k 2 + O ( k ) ,
which is analytic at k = 0 only for the isotropic case of N / / = 1 2 N . Note that N / / and N can be written in terms of the spectral function as
N / / + 1 2 N = Φ θ ( k / / = 0 , k 0 ) = 1 ( 2 π ) 3 Lim k 0 e j k r θ θ d r ,
and
N = Φ θ ( k = 0 , k / / 0 ) = 1 ( 2 π ) 3 Lim k / / 0 e j k / / r / / θ θ d r .
Thus, we see that N / / and N are closely related to the integral Lim V   V θ θ d r , which we shall shortly see is an invariant of scalar mixing.
Let us now turn to the Corrsin integral
L θ = Lim V   V θ θ d r ,
which is central to scalar mixing analysis. Since θ θ = O ( r 3 ) , L θ is, in general, only conditionally convergent, and so it is important to determine the precise form of θ θ in the far-field. Here, the analogy to Φ i i ( k ) is again helpful. Based on (5), we have
u u = π 2 2 M / / M r 2 3 r / / 2 r 5 ,
and so the structural similarity between (4) and (28) tells us that
θ θ = π 2 2 N / / N r 2 3 r / / 2 r 5 .
Note that the long-range correlation, θ θ = O ( r 3 ) , vanishes only in the isotropic case.
The key point about (34) is that its right-hand side may be written as a divergence of a vector of magnitude r 2 ,
θ θ = π 2 2 N / / N   r / / 1 r e ^ z .
Thus, if we integrate θ θ ( r ) over the volume bounded by two large concentric surfaces of the same shape but different sizes, the net result is zero. Moreover, it remains zero as the outer surface recedes to infinity. This is important because it means that L θ is convergent, if only conditionally convergent (in the sense defined above). For example, in the case of a large cylindrical control volume of radius R and height 2H, (30) and (31) suggest the following:
θ θ d r = ( 2 π ) 3 N / / + 1 2 N   ,                   ( R / H 0 ,   long cylinder ) ,
θ θ d r = ( 2 π ) 3 N   ,                           ( R / H ,   thin disc ) .
In summary, in the isotropic case of N / / = 1 2 N , there are no long-range correlations of the form θ θ = O ( r 3 ) , and so L θ is absolutely convergent and given by either (30) or (31):
L θ = θ θ d r = ( 2 π ) 3 N / / + 1 2 N = ( 2 π ) 3 N .
However, in the anisotropic case of N / / 1 2 N , we have θ θ = O ( r 3 ) and L θ is only conditionally convergent, and its value depends on the shape of V.

4. An Invariant for Passive Scalar Mixing in Axisymmetric Saffman Turbulence

Crucially, integral L θ is an invariant of freely decaying Saffman turbulence. This may be shown as follows. From (26), we have the two-point evolution equation:
t θ θ = r i ( δ u i ) θ θ + 2 α 2 θ θ ,
where δ u = u u is a velocity increment. By integrating over a volume V in r-space, and noting that θ θ O ( r 3 ) , we obtain
d d t L θ = S δ u θ θ d S ,
where S is the surface of V as that surface recedes to infinity. However, in Saffman turbulence, the long-range correlations are no stronger than order r−3, and so δ u θ θ O ( r 3 ) . The surface integral in (40) now vanishes, yielding
L θ = Lim V   V θ θ d r = constant .
The physical interpretation of (40) and (41) is clear. Homogeneity allows us to rewrite L θ as
L θ = Lim V   V θ θ d r = Lim V 1 V V θ d V 2 ,
and so L θ is conserved if the scalar flux across the bounding surface of V, say S, is less than order V 1 / 2 (or equivalently S 3 / 4 ) as V → ∞. However, the central limit theorem suggests that this flux should be no larger that S 1 / 2 , consistent with L θ being an invariant. In fact, (40) is nothing more than a concise form of the scalar budget equation
d d t L θ = d d t 1 V V θ d V 2 = 2 V V θ d V S θ u d S = 2 V V S θ θ u d S d V ,
and the vanishing of the surface integral in (40) amounts to the statement that there is negligible convection of the passive scalar in or out of V in the limit of V .
Perhaps some comments are in order at this point. First, since N / / and N are related to L θ through (30) and (31), we would expect the invariance of L θ to translate to the statement that N / / and N are both invariants. This is indeed the case (Yoshimatsu and Kaneda, 2018 [14]). Second, if l / / θ and l θ are integral scales for the passive scalar, say
l / / θ = 1 θ 2 θ ( x ) θ ( x + r e ^ z ) d r ,
l θ = 1 θ 2 θ ( x ) θ ( x + r e ^ x ) d r ,
then, since L θ , N / / , and N are all dominated by large scales, we have
L θ ~ N / / ~ N ~ θ 2 l / / θ l θ 2 .
Moreover, when the large scales in the passive scalar field evolve in a self-similar way, i.e., when θ θ / θ 2 is a function of r / / / l / / θ and r / l θ only, then the invariance of N / / and N demands
θ 2 l / / θ l θ 2 = constant ,
which places a strong constraint on the evolution of the scalar variance θ 2 .
Third, the constraint (47) is similar to that imposed on the velocity field as a result of (15)→(17). That is to say, M / / and M are dominated by large scales and so
M / / ~ u / / 2 l / / u l u 2 ,   M ~ u 2 l / / u l u 2 ,
where l / / u and l u are suitable integral scales for the velocity field (see Section 5). Moreover, when the large-scale velocity correlations are self-similar, as is typical of fully developed turbulence, (15) →(17) demand
u / / 2 l / / u l u 2 = constant ,   u 2 l / / u l u 2 = constant .
The existence of invariants M / / and M (or, equivalently, L / / and L ), as well as the associated conservation relationships (49), is confirmed in DOK12.
Finally, we combine δ u θ θ O ( r 3 ) with (39). This tells us that, if θ θ = O ( r 3 ) at t = 0, then we have θ θ = O ( r 3 ) for all t > 0. Thus, the initial Saffman-like spectrum E θ ( k 0 ) ~ k 2 is preserved for all time, as noted above. Indeed, this has to be the case since
E θ ( k 0 ) = 4 3 π N α α k 2 ,
and N i i = N / / + N is an invariant. Conversely, if the initial spectrum takes an alternative form, say E θ ( k 0 ) ~ k 4 , then the Saffman spectrum E θ ( k 0 ) ~ k 2 remains kinematically inaccessible to the passive scalar.

5. The Predicted Decay Rate of Energy and Scalar Variance

We shall now use the invariants M / / , M , N / / , and N , or equivalently L / / , L , and L θ , combined with the assumption of self-similarity of large scales, to predict the evolution of kinetic energy and scalar variance in freely evolving turbulence. The approach is the same as that used by Saffman [12] and DOK12, where the invariance of M / / and M is used to predict the decay rate of u / / 2 and u 2 .

5.1. The Rate of Energy Decay

Before discussing the decay of scalar variance, we should note how the kinetic energy and integral length scales l / / u and l u evolve, as the passive scalar is slave to the velocity field.
Let us start by introducing some notation. We define the velocity integral scales as follows:
l / / u = 1 u / / 2 u / / ( x ) u / / ( x + r e ^ z ) d r ,   u / / = u / / 2 1 / 2 ,
l u = 1 2 u 2 u x ( x ) u x ( x + r e ^ x ) d r + u y ( x ) u y ( x + r e ^ y ) d r ,   u = 1 2 u 2 1 / 2 .
Next, we assume that, once the turbulence is fully developed, the large scales in the velocity field (scales of order l / / u and l u ) are self-similar, i.e., u / / u / / / u / / 2 and u u / u 2 are functions of r / / / l / / u and r / l u only. This is certainly observed in the numerical experiments of DOK12. We shall refer to this as partial self-similarity, as it does not extend down to Kolmogorov scales, nor perhaps to the far-field correlations u i u j given by (5). Now, M / / and M are absolutely convergent and so are independent of the far-field contribution u i u j . (The long-wavelength oscillation associated with e j k r and e j k / / r / / kills off the far-field contributions to these integrals, leaving (15)→(17) well behaved.) It follows that partial self-similarity may be applied to (15)→(17), which then demands that, in fully developed turbulence,
M = c u 2 l u 2 l / / u = constant ,   M / / = c / / u / / 2 l u 2 l / / u = constant ,
where c and c / / are dimensionless constants. In such a situation, we have
u 2 / u / / 2 = constant .
This suggests that, rather surprisingly, any velocity anisotropy which is present at the beginning of the fully developed stage will be locked into the turbulence thereafter. The persistence of anisotropy in the form of (54) is confirmed in the numerical simulations of DOK12.
We now estimate the decay rate of energy under the assumption that partial self-similarity holds and that, as usual, the flux of energy to small scales is controlled by large scales. Since u / / 2 / u 2 = constant , we only need to consider one component of energy, say u 2 , and dimensional analysis then gives us
d u 2 d t = A u 3 l u   ,                     A = A l / / u / l u , [ fully developed turbulence only ]
where A is a dimensionless function of l / / u / l u . Now, certain types of axisymmetric turbulence satisfy u / l u ~ u / / / l / / u as a result of continuity. In such cases, the constraint u / / 2 / u 2 = constant demands that l / / u / l u = constant , and this is certainly observed in the fully developed state in the numerical experiments of DOK12. Equations (53) and (55) then simplify to
u 2 l u 3 = constant ,   d u 2 d t = A u 3 l u ,
which results in the following decay laws:
u 2 ~ u / / 2 ~ t 6 / 5 ,   l u ~ l / / u ~ t 2 / 5 , [ fully developed turbulence only ]
as suggested by Saffman [12], and as observed in the direct numerical simulations of DOK12.

5.2. The Decay of Scalar Variance

Given (57), we can now predict the decay rate of scalar variance in fully developed turbulence. The analysis differs somewhat depending on whether l / / u and l u are greater or smaller than the scalar integral scales l / / θ and l θ . Consider first the case where the initial conditions have l min θ >> l max u . (Here, the subscripts indicate whether the minimum or maximum of the two possible length scales is taken.) Then, mixing occurs on the scale of l max u , and dimensional analysis, combined with u / / 2 / u 2 = constant and l / / u / l u = constant , gives us
d θ 2 d t = C u l u θ 2 ,
where C is a constant. Given (57), and the fact that mixing occurs on a scale much smaller than l min θ , we expect θ 2 to decay as a power law, while l / / θ   and   l θ   will be more or less constant, and this will continue for as long as l min θ >> l max u . Note that the scalar field is not self-similar during this stage, so that (47) does not apply. Eventually, l max u will grow to be of the same order as l min θ , and thereafter l / / θ and l θ will grow along with l / / u and l u , as t 2 / 5 . If the scalar spectrum exhibits partial self-similarity during this growth, then invariant (47) combined with (57) yields
θ 2 ~ t 6 / 5 ,   l θ ~ l / / θ ~ t 2 / 5 . [ fully   developed   turbulence ,   l θ   ~   l u ] .
Now consider the initial conditions in which l max θ << l min u . In this case, mixing is no longer accomplished by the integral-scale eddies, but rather by inertial-range eddies of scale l max θ . In such a situation, we have
d θ 2 d t ~   C       θ 2 τ θ   ,             C = C l / / θ / l θ ,
where C is a dimensionless function of l / / θ / l θ and τ θ is the eddy turnover time of scales of order l max θ . Now, τ θ can be estimated from Kolmogorov’s 2/3 law as l max θ 2 / 3 / ε 1 / 3 . (Here, ε is the rate of dissipation of kinetic energy.) Since the anisotropy of the velocity integral scales is fixed in fully developed turbulence, i.e., u / / 2 / u 2 = constant and l / / u / l u = constant , we would expect the velocity anisotropy at scale l max θ to also be constant. This in turn suggests that l / / θ / l θ will be constant, or equivalently l / / θ ~ l θ ~ l θ , provided that the velocity anisotropy is not too large at scale l max θ . In such cases, (47) and (60) simplify to
d θ 2 d t ~ C   ε 1 / 3 θ 2 l θ   2 / 3       ,                 θ 2 l θ   3 = constant ,
from which
d θ 2 d t ~   ε 1 / 3 θ 2 11 / 9 ,
where (57) has   ε ~ t 11 / 5 . Once again, we conclude that
θ 2 ~ t 6 / 5 ,   l θ ~ l / / θ ~ t 2 / 5 . [ fully   developed   turbulence ,   l θ   < <   l u ]
In summary, we would expect the scalar variance to decay as θ 2 ~ t 6 / 5 , and the scalar integral scales to grow as l θ ~ l / / θ ~ t 2 / 5 , whenever the turbulence is fully developed and l θ O( l u ). These various predictions are confirmed in the direct numerical simulations of Yoshimatsu and Kaneda [7].

6. Discussion

Our understanding of the large-scale dynamics of homogeneous turbulence has come a long way since those early discussions at NCAR in 2002, and this author has greatly enjoyed the journey, particularly those collaborations with Jack Herring, Yukio Kaneda, and Per-Age Krogstad. We now have a reasonably complete understanding of large-scale dynamics, both in terms of kinetic energy and passive scalar mixing, and this is true of both Batchelor E ( k 0 ) ~ k 4 and Saffman E ( k 0 ) ~ k 2 turbulence.
As noted by Kolmogorov, in Batchelor turbulence, the kinetic energy should decay as u 2 ~ t 10 / 7 , provided that Loitsyansky’s integral is conserved. This decay law was confirmed in the direct numerical simulations of Ishida, Davidson, and Kaneda [2], but only for fully developed turbulence. (Loitsyansky’s integral is not conserved during an initial transient, before the turbulence becomes fully developed.) Similarly, for isotropic Saffman turbulence, the prediction of u 2 ~ t 6 / 5 is observed in the wind tunnel experiments, as stipulated by Krogstad and Davidson [4]. Finally, in anisotropic (but axisymmetric) Saffman turbulence, the theoretical predictions of
u 2 ~ u / / 2 ~ t 6 / 5 ,   l u ~ l / / u ~ t 2 / 5 ,
are well supported by the direct numerical simulations of DOK12, while the estimates
θ 2 ~ t 6 / 5 ,   l θ ~ l / / θ ~ t 2 / 5 ,
are consistent with the direct numerical simulations of Yoshimatsu and Kaneda [7], as discussed above.
Ironically, the one question which has proved to be most resistant to theoretical attack turns out to be the very subject of those early discussions with Jack Herring back in 2002, i.e., the observation that Loitsyansky’s integral is effectively constant in fully developed Batchelor turbulence. The point is this: Batchelor and Proudman [13] noted that the long-range velocity correlations in E ( k 0 ) ~ k 4 turbulence fall off sufficiently slowly with r so as to invalidate the usual proof that Loitsyansky’s integral is an invariant. However, the direct numerical simulations show that, once the turbulence is fully developed, but not before, Loitsyansky’s integral is indeed constant (Herring et al. [1]; Ishida, Davidson, and Kaneda, [2]). It seems that there is something particular about the structure of fully developed turbulence that suppresses the long-range interactions of Batchelor and Proudman [13], and we still do not understand what that might be. Evidently, there is still much to accomplish, and so this is a story without an ending.

7. Conclusions

We reviewed the integral invariant theory of Saffman turbulence, showing that the various statistical invariants, both the Saffman and Corrsin integrals, have a simple physical interpretation. We have also shown that, under a wide range of conditions, the integral scales of a passive scalar evolve as
θ 2 ~ t 6 / 5 ,   l θ ~ l / / θ ~ t 2 / 5 .
The primary point of novelty of our approach lies in its use of various integral invariants to deduce the evolution of integral scales.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The author would like to thank the Leverhulme Trust for its generous support through the grant RPG-2015-195/RG77943, as well as Yukio Kaneda for their fruitful discussions on passive scalar mixing.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Herring, J.R.; Kimura, Y.; James, R.; Davidson, P.A. Statistical and dynamical aspects of isotropic and stably stratified turbulence. In Mathematical and Physical Theory of Turbulence; Cannnon, J., Shivamoggi, B.K., Eds.; Chapman & Hall: London, UK, 2006. [Google Scholar]
  2. Ishida, T.; Davidson, P.A.; Kaneda, Y. On the decay of isotropic turbulence. J. Fluid Mech. 2006, 564, 455–475. [Google Scholar] [CrossRef]
  3. Davidson, P.A.; Okamoto, N.; Kaneda, Y. On freely-decaying, anisotropic, axisymmetric, Saffman turbulence. J. Fluid Mech. 2012, 706, 150–172. [Google Scholar] [CrossRef]
  4. Krogstad, P.-A.; Davidson, P.A. Is grid turbulence Saffman turbulence? J. Fluid Mech. 2010, 642, 373–394. [Google Scholar] [CrossRef]
  5. Davidson, P.A. On the Decay of Saffman Turbulence Subject to Rotation, Stratification or an Imposed Magnetic Field. J. Fluid Mech. 2010, 633, 268–292. [Google Scholar] [CrossRef]
  6. Davidson, P.A. Turbulence in Rotating, Stratified and Electrically Conducting Fluids; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  7. Yoshimatsu, K.; Kaneda, Y. Anisotropy freezing of passive scalar fields in anisotropy growing homogeneous turbulence. Phys. Rev. Fluids 2020, 5, 014604. [Google Scholar] [CrossRef]
  8. Richardson, L.F. Atmospheric diffusion shown on a distance-neighbour graph. Proc. R. Soc. A 1926, 110, 709–737. [Google Scholar]
  9. Monin, A.S.; Yaglom, A.M. Statistical Fluid Mechanics II; MIT Press: Cambridge, MA, USA, 1975. [Google Scholar]
  10. Corrsin, S. The decay of isotropic temperature fluctuations in an isotropic turbulence. J. Aeronaut. Sci. 1951, 18, 417–423. [Google Scholar] [CrossRef]
  11. Saffman, P.G. The large-scale structure of homogeneous turbulence. J. Fluid Mech. 1967, 27, 581–593. [Google Scholar] [CrossRef] [Green Version]
  12. Saffman, P.G. Note on decay of homogenous turbulence. Phys. Fluids 1967, 10, 1349. [Google Scholar] [CrossRef] [Green Version]
  13. Batchelor, G.K.; Proudman, I. The large-scale structure of homogeneous turbulence. Philos. Trans. R. Soc. Lond. A 1956, 248, 369–405. [Google Scholar]
  14. Yoshimatsu, K.; Kaneda, Y. Large-scale structure of velocity and passive scalar fields in freely decaying homogeneous anisotropic turbulence. Phys. Rev. Fluids 2018, 3, 104601. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Davidson, P.A. The Decay of Energy and Scalar Variance in Axisymmetric Turbulence. Atmosphere 2023, 14, 1019. https://doi.org/10.3390/atmos14061019

AMA Style

Davidson PA. The Decay of Energy and Scalar Variance in Axisymmetric Turbulence. Atmosphere. 2023; 14(6):1019. https://doi.org/10.3390/atmos14061019

Chicago/Turabian Style

Davidson, Peter A. 2023. "The Decay of Energy and Scalar Variance in Axisymmetric Turbulence" Atmosphere 14, no. 6: 1019. https://doi.org/10.3390/atmos14061019

APA Style

Davidson, P. A. (2023). The Decay of Energy and Scalar Variance in Axisymmetric Turbulence. Atmosphere, 14(6), 1019. https://doi.org/10.3390/atmos14061019

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop