Markov property of Lagrangian turbulence

Based on direct numerical simulations with point-like inertial particles, with Stokes numbers, $\textrm{St}=0, 0.5$, $3$, and $6$, transported by homogeneous and isotropic turbulent flows, we present in this letter for the first time evidence for the existence of Markov property in Lagrangian turbulence. We show that the Markov property is valid for a finite step size larger than a Stokes number-dependent Einstein-Markov coherence time scale. This enables the description of multi-scale statistics of Lagrangian particles by Fokker-Planck equations, which can be embedded in an interdisciplinary approach linking the statistical description of turbulence with fluctuation theorems of non-equilibrium stochastic thermodynamics and local flow structures. The formalism allows estimation of the stochastic thermodynamics entropy exchange associated with the particles' Lagrangian trajectories. Entropy consuming trajectories of the particles are related to specific evolution of velocity increments through scales and may be seen as intermittent structures. Statistical features of Lagrangian paths and entropy values are thus fixed by the fluctuation theorems.

Abstract -Based on direct numerical simulations with point-like inertial particles, with Stokes numbers, St = 0, 0.5, 3, and 6, transported by homogeneous and isotropic turbulent flows, we present in this letter for the first time evidence for the existence of Markov property in Lagrangian turbulence. We show that the Markov property is valid for a finite step size larger than a Stokes number-dependent Einstein-Markov coherence time scale. This enables the description of multi-scale statistics of Lagrangian particles by Fokker-Planck equations, which can be embedded in an interdisciplinary approach linking the statistical description of turbulence with fluctuation theorems of non-equilibrium stochastic thermodynamics and local flow structures. The formalism allows estimation of the stochastic thermodynamics entropy exchange associated with the particles' Lagrangian trajectories. Entropy consuming trajectories of the particles are related to specific evolution of velocity increments through scales and may be seen as intermittent structures. Statistical features of Lagrangian paths and entropy values are thus fixed by the fluctuation theorems.
Introduction.-The physics of particles submerged in fluids has played a central role in the development of statistical mechanics, and in our current understanding of out-of-equilibrium systems.
Inertial particles are inclusions in the flow which are denser or lighter than the fluid, and have a size smaller or larger than the smallest relevant flow scale (the scale of the smallest eddies, also called the flow dissipative scale). Such particles are carried by the fluid, but they also have their own inertia, and therefore the fluid and the particles have different dynamics. In the limit of point-wise particles with negligible inertia, the particles become Lagrangian tracers, which perfectly follow the fluid elements. In all these two-phase systems the mechanisms that explain how turbulence affects the motion of the particles are not completely clear. As an example, turbulence can both enhance or hinder the settling velocity of inertial particles [1]. For heavy particles, an initially homogeneous distribution of particles may, after interacting with a turbulent flow, regroup into clusters forming dense areas and voids, in a phenomenon called preferential concentration where turbulence somehow unmixes the particles [2]. Thus, the velocity and statistics of tracers and inertial particles is also different: while tracers sample the flow homogeneously, inertial particles only see specific regions of the flow, with either low acceleration or vorticity, and their velocity can differ significantly from the fluid velocity. It is only in the last decades when sufficient time and spatial resolution have been achieved in experiments and numerical studies to allow the analysis of these phenomena. Frequently new numerical and experimental data has been in contradiction with theoretical models, and previous knowledge on fluid-particle interactions had to be reconsidered even in simplified cases [3]. Furthermore, there are many open questions concerning inhomogeneous flows [4], finite-size [5,6] and non-spherical particles [7], among others. More recently, colloidal particles were used in the first experiments to verify fluctuation relations such as the Jarzynski equality [8], which links the statistics of fluctuating quantities in a non-equilibrium process with equilibrium quantities. Colloidal particles were also used to verify the thermodynamic cost of information processing [9] proposed by Landauer. Fluctuation theorems for more general p-1 problems, in which a physical system is coupled with (and driven by) another out-of-equilibrium system, would open potential applications in other complex coupled systems as in soft and active matter. But what fundamental statistical relations are satisfied by particles that interact with a complex and out-of-equilibrium turbulent flow? Brownian motion set the path for the study of diffusion and random processes. It becomes more complicated if the particle motion is driven by turbulence; even a simple point-wise passive particle provides a challenge in the modeling of many aspects of its dynamics [10,11].
Friedrich-Peinke approach to turbulence.-In this letter we show that the Markov property is valid for the dynamics of dense sub-Kolmogorov particles coupled to a turbulent velocity field of three-dimensional homogeneous and isotropic turbulence (HIT) at inertial range time scales. To this end we use pseudo-spectral direct numerical simulations (DNSs) solving the incompressible Navier-Stokes equation with a simple point particle model [12]. To characterize HIT it is customary to study the statistics of Eulerian velocity increments u r = [u(x + r) − u(x)] · r/|r| for scales r = |r|. However, two-point statistics of velocity increments do not fully characterize small-scale turbulence [13], and many attempts at dealing with multi-point statistics have been considered. One way to do this is to use the Friedrich-Peinke approach [14], in which the stochastic dynamics of velocity increments u r are considered as they go through the cascade from large to small scales r. A central assumption is that the evolution of the stochastic variable u r possesses a Markov process "evolving" in r. Previous studies showed that u r can be considered as Markovian to a reasonable approximation [14][15][16], at least down to a scale close to the Taylor scale [16,17]. Furthermore, u r satisfies a diffusion process [18] −∂ r u r = D (1) (u r , r) + [D (2) where a linear dependence on the value of the increment for the drift coefficient D (1) (u r , r), and a quadratic dependence for the diffusion coefficient D (2) (u r , r), was found [13,16,[19][20][21]. Further details on the noise term Γ(r) and the empirical estimation of the drift and diffusion coefficients D (1,2) (u r , r) are given in [22] (Appendix A) .
Up to now, this approach has been applied to an Eulerian description of turbulent flows. Based on the general interest in the relation between Eulerian and Lagrangian properties of turbulence the question of the application of this approach and, in particular, the verification of the necessary condition of Markov property for Lagrangian turbulence and inertial particles arises.
DNS of inertial particles in turbulence.-To this end we use DNSs of forced HIT. A large-scale external mechanical forcing is given by a superposition of modes with slowly evolving random phases, following standard practices for its temporal integration and de-aliasing procedures. An adequate spatial resolution of the smallest scales, i.e., κη 1 is chosen [23]. Here, η is the Kolmogorov or dissipation length scale, η = (ν 3 /ε) 1/4 (where ε is the kinetic energy dissipation rate, and ν the kinematic viscosity of the fluid), and κ = N DNS /3 is the maximum resolved wavenumber in Fourier space (with N DNS = 512 the linear spatial resolution in each direction). The DNSs are characterized by a Reynolds number based on the Taylor microscale of Re λ ≈ 240 (see [2] for details). Inertial particles were modeled using the Maxey-Riley-Gatignol equation in the limit of point heavy particles, which for a particle with velocity v in the position x p submerged in a flow with velocity u, readṡ with τ p the particle Stokes time. For tracers, τ p = 0.
There is no particle-particle or particle-fluid interaction (i.e., we use a one-way coupling approximation). The distinction between particles is quantified by the Stokes number St = τ p /τ η (the ratio of the particle relaxation time to the Kolmogorov time, defined as τ η = [ν/ǫ] 1/2 ).
Cascade trajectories.-To study the Markov property of particles trajectories, we study v(t) (the Lagrangian velocity of the particle at time t) for St > 0, and for Lagrangian tracers with St = 0 the analysis is performed using v(t) = u(x p , t) (the velocity of the fluid element where the tracer is located). Building now velocity increments in time instead of the commonly used spatial scales r (e i corresponds to the unitary vector in the direction x, y, or z specified in the figures) we define for every particle trajectory a "cascade trajectory" (or "sequence" of velocity increments at decreasing scales) [u(·)] = {u T , . . . , u τ f } for different time-separations or time scales τ , from the initial T to the final time scale τ f with T > τ f . The notation [u(·)] indicates the entire path through the hierarchy of time scales from large to small scales instead of a distinct value u τ ; T ≈ 120τ η and τ f ≈ 0.2τ η were chosen as references.
The advantage of a Lagrangian perspective compared to a spatial (Eulerian) perspective is that Lagrangian particle trajectories have a clear physical meaning: each particle follows a path according to the fluid dynamics and its own inertia. In other words, the time increments have a direct interpretation as values corresponding to changes in the velocity along trajectories as the system evolves. In the Eulerian case, the spatial distances (or the "evolution" in r) have an indirect interpretation in the context of a Fokker-Planck formulation of the dynamics.
Results.-To verify the existence of the Markov property for Lagrangian turbulence, we integrate three sets of 1.3 × 10 5 inertial particles, respectively with St = 0, 0.5, 3, and 6.
In Fig. 1 a qualitative validation of the Markov property based on the alignment of the single conditioned p (u τ2 |u τ1 ) and double conditioned p (u τ2 |u τ1 , u τ0 ) probability density functions (PDFs) of datasets of increments for a chosen set of three different time scales τ 0 > τ 1 > τ 2 is shown as contour plots. Each scale is separated by ∆τ ≈ 13τ η , which is equal to the Einstein-Markov coherence time scale ∆ EM for St = 3. As described below and in the discussion accompanying Fig. 2, ∆ EM is the smallest time-scale for which the Markovian assumption can be considered valid. For Markovian processes the relation p (u τ2 |u τ1 ) = p (u τ2 |..., u τ1 , u τ0 ) holds. For finite datasets, is commonly assumed to be a sufficient condition. The close alignment between the PDFs (black and red solid lines in Fig. 1) confirms the validity of the Markov property. Note that while Fig. 1 shows results for St = 3, all our datasets give similar results. In Fig. 2 the Markovian approximation is checked systematically for different values of ∆τ , using the Wilcoxon test, which is a quantitative and parameter-free test that determines the Einstein-Markov coherence time scale ∆ EM (see [16,17] and [22], Appendix G, for details and the description of essential notions). This test is a reliable procedure to validate Eq. (4) [16,17]. Based on this analysis, the Markovian approximation is valid for time-scales larger than or equal to this St-dependent critical time separation ∆ EM ≈ 7, 10, 13 and 16τ η for St = 0, 0.5, 3, and 6 respectively (see Fig. 2(a)). Accordingly, the complexity of the dynamics of inertial particles in turbulent flows, expressed by the evolution of the stochastic variables u τ at decreasing time scales, can be treated as a Markov process, with an Einstein-Markov coherence time scale that increases with St. Our analysis presented in [22] (Appendix C and D) demonstrates that on the scale that the Markov property holds the non-Gaussian behavior is still present and that our results are well linked to the range containing Lagrangian small-scale intermittency. Furthermore, we see that ∆ EM is the same (at fixed St) for all velocity components (see Fig. 2(b)). We performed a similar analysis for the fluid velocity at the particle position u(x p , t) shown in [22] (Appendix F). Interestingly, the Einstein-Markov time-scale decreases with increasing St, unlike the behavior observed for v(t). While this remarkable behavior and the effect this has on the fluid velocity observed by the particles requires further study, it may be related to the fact that a larger St number implies a larger drag, and therefore a tendency of the system to dissipate energy more efficiently and faster. This would result in a reduction of the associated Einstein-Markov time of u(x p , t). . Note, that the normalized expectation value t(τ, ∆τ ) is supposed to be close to 1 if the Markovian assumption is valid (see [16,17] and [22] (Appendix G) for details and the description of essential notions).
The results shown above can be linked to stochastic thermodynamics [20,24] (also called stochastic energetics [25]), developed for many different physical systems [26][27][28][29], and to integral and detailed fluctuation theorems. Experimental studies of stochastic thermodynamics mainly focus on nanoscale or quantum systems, or when dealing with classical systems, on biological systems [30,31], which are assumed to be well off the thermodynamic limit so that the probabilistic nature of balance relations becomes clearer. We test these theorems starting with the integral case as it should be more robust.
For the case of passive particles, it is the turbulent surroundings that impose the non-equilibrium condition. Turbulence can be understood, and especially its energy cascade, as a process leading a fluid under specific conditions and parameters (the Reynolds number) from a nonequilibrium state into another non-equilibrium state, by the combination of energy injection and dissipation. In the spirit of non-equilibrium stochastic thermodynamics [32], by using the Fokker-Planck equation, it is possible to assop-3 ciate with every individual trajectory [u(·)] a total entropy variation [20,24,25,29,32] given by the sum of two terms A thermodynamic interpretation of this quantity based on the relation between heat, work, and inner energy, can be given [25,29,32]. ∆S sys is the change in entropy associated with the change in state of the system for a single realization of a process that takes some initial probability p (u T , T ) and changes to a different final probability p u τ f , τ f . It is simply the logarithmic ratio of the probabilities of the stochastic process, so that its variation is The other term, the entropy exchanged with the surrounding medium ∆S med from the initial to the final time scale, measures the irreversibility of the trajectories: Then, the dynamics of inertial particles coupled to a turbulent flow are characterized by the empirically estimated drift and diffusion coefficients D (1,2) (see [22] for a presentation of D (1,2) for all Stokes numbers considered). The knowledge of drift and diffusion coefficient allows to determine for each [u(·)] an entropy value, for which the thermodynamic fluctuations can be investigated. The dashed line in Fig. 3(a)  . We find that for the four Stokes numbers the results are in good agreement with the IFT, which is a fundamental entropy law of non-equilibrium stochastic thermodynamics [29,32]. This is not trivial, as the IFT is extremely sensitive to errors in the empirical estimation of D (1,2) defining the Fokker-Planck equation. For this reason, our empirical observation on the validity of the IFT for the inertial particles provides another verification of our approach. Furthermore, inertial particles (i.e., for sufficiently large St) sample preferentially some specific regions of the flow [2]. The remarkable observation here is that they do it always in agreement with this theorem. Figure 3(b) shows that in addition to the IFT also the detailed fluctuation theorem (DFT) holds for our data coming from a turbulent flow with finite Reynols number, which is expressed as Thus, in addition to the IFT, the DFT expresses the balance (explicit exponential symmetry constraint) between entropy-consuming (∆S tot < 0) and entropy-producing (∆S tot > 0) trajectories. See [22] for the illustration of p (∆S tot ) for St = 0.5, 3 and 6. Finally we ask whether these findings of the entropy fluctuation are just a statistical result or if they have some relation to flow structures. In particular, it is of interest if the negative entropy events have some special features. Therefore we study the mean absolute velocity increment trajectories conditioned on a specific total entropy variation |u τ | ∆Stot . Figures 4(a) and (c) show that the total entropy variation is linked to distinct trajectories. Entropy-consuming trajectories are characterized by an increase in the averaged absolute values of the increments |u τ | ∆Stot=−2 with decreasing time scale τ , while trajectories marked by entropy-production |u τ | ∆Stot=3 smoothly decrease their absolute increments with decreasing τ . This is further highlighted by Fig. 4(b) and (d). In these three-dimensional representations the absolute velocity increment at initial and final scale, T and ∆ EM , conditioned on the entropy for all trajectories for which the conditional averaging was performed, are shown. In [22] (Appendix E) we show that the majority of large and intermittent small scale increments (heavy tailed statistics of the small-scale increment PDF) are contained in the negative entropy trajectories.
Conclusion.-In this letter we showed the existence of Markov property for Lagrangian turbulence. The dynamics of inertial particles coupled to a turbulent flow can be treated as a Markov process for a finite step size larger than a Stokes number-dependent Einstein-Markov coherence time scale. This opens an interesting interpretation of the St number in terms of the Markov memory of the particles trajectories, and is compatible with the picture that particles with more inertia filter fast and small scale fluctuations of the carrying flow: for particles with larger St (and thus larger particle response times) the Markovianization of trajectories by the turbulence takes place at longer times. This can be used to quantify effective parameters in cases where the particles' inertia is not clearly defined, like finite-size, or non-spherical particles.
Moreover, based on the theory of Markov processes it is possible to derive an explicit Fokker-Planck equation for the description of multi-scale statistics of Lagrangian particles. The Markov property corresponds to a threepoint (two-scale) closure of the general joint multi-scale statistics [33], representing a step towards a more comprehensive understanding of turbulence, which is characterized by a challenging complexity with many high order correlations. This results in a tremendous reduction of the degrees of freedom, and shows that the long-standing problem in turbulence of the closure of the hierarchy of equations can be approached in a Lagrangian framework.
The interpretation of the stochastic process as an analogue of a non-equilibrium thermodynamic process [20,24,33] when applied to Lagrangian turbulence allowed us to show that inertial particles satisfy both in-tegral and detailed fluctuation theorems in a strict sense. The verification of these theorems, which have contributed significantly to the understanding of open and out-ofequilibrium systems and to the development of statistical physics, adds strength to the validity of the Markov property and the description by a Fokker-Planck equation. Previous observations of the irreversibility of particles in turbulence considered a Jarzynski-like equality for the particles energetics [34], but the slope was not one as expected for the DFT. Reproductions of such exact results are found very rarely in turbulence research, so that the result provides a novel constraint for the particles' evolution. The fluctuation theorems may be taken also as a kind of new conservations laws.
These fluctuation theorems also shed light on asymmetries in particle trajectories, as reported before in [34], and as shown here for individual trajectories depending on their entropy evolution. A connection can be made between entropy-producing and -consuming trajectories with specific evolution of velocity increments with decreasing time scale. Interestingly, our study also shows that particle trajectories are out-of-equilibrium while keeping ergodicity. These abnormal behaviors (or structures) in the Lagrangian paths are statistically embedded in the fulfillment of the two fluctuation theorems. Thus, the presented results may shed new light on the long-lasting debate of the statistical approach to turbulence versus the study of flow structures.
The case of out-of-equilibrium systems like turbulence coupled to dense particles has applications in geophysics, astrophysics, and aerosol dispersion. In the cases studied here we have non-trivially coupled systems, whose particles not always randomly sample the flow topology [2]. Note that preferential concentration for St = 0 implies that particles are not sampling the flow homogeneously, and a non-random sub-sampling of the Eulerian flow velocity, plus the particles own inertia, could break down the IFT and DFT theorems. In spite of this, the theorems are satisfied for particles with different response times (albeit, in the limit St → ∞, the relations should not hold as particles decouple from the flow). Also, the formalism proposed here only requires measuring the trajectory and velocity of particles. Such study is accessible with stateof-the-art experimental techniques, even at large Reynolds numbers, while most models to describe the dynamics of these systems require much more detail, such as measurements also of particles' acceleration. At last, we point out that the entropy analysis may also serve to compare different turbulent situations, as it was successfully done by the statistics p(∆S tot ) for turbulent wave situations [35]. * * * This work has been partially supported by the ECOS project A18ST04, by the Volkswagen Foundation (96528) and by the Laboratoire d'Excellence LANEF in Grenoble (ANR-10-LABX-51-01). p-5