Advances in Mathematical Physics

Volume 2018, Article ID 2365602, 12 pages

https://doi.org/10.1155/2018/2365602

## Statistics of Finite Scale Local Lyapunov Exponents in Fully Developed Homogeneous Isotropic Turbulence

“La Sapienza” University, Dipartimento di Ingegneria Meccanica e Aerospaziale, Via Eudossiana, 18, 00184 Rome, Italy

Correspondence should be addressed to Nicola de Divitiis; moc.liamg@siitivided.n

Received 21 February 2018; Accepted 2 May 2018; Published 7 June 2018

Academic Editor: Xavier Leoncini

Copyright © 2018 Nicola de Divitiis. 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.

#### Abstract

The present work analyzes the statistics of finite scale local Lyapunov exponents of pairs of fluid particles trajectories in fully developed incompressible homogeneous isotropic turbulence. According to the hypothesis of fully developed chaos, this statistics is here analyzed assuming that the entropy associated with the fluid kinematic state is maximum. The distribution of the local Lyapunov exponents results in an unsymmetrical uniform function in a proper interval of variation. From this PDF, we determine the relationship between average and maximum Lyapunov exponents and the longitudinal velocity correlation function. This link, which in turn leads to the closure of von Kármán–Howarth and Corrsin equations, agrees with results of previous works, supporting the proposed PDF calculation, at least for the purposes of the energy cascade main effect estimation. Furthermore, through the property that the Lyapunov vectors tend to align the direction of the maximum growth rate of trajectories distance, we obtain the link between maximum and average Lyapunov exponents in line with the previous results. To validate the proposed theoretical results, we present different numerical simulations whose results justify the hypotheses of the present analysis.

#### 1. Introduction

In fully developed turbulence, the finite scale Lyapunov exponents of the fluid kinematic field are of paramount importance because they (i) describe the turbulent energy cascade phenomenon and (ii) give the fluid viscous dissipation when the length scale goes to zero. One of the characteristics of these exponents in turbulence is that their statistics is related to the instantaneous velocity field, whereas they do not depend directly on the time variations of this latter. Such characteristic, which represents the crucial point of this work, is consequence of the fact that the times of variations of the velocity field, which changes according to the Navier–Stokes equations, are much greater than those of the fluid displacements which follow the fluid kinematics [1, 2]. Specifically, whereas the velocity field is a function of slow growth of the time, the distance between particles and the local fluid deformation are unbounded quantities which rise exponentially with the time [1]. Accordingly, the fluid particles displacements statistics, being not directly linked to the time variations of the velocity field, can be considered to be the result of the current fluid act of motion.

It is worth remarking that the adoption of the finite scale Lyapunov exponents in place of the classical ones is justified by the fact that the turbulence is a complex phenomenon involving numerous length scales; each of them is characterized by different properties. The several perturbations of finite size will vary following nonlinear differential equations out of tangent space; thus the sole use of classical Lyapunov exponents is not adequate to describe the perturbations behavior associated with the different length scales and therefore the energy cascade.

Next, as a result of item (i), the knowledge of the Lyapunov exponents statistics would lead to the closure of the equations of velocity and temperature correlations and therefore to the determination of kinetic energy and temperature spectra. Therefore, the distribution of such exponents is very important for quantifying the effects of the turbulent energy cascade.

Several works dealing with the closure of the correlations equations are present in the literature, for instance, [3–9]. Nevertheless, to the author’s knowledge, the influence of the finite scale Lyapunov exponents statistics on turbulent energy cascade and on the closure of the correlations equations has not received due attention. Hence, the objective of the present work is to develop a theoretical analysis based on the aforementioned properties which leads to determining the statistics of the local finite scale Lyapunov exponent in fully developed homogeneous isotropic turbulence for incompressible fluids. This local exponent defined as provides the instantaneous growth rate of the distance between two fluid particles trajectories and , where is the separation vector (finite scale Lyapunov vector). This exponent is linked to the so called finite size Lyapunov exponent, FSLE, defined in [10] in a more general framework of dynamic systems. Specifically, FSLE can be obtained in terms of as a proper average of over a given interval, say , where is the subset in which rises from to , is an assigned threshold slightly greater than the unity to avoid interferences between the scales, and , , are rescaled along the direction .

Because of nonsmooth spatial variations of the velocity field, can exhibit fluctuations of sizable amplitude with respect to its average value; thus plays the role of a stochastic variable and will be distributed according to a certain PDF. To define the magnitude of these oscillations, average and maximum finite scale Lyapunov exponents, and respectively, are defined aswhere and denote the average over the entire ensemble of and the average calculated on the ensemble where .

As a consequence of the aforementioned properties, the present analysis assumes that the kinematics of a pair of fluid particles, characterized by , is much faster and statistically independent with respect to the time variations of velocity field. This property, just discussed in [1, 2] for what concerns the closure of von Kármán–Howarth and Corrsin equations [11–13], was previously supported by the arguments presented in [14] (and references therein), where the author observes the following: (i) the velocity fields produce chaotic trajectories also for relatively simple mathematical structure of ; (ii) the flows given by stretch and fold continuously and rapidly causing an effective mixing of the particles trajectories.

Through the hypotheses of fully developed chaos and fluid incompressibility, we first estimate the interval of variation of and thereafter determine the distribution function of by maximizing the entropy associated with the fluid kinematic state. The maximization of such entropy is justified by the fact that the regime of fully developed turbulence corresponds to a situation of maximum chaos where the bifurcations cause a total loss of the initial condition data of and . As a consequence, we show that will be uniformly distributed in such interval; in particular, we determine the relationship between and , resulting in .

A further confirmation of such link is obtained by exploiting the alignment property of , following which the Lyapunov vectors tend to align the direction of the maximum growth rate of [15].

In order to compare the proposed statistics with the average energy cascade effects, and are expressed in terms of the longitudinal velocity correlation function , where , , and . This relationship leads to the closure formulas of von Kármán–Howarth and Corrsin equations and coincides with that just presented in [1], where the author adopts only and without considering the distribution of the local exponents. There, the proposed closure formula leads to a value of the skewness of equal to −3/7, in agreement with those obtained by the several authors with direct numerical simulation of the Navier–Stokes equations (DNS) [16–18] and Large–eddy simulations (LES) [19–21]. Therefore, the here proposed statistics should be adequate for describing the distribution of , at least for what concerns the main properties of the energy cascade phenomenon. The novelty of this work with respect to the literature and in particular with respect to [1] is represented by the statistics of and its distribution function. This latter provides much more detailed statistical information about the energy cascade and viscous dissipation than the analysis of [1] which, adopting average exponents, describes only the mean effects. The present analysis provides, by means of such PDF, the estimation of all the Lyapunov exponent statistical moments. Accordingly, this study recovers, among the other things, the results of [1] and, in addition, gives all the statistical properties arising from this specific PDF, in particular, the deviations of such effects of dissipation and energy cascade with respect to their average values.

Furthermore, to justify the hypotheses of the proposed statistics, the theoretical results of this latter are compared with numerical simulations of a proper differential system representing the incompressible fluid kinematics.

#### 2. Interval of Variation of in Incompressible Turbulence

This section proposes an analysis for estimating the range of variations of based on the hypotheses of fluid incompressibility and fully developed chaos.

The present analysis starts from the consideration that the turbulent energy cascade is related to the fluid particles trajectories divergence; therefore the relative fluid kinematics plays an important role in the estimation of the properties of such energy cascade [1, 2]. The relative kinematics is expressed by finite scale Lyapunov vector which satisfies where varies according to the Navier–Stokes equations and and are two fluid particles trajectories. The local divergence between and is quantified by which, due to the bifurcations of the kinematic field (kinematic bifurcations), exhibits oscillations of sizable amplitude with respect to its average value. These bifurcations continuously happen in those points of the space where Observe that the kinematic bifurcations defined by (4) are not the Navier–Stokes equations bifurcations (dynamic bifurcation) but arise from these latter [2]. In fact, the Navier–Stokes bifurcations frequently occur determining continuously nonsmooth spatial variations of which in turn lead to condition (4) in the several points of the space.

The definition of given by (1) implies that can be locally expressed as where plays the role of the stochastic variable and is a proper orthogonal matrix providing the orientation of with respect to the inertial frame . Accordingly, is in which defines the angular velocity of with respect to .

For an assigned velocity field, finite scale Lyapunov exponents and vectors are formally calculated with (1) through the following orthogonalization procedure: (a) the maximal local Lyapunov exponent, say , is first obtained by choosing the direction which maximizes in (1), (b) the second exponent is calculated by selecting the direction in the subspace (plane) orthogonal to which maximizes , and (c) finally, the third one corresponds to the direction normal to both and , where . The so obtained vectors system defines a rigid space which moves with respect to with a given angular velocity depending on the local fluid motion. Therefore, and , , are locally expressed byThe classical local Lyapunov exponents are defined for , .

Now, in order to estimate the set of variations of , observe that, due to fluid incompressibility, and the classical local exponents obey to the following condition:In general, (8) does not hold for finite scale Lyapunov vectors. Nevertheless, (8) is valid for those finite scale vectors for which the volume is locally preserved: therefore, without lack of generality, these exponents can be written in the form where and are variables depending on the current act of motion and Now, following the hypothesis of fully developed chaos, it is reasonable that ranges in the set , where assumes its maximum value compatible with (9) and is consequently calculated. This implies that , that is, It is worth remarking that depends on the instantaneous velocity field which in turn is the result of time evolution of the fluid motion starting from the initial condition following the Navier–Stokes equations. Hence, at the current time, assumes values related to viscous dissipation and kinetic energy both consequence of the fluid motion.

#### 3. Incompressible Fully Developed Turbulence

This section studies the distribution functions and in incompressible fully developed turbulence, where and are, respectively, the distribution functions of fluid particles position and Lyapunov vector and of local Lyapunov exponent. In particular, we will show that the proposed statistics leads to the following relations: and to the link between and the longitudinal velocity correlation function.

According to (1), depends on and is expressed in terms of this latter through the Frobenius–Perron equation where is the Dirac’s delta. Hence, the distribution function is first studied. This PDF changes with the time according to the Liouville theorem [22] associated with (3). This theorem, arising from the relation and from (3), provides the evolution equation of [22]where and denote the divergence of defined in the spaces and respectively and and are the elemental volumes in the corresponding spaces.

Taking into account (14) and that the homogeneous isotropic turbulence is defined for unbounded fluid domains, will satisfy the following boundary condition: Accordingly, the statistical average of an integrable function of and , say , is calculated in terms of In particular, average and maximum finite scale Lyapunov exponents are

*Remark 1. *It is very important to observe that the Liouville theorem in the form of (15) holds also when the variable velocity field is replaced with the same field frozen at the current time. This is due to the hypothesis that the times of variations of the velocity field, which changes according to the Navier–Stokes equations, are much larger than those associated with and whose times of variations are of the order of . Following such hypothesis, the statistics of and does not depend on the time variations of and can be considered to be the result of the instantaneous fluid act of motion.

##### 3.1. Distribution Function of and

For the sake of our convenience, we introduce the quantity which defines the relative kinematics, being ; thus the Liouville theorem reads as where Now, the entropy associated with is defined as Due to fluid incompressibility, ; therefore from the Liouville theorem the entropy rate identically vanishes In order to estimate the steady distribution of , observe that the fully developed turbulence corresponds to a situation of maximum chaos where the kinematic bifurcations cause a total loss of the initial condition data of . Accordingly, it is reasonable that assumes its maximum value compatible with (14) and (20) where . This corresponds to the following variational problem: in which is the Lagrangian of the problem and and are the Lagrange multipliers associated with conditions (14) and (20), respectively. The maximum of is then obtained as steady condition for () through the variational calculus and this leads to the Euler–Lagrange equation whose solutions are searched for substituting (25) into (26) where is the normalization constant related to whose value is calculated with (14), whereas is obtained through the Liouville equation. This latter gives in which denotes the dyadic product between vectors. Integrating (28), we obtain where is an arbitrary solenoidal vector field. Now, the determinant of identically vanishes; thus (29) admits solutions only for . Moreover, as exhibits minimum rank, (29) has only the trivial solution i.e., is uniformly distributed on , being in which indicates the measure of the set .

##### 3.2. Lyapunov Exponents Distribution Function

In order to estimate , the following should be noted.

If is given, (1) corresponds to a hypersurface of the space whose equation reads asWhen varies, and its points will change according to (32) where and are, respectively, local normal and tangent unit vectors to both calculated in , being not determined. On the other hand, can be determined considering that does really not depend on , whereas the variations of are related to those of by means of (5) and (6). Thus, is locally expressed as As is solenoidal, the surface integrals of over arbitrary hypersurfaces (flow of through ) assume the same value which does not depend on , i.e., Therefore, from (33), we obtain

Now, the distribution function of is calculated substituting (31) into the Frobenius–Perron equation (13) and considering that is uniformly distributed on where the surface integral over , described by , is formally calculated according to the Minkowski measure theory [23]. Taking into account (36) and that , results to be uniformly distributed in , being Hence, and are calculated in terms of with (18) and (19) Next, it is useful to calculate the mean square According to (40), the mean square of equals the square of the average of calculated for , and the standard deviation is proportional to

It is worth remarking that the obtained distribution (38) is the result of three effects, of which two are in competition with each other. The first one of these, due to the Lyapunov vectors tendency to align the maximum growth rate direction of , is responsible for variation intervals in which , producing trajectories instability. The second one, related to the fluid incompressibility, acts in opposite sense preserving the volume, determining regions where . The third element is the chaotic regime, here given by imposing , which causes a continuous distribution of .

#### 4. Analysis through Alignment Property of

One reasonable confirmation of the previous results is here given by exploiting the alignment property of , following which tends to align with the direction of the maximum growth rate of [15] and the fluid incompressibility. This leads to an alternative way to achieve (12). Accordingly, is now calculated adopting a proper PDF obtained as projection of at the time , where is considered to be constant and is the Lyapunov time defined by After the time , the alignment tendency of provides the fact that mostly all the Lyapunov vectors calculated at = const will be such that . The vectors lying in subspaces of orthogonal to the maximum rising rate direction of which do not follow such alignment form a null measure set in ; thus is a distribution function which satisfies and is expressed in function of and where and . Neglecting the higher order terms, is calculated as This PDF identically satisfies (14). In fact, the integral over of the first term at the RHS of (45) is equal to one, whereas the integral of the second one can be reduced to a proper surface integral of over , where through Green’s identity; thus this identically vanishes. Furthermore exhibits the same entropy of at least of higher order terms, which in turn does not vary with the time due to the fluid incompressibility The second addend at the RHS of (46) identically vanishes as it can be reduced to be a surface integral of over , where . Therefore At the least higher order terms, maintains the same level of information of ; thus it is adequate to estimate . Accordingly, this latter is calculated as Substituting (45) and (42) into (48), we have Integrating by parts the second addend and taking into account the fluid incompressibility and the boundary conditions (16), we obtain Hence in agreement with previous results.

#### 5. Lyapunov Exponents in Terms of Longitudinal Velocity Correlation

The results of the previous analysis allow achieving the link between and the longitudinal velocity correlation function. In fact, the standard deviation of longitudinal velocity difference directly depends on and according to On the other hand, the Lyapunov theory gives the longitudinal velocity difference in terms of Taking into account the isotropy and (40), reads as thus, and are expressed in function of the finite scale by means of the longitudinal velocity correlation Equation (55) coincides with that proposed in [1] which leads to the closure formulas of the von Kármán and Corrsin equations. Unlike [1], (55) is here achieved exploiting the shape of the distribution (38). Reference [1] shows that (55) provides a value of the skewness of equal to (see the appendix), in good agreement with the results obtained by the several authors with direct numerical simulation of the Navier–Stokes equations (DNS) [16–18] () and Large–eddy simulations (LES) [19–21] (). In detail, Table 1 recalls the comparison between the value of the skewness calculated with the proposed analysis and those obtained by the aforementioned works. As a result the maximum absolute difference between the proposed value and the other ones is less than 10%. Furthermore, other studies [24–26] have shown that the closure formulas referable to this analysis provide kinetic energy and temperature spectra which exhibit scaling laws in agreement with the theoretical arguments of Kolmogorov, Obukov–Corrsin, and Batchelor [27–29]. Therefore, the adopted hypothesis and the consequent distribution (38) seem to be adequate assumptions at least for what concerns the estimation of turbulent energy cascade main effects.