Next Article in Journal
As Regular as Clockwork: Alexander von Humboldt, Robert de Lamanon and the Beginning of the Scientific Investigation of the Tidal Barometric Oscillation
Previous Article in Journal
Numerical Simulation of the Effects of Surface Roughness on Light Scattering by Hexagonal Ice Plates
Previous Article in Special Issue
Theory and Modelling of Isotropic Turbulence: From Incompressible through Increasingly Compressible Flows
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Chaotic Measures as an Alternative to Spectral Measures for Analysing Turbulent Flow

1
Njord Centre, Physics Department, University of Oslo, 0316 Oslo, Norway
2
School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3FD, UK
3
Thornton Tomasetti, 2 St Davids Drive, Dunfermline KY11 9PF, UK
*
Author to whom correspondence should be addressed.
Atmosphere 2024, 15(9), 1053; https://doi.org/10.3390/atmos15091053
Submission received: 20 June 2024 / Revised: 31 July 2024 / Accepted: 27 August 2024 / Published: 30 August 2024
(This article belongs to the Special Issue Isotropic Turbulence: Recent Advances and Current Challenges)

Abstract

:
Turbulence has associated chaotic features. In the past couple of decades, there has been growing interest in the study of these features as an alternative means of understanding turbulent systems. Our own input to this effort is in contributing to the initial studies of chaos in Eulerian flow using direct numerical simulation (DNS). In this review, we discuss the progress achieved in the turbulence community in understanding chaotic measures including our own work. A central relation between turbulence and chaos is one by Ruelle that connects the maximum Lyapunov exponent and the Reynolds number. The first DNS studies, ours amongst them, in obtaining this relation have shown the viability of chaotic simulation studies of Eulerian flow. Such chaotic measures and associated simulation methodology provides an alternative means to probe turbulent flow. Building on this, we analyze the finite-time Lyapunov exponent (FTLE) and study its fluctuations; we find that chaotic measures could be quantified accurately even at small simulation box sizes where for comparative sizes spectral measures would be inconclusive. We further highlight applications of chaotic measures in analyzing phase transition behavior in turbulent flow and two-dimensional thin-layer turbulent systems. This work shows that chaotic measures are an excellent tool that can be used alongside spectral measures in studying turbulent flow.

1. Introduction

Turbulence is a descriptor for fluid flow with highly complex structure, where velocity varies rapidly with space and time [1]. Although the flow is seemingly random, its statistical quantities and their dependence on length scales are reproducible. These statistics strongly depend on the boundaries and other macroscopic restrictions. Along with global statistics such as dissipation and turbulent energy, structure functions can characterize turbulent flows. For homogeneous isotropic turbulence (HIT), the structure functions can be related to the spectral quantities, which are statistics that depend on the wave number [2]. These are relevant in statistical closure approximations, where higher-order structure functions are related to lower-order ones [3,4].
However, chaotic measures are another set of statistics which are comparatively less well established in the turbulence literature even though there are several advantages of using them over spectral measures. We believe that the analysis of chaotic properties of turbulence presents an untapped potential for a better understanding of turbulence. Recently, new advances have been made in quantifying chaotic measures in turbulence, and this review provides context and framework for these studies, tracks the evolution of the main ideas, as well as summarizes their key implications.
In discussing the advantages of chaotic measures, we very briefly discuss the disadvantages of spectral measures in the context of direct numerical simulation (DNS). Within the K41 theory, there is the large length scale where energy is inputted, and then the dissipation at small-length scale occurs. These scales are decoupled, and between them an inertial range should exist where energy transfer occurs from larger to smaller scales. Prediction of the wave number scaling of this range as well as how the longitudinal structure functions relate to the dissipation can then be made [5,6,7]. Intermittency effects can be said to account for the degree that these predictions agree or do not agree with simulations and experiment [8]. However, with DNS, it is often hard to generate a well-defined inertial range. To achieve this, we need a large separation of scales, which requires a large simulation and thus large computational resources. Unlike the scaling of any inertial range, many chaotic quantities are more like the turbulent energy or the local dissipation rate—they exist regardless of the existence of a well-defined inertial range. They thus can be calculated with much more modest simulation sizes. Laws for scaling such as S n ( r ) r ζ n for structure functions, where ζ n is the scaling exponent, require very large Reynolds numbers and a well-defined inertial range as shown by Qian [9,10,11]. The fact that these conditions are not required for the chaotic measures adds a strong argument for their use. It is notable that even in 1986, as shown by Deissler [12], robust chaos was observed where an inertial range would definitely not have been found.

Advantages of Chaotic Measures Applied to Turbulence

Chaos captures the idea that a small initial difference between two systems can grow exponentially in time. We can study an evolving fluid in two ways: we may either imagine following one fluid parcel in the flow and see how its position evolves in time, which is the Lagrangian approach, or we may also imagine standing in one position and seeing what the fluid does at that position, which is the Eulerian approach. These approaches are inextricably linked. If we knew the precise details of the flow at every point at all times, we could determine how each fluid parcel behaves. Similarly, if we knew the behavior of each parcel, we could recreate the evolution at each position.
Turbulence displays chaotic dynamics [13] and ideas from chaos theory find many different applications in turbulence including the dispersion of pairs of particles [14,15,16,17], the presence of Lagrangian coherent structures [18], turbulent mixing [19], turbulent transitions [20], and predictability [21,22,23,24,25,26]. Whilst the Lagrangian aspects of chaos in turbulence have been heavily studied, the Eulerian picture has been comparatively neglected up until the past few years when advances in computing power and theory have allowed novel insights to be found in long-pondered questions. These recent developments are the main focus of this review, including the contributions we make to the overall effort.
One of the reasons we wish to apply chaotic measures to turbulence is in the name of predictability. We may believe that given more and more accurate information of a system and its governing equations, we are able to predict the future of that system with more and more precision. However, chaos means that a small difference from the true evolution diverges exponentially. In many short introductions of turbulence, the fact that turbulence is chaotic is mentioned and then not further explored. Sometimes a calculation is made of the predictability time, whereby, because the chaos is very high, the predictability is declared to be very low. This, however, goes against our experience of the atmosphere, which is highly turbulent and chaotic, yet weather prediction is accurate for a non-trivial time frame.
Another reason to apply chaotic measures is to understand the nature of turbulence itself. Different hypotheses of the nature of turbulence lead to different predictions of the implications for chaos in turbulent dynamics. These implications can be tested and hypotheses strengthened or discarded. Beyond theoretical considerations, there are practical ones as well. Turbulence is frequently studied computationally. In this review, one of the themes we highlight is that chaotic measures can be used to quantify the resolution of simulations. In structured flows, this allows us determination of where a simulation needs more resources, thus allowing more efficient computation. Whilst this review is mainly concerned with chaotic measures in HIT, which is an idealized case of primarily theoretical interest, it acts as a testbed to develop the method for more practical use. For example, the application of chaotic measures could perhaps be generalized to wave turbulence, which describes a broad category of weakly nonlinear dispersive waves [27], including gravity–capillary wave turbulence [28,29,30], magneto-hydrodynamic free-surface wave turbulence [31,32], and acoustic turbulence [33,34]. Researchers in this field may find the ideas presented in this review relevant to their own field of study.

2. Chaotic Measures of Turbulence

2.1. Lyapunov Exponents

The nature of the chaos in a system can be quantified by its Lyapunov spectrum. We can define a dynamical system as an N-dimensional system of first-order differential equations specifying the time evolution given by trajectory x ( t ) of initial condition x ( 0 ) . For an explicit example, we consider
d x ( t ) d t = F [ x ( t ) ] ,
such as Navier–Stokes equations. For an infinitesimally small perturbation to x ( 0 ) given by δ ( 0 ) , at later times we have
x 0 ( t ) = x ( t ) + δ ( t )
where δ ( t ) denotes the separation between the two trajectories. The maximum Lyapunov exponent, λ 1 , has the greatest influence on the predictability of the system. It offers a relation for the separation between trajectories δ ( t ) for a finite time interval
| δ ( t ) | | δ ( 0 ) | e λ 1 t .
Therefore, if λ 1 is positive, we have an exponential divergence between x ( t ) and x 0 ( t ) , and the system is said to exhibit deterministic chaos. This rapid divergence between two initially close states represents an extreme dependence on initial conditions.
As long as δ ( t ) = x 0 ( t ) x ( t ) is small, it can be considered a tangent vector. We can then find the time evolution of δ ( t ) from the linearized equation
d δ i ( t ) d t = j = 1 d F i x j | x ( t ) δ j ( t ) .
Then there exists a basis e i in this tangent space such that for large t [35],
δ ( t ) = i = 1 N c i e λ i t e i .
The λ i are known as the Lyapunov exponents of the system and are ordered by convention as λ 1 > λ 2 > > λ N . This collection of exponents is called the Lyapunov spectrum of the system [36]. In a conservative system, the sum of Lyapunov exponents must be zero [37]. As some trajectories diverge exponentially (they have positive Lyapunov exponent), others converge. A dissipative system with at least one λ i > 0 is said to be chaotic. We look at this spectrum further in the context of the Kolmogorov–Sinai entropy below.
Under strict mathematical requirements, the exact Lyapunov spectrum requires an infinite amount of time to calculate, which is infeasible. For estimating the maximum Lyapunov exponent, one can adopt the naive method, where a random small perturbation is made, and the divergence tracked, with a fit made to the subsequent exponential growth. However, this takes a long time for each Lyapunov exponent calculated and could introduce biases depending on the form of the perturbation. As such, different generalizations have been made to calculate the maximum Lyapunov exponent. Two commonly used methods are the finite-size and finite-time Lyapunov exponents, denoted FSLE and FTLE, respectively.
For a perturbation of an initial size, δ , one may wish to know, on average, how long it takes to grow to a certain tolerance, Δ . This time is denoted T ( δ , Δ ) . We can then define the finite-size Lyapunov exponent (FSLE) as [23]
λ ( δ , Δ ) = 1 T ( δ , Δ ) ln Δ δ ,
where the angled brackets represent an average and which recovers the usual definition of the maximum Lyapunov exponent when δ and Δ are infinitesimal. If disturbances grow at a rate equal to the local eddy turnover time, then λ ( δ ) δ 2 [22,38].
We can change the definition of the FSLE in Equation (1) so that, instead of fixing an initial tolerance Δ and then measuring T ( δ , Δ ) , we fix a time T and measure the error Δ ( δ , T ) after that time. This defines the finite-time Lyapunov exponent (FTLE) as
λ ( δ , T ) = 1 T ln Δ ( δ , T ) δ ,
where we do not currently worry about whether this averaging procedure takes place inside or outside the fraction or the natural logarithm. We can also define a definition for λ which is spatially dependent for chaos in the Lagrangian sense. These can be used to characterize the presence of Lagrangian coherent structures (LCSs) [18]. These are surfaces on which the infinitesimal deformation of the flow is maximized, represented by extremizing surfaces of λ .
In the Eulerian sense, we can implement the FTLE with the following procedure: we take the original field u 1 and create another field u 2 = u 1 + δ u , with some randomisation procedure. Then, we evolve the system over a finite-time T, measure the field difference, and then introduce a successive perturbation, starting the process again. This subsequent perturbation is given by
u 2 ( x , T ) = u 1 ( x , T ) + δ Δ ( δ , T ) δ u ( x , T ) .
The maximum Lyapunov exponent, which fluctuates in time, λ ˜ , can then be calculated as
λ ˜ = 1 T ln Δ δ ,
and this procedure is repeated, obtaining a set of values for λ ˜ .
One advantage of the FTLE method is that it produces a sample of maximum Lyapunov exponents that allows for us to perform a statistical analysis of the measurements. Compared to the naive method of direct measurement of | δ ( t ) | , using the FTLE method, it is possible to obtain a large number of maximum Lyapunov exponents which can be used to analyze the distribution for a given steady state. This procedure shows that time fluctuations in the maximum Lyapunov exponent is quite significant, which is difficult to show using the direct method. However, we also lose the ability to directly control the form of the perturbation, which can provide us with important physical information about turbulence.

2.2. KS Entropy

The unpredictability of deterministically chaotic systems can also be viewed from an information theoretical viewpoint. To this end, the concept of entropy can be extended to such systems through the Kolmogorov–Sinai (KS) entropy which provides a measure of the information production rate and is an invariant property of the system.
To define the KS entropy, which is a generalization of the classical Shannon entropy to dynamical systems, we begin by performing a discretization of the d-dimensional state space of the system, i.e., a space where each point represents a unique state of the system. After discretization, the space is composed of cells of side length l and thus volume l d . Our space is equipped with an invariant probability density function, ρ ( x ) ; hence, we can write the probability of the system being in a given cell as
p i = P i ρ ( x ) , d V i .
These probabilities could be used to define a Shannon entropy for the system. However, as it stands, the above definition is not independent of how the space is partitioned, which is required of a system invariant such as the KS entropy [13,39]. Moreover, the KS entropy provides the rate of information production; thu, its definition should include time dependence.
We now extend our probabilistic definition to consider trajectories through state space. By looking at these trajectories, we are using the Eulerian framework for chaos. As we have a measure of the probability that the system is located in a given cell, we can define the probability of a trajectory comprising the following: being in cell i 1 at t = τ , then i 2 at t = τ + δ t , and so on until at t = τ + n δ t , the system is in cell i n , which we denote as p i 1 i 2 i n . Then, still motivated by the Shannon entropy, we define the entropy of such a trajectory as
H n = i 1 , , i n p i 1 i 2 i n log 2 p i 1 i 2 i n .
The difference in entropy between a trajectory at step n and the same trajectory at step n + 1 provides a measure of information production. Then, by taking appropriate limits to remove dependence on trajectory length, sampling frequency and how we choose to discretize space, the KS entropy, h K S , is defined [40]:
h K S = lim δ t 0 lim l 0 lim N 1 N τ n = 0 N 1 ( H n + 1 H n ) .
Systems which exhibit deterministic chaos have positive, but finite, KS entropy whilst purely random systems have infinite entropy.
A more practical definition of the KS entropy can be found by making a connection with the Lyapunov spectrum of the system through Pesin’s theorem,
h K S = λ i > 0 λ i .
Considering once more a discretized state space, if two trajectories are initially located inside the same cell, then they are effectively identical. If the system exhibits deterministic chaos, then these trajectories diverge over time, with the rate of divergence along a given eigenvector of the system given by the corresponding Lyapunov exponent. Eventually, this divergence leads to the two trajectories occupying different cells and thus becoming distinguishable, which can be viewed as producing information. The above definition of the KS entropy in terms of the sum of positive Lyapunov exponents, i.e., the cumulative divergence rate of the trajectories, then makes clear the link between chaos and information theory.
In contrast with conservative systems, where Liouville’s theorem applies, for dissipative systems such as fluid turbulence, phase space volumes are not constant along trajectories and are instead found to contract as they move along the trajectory in phase space. However, due to the presence of a subset of positive Lyapunov exponents, there can exist sub-volumes, i.e., spaces of lower dimensionality than the overall state space, which do preserve their volume along trajectories. This leads to the idea of an attractor, which is a subset of the state space to which the system is said to evolve towards. The dimension of this attractor can also be related to the Lyapunov spectrum via the Kaplan–Yorke conjecture [41,42] as
dim ( A ) = j + i = 1 j λ i | λ j + 1 | ,
where j is determined by
i = 1 j λ i > 0 and i = 1 j + 1 λ i < 0 .
For chaotic systems, dim ( A ) is expected to be non-integer and the resulting attractor is referred to as a strange attractor.

3. Theoretical Predictions

Multiple theoretical predictions have been made about the Lyapunov exponents and other chaotic measures in turbulence and specifically how they relate to other measures of turbulence. We review some of these below.

3.1. Theory of Lorenz

Since some of the earliest observations of chaos were found in atmospheric models by Lorenz [21], and as atmospheric dynamics are strongly influenced by turbulence, it is only fitting that predictions about the level of predictability of turbulence immediately arose from these papers.
Lorenz considered a fluid with energy spectra
E ( k ) = A k n .
At each scale k, he posited a local time scale τ ( k ) , which is proportional to the turnover time of those eddies which have a wave length 2 π / k . On dimensional grounds, this time scale is [22]
τ ( k ) = A 1 / 2 k ( n 3 ) / 2 ,
and Lorenz then suggests that this is the time it takes for a perturbation to grow from perturbation scale k to 2 k . In sum, this implies [13]
τ ( k ) 1 k v k , v k 2 k 2 k E ( k ) d k ,
with v k being the typical velocity difference. If we consider that τ ( k ) is minimized when k v k is maximized, which is also when k 2 v k 2 is maximized, the τ ( k ) should be dominated by the peak of E ( k ) k 3 . This is different to the theory proposed by Ruelle below that the error is dominated by the Kolmogorov scale eddies.
Because this was derived within the context of atmospheric prediction, it was immediately recognized that the degrees of freedom in the system could propagate errors in a way that was more important than the divergence of close trajectories. Specifically, there could be different timescales that would grow the error at different rates and interact non-linearly. Because of these, it was expected that the predictability of the entire system did not have the maximum Lyapunov exponent as the key determiner of predictability. Indeed, it was found that the largest length scales had predictability times independent of the Reynolds number [13].

3.2. Theory of Ruelle

An alternative approach is that of Ruelle [43]. He defines a characteristic exponent for the exponential growth of a disturbance Δ
λ = lim t 1 t ln Δ ( t ) ln Δ ( 0 ) .
He then reasons that this should depend on eddy size l. Using dimensional considerations, he derives
λ ( l ) ε 1 / 3 l 2 / 3 ,
where ε is the dissipation rate. He then states that the largest characteristic exponent should be associated with the smallest eddies, which are the Kolmogorov scale eddies η ( ν 3 / ε ) 1 / 4 , where ν is viscosity. Overall, this implies
λ ( l ) ε 1 / 3 η 2 / 3 ε 1 / 2 ν 1 / 2 1 τ ,
which in essence implies that the chaos is dominated by the smallest time scale, which is the Kolmogorov time scale. This theory is fundamentally based on the Eulerian sense of chaos, but has often been used in a Lagrangian sense. Comment is made upon this discrepancy in [44].

3.3. Multifractal Corrections

The above Ruelle theory implicitly relies on the K41 theory. An attempted correction is made in [45] to account for intermittency effects. Using the Kolmogorov theory, we have the prediction that
η L Re 3 / 4 ,
where L is the integral length scale, and this implies
λ 1 T 0 Re 1 / 2 .
However, the turnover time of an eddy of size l should be
τ ( l ) v ( l ) l T 0 l 1 h ,
where v ( l ) is the velocity of the eddy and h is the Hölder exponent of the velocity difference
δ v ( l ) = | v ( x + r ) v ( x ) | V l h ,
and l = r / L is some scaling parameter. This implies that
η L Re 1 / ( 1 + h ) ,
and so,
λ 1 τ 1 T 0 R e α , α = 1 h 1 + h .
Prediction α = 1 / 2 is then a consequence of the Kolmogorov theory. If ζ n is the scaling exponent for structure function S n ( r ) r ζ n , then in the Kolmogorov theory, ζ n = n / 3 ; hence, h = 1 / 3 , and so α is predicted to be 1 / 2 . Differences from this value can then come from intermittency effects, where it is shown with a shell model of turbulence that α < 1 / 2 . However, in [44], Fouxon et al. state that this derivation is wrong, since it relies on a questionable averaging procedure. In summary, the derivation relies on the regions with reduced chaos becoming larger due to intermittency, but the chaos is actually dominated by the remaining regions which have increased chaos and contribute more to λ . Such structural behavior is found in [46].

3.4. Predictions for Attractor Dimension and KS Entropy

Within the literature, there exist a number of complex and mathematically sophisticated predictions for the scaling behavior of the KS entropy and the attractor dimension for isotropic turbulence. However, simplified expressions can be obtained by considering the K41 theory and thus neglecting potential intermittency corrections.
Considering first the attractor dimension, a dimensional argument was provided by Landau and Lifshitz [47] which estimated the degrees of freedom of a turbulent flow as the ratio of the largest to the smallest length scales. This is expressed as
dim ( A ) L η 3 R e 9 / 4 .
In [48], Ruelle derived a more general result which allowed for a spatially intermittent dissipation rate; however, this requires detailed knowledge of the spatial probability distribution for the dissipation rate which is to date unknown. Instead, if K41 is assumed, this result is consistent with the simple dimensional arguments. Further estimates of scaling of the attractor dimension can be found in [49,50]. Finally, the multi-fractal model can also be used to derive a scaling prediction with intermittency corrections [51]. Here, the scaling is still in terms of the Reynolds number; however, the exponent is lowered compared to the K41 case.
For KS entropy, fewer predictions have been made for isotropic turbulence. Ruelle [48] again provided a result that can account for spatial intermittency, provided statistical information about the dissipation rate. Once more, if K41 is assumed, then this result can be simplified and predicts
h K S 1 τ η L η 3 1 T R e 11 / 4 .
Interestingly, this result is simply the product of the scaling predictions for the maximum Lyapunov exponent and the attractor dimension.

3.5. Two Dimensional Turbulence

In the above sections, we focused on theoretical predictions for three-dimensional turbulence; however, the same ideas can also be applied to the two-dimensional case. The inverse cascade of energy in two-dimensional turbulence [52,53,54] has been found to lead to the formation of long-lived and predictable coherent structures which has profound effects for predictability.
For the maximum Lyapunov exponent, just as in three dimensions, we follow the dimensional argument put forward by Ruelle that suggests it should be inversely proportional to the fastest timescale of the flow. We begin by assuming that in the direct cascade of enstrophy, the energy spectrum takes the form E ( k ) k 3 , ignoring, at this stage, any logarithmic corrections. Importantly, this form suggests that there is only a single timescale, τ , associated with the entire forward cascade which is set by the enstrophy dissipation rate, η , such that τ η 1 / 3 . This then implies
λ 1 1 τ η 1 / 3 .
The enstrophy dissipation rate is controlled by the enstrophy injection scale at the top of the cascade. Interestingly, in contrast with three-dimensional turbulence, this suggests that in two-dimensional turbulence the smallest scales of the flow contain information about the large length scales that generated them.
Applying similar dimensional arguments to the KS entropy, we can also estimate its scaling behavior. In line with the three-dimensional case, we expect the entropy scaling to be given as the product of the inverse of the fastest timescale and the total number of excited modes. The former we know from the maximal exponent and the latter offers the scaling of the attractor dimension. We consider the ratio of the largest and smallest scales of the flow and use this as an estimate of the scaling of the attractor dimension. In doing so, we require the dissipation length scale, χ = ( ν 3 η ) 1 / 6 , where ν is viscosity. The ratio then implies
dim ( A ) L χ 2 Re ,
and hence the KS entropy should scale as
h K S 1 τ Re = η 1 / 3 Re .
The above picture is incomplete since it ignores the logarithmic corrections to the energy spectrum that were suggested by Kraichnan [55]. These corrections are necessary for a constant enstrophy flux in the direct cascade. The influence of this logarithmic correction on the chaotic properties of two-dimensional turbulence were investigated by Ohkitani [56] and modify the above scaling predictions,
λ ( η log ( Re ) ) 1 / 3 ,
dim ( A ) Re ( log ( Re ) ) 1 / 3 ,
and
h K S η 1 / 3 Re ( log ( Re ) ) 2 / 3 .
Ruelle’s predictions in [48] were modified for the two-dimensional case in [57] by Lieb. Whilst it is suggested that during the inverse cascade there is intermittency in the energy dissipation rate, in the absence of this effect, the scaling for entropy is
h K S ε ν 2 V ,
where V is the volume of the system, whereas for the attractor dimension it was found that
dim ( A ) ε ν 3 V .
These are consistent with the picture suggested by the dimensional argument in which the chaotic properties are influenced by all scales in the system.

4. Results from Simulations

The above theory has been tested using numerical simulations. In discussing the results of simulations, we split the discussion into various categories. However, they interlink, and it is impossible to create a totally clear division between them. The main ideas that come up again and again are the following: how the chaos depends on the turbulent statistics, the different length scales that exist in the way that chaos presents in turbulence, and how this can allow for more efficient simulation.

4.1. Assessing Theoretical Predictions

4.1.1. Ruelle’s Prediction

Multiple simulations into the Eulerian characteristics of chaos in turbulence, including those by us, have shown that Ruelle’s prediction of constant λ τ does not hold, and instead that it increases with Re [38,46,58]. These simulations reached as high as 2048 3 collocation points [58] by us and Re = 8224 [38] by Boffetta and Musacchio. These three papers had slightly different emphases. In [46], Mohan et al. explicitly measured the dependence of λ τ on Re λ with dependence λ τ Re λ β . This exponent was estimated with Baesian inference finding that the most likely value had 1 / 4 < β < 1 / 3 , “with the possibility that the value is zero essentially precluded”. In our work [58], the relation of λ T 0 to Re was measured and shown in Figure 1, where we found λ T 0 C Re α , and α = 0.53 ± 0.03 . This fit was performed to directly compare it to previous simulations involving shell models and predictions of intermittency using the Hölder exponent, which predicted that the value for α would decrease below 1/2 [59]. In Boffetta and Musacchio [38], the relation of λ T E 0 to Re was measured, where T E 0 = E / ε , where it was found that λ T E 0 Re β , with β = 0.64 ± 0.05 . They found that λ τ being constant is a poor assumption. This exponent was significantly different to that in [58], but this discrepancy was resolved in [60], which concluded that the difference between the value of 0.64 and 0.53 could be accounted for by the different timescales used, namely E / ε for β and L / U for α .
Interestingly, an even earlier paper by Mukherjee et al. looking at predictability in dry convective boundary layers using large eddy simulation (LES) found a scaling equivalent to the prediction of λ Re 1 / 2 with an exponent possibly even greater than 1 / 2 [61]. These simulations were achieved using a simplified LES subgrid, which is similar to DNS. However, the direct Ruelle relation of λ 1 / τ was not directly measured. These simulations also showed a growth rate independent of perturbation length scale, and a self-similar spectral shape, also discussed in depth below.
Whilst the simulations looking directly at the Ruelle prediction all showed non-constant λ τ , in [62] Budanur and Kantz used simulations with very high resolution (meaning that k m a x η was much larger than unity) and stated that the reason for an increase in λ τ with Re could be due to the unresolved scales in the previous simulations, since they found no increase in λ τ with Re. It was also stated that the artificial forcing used in these simulations could have resulted in non-physical Lyapunov scaling. However, explicit tests by Mukherjee et al. using identical Reynolds numbers found no dependence on resolution in error growth rates [61].
Since the inverse of the maximum Lyapunov exponent can be associated with a timescale, a continuing increase in λ τ with Re eventually leads to λ τ > 1 , which implies that there are timescales in the system faster than the Kolmogorov microscale time. Whilst it would be interesting to discover what happens at such scales in DNS, the observed very slow scaling of λ τ with Re means that this would require Re on the order of 10 9 , which is outside the realm of possibility for simulation in the foreseeable future. If such a scaling did hold, it might have implied that the instabilities act at even smaller scales than the Kolmogorov microscale, which is below the scale thought to be the smallest relevant scale in turbulence. Sub-Kolmogorov scales would then be of importance, and this should be tested as suggested in [62] by Budanur and Kantz. Although we believe the simulations performed in [38,46,58] are robust, this does not preclude the fact that new and interesting results can be found by testing this relation further. Indeed, it has been suggested that thermal noises, which happen at smaller scales than the Kolmogorov microscale, may strongly affect predictability [63]. Interestingly, in thermalized fluids, a link between the maximum Lyapunov exponents and the temperature has already been seen in many body chaos, suggesting an underlying universality [64].
When observing the fit to data, the fit of λ T E 0 to Re α is excellent [38,60], with scaling that does not vary greatly over the entire parameter range. This begs the question: Perhaps this relation to Re is the real one? If it were, this would imply that λ τ R e α 1 / 2 C ( R e ) , where C ( R e ) = ε L / U 3 is the dimensionless dissipation rate. In principle, at low Re, this may depend on Re [65]. From the numbers found in [38], and since Re λ Re 1 / 2 , this would imply that λ τ Re λ 0.28 , which is within the range calculated in [46], despite these being independently measured. The argument that τ is the relevant timescale of chaos is based on a simple physical argument, namely that chaos should be dominated by the smallest timescale of the system [43], although care needs to be taken in this interpretation, as noted by Fouxon et al. in [44]. But we know that chaos is a multiscale phenomenon, and it stands to reason that chaos in turbulence may also be a multiscale phenomenon. This is supported by the fact that the predictability of the large length scales acts fundamentally differently to that of the smallest scales. As such, it could be thought unlikely that the level of chaos in turbulence depends only on a small length scale quantity.
If the chaos is not totally derived from the smallest turbulence scale, the question is raised about where it does come from. A possible resolution to this question is provided by Mohan et al. in [46] itself, where it is discussed what the spatial contributions to the error were. The discussion there is particularly enlightening and should be read in full by the interested reader. Here, it was found, by looking at times where the field displayed extreme values of FTLE, that the difference between the two fields, i.e., where the error was located spatially, was dominated by local disturbances that were long-lived. In effect, the error growth can be dominated globally by large coherent structures in the error itself. The average growth is dominated by these extreme values and so they are particularly important. Within this region of high error, “the base field exhibits a pair of corotating vortex tubes” [46], which suggests that the growth of error is fundamentally linked to vortices in the flow, and is thus affected by the rate of vortex stretching. This idea is further explored later in the review. This effect of error growth being dominated by the vortices has also been seen in two-dimensional decaying turbulence [66,67].
The interpretation of λ τ increasing with Re is that it might be more realistic to expect that, as the flow becomes more extreme, there are more extreme regions that affect the maximum Lyapunov exponent, and since the extremes are non-Gaussian, and this extremity increases with Re, then the value should increase with the Reynolds number. This does not require there to be timescales that act, on average, faster than the Kolmogorov microscale time.
Because of the impossibility of directly testing Re where we expect that λ τ > 1 , we need to turn to theory, which can also provide exciting new directions for numerical studies of turbulence. An example of such an intriguing theoretical contribution is found in [44]. Fouxon et al. note that the Ruelle relation is reliant on a lack of intermittency, whereas the Kolmogorov microscale time can change depending on location. Applied to fluid particles, this implies that there is a dependence on the location for fluid particles and so the two relations are not necessarily the same. Furthermore, the Ruelle relation was derived for turbulence, but has seen successful application to fluid particles embedded in turbulence. The authors show how one prediction can apply to the other. The statistics of a field and the particles in that field are not necessarily the same and are dependent on how statistics are calculated, and thus there needs to be a careful consideration of statistics, which can prove fruitful theoretically and in simulations. This would have implications on any applicability of λ τ Re α . These are promising avenues to pursuit. The authors also clarify the nature of intermittency in turbulence, stating that there are in fact at least two sources. The first is due to random fluctuations in the field that exist even with Gaussian fields. The second is due to the intermittent nature of turbulence, which they define as “atypically large persistent velocity gradients” [44].

4.1.2. KS Entropy and the Lyapunov Spectrum

Extending beyond the maximum Lyapunov exponent, a robust use of chaotic measures requires computing the spectrum of Lyapunov exponents, which is an alternative complete description of turbulence in distinction to spectral measures. This is a numerically intensive problem. The first attempt at such simulations was performed by us [68] and independently in [69,70] by Hassanaly and Raman.
In [68], we directly computed the Lyapunov spectrum for three-dimensional HIT using DNS and reproduce it here in Figure 2. This allowed numerically rigorous tests of the mathematical relations and conjectures in the literature. Since a very large spectrum was obtained, it allowed calculation of the KS entropy and thus the rate of information production. The dimension is shown to be very large, even at low Re, being as high as 5704 for Re = 212. Thus, low-dimensional approximations offer poor estimates. The results allow a direct comparison to the results of Ruelle, intermittency, and finite Reynolds number effects, given the model independent results. The calculations are achieved using the Benettin algorithm, which is very costly, involving M FTLEs, one for each calculated exponent, and a Gram–Schmidt algorithm to orthogonalize the measurements such that a well-ordered spectrum can be obtained, which means that Re needs to be modest, ranging from Re 50 to 212, and M has to scale very quickly to capture all positive exponents. We find that h K S T = C Re α with C = 0.0008 ± 0.0002 and α = 2.65 ± 0.06 , with commentary on agreement with predictions and problems with using shell models. However, the low values of Re tested place a large uncertainty in the scaling exponent. The shape of the spectrum is then compared to theory, and here, a divergence density of exponents around λ 0 predicted by Ruelle [48] is not observed. An Re-independent shape for the spectra is found consistent with what was seen in Poiseuille flow by Keefe et al. [71].
Finally, an estimate of the attractor dimension for a given Re is made. For the lowest Re cases, this is calculated directly using the Kaplan–Yorke conjecture, whereas for the higher Re cases, this relies on the shape similarity property of the spectrum to enable a reliable prediction of the crossing point for the sum of exponents. The dimension is extremely high. Landau and Lifshitz offer a good approximation, with others compared for the number of positive exponents, finding N = b Re γ with b = 0.008 ± 0.002 and γ = 2.35 ± 0.05 .
We performed a similar set of simulations for two-dimensional turbulence in [72], which extended the results in a systematic way, pointing out its greater applicability to thin layer turbulence such as in the atmosphere. There are different timescales in two and three dimensions, which results in different Re scaling. The simulations use large length scale dissipation to prevent the formation of a condensate at low k and a stochastic forcing. We found very good agreement for λ 1 = 0.42 η 1 / 3 , which agrees with prediction, although the fit was performed by varying the prefactor since varying the exponent in this case would not be dimensionally consistent. We then showed that some correction to λ τ is required, resulting in some Re dependence, as suggested by Kraichnan.
In two dimensions, it was found that in addition, there is an influence on the forcing length scale and system size, suggesting non-universal behavior in addition to what was predicted from the Kraichnan theory. The attractor dimension grew with the inverse cascade inertial range, but had a contribution from enstrophy. As in three dimensions, no divergence around λ = 0 in the spectrum was seen. It was further found that dimensional arguments were not sufficient to account for predictions, with corrections based on forcing and length scale, due to the non-local effects of coherent vorticies. This again shows how important vortices are for error growth.
A further study was performed in [69]. Hassanaly and Raman note how a dynamical systems understanding of turbulence can act as a complement to the statistical view of averages and deviations from the average. Whilst there are statistical quantities to capture these effects, their influence on the dynamics of the system can be illuminated with a dynamical systems approach. In the study, they calculate the Lyapunov spectrum for three-dimensional HIT using low-mach-number incompressible flows with constant density. They calculate the Kaplan–Yorke dimension as D K Y = ( L / η ) 2.8 ± 0.095 . From their plots of the Lyapunov vectors, the authors conclude that the vectors “themselves are disorganized fields, containing structures similar to the original velocity field”. After performing a statistical analysis between the velocity field and the LV field, they find results similar to [46], where the contribution to the error growth is contained in very small areas, furthermore finding that this fraction decreases with increasing Re. Hassanaly and Raman find “a strong correlation between the chaotic parts of the domain and the velocity gradients at low Reynolds number. At higher Reynolds numbers, enstrophy becomes a good marker of chaoticity”. Lyapunov exponent and Lyapunov vector analysis has also been applied in [70] by Hassanaly and Raman, and to magnetohydrodynamic (MHD) turbulence on a Shannon information entropy level by Vasey et al. in [73].
The full spectrum of Lyapunov exponents would be an alternative complete description of the fluid flow. Our studies showing the stability of such quantities for the maximal exponent, we expect, also hold for the other exponents. Thus, computing these exponents may prove to have considerable numerical advantages. We performed an initial examination of such a calculation in [68], and it seems beneficial to develop such studies in the future.

4.2. Predictability

4.2.1. Error Spectrum

By looking directly at the error spectrum, we can see the spatial distribution of the error in a similar way to looking at it visually. By following the direct approach, we can see the different scale behavior as opposed to FTLE [58]. It is harder to see the different scale behaviors in FTLEs since those saturate at exponential growth regimes, the spectral information is lost, and since different scales are important for different problems. The “naive” approach has validity, especially since we can directly probe different scales (see the eddy damped quasi-normal Markovian (EDQNM) approximation for an important example). We looked at the difference spectrum in DNS in [58], reproduced here in Figure 3. There is a power law dependence in the small wavenunmber region where the exponent grows as three or four, and then drops off exponentially beyond the peak, which occurs at roughly k η , the wavenumber of the Taylor microscale. As Re increases, this exponent becomes increasingly flat (a white spectrum) and indicates that below that scale, the noise is distributed at all scales equally. Further exploration shows that using the error spectrum, we see dynamical changes that depend on the saturation of the error spectrum. The transition from exponential to linear growth happens when the error saturation occurs at k η [38], with the FSLEs also offering important structural information, finding “the time for the perturbation to affect a wave number k in the inertial range is proportional to ε 1 / 3 k 2 / 3 ”.
In [74], Yoshimatsu and Ariki analyze an initially high-k error in three-dimensional turbulence and find similar results. They use self-similarity arguments to find an integral length scale of the error field and show how the normalized error spectra collapse, and in the low wavenumber range find a k 4 dependence. They use a spectral cutoff filter to split the error production into high and low wavenumbers dependent on a critical wavenumber, finding that the interaction of the two contributes to the error scaling at low wavenumber, whereas the growth of the absolute value of the field is mainly provided by the low wavenumbers. This goes against the common idea that the error growth is driven mainly by the smallest scale features which propagate upwards. Instead, we have the large length scale features being important as well. One of the advantages of DNS is being able to probe such realizations and understand turbulence better, and chaotic measures are vital for this understanding.

4.2.2. Cascades of Error

In [22], Lorenz was primarily concerned with the way that an error in the small length scale could propagate upwards in scale, so that a lack of knowledge at small length scales would imply limits on predictions of the large length scales. However, it has been pointed out that when applying the analysis to a turbulent fluid, there are drawbacks since the multiscale error interacts not just in an inverse cascade, but forwards as well. This is partially analogous to the situation of energy cascades in turbulence.
In Fourier space, Waleffe analyzed turbulence by its fundamental triadic interactions. Under the assumption of statistical homogeneity, the nonlinearity of the Navier–Stokes equations can be understood in the Fourier space as the interactions of triads of wavevectors k , p , q , such that k + p + q = 0 . Formally, each such triad of wavevectors defines a dynamical system describing the interaction of the corresponding Fourier modes, which makes it possible to investigate the behavior of different interactions between helical Fourier modes. Most attention was historically spent on analyzing the forward cascade of energy; however, certain triadic interactions of turbulence actually sum to an inverse cascade of energy [75], as has been shown in a numerical simulation [76]. Such triadic decompositions in MHD also display intriguing behavior [77], which is relevant because MHD also displays novel chaotic behavior compared to fluid turbulence as discussed below. Similar insights could be gained by considering error through triadic interactions. In [78], Rotunno and Snyder generalized the Lorenz model of propagation of errors upwards to “surface quasigeostrophic equations”. Here, it was found that the downscale error propagation is of vital importance for understanding the predictability dynamics, and that the importance of the interaction of large and small length scale features for the error has also been noted as long ago as 1957 [79]. It has also been shown that in weather forecasting, it is not only errors from the small scales that affect the predictability, but that large scale errors may be a limiting factor [80].
We looked at the cascade of error in [58] implicitly by testing perturbations made at different scales. The perturbations were spatially isolated to largest scales for the main results, but tests were performed with high wavenumber (large length scale) perturbations. The simulations showed that perturbations made at large length scales lead to an initial convergence of the two fields, suggesting that a portion of the error must cascade forwards to smaller scales in order to grow them, and that also there is a diminution of error. This was not the case for the perturbations made at small length scales (low wavenumbers). In addition, exponential growth in the low wave number region was seen if forcing was performed at intermediate wavenumbers, suggesting that the cascade of error is not necessarily dependent on the forwards cascade of energy, even if it does move forward. Further insights on the scale-of-error growth were made in [81] and examined more explicitly by Ge et al. in [67].
The fact that there is scale-dependent error growth has important implications for predictability and simulations of LES [62]. The scale dependence was one of the original motivators of the FSLE. In LES, if there is a scale-dependent error growth, the prediction horizon depends on the grid spacing. Numerical simulation on LES shows that as the resolution is increased, there is a trend towards the DNS value. LES simulations show that the smaller the resolution, the faster the perturbation/error can grow [62]. This tradeoff between predictability and error growth allows one to determine when “the point of diminishing returns is reached”, and Budanur and Kantz find that “the forecast horizon cannot be extended by simply improving resolution alone”.

4.2.3. Linear Growth of Error

After the initial exponential growth phase, linear growth of error was seen in multiple simulations including ours [38,58]. It was found to depend linearly on dissipation, and this was linked to KS entropy. However, a possibility was mentioned that this could be forcing-dependent [38,67]. At late stages, the predictability is dominated by these late stages, which allows for the long time predictability in atmospheric systems. In MHD, as mentioned below, we further found that the linear growth rate depended on the dissipation rate of the particular field [82]. Since in cosmological scale magnetic fields, the magnetic dissipation can be extremely low, this allows, in principle, extremely long prediction horizons even in turbulent electrically conducting fluids.

4.2.4. FTLE Statistics

The FTLE are themselves a statistical quantity and can be analyzed in a similar way to fluctuations of the energy and dissipation, and looking at the instantaneous values of the FTLE reveals important details for the flows at that point in time [46]. FTLE have been used to characterize predictability of weather events [83] and extreme values in geophysics models [84].
A more systematic look at the FTLE was performed by us in [60]. It was noted that the statistics of the FTLE are in some way better than using the kinetic energy and many other statistics since their autocorrelation time was much shorter than for any other analyzed statistic or energy at any scale or dissipation, a finding similar to [46]. The Eulerian maximum Lyapunov exponent is a robust measure of level of chaos, and it provides better statistics than the Reynolds number or energy, which means it may be a better characterization of the system. We believe this is because the deviation of trajectories is sampled at every grid point in the simulation, the number of which is always greater than the number of grid points contributing to any spectral quantity. FTLEs are stable even for short steptimes for the FTLE, as we found in this paper, thus reducing the computational cost needed to obtain large sets and improve statistics (the energy varies only as the large eddy turnover time, for instance). Whilst the direct method offers better knowledge of the spectrum, we receive less statistics for the maximum Lyapunov exponent itself and its fluctuation. The simulations also showed that the FTLE offered steady statistics at a much shorter time after initialization than Re and E, which suggested that meaningful results may be gained from them much more quickly than these other quantities more affected by hysteresis. Allowing less relaxation time for the simulation is also useful for efficient simulation, which as a topic is explored further below. Contrary to the findings of Budanur and Kantz [62], we found in these DNS simulations that the maximum Lyapunov exponent did not depend on the lattice size at all, even when the simulation was super well-resolved.
Statistical analysis showed that the skewness of total energy, Re, and dissipation (which is a small length scale quantity) are positive, whereas the skewness is negative for λ . It was also found that the fluctuations of the maximum Lyapunov exponent were significant, and calculating them was not captured by a Gaussian propagation of errors. “ σ λ / λ /fluctuations are significant and they cannot be simply estimated by measuring the values of Re, σ R e , T 0 , and σ T 0 ” [60]. This suggests that the extreme values of the FTLE, which can drive growth of error in the long term, cannot be determined simply. The fluctuations of the Laypunov exponent increased with Re as a power law. The rate of these fluctuations allows one to estimate the time needed to run a simulation for good statistics, and is thus an explicit way of attaining efficient simulation. Specifically, they would allow us to predetermine over how many large eddy turnover times statistics need to be taken to be representative. The fact that FTLEs are the most robust statistical measure we could find offers impetus to the use of it as a standard measure in DNS of HIT, even working in MHD.

4.3. Efficient Simulation

Chaos can be seen in turbulence even at very small box sizes, as was seen already using 32 3 collocation points by Deissler in 1986 [12]. Whilst some chaotic properties may require heavy numerical work, using chaos as a measure can improve the efficiency of simulation by directing resources where they are most appropriately needed. For instance, if only a certain tolerance is needed on predictability, it may be computationally wasteful to exceed that tolerance by an order of magnitude. Understanding the chaotic measures of turbulence allows us to determine this point explicitly.
The fact that the predictability can be quantified was used to achieve more efficient simulations in [85], where Wang and Zaki used it to argue about the fact that the small length scale errors cannot contaminate the large length scales by the time it reaches downstream. Whilst the analysis of the Lyapunov spectrum in turbulence may seem of mainly theoretical interest, the application of the method to practical problems, such as analyzing separated flow around an airfoil [86], shows that it can be a useful tool even in practical situations. In this study, Fernandez and Wang found a dependence on the mesh refinement, as seen elsewhere.
Large-eddy simulations work by a filtering of the Navier–Stokes equations [87], with a model for subgrid scales. Within the context of LES, authors often use the term “error” to refer to the difference between a simulated value and the statistics that would have been obtained from a full simulation. The maximum Lyapunov exponent has been used to characterize the growth of numerical error in LES simulations by Bae and Lozano-Duran [88].
In [89], Nastac et al. proposed using Lyapunov exponents to evaluate LES quality by analyzing the specific problem of a turbulent jet flame and turbulent combustion and their effects on predictability. They used the Eulerian field, finding that whilst in the DNS regime there is no effect of grid size on predictability, here, there is, but only for inert flows. When they looked at reaction flows, they found a different effect and further studied how chaotic analysis can improve simulations as a metric for determining what are the relevant scales of a specific problem. The Lyapunov exponent can determine regions of high error growth, and thus the regions where grid resolution should be increased. The exponents can be used “to estimate the time required for upstream turbulence to propagate downstream”. The authors further noted that for an inert jet, this is trivial, but can be useful for complex flows and geometries generating turbulent structures. Chaotic measures analysis is a general method for LES and does not depend on the subgrid scale modeling. Once the limit of the maximum Lyapunov exponent is reached, it should approximate the chaotic dynamics accurately. Naive application of the method showed the physics of the inert jet, namely that the most important region was close to the nozzle. There was then more specific comment on how it applied to the physics of turbulent jets and showed the strength of the application, especially to how combustion can make the dynamics more predictable, partly due to flame relaminarization. This sort of analysis may be especially promising since measuring λ is fairly simple.
A different analysis of LES quality was performed in [90]. Wu et al. pointed out how Lyapunov exponents are useful in transient simulations, where statistically stationary flow is not achieved. This is also useful for capturing rare events, typical in combustion, happening at short time scales. The maximum Lyapunov exponent can be found for arbitrary geometries. They found that once the dynamical scales are fully resolved, the maximum Lyapunov exponent saturates, and thus we know the maximum resolution needed for efficient simulation. However, they warned that this cannot be known beforehand and requires repeated simulation to discover. Even so, the authors concluded that these metrics are “capable of identifying quantity-specific sensitivities with respect to the numerical resolution while requiring significantly less computational resources than acquiring profiles of conventional turbulent statistics”. This analysis was extended in Engelmann et al. [91], using the Shannon entropy in LES flows and its specific application to wall and shear flows. Here, they also found that Shannon entropy can be “a more stringent quality criterion than resolved energy” [91]. In effect, information entropy can be a useful tool for assessing LES quality, not just statistics of energy.
This idea that we can use the entropy production (or other dynamical systems quantities) to analyze the suitability of simulations can extend to other non-linear systems on unstructured grids (e.g., for biological simulations), and it can provide a different metric for determining where computational resources should be spent efficiently. Of course, resolvability criteria are important, but this is another tool.
In engineering devices, one common problem encountered is the need to optimize some time averaged quantity. When performed using least squares shadowing, this depends on the integration of a number of equations equal to the number of positive Lyapunov exponents. Thus, understanding the way that the Lyapunov spectrum changes with parameters such as the Reynolds number is of key importance. Various strategies have been proposed to surmount this problem [92]. A recent theoretical framework of data assimilation approaches, depending on the chaotic properties of turbulence, allows one “to infer small-scale turbulent structures based on observational data from large-scale ones” [93]. Lyapunov analysis has also been performed to provide predictions of combustion events in the context of sequential data assimilation [94].

4.4. Insights

4.4.1. Theoretical Extensions

The analysis of turbulence using DNS allows different mathematical theories to be tested numerically, such as the idea of superfast amplification and nonlinear saturation of perturbations [95]. Lyapunov analysis of turbulence has also aided thoughts about stochasticity, such as in renomalization group approaches by Eyink and Bandak [96]. This paper understands turbulence from a more dynamical systems approach, and the understanding of turbulence through the lens of chaos is intimately linked to this approach. The question of how the “forgetting” is achieved can be understood by the drift over time, D, where
D ( t ) = 1 2 ( u ( t ) u ( t 0 ) ) 2 d V .
The behavior for a sample run is shown in Figure 4 [97]. In additionl, we can define the phase space speed, s ( t ) , where
s ( t ) 2 = ( u ( t + Δ t ) u ( t ) ) 2 d V Δ t 2 ,
with Δ t being the simulation time step. The dependence on Re is shown in Figure 5, which only depends on the large length scale quantities indicating that the majority of the speed in phase space is dominated by the large length scales. This can be expected due to the fact that these scales contain the majority of the energy. That such easily obtainable quantities may be of use to theorists highlights that simulations of turbulence still contain rich veins to be drawn from and that a quest for ever larger box sizes is not necessarily the best use of computational resources.
Further, analysis of the maximum Lyapunov exponent on active turbulence may reveal interesting things about predictability of active matter [98], which is important since these are key biological processes. This analysis was performed in [99], where Mukherjee et al. found different behaviors at different levels of activity.

4.4.2. Changes in System Dimension

Changes in system dimension can change the chaotic properties of turbulence. We were motivated [100] by studies in critical phenomenon by Wilson and Fisher looking for a critical spatial dimension above which the critical exponent is given by mean field theory [101]. We wanted to discover if similar effects held in fluid turbulence. It has been suggested that K41 is in fact a mean field theory, exact only above a critical dimension [102,103]. In our previous studies, we found that chaotic properties are highly dependent on the spatial dimensions for two and three spatial dimensions [58,72,104]. We thus utilized such measures to examine higher spatial dimensions using EDQNM and some DNS [81].
In [104], we studied the transition from two-dimensional to three-dimensional turbulence using chaotic measures. These are very sharp measures as compared to looking at turbulence spectra, which are much more affected by noisy data. The study uses a box with decreasing height to study the transition, with a stochastic external forcing. As the box height decreases, the phenomenology of the system transitions from purely three-dimensional to an intermediate state of mixed two- and three-dimensional behavior and finally on to purely two-dimensional [105,106,107,108]. The very different behavior of vorticity between two and three dimensions is likely part of the explanation for the sharpness of chaotic measures in this case.
We observed a discontinuous jump in λ when properly scaled as shown in Figure 6 [104]. This was in contrast to more standard statistics which exhibited a smooth transition. Additionally the time history of the FTLE becomes more intermittent around the transition point showing large fluctuations; hence, these chaotic measures really are capturing information about the flow. The observed discontinuity in the maximum Lyapunov exponent as a function of box height as well as the enhanced fluctuations around this point are reminiscent of a phase transition. An example of these fluctuations is shown in Figure 7. This suggests that chaotic measures can reveal additional information that is not accessible to standard statistical measurement. The use of chaos to study phase transitions has only seen limited discussion in the literature [109,110,111]; however, our findings in [104] suggest this approach should be examined more carefully. We also saw an interesting response to system dimension in EDQNM models [81], which show strong dependence of the chaotic measures on dimension, discussed below.

4.4.3. Importance of Vortex Stretching

One of the themes that emerges from our and other studies into chaotic measures in turbulence is the importance of the vortical structures to the Eulerian chaos. We emphasized this in [81], which followed on from our previous work in [112] and looked specifically at the effect of spatial dimension on EQDNM models of turbulence. The introduction of this paper contains a succinct overview of previous uses of EDQNM to study dimensionality effects and further references can be found therein. Briefly, the EDQNM model started by using the Navier–Stokes equations in the Fourier space to derive an equation for the evolution of the third-order velocity field moment. A quasi-normal approximation was used to close the infinite moment hierarchy and further adapted, resulting in the eddy damped quasi-normal Markovian (EDQNM) approximation, as shown by Orszag [113]. In this study, we then derived an equation for the correlation spectrum between two different velocity fields under this closure which can then be used to measure the maximal Lyapunov exponent within the EDQNM closure.
The study was motivated by wanting to better understand criticality and scale invariance. Notably, EDQNM models provide an algebraic method by which to vary the system dimension resulting in little additional computational expense. Physically, it was found that the enstrophy production peaked near six dimensions, resulting in changed dynamical behavior.
Whilst the equations themselves are not chaotic, it is possible to deduce the error growth within the model [25,113]. We found that DNS and EDQNM offer a similar order of results [81], where DNS and EDQNM λ decrease with dimension. The EDQNM runs showed that the Ruelle relation breaks down even in EDQNM in a way that is similar to DNS because of a lack of “scaling range and hence no separation between the viscous and large-eddy time scales”. Most interestingly, a critical dimension was found for error growth. From the analysis of the error spectrum, we looked at the predictability time and found good agreement with DNS, but we obtained some late time saturation even in 3D. In terms of error cascades, there was a “competition between the generation of uncorrelated energy and its sweeping out by the cascade” [81]. In the EDQNM runs, we saw that this forward cascade saturates at a higher wavenumber. There was a cascade of correlated energy which is more effective with higher dimension. It was speculated that the diminution of error is then caused by the change in the vortex stretching term as we increase dimensions.
In [67], Ge et al. ran with these and other ideas to more explicitly analyze the effect of vortex stretching and compression on error growth, finding competition between them. The authors derived the evolution equation for the error for HIT, finding “three different mechanisms: internal production resulting from interactions between the strain rate and the velocity-difference field, dissipation of the velocity-difference field and external force input”, finding that, indeed, the production of error “can be related to the strain self-amplification and vortex-stretching”. They then clarified the effects that the forcing can have on the production of error. They found that “spontaneous decorrelation of a flow from its perturbed flow in the absence of external inputs of uncertainty can only occur through local compressions”, and a self-similar uncertainty spectrum similar to [74]. This spectrum implies a characteristic length of uncertainty that is defined and tracked. The authors associated it with possibly being related to the Taylor length scale, and linked this to the finding by Goto and Vassilicos in [114], where Ge et al. asserted the following: “the mean distance between stagnation points in a HIT and therefore tends to represent the average size of turbulent eddies which is highly weighted towards the more numerous smallest ones”. They proposed a large length scale random sweeping effect as being a possible driver of error growth, which is in agreement with the highly localized nature of the error found in [46] and is motivated by their finding that the chaos is driven by events with extremely high kurtosis and positive skewness. Finally, the authors concluded that “the identification of local compression and stretching events as key to the development of uncertainty means that future prediction methods may benefit from strategies for early detection of such events so as to concentrate maximum accuracy on the compression events and less accuracy on the stretching events”. This again proves the necessity of the method for efficient simulation, and how different strands of the theory and results combine to mutually advance each other.
Whilst we emphasised chaotic measures as an alternative for analyzing turbulent flow in this review, we do not believe that they are exclusive. We have already seen how the error spectrum is a useful way of understanding the error; it would be useful to also analyze the error transfer, as derived in [67]. This can be spectrally decomposed and provide a measure of the transfer of error across scales in an analogous way to the transfer of energy in HIT.

4.4.4. MHD

We extended the study of Eulerian chaotic measures by tracking trajectories in phase space in turbulence to MHD in [60,82], whilst a theoretical study by Fouxon et al. was conducted in [44]. Shell models of turbulence applied to MHD with Prandtl number ∼ 1 found λ ν 1 / 2 [115]; however, since shell models may not be appropriate for the study, we used DNS simulations.
In MHD, a magnetic field is coupled to an electrically conducting fluid, with dynamics given by the MHD equations [116,117]. The equivalent of the Ruelle relation was complicated by the fact that there are multiple microscales in this system. The non-linearity of the equation is direct for the velocity field, but for the magnetic field the non-linearity arises indirectly through coupling with the velocity field. Thus, we may expect that the chaos should be dominated by the velocity field and not the magnetic field, which was indeed found to be the case. A test of removing the Lorentz force confirmed that it really is the velocity that drives the chaos [82]. No stable trend was associated with the Prandtl number, though using λ 1 / τ was found to fit better due to the noise from the Prandtl number and the fact that the relationship between Re and τ is not as simple in MHD as for HD. For relationship λ T 0 Re α , it was found that α = 0.43 ± 0.09 .
Interestingly, we found that magnetic helicity can destroy the chaos, which was not just due to inverse transfer of energy due to chaos as seen in two-dimensional turbulence. Magnetic helicity effectively increases the predictability of the system. In other simulations, we also found a linear growth of error, as for the pure hydrodynamic case. However, this linear growth rate was separate for both fields, depending on the corresponding field’s dissipation rate. Since magnetic dissipation can be much smaller than the kinetic dissipation, the magnetic field has much greater predictability than the velocity field. This suggests that the velocity details might vary greatly whilst the magnetic field is very similar. Both of these effects suggest that galactic scale magnetic fields, if they have high magnetic helicity or Prandtl number, though turbulent, may have very long predictability times.

5. Discussion and Conclusions

In this review, we tried to highlight the contributions of many different researchers in advancing the field of knowledge of turbulence and specifically chaotic measures as a means to understanding it. There are many different strands to this thought, but we pick out a few key ones below:
  • Turbulence is a multiscale phenomena both spatially and temporally. This is true even in chaos, where the different timescales interact to lead to a long-term large length scale predictability co-existing with increasingly rapid error growth at the smallest scales.
  • The Eulerian chaos and vortex stretching appear to be fundamentally linked, and may contribute to the changes seen in different dimensions.
  • The Lyapunov spectrum may emerge as a practical measure for a complete description of the fluid flow.
  • Chaotic measures, especially the maximum Lyapunov exponent, are more stable than spectral quantities and show spatial features that are hard to discern from spectral quantities alone.
  • By using chaotic measures, we can quantify simulations in a model-free way that does not rely on a priori unknown spectral quantities and thus improve computational efficiency.
  • Chaotic measures can be used to detect phase transitions.
These recent discoveries are not only of theoretical but also of practical importance, since these measures are easily implementable in simulations of real-world systems. In addition, we hope that this review will be useful for researchers exploring turbulence in other contexts, such as in spontaneous stochasticity [63,96] and molecular-gas-dynamics simulations [118], where chaotic measures can also provide an alternative measure to spectral quantities to understand the difference between these and the Navier–Stokes equations.

Author Contributions

Conceptualization, R.D.J.G.H., D.C. and A.B.; writing—original draft, R.D.J.G.H., D.C. and A.B.; writing—review and editing, R.D.J.G.H., D.C. and A.B. All authors have read and agreed to the published version of the manuscript.

Funding

R.H. acknowledges partial funding by the University of Oslo, UiO: Life Science through the convergence environment ITOM. A.B. acknowledges partial funding by STFC.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created for this review.

Acknowledgments

For the purpose of open access, the authors applied a Creative Commons Attribution (CC BY) licence to any Author Accepted Manuscript version arising from this submission.

Conflicts of Interest

Daniel Clark is an employee of Thornton Tomasetti. This paper reflects the views of the scientist and not the company.

References

  1. McComb, W.D. The Physics of Fluid Turbulence; Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  2. McComb, W.D. Homogeneous, Isotropic Turbulence: Phenomenology, Renormalization and Statistical Closures; OUP Oxford: Oxford, UK, 2014; Volume 162. [Google Scholar]
  3. Kraichnan, R.H. The structure of isotropic turbulence at very high Reynolds numbers. J. Fluid Mech. 1959, 5, 497–543. [Google Scholar] [CrossRef]
  4. Zhou, Y. Turbulence theories and statistical closure approaches. Phys. Rep. 2021, 935, 1–117. [Google Scholar]
  5. Taylor, G.I. The spectrum of turbulence. Proc. R. Soc. London. Ser. Math. Phys. Sci. 1938, 164, 476–490. [Google Scholar] [CrossRef]
  6. Kolmogorov, A.N. The local structure of turbulence in incompressible viscous fluid for very large Reynolds. Dokl. Akad. Nauk. SSSR 1941, 30, 301. [Google Scholar]
  7. Kolmogorov, A.N. A refinement of previous hypotheses concerning the local structure of turbulence in a viscous incompressible fluid at high Reynolds number. J. Fluid Mech. 1962, 13, 82–85. [Google Scholar] [CrossRef]
  8. Frisch, U. Turbulence: The Legacy of an Kolmogorov; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  9. Qian, J. Inertial range and the finite Reynolds number effect of turbulence. Phys. Rev. E 1997, 55, 337. [Google Scholar] [CrossRef]
  10. Qian, J. Slow decay of the finite Reynolds number effect of turbulence. Phys. Rev. E 1999, 60, 3409. [Google Scholar] [CrossRef] [PubMed]
  11. Sagaut, P.; Cambon, C. Homogeneous Turbulence Dynamics; Springer: Berlin/Heidelberg, Germany, 2008; Volume 15. [Google Scholar]
  12. Deissler, R. Is Navier–Stokes turbulence chaotic? Phys. Fluids 1986, 29, 1453–1457. [Google Scholar] [CrossRef]
  13. Bohr, T.; Jensen, M.H.; Paladin, G.; Vulpiani, A. Dynamical Systems Approach to Turbulence; Cambridge University Press & Assessment: Cambridge, UK, 1998. [Google Scholar]
  14. Taylor, G.I. Diffusion by continuous movements. Proc. Lond. Math. Soc. 1922, 2, 196–212. [Google Scholar] [CrossRef]
  15. Richardson, L.F. Atmospheric diffusion shown on a distance-neighbour graph. Proc. R. Soc. London. Ser. Contain. Pap. Math. Phys. Character 1926, 110, 709–737. [Google Scholar]
  16. Salazar, J.P.; Collins, L.R. Two-particle dispersion in isotropic turbulent flows. Annu. Rev. Fluid Mech. 2009, 41, 405–432. [Google Scholar] [CrossRef]
  17. Biferale, L.; Boffetta, G.; Celani, A.; Devenish, B.; Lanotte, A.; Toschi, F. Lagrangian statistics of particle pairs in homogeneous isotropic turbulence. Phys. Fluids 2005, 17, 115101. [Google Scholar] [CrossRef]
  18. Haller, G. Lagrangian coherent structures. Annu. Rev. Fluid Mech. 2015, 47, 137–162. [Google Scholar] [CrossRef]
  19. Ottino, J.M. Mixing, chaotic advection, and turbulence. Annu. Rev. Fluid Mech. 1990, 22, 207–254. [Google Scholar] [CrossRef]
  20. Eckhardt, B.; Schneider, T.M.; Hof, B.; Westerweel, J. Turbulence transition in pipe flow. Annu. Rev. Fluid Mech. 2007, 39, 447–468. [Google Scholar] [CrossRef]
  21. Lorenz, E.N. Deterministic nonperiodic flow. J. Atmos. Sci. 1963, 20, 130–141. [Google Scholar] [CrossRef]
  22. Lorenz, E.N. The predictability of a flow which possesses many scales of motion. Tellus 1969, 21, 289–307. [Google Scholar] [CrossRef]
  23. Aurell, E.; Boffetta, G.; Crisanti, A.; Paladin, G.; Vulpiani, A. Predictability in the large: An extension of the concept of Lyapunov exponent. J. Phys. Math. Gen. 1997, 30, 1. [Google Scholar] [CrossRef]
  24. Leith, C. Atmospheric predictability and two-dimensional turbulence. J. Atmos. Sci. 1971, 28, 145–161. [Google Scholar] [CrossRef]
  25. Leith, C.; Kraichnan, R. Predictability of turbulent flows. J. Atmos. Sci. 1972, 29, 1041–1058. [Google Scholar] [CrossRef]
  26. Boffetta, G.; Cencini, M.; Falcioni, M.; Vulpiani, A. Predictability: A way to characterize complexity. Phys. Rep. 2002, 356, 367–474. [Google Scholar] [CrossRef]
  27. Nazarenko, S. Wave Turbulence; Springer: Berlin/Heidelberg, Germany, 2011; Volume 825. [Google Scholar]
  28. Korotkevich, A.O. Inverse cascade spectrum of gravity waves in the presence of a condensate: A direct numerical simulation. Phys. Rev. Lett. 2023, 130, 264002. [Google Scholar] [CrossRef] [PubMed]
  29. Falcon, E.; Mordant, N. Experiments in surface gravity–capillary wave turbulence. Annu. Rev. Fluid Mech. 2022, 54, 1–25. [Google Scholar] [CrossRef]
  30. Pan, Y.; Yue, D.K. Direct numerical investigation of turbulence of capillary waves. Phys. Rev. Lett. 2014, 113, 094501. [Google Scholar] [CrossRef]
  31. Kochurin, E.; Ricard, G.; Zubarev, N.; Falcon, E. Three-dimensional direct numerical simulation of free-surface magnetohydrodynamic wave turbulence. Phys. Rev. E 2022, 105, L063101. [Google Scholar] [CrossRef] [PubMed]
  32. Ricard, G.; Falcon, E. Transition from wave turbulence to acousticlike shock-wave regime. Phys. Rev. Fluids 2023, 8, 014804. [Google Scholar] [CrossRef]
  33. Griffin, A.; Krstulovic, G.; L’vov, V.S.; Nazarenko, S. Energy spectrum of two-dimensional acoustic turbulence. Phys. Rev. Lett. 2022, 128, 224501. [Google Scholar] [CrossRef]
  34. Kochurin, E.A.; Kuznetsov, E.A. Direct numerical simulation of acoustic turbulence: Zakharov–Sagdeev spectrum. JETP Lett. 2022, 116, 863–868. [Google Scholar] [CrossRef]
  35. Oseledec, V.I. A multiplicative ergodic theorem, Lyapunov characteristic numbers for dynamical systems. Trans. Mosc. Math. Soc. 1968, 19, 197–231. [Google Scholar]
  36. Wolf, A.; Swift, J.B.; Swinney, H.L.; Vastano, J.A. Determining Lyapunov exponents from a time series. Phys. Nonlinear Phenom. 1985, 16, 285–317. [Google Scholar] [CrossRef]
  37. Dettmann, C.P.; Morriss, G. Proof of Lyapunov exponent pairing for systems at constant kinetic energy. Phys. Rev. E 1996, 53, R5545. [Google Scholar] [CrossRef]
  38. Boffetta, G.; Musacchio, S. Chaos and predictability of homogeneous-isotropic turbulence. Phys. Rev. Lett. 2017, 119, 054102. [Google Scholar] [CrossRef] [PubMed]
  39. Eckmann, J.P.; Ruelle, D. Ergodic theory of chaos and strange attractors. Rev. Mod. Phys. 1985, 57, 617. [Google Scholar] [CrossRef]
  40. Schuster, H.G.; Just, W. Deterministic Chaos: An Introduction; John Wiley & Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  41. Kaplan, J.; Yorke, J. Functional differential equations and approximation of fixed points. Lect. Notes Math. 1979, 730, 204–227. [Google Scholar]
  42. Frederickson, P.; Kaplan, J.L.; Yorke, E.D.; Yorke, J.A. The Liapunov dimension of strange attractors. J. Differ. Equ. 1983, 49, 185–207. [Google Scholar] [CrossRef]
  43. Ruelle, D. Microscopic fluctuations and turbulence. Phys. Lett. A 1979, 72, 81–82. [Google Scholar] [CrossRef]
  44. Fouxon, I.; Feinberg, J.; Käpylä, P.; Mond, M. Reynolds number dependence of Lyapunov exponents of turbulence and fluid particles. Phys. Rev. E 2021, 103, 033110. [Google Scholar] [CrossRef] [PubMed]
  45. Crisanti, A.; Jensen, M.; Vulpiani, A.; Paladin, G. Intermittency and predictability in turbulence. Phys. Rev. Lett. 1993, 70, 166. [Google Scholar] [CrossRef]
  46. Mohan, P.; Fitzsimmons, N.; Moser, R.D. Scaling of Lyapunov exponents in homogeneous isotropic turbulence. Phys. Rev. Fluids 2017, 2, 114606. [Google Scholar] [CrossRef]
  47. Landau, L. Em lifshitz, fluid mechanics. In Course of Theoretical Physics; Elsevier: Amsterdam, The Netherlands, 1959; Volume 6. [Google Scholar]
  48. Ruelle, D. Large volume limit of the distribution of characteristic exponents in turbulence. Commun. Math. Phys. 1982, 87, 287–302. [Google Scholar] [CrossRef]
  49. Constantin, P.; Foias, C.; Manley, O.P.; Temam, R. Determining modes and fractal dimension of turbulent flows. J. Fluid Mech. 1985, 150, 427–440. [Google Scholar] [CrossRef]
  50. Gibbon, J.D.; Titi, E. Attractor dimension and small length scale estimates for the three-dimensional Navier-Stokes equations. Nonlinearity 1997, 10, 109. [Google Scholar] [CrossRef]
  51. Meneveau, C.; Nelkin, M. Attractor size in intermittent turbulence. Phys. Rev. A 1989, 39, 3732. [Google Scholar] [CrossRef] [PubMed]
  52. Batchelor, G.K. Computation of the energy spectrum in homogeneous two-dimensional turbulence. Phys. Fluids 1969, 12, II-233. [Google Scholar] [CrossRef]
  53. Kraichnan, R.H. Inertial ranges in two-dimensional turbulence. Phys. Fluids 1967, 10, 1417. [Google Scholar] [CrossRef]
  54. Leith, C.E. Diffusion approximation for two-dimensional turbulence. Phys. Fluids 1968, 11, 671–672. [Google Scholar] [CrossRef]
  55. Kraichnan, R.H. Inertial-range transfer in two-and three-dimensional turbulence. J. Fluid Mech. 1971, 47, 525–535. [Google Scholar] [CrossRef]
  56. Ohkitani, K. Log-corrected energy spectrum and dimension of attractor in two-dimensional turbulence. Phys. Fluids Fluid Dyn. 1989, 1, 451–452. [Google Scholar] [CrossRef]
  57. Lieb, E.H. On characteristic exponents in turbulence. Commun. Math. Phys. 1984, 92, 473–480. [Google Scholar] [CrossRef]
  58. Berera, A.; Ho, R.D. Chaotic properties of a turbulent isotropic fluid. Phys. Rev. Lett. 2018, 120, 024101. [Google Scholar] [CrossRef]
  59. Aurell, E.; Boffetta, G.; Crisanti, A.; Paladin, G.; Vulpiani, A. Growth of noninfinitesimal perturbations in turbulence. Phys. Rev. Lett. 1996, 77, 1262. [Google Scholar] [CrossRef]
  60. Ho, R.D.; Armua, A.; Berera, A. Fluctuations of Lyapunov exponents in homogeneous and isotropic turbulence. Phys. Rev. Fluids 2020, 5, 024602. [Google Scholar] [CrossRef]
  61. Mukherjee, S.; Schalkwijk, J.; Jonker, H.J. Predictability of dry convective boundary layers: An LES study. J. Atmos. Sci. 2016, 73, 2715–2727. [Google Scholar] [CrossRef]
  62. Budanur, N.B.; Kantz, H. Scale-dependent error growth in Navier-Stokes simulations. Phys. Rev. E 2022, 106, 045102. [Google Scholar] [CrossRef] [PubMed]
  63. Bandak, D.; Mailybaev, A.A.; Eyink, G.L.; Goldenfeld, N. Spontaneous stochasticity amplifies even thermal noise to the largest scales of turbulence in a few eddy turnover times. Phys. Rev. Lett. 2024, 132, 104002. [Google Scholar] [CrossRef]
  64. Murugan, S.D.; Kumar, D.; Bhattacharjee, S.; Ray, S.S. Many-body chaos in thermalized fluids. Phys. Rev. Lett. 2021, 127, 124501. [Google Scholar] [CrossRef]
  65. McComb, W.; Berera, A.; Yoffe, S.; Linkmann, M. Energy transfer and dissipation in forced isotropic turbulence. Phys. Rev. E 2015, 91, 043013. [Google Scholar] [CrossRef]
  66. Boffetta, G.; Celani, A.; Crisanti, A.; Vulpiani, A. Predictability in two-dimensional decaying turbulence. Phys. Fluids 1997, 9, 724–734. [Google Scholar] [CrossRef]
  67. Ge, J.; Rolland, J.; Vassilicos, J.C. The production of uncertainty in three-dimensional Navier–Stokes turbulence. J. Fluid Mech. 2023, 977, A17. [Google Scholar] [CrossRef]
  68. Berera, A.; Clark, D. Information production in homogeneous isotropic turbulence. Phys. Rev. E 2019, 100, 041101. [Google Scholar] [CrossRef]
  69. Hassanaly, M.; Raman, V. Lyapunov spectrum of forced homogeneous isotropic turbulent flows. Phys. Rev. Fluids 2019, 4, 114608. [Google Scholar] [CrossRef]
  70. Hassanaly, M.; Raman, V. Numerical convergence of the Lyapunov spectrum computed using low Mach number solvers. J. Comput. Phys. 2019, 386, 467–485. [Google Scholar] [CrossRef]
  71. Keefe, L.; Moin, P.; Kim, J. The dimension of attractors underlying periodic turbulent Poiseuille flow. J. Fluid Mech. 1992, 242, 1–29. [Google Scholar] [CrossRef]
  72. Clark, D.; Tarra, L.; Berera, A. Chaos and information in two-dimensional turbulence. Phys. Rev. Fluids 2020, 5, 064608. [Google Scholar] [CrossRef]
  73. Vasey, G.; Messenger, D.; Bortz, D.; Christlieb, A.; O’Shea, B. Influence of initial conditions on data-driven model identification and information entropy for ideal mhd problems. arXiv 2023, arXiv:2312.05339. [Google Scholar]
  74. Yoshimatsu, K.; Ariki, T. Error growth in three-dimensional homogeneous turbulence. J. Phys. Soc. Jpn. 2019, 88, 124401. [Google Scholar] [CrossRef]
  75. Waleffe, F. The nature of triad interactions in homogeneous turbulence. Phys. Fluids Fluid Dyn. 1992, 4, 350–363. [Google Scholar] [CrossRef]
  76. Biferale, L.; Musacchio, S.; Toschi, F. Inverse energy cascade in three-dimensional isotropic turbulence. Phys. Rev. Lett. 2012, 108, 164501. [Google Scholar] [CrossRef]
  77. Linkmann, M.; Berera, A.; McKay, M.; Jäger, J. Helical mode interactions and spectral transfer processes in magnetohydrodynamic turbulence. J. Fluid Mech. 2016, 791, 61–96. [Google Scholar] [CrossRef]
  78. Rotunno, R.; Snyder, C. A generalization of Lorenz’s model for the predictability of flows with many scales of motion. J. Atmos. Sci. 2008, 65, 1063–1076. [Google Scholar] [CrossRef]
  79. Thompson, P.D. Uncertainty of initial state as a factor in the predictability of large scale atmospheric flow patterns. Tellus 1957, 9, 275–295. [Google Scholar] [CrossRef]
  80. Durran, D.R.; Gingrich, M. Atmospheric predictability: Why butterflies are not of practical importance. J. Atmos. Sci. 2014, 71, 2476–2488. [Google Scholar] [CrossRef]
  81. Clark, D.; Armua, A.; Ho, R.D.; Berera, A. Critical transition to a non-chaotic regime in isotropic turbulence. J. Fluid Mech. 2022, 930, A17. [Google Scholar] [CrossRef]
  82. Ho, R.D.; Berera, A.; Clark, D. Chaotic behavior of Eulerian magnetohydrodynamic turbulence. Phys. Plasmas 2019, 26, 042303. [Google Scholar] [CrossRef]
  83. Pandey, L.K.; Dwivedi, S.; Martin, M. Short-term predictability of the Bay of Bengal region using a high-resolution Indian ocean model. Mar. Geod. 2021, 44, 215–237. [Google Scholar] [CrossRef]
  84. Sterk, A.; Holland, M.; Rabassa, P.; Broer, H.; Vitolo, R. Predictability of extreme values in geophysical models. Nonlinear Process. Geophys. 2012, 19, 529–539. [Google Scholar] [CrossRef]
  85. Wang, M.; Zaki, T.A. Synchronization of turbulence in channel flow. J. Fluid Mech. 2022, 943, A4. [Google Scholar] [CrossRef]
  86. Fernandez, P.; Wang, Q. Lyapunov spectrum of the separated flow around the NACA 0012 airfoil and its dependence on numerical discretization. J. Comput. Phys. 2017, 350, 453–469. [Google Scholar] [CrossRef]
  87. Leonard, A. Energy cascade in large-eddy simulations of turbulent fluid flows. In Advances in Geophysics; Elsevier: Amsterdam, The Netherlands, 1975; Volume 18, pp. 237–248. [Google Scholar]
  88. Bae, H.J.; Lozano-Duran, A. Numerical and modeling error assessment of large-eddy simulation using direct-numerical-simulation-aided large-eddy simulation. arXiv 2022, arXiv:2208.02354. [Google Scholar]
  89. Nastac, G.; Labahn, J.W.; Magri, L.; Ihme, M. Lyapunov exponent as a metric for assessing the dynamic content and predictability of large-eddy simulations. Phys. Rev. Fluids 2017, 2, 094606. [Google Scholar] [CrossRef]
  90. Wu, H.; Ma, P.C.; Lv, Y.; Ihme, M. Lyapunov exponent and Wasserstein metric as validation tools for assessing short-time dynamics and quantitative model evaluation of large-eddy simulation. In Proceedings of the 2018 AIAA Aerospace Sciences Meeting, Kissimmee, FL, USA, 8–12 January 2018; p. 0440. [Google Scholar]
  91. Engelmann, L.; Ihme, M.; Wlokas, I.; Kempf, A. Towards the suitability of information entropy as an LES quality indicator. Flow Turbul. Combust. 2022, 108, 353–385. [Google Scholar] [CrossRef]
  92. Kantarakias, K.D.; Papadakis, G. Sensitivity analysis of chaotic systems using a frequency-domain shadowing approach. J. Comput. Phys. 2023, 474, 111757. [Google Scholar] [CrossRef]
  93. Inubushi, M.; Saiki, Y.; Kobayashi, M.U.; Goto, S. Characterizing Small-Scale Dynamics of Navier-Stokes Turbulence with Transverse Lyapunov Exponents: A Data Assimilation Approach. Phys. Rev. Lett. 2023, 131, 254001. [Google Scholar] [CrossRef] [PubMed]
  94. Magri, L.; Doan, N.A.K. Physics-informed data-driven prediction of turbulent reacting flows with Lyapunov analysis and sequential data assimilation. In Data Analysis for Direct Numerical Simulations of Turbulent Combustion: From Equation-Based Analysis to Machine Learning; Springer: Berlin/Heidelberg, Germany, 2020; pp. 177–196. [Google Scholar]
  95. Li, Y.C.; Ho, R.D.; Berera, A.; Feng, Z. Superfast amplification and superfast nonlinear saturation of perturbations as a mechanism of turbulence. J. Fluid Mech. 2020, 904, A27. [Google Scholar] [CrossRef]
  96. Eyink, G.L.; Bandak, D. Renormalization group approach to spontaneous stochasticity. Phys. Rev. Res. 2020, 2, 043161. [Google Scholar] [CrossRef]
  97. Ho, R.D. High dimensional chaotic systems which behave like random walks in state space. arXiv 2019, arXiv:1908.05989. [Google Scholar]
  98. Singh, R.K.; Mukherjee, S.; Ray, S.S. Lagrangian manifestation of anomalies in active turbulence. Phys. Rev. Fluids 2022, 7, 033101. [Google Scholar] [CrossRef]
  99. Mukherjee, S.; Singh, R.K.; James, M.; Ray, S.S. Intermittency, fluctuations and maximal chaos in an emergent universal state of active turbulence. Nat. Phys. 2023, 19, 891–897. [Google Scholar] [CrossRef]
  100. Berera, A.; Ho, R.D.; Clark, D. Homogeneous isotropic turbulence in four spatial dimensions. Phys. Fluids 2020, 32, 085107. [Google Scholar] [CrossRef]
  101. Wilson, K.G.; Fisher, M.E. Critical exponents in 3.99 dimensions. Phys. Rev. Lett. 1972, 28, 240. [Google Scholar] [CrossRef]
  102. Siggia, E.D. Origin of intermittency in fully developed turbulence. Phys. Rev. A 1977, 15, 1730. [Google Scholar] [CrossRef]
  103. Bell, T.L.; Nelkin, M. Time-dependent scaling relations and a cascade model of turbulence. J. Fluid Mech. 1978, 88, 369–391. [Google Scholar] [CrossRef]
  104. Clark, D.; Armua, A.; Freeman, C.; Brener, D.J.; Berera, A. Chaotic measure of the transition between two-and three-dimensional turbulence. Phys. Rev. Fluids 2021, 6, 054612. [Google Scholar] [CrossRef]
  105. Benavides, S.J.; Alexakis, A. Critical transitions in thin layer turbulence. J. Fluid Mech. 2017, 822, 364–385. [Google Scholar] [CrossRef]
  106. Ecke, R.E. From 2D to 3D in fluid turbulence: Unexpected critical transitions. J. Fluid Mech. 2017, 828, 1–4. [Google Scholar] [CrossRef]
  107. van Kan, A.; Alexakis, A. Condensates in thin-layer turbulence. J. Fluid Mech. 2019, 864, 490–518. [Google Scholar] [CrossRef]
  108. Musacchio, S.; Boffetta, G. Split energy cascade in turbulent thin fluid layers. Phys. Fluids 2017, 29. [Google Scholar] [CrossRef]
  109. Butera, P.; Caravati, G. Phase transitions and Lyapunov characteristic exponents. Phys. Rev. A 1987, 36, 962. [Google Scholar] [CrossRef] [PubMed]
  110. Caiani, L.; Casetti, L.; Clementi, C.; Pettini, M. Geometry of dynamics, Lyapunov exponents, and phase transitions. Phys. Rev. Lett. 1997, 79, 4361. [Google Scholar] [CrossRef]
  111. Barré, J.; Dauxois, T. Lyapunov exponents as a dynamical indicator of a phase transition. Europhys. Lett. 2001, 55, 164. [Google Scholar] [CrossRef]
  112. Clark, D.; Ho, R.D.; Berera, A. Effect of spatial dimension on a model of fluid turbulence. J. Fluid Mech. 2021, 912, A40. [Google Scholar] [CrossRef]
  113. Orszag, S.A. Analytical theories of turbulence. J. Fluid Mech. 1970, 41, 363–386. [Google Scholar] [CrossRef]
  114. Goto, S.; Vassilicos, J. The dissipation rate coefficient of turbulence is not universal and depends on the internal stagnation point structure. Phys. Fluids 2009, 21, 035104. [Google Scholar] [CrossRef]
  115. Grappin, R.; Leorat, J.; Pouquet, A. Computation of the dimension of a model of fully developed turbulence. J. Phys. 1986, 47, 1127–1136. [Google Scholar] [CrossRef]
  116. Elsasser, W.M. The hydromagnetic equations. Phys. Rev. 1950, 79, 183. [Google Scholar] [CrossRef]
  117. Cramer, K.R.; Pai, S.I. Magnetofluid Dynamics for Engineers and Applied Physicists; International Atomic Energy Agency: Vienna, Austria, 1973. [Google Scholar]
  118. McMullen, R.M.; Krygier, M.C.; Torczynski, J.R.; Gallis, M.A. Navier-Stokes equations do not describe the smallest scales of turbulence in gases. Phys. Rev. Lett. 2022, 128, 114501. [Google Scholar] [CrossRef]
Figure 1. Re against λ T 0 with dashed red line showing fit of α = 0.53 whilst solid black line shows α = 0.5 . Inset shows a direct test of constant λ τ . Figure taken from [58].
Figure 1. Re against λ T 0 with dashed red line showing fit of α = 0.53 whilst solid black line shows α = 0.5 . Inset shows a direct test of constant λ τ . Figure taken from [58].
Atmosphere 15 01053 g001
Figure 2. Lyapunov spectrum for Re from 90 to 212, rescaled by h K S and number of positive exponents. Scaling is indicated by the collapse onto the same curve. Figure taken from [68].
Figure 2. Lyapunov spectrum for Re from 90 to 212, rescaled by h K S and number of positive exponents. Scaling is indicated by the collapse onto the same curve. Figure taken from [68].
Atmosphere 15 01053 g002
Figure 3. Characteristic error spectrum E d ( k ) in black, shown for Re = 2500 with log-log main figure and semi-log inset. Main figure has dashed line k 3 , inset has exponential , both in red. Figure taken from [58].
Figure 3. Characteristic error spectrum E d ( k ) in black, shown for Re = 2500 with log-log main figure and semi-log inset. Main figure has dashed line k 3 , inset has exponential , both in red. Figure taken from [58].
Atmosphere 15 01053 g003
Figure 4. Forgetting rate with time shown as change in drift energy D. The dashed line in main figure is the prediction D = 2 E ( 1 exp ( t / τ ) and in the inset it is t 2 . Figure taken from [97].
Figure 4. Forgetting rate with time shown as change in drift energy D. The dashed line in main figure is the prediction D = 2 E ( 1 exp ( t / τ ) and in the inset it is t 2 . Figure taken from [97].
Atmosphere 15 01053 g004
Figure 5. Normalized phase space speed and Re. Interestingly, this uses only large length scale quantities: T, the large eddy turnover time; and u, the rms velocity. Dashed line shows a least squares fit.
Figure 5. Normalized phase space speed and Re. Interestingly, this uses only large length scale quantities: T, the large eddy turnover time; and u, the rms velocity. Dashed line shows a least squares fit.
Atmosphere 15 01053 g005
Figure 6. Reynolds number scaling of the maximum Lyapunov exponent in thin-layer turbulence. A discontinuous jump is found as a function of the aspect ratio of the system. Figure taken from [104].
Figure 6. Reynolds number scaling of the maximum Lyapunov exponent in thin-layer turbulence. A discontinuous jump is found as a function of the aspect ratio of the system. Figure taken from [104].
Atmosphere 15 01053 g006
Figure 7. Time history of the maximal FTLE computed near to the transition point (black) and far from the transition (blue). The statistics of the point near the transition displays large fluctuations. Figure taken from [104].
Figure 7. Time history of the maximal FTLE computed near to the transition point (black) and far from the transition (blue). The statistics of the point near the transition displays large fluctuations. Figure taken from [104].
Atmosphere 15 01053 g007
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

Ho, R.D.J.G.; Clark, D.; Berera, A. Chaotic Measures as an Alternative to Spectral Measures for Analysing Turbulent Flow. Atmosphere 2024, 15, 1053. https://doi.org/10.3390/atmos15091053

AMA Style

Ho RDJG, Clark D, Berera A. Chaotic Measures as an Alternative to Spectral Measures for Analysing Turbulent Flow. Atmosphere. 2024; 15(9):1053. https://doi.org/10.3390/atmos15091053

Chicago/Turabian Style

Ho, Richard D. J. G., Daniel Clark, and Arjun Berera. 2024. "Chaotic Measures as an Alternative to Spectral Measures for Analysing Turbulent Flow" Atmosphere 15, no. 9: 1053. https://doi.org/10.3390/atmos15091053

APA Style

Ho, R. D. J. G., Clark, D., & Berera, A. (2024). Chaotic Measures as an Alternative to Spectral Measures for Analysing Turbulent Flow. Atmosphere, 15(9), 1053. https://doi.org/10.3390/atmos15091053

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