How tracer particles sample the complexity of turbulence

On their roller coaster ride through turbulence, tracer particles sample the fluctuations of the underlying fields in space and time. Quantitatively relating particle and field statistics remains a fundamental challenge in a large variety of turbulent flows. We quantify how tracer particles sample turbulence by expressing their temporal velocity fluctuations in terms of an effective probabilistic sampling of spatial velocity field fluctuations. To corroborate our theory, we investigate an extensive suite of direct numerical simulations of hydrodynamic turbulence covering a Taylor-scale Reynolds number range from 150 to 430. Our approach allows the assessment of particle statistics from the knowledge of flow field statistics only, therefore opening avenues to a new generation of models for transport in complex flows.


Introduction
Turbulence, as frequently encountered in our atmosphere and oceans as well as in a plethora of astrophysical phenomena, exhibits a stunning complexity in space and time. This complexity can be perceived from two complementary perspectives: from the Eulerian point of view turbulence is observed in the laboratory frame by probing the flow at one or more fixed spatial points; in the Lagrangian frame, the flow is probed by tracer particles which follow the velocity field (see movie SM1 which is available online at http://stacks.iop.org/NJP/20/013001/mmedia). This perspective is of particular importance for characterizing turbulent mixing and transport in particle-laden flows.
One of the most challenging phenomena of turbulence is the breaking of statistical self-similarity: probing velocity fluctuations across scales reveals a continuous change in shape of the corresponding distributions from the largest, energy-containing scales of the flow down to the smallest, dissipative scales, on which extreme events occur with astonishing frequency [1,2]. This phenomenon, known as intermittency, can be cleanly isolated in homogeneous isotropic turbulence, which shares universal small-scale properties with a wide class of hydrodynamic flows even at moderate Reynolds numbers [3]. Experimental [4][5][6][7][8] and numerical [9][10][11][12][13][14][15][16] studies of isotropic turbulence have traced the signatures of Eulerian statistics back to the intricate spatial structure of turbulence.
As tracer particles sample the fine-scale structure of turbulence, they encounter violent velocity and acceleration fluctuations on the smallest temporal scales [17][18][19][20]. In analogy to the Eulerian observations, Lagrangian velocity increment distributions are also nonself-similar [17] when probing the velocity fluctuations across temporal scales. However somewhat puzzlingly, extreme events are even more frequent in the Lagrangian than in the Eulerian frame.
This leads to the long-standing riddle of relating the observations in the two complementary frames, and in particular to the question: how do tracer particles sample the spatio-temporal velocity fluctuations of turbulence? Here we provide an answer by combining an exact probabilistic framework first introduced in [21,22] with the physics of Lagrangian particle transport and Eulerian temporal decorrelation into a predictive theory relating Lagrangian to Eulerian statistics. We show that, remarkably, temporal velocity fluctuations along Lagrangian tracers can be statistically predicted from instantaneous spatial velocity fluctuations by properly mixing Eulerian statistics of various scales. We demonstrate that this probabilistic mixing can be captured in terms of the probability density function (PDF) of an effective Lagrangian dispersion process, whose properties we determine. The investigation of an extensive suite of direct numerical simulations of homogeneous isotropic turbulence furthermore reveals its universality with respect to a range of Reynolds numbers.

Bridging Eulerian and Lagrangian statistics
The starting point for the development of our theory is that the Lagrangian velocity increment v = u (X(x 0 , t 0 + τ ), t 0 + τ ) − u(x 0 , t 0 ) can be perceived as an Eulerian twopoint-two-time velocity increment between two space-time points which are connected by a Lagrangian trajectory X(x 0 , t 0 + τ ) starting from x 0 at t 0 . The Lagrangian increment statistics can therefore be related to the Eulerian statistics in an exact manner (see Appendix A.1 for a detailed derivation): Here, f L (v; τ ) is the PDF of the Lagrangian velocity increment, f E (v; R, τ ) is the Eulerian two-point-two-time increment PDF, and P (R|v, τ ) is the conditional Lagrangian dispersion PDF of the particle displacement R = X(x 0 , t 0 + τ ) − x 0 , which describes how particles spread over time lag τ depending on the velocity increment v.
To establish a bridging relation between temporal Lagrangian and instantaneous spatial Eulerian increment statistics, we need to develop an understanding of the Eulerian temporal decorrelation (contained in the Eulerian two-point-two-time increment PDF) as well as of the Lagrangian single-particle dispersion. For long times, Lagrangian tracer particles have typically travelled far from their point of origin. As the particles are advected with the local instantaneous velocity, all correlations with the velocity increment taken over the travelled distance have decayed, an assumption consistent with Corrsin's hypothesis [23]. As a consequence Figure 1. Dispersion of Lagrangian trajectories with fixed initial velocities (twice the rms velocity u) from a direct numerical simulation at R λ = 316. For this visualization, particles with identical initial velocity magnitudes have been selected from the flow domain and shifted to a common origin and orientation. (a): Initially, their motion is almost ballistic and remains close to the mean Lagrangian trajectory shown in magenta. (b)-(d): Over time, the particles disperse and randomly sample the region surrounding the mean trajectory. The black curve shows the ballistic approximation to the mean Lagrangian trajectory. The distance between the mean Lagrangian trajectory and the individual trajectories (as indicated for one sample trajectory by the arrow) is the effective Lagrangian displacement. (See movie SM2 for an animated version of this figure which can be found online at http://stacks.iop.org/NJP/00/000000/ mmedia.) P (R|v, τ ) becomes independent of v and simplifies to the unconditional standard Lagrangian dispersion PDF. The Lagrangian dispersion PDF is close to Gaussian in homogeneous isotropic turbulence, and therefore fully characterized by its mean squared displacement, which grows diffusively for long times, i.e. R 2 ∼ τ . While the particles disperse, the flow also decorrelates at their point of origin. The Eulerian temporal decorrelation is primarily governed by the random advection of velocity fluctuations by larger-scale eddies past the fixed Eulerian position, so-called random sweeping [24,25]. Based on this hypothesis the Eulerian two-point-two-time increment PDF is obtained by blurring out the instantaneous Eulerian increment PDF with a Gaussian "random sweeping decorrelation PDF" (see Appendix A.2 for a detailed derivation), whose variance also grows diffusively for long times. The long-time behavior therefore can be captured by combining the two classical ideas of Corrsin's statistical independence hypothesis for the Lagrangian particle dispersion and the Tennekes-Kraichnan random sweeping hypothesis for the Eulerian decorrelation in one effective PDF; we obtain is close to Gaussian with a variance of R 2 ∼ τ containing contributions from both Lagrangian single-particle dispersion as well as random sweeping decorrelation. In comparison to (1) this relation features Eulerian single-time information only and shows that particle dispersion and random sweeping introduce a probabilistic sampling of Eulerian velocity fluctuations over a range of spatial scales.
For short times, Corrsin's approximation as well as the assumption that random sweeping decorrelation and particle dispersion are independent are clearly violated. In fact, the latter two effects are induced both by the same advection velocity and happen therefore coherently. In the following, this will be accounted for by removing spurious sweeping effects in the spirit of Kolmogorov's original works [26,27]. Furthermore, statistical correlations between particle dispersion and velocity fluctuations play a crucial role at short times.
To determine the short-time behavior, we consider Lagrangian trajectories which have a common initial velocity, visualized in figure 1. Although the particles start with the same initial velocity, the individual trajectories deviate very soon. This virtual blob rapidly deforms and spreads across the volume due to turbulent mixing. The center of mass of this blob of particles, however, drifts in very good approximation ballistically for short times. To see that, we introduce the mean Lagrangian displacement X(x 0 , t)− x 0 = X(x 0 , t) − x 0 |u 0 where the average is performed over all Lagrangian particles with fixed initial velocity u 0 . For short times, the mean Lagrangian displacement can be expanded in a Taylor series: Here, a 0 |u 0 denotes the average acceleration at the initial Lagrangian point conditional on the initial velocity. Because the average is taken over many different flow configurations with fixed initial velocity, we obtain a 0 |u 0 ≈ 0 due to a weak directional correlation between the large-scale velocity and the small-scale acceleration field. The same argument applies to the higher-order terms. Therefore, the short-time behavior is given by a ballistic flight. In fact, this Taylor expansion is expected to hold for much longer than the Taylor expansion of an individual Lagrangian trajectory -there higherorder terms cannot be neglected even on the order of the Kolmogorov time scale, which characterizes the small time scales of the flow. As the mean Lagrangian trajectory X(x 0 , t) moves ballistically, it approximately maintains its initial velocity. As a result the Lagrangian velocity increment at small time lags can be expressed as v ≈ u (X(x 0 , t 0 + τ ), t 0 + τ )−u(X(x 0 , t 0 +τ ), t 0 +τ ) such that velocity increment is now taken at the same time t 0 + τ over the "effective Lagrangian displacement" R = X(x 0 , t 0 + τ ) − X(x 0 , t 0 + τ ). The instantaneous Eulerian field is therefore sampled on a scale smaller than naïvely expected based on single-particle dispersion. As a result we obtain a bridging relation very similar to (1) (see Appendix A.3 for a detailed derivation): This relation again features only instantaneous Eulerian information, now mixed by an effective Lagrangian dispersion process which removes sweeping effects. Comparing the two bridging relations for short and long times, we arrive at the remarkable conclusion that they can be unified by identifying P (R|v, τ ) as the "effective Lagrangian dispersion PDF" that combines the dispersion of tracer particles for short times and the combined effects of particle dispersion and random sweeping for long times. We have already established its shape for large time lags as Gaussian with a diffusively increasing variance. For short times, the effective Lagrangian dispersion PDF in comparison has to be sharply localized. Here, the effective Lagrangian displacement is in leading order given by the initial acceleration, X(x 0 , t 0 +τ )−X(x 0 , t 0 +τ ) ≈ a 0 τ 2 /2. For such small time lags, the acceleration can be expressed in terms of the velocity increment, a 0 ≈ v/τ . This exact result, also discussed in [28], is reminiscent of the phenomenological bridging relation used in the multifractal framework [29][30][31][32]. The two limiting cases (short and long time behavior) can be unified into a model for the effective Lagrangian dispersion PDF taking the form: where for short times σ(τ ) approaches zero such that the Gaussian asymptotically converges to a delta function while µ(τ ) behaves like τ /2. For long times, however, µ(τ ) goes to zero as the correlations with the velocity increment decay, whereas σ(τ ) reaches a diffusive regime.

Results from direct numerical simulations
To test our theory, we ran an extensive suite of direct numerical simulations (DNS) of homogeneous isotropic turbulence covering a Reynolds number range from R λ ≈ 150  to R λ ≈ 430. The data was generated with a standard pseudo-spectral code solving the Navier-Stokes equation in the vorticity formulation. DNS setup and algorithms are described in more detail in Appendix A.5.
To further simplify our theoretical predictions we make use of isotropy. As one implication, the full statistical information of the Lagrangian velocity increment vector is contained in the velocity increment magnitude, and the corresponding bridging relation is readily derived (see Appendix A.4 for details): For P (R|v, τ ) we take the angle-integrated version of (4), which takes the form where sinhc(x) = sinh(x)/x. The angle-integrated effective Lagrangian dispersion PDF simplifies to a Maxwellian for µ(τ ) = 0, but otherwise displays a non-trivial dependence on the velocity increment. Because the proposed model PDF mixes different Eulerian scales for different velocity increments, we also have to include the factor N (τ ) to keep the model Lagrangian PDFs normalized. For the data sets under consideration N (τ ) is of order unity, as expected. To compare our model to the DNS data we combine (6) with (5) and determine µ(τ ) and σ(τ ) by a least square fitting of the structure functions (moments of the velocity increment PDF) up to order eight. The detailed fitting procedure is explained in Appendix A.6. The results are shown in figure 2. The width of the effective Lagrangian dispersion PDF is controlled by σ(τ ), which is shown in panels (a) and (b), non-dimensionalized by Kolmogorov units and integral scale units, respectively. Remarkably, non-dimensionalization with the Kolmogorov units leads to the collapse of the curves obtained from various Reynolds number data sets at small scales, whereas non-dimensionalization with integral scales leads to a collapse at large scales. Closer inspection of the results suggests that σ(τ ) first increases linearly and then continues to grow approximately like τ 3/2 . The effective Lagrangian dispersion takes out sweeping effects and describes the two-particle separation between real and the virtual mean Lagrangian trajectories. It is therefore tempting to assume that the τ 3/2growth is due to a Richardson dispersion-like behavior. On theoretical grounds, σ(τ ) is expected to reach the long-time limit ∼ uT 1/2 τ 1/2 , which is in reasonable agreement with σ(τ ) obtained from the DNS data. Deviations at the largest time lags are likely rooted in the fact that the fitting procedure becomes less sensitive to changes in σ(τ ) as the increment PDFs reach their asymptotic close-to-Gaussian shape. The cross-over between the regimes depends on the Reynolds number. In fact, we find that this behavior can be captured for all Reynolds numbers in the empirical formula Interestingly, our investigations identify the Eulerian integral time scale as the one controlling the transition at large scales. This could be rooted in the fact that the longtime limit is dominated by random sweeping decorrelation along with the expectation that the Eulerian and Lagrangian integral time scale are roughly of the same order [33].
The shift of the effective Lagrangian dispersion PDF is controlled by µ(τ ), which is shown in panel (c) of figure 2. For short times the estimate obtained from the data is in good agreement with our analytical prediction of µ(τ ) ≈ τ /2, whereas µ(τ ) drops  on the order of the integral time scale. When scaled by the integral time scale, the fits for various Reynolds numbers collapse. For the Reynolds numbers under consideration, we find that captures the numerical observations well. It has to be mentioned though that the exponential decay is delicate to estimate from the data: as µ becomes less significant, its accurate estimation becomes very difficult. Figure 3 shows the effective Lagrangian dispersion PDF (6) for various time lags τ with σ(τ ) and µ(τ ) given by (7) and (8), respectively. For small to moderate time lags, the PDF displays a pronounced dependence on the velocity increment, as shown in panels (a) and (b): for larger velocity increments, relatively larger Eulerian scales with markedly larger variances are mixed into the Lagrangian statistics. This provides a mechanism to generate heavier tails of the Lagrangian PDFs compared to the Eulerian. This effect weakens with increasing time lag. On the order of the integral time scale, Corrsin's approximation begins to hold and the effective Lagrangian dispersion PDF becomes independent of the velocity increment; as a result it simplifies to a standard Maxwellian, as demonstrated in panel (c). In essence, the effect of the bridging relation is two-fold: Eulerian statistics of various scales are mixed to obtain Lagrangian statistics. Additionally, the scale-mixing depends on the magnitude of the velocity fluctuations. This allows for a different statistical behavior of higher-order statistics, dominated by the tails of the PDF, compared to lower orders, dominated by its core.
To test the fidelity of the bridging relation, we obtain Eulerian velocity increment magnitude PDFs from the DNS and compute Lagrangian increment magnitude PDFs using the angle-integrated version of the effective Lagrangian dispersion PDF along with the formulas for σ(τ ) and µ(τ ). As it is more common to consider PDFs of the velocity   increment components rather than the PDFs of magnitudes, we compute the standard component velocity increment PDFs from

Conclusions and Outlook
In summary, we have introduced a novel bridging relation that relates Lagrangian temporal velocity fluctuations to instantaneous Eulerian spatial velocity fluctuations. We showed that Lagrangian increment statistics can be perceived as a scale-mixing of Eulerian increment statistics by an effective Lagrangian dispersion. Based on the physics of single-particle transport as well as random sweeping effects, we derived the effective Lagrangian dispersion PDF and determined its asymptotic behavior. The predictions for the Lagrangian velocity increment PDFs were found to be in excellent agreement with results from high-resolution direct numerical simulations of fully developed turbulence. Supported by the extensive numerical investigations, we conclude that for the Reynolds number range under consideration Lagrangian intermittency can be captured with the effective Lagrangian dispersion PDF, whose shape is independent of Reynolds number. The Reynolds number dependence enters through the Eulerian statistics and through the width σ(τ ). The fact that the effective Lagrangian dispersion mixes Eulerian statistics across a wide range of scales indicates that Lagrangian velocity increment statistics typically contain both dissipative as well as integral-scale contributions (see also the related discussion in [21]). This may provide an explanation of a reduced scaling range of Lagrangian structure functions compared to their Eulerian counterparts; for instance, it was argued that even for the second-order (i.e. the lowest non-trivial order) Lagrangian structure function a clear scaling range is only expected starting at R λ ≈ 5 · 10 3 [34] or even beyond R λ ≈ 3 · 10 4 [35], where Eulerian scaling is already well established.
The presented approach is very general and therefore opens avenues to a new generation of models for transport in complex particle-laden flows by allowing an assessment of particle statistics from the knowledge of the flow field statistics only. For instance, an exciting future direction would be the extension of the theory to inertial and finite-sized particles. Here, the effect of preferential concentration will imply a biased sampling of the flow leading to modifications of the effective Lagrangian dispersion PDF. It also appears worthwhile to pursue applications to atmospheric flows and turbulent convection by including characteristic large-scale features of these flows. Finally, it will be most interesting to go beyond hydrodynamic flows and investigate the bridging induced by the effective Lagrangian dispersion PDF, for example, in magnetohydrodynamic flows, which display a fine-scale structure markedly different from hydrodynamic turbulence [22,36].
Planck Computing and Data Facility. This work was supported by the Max Planck Society.

Appendix A.
Appendix A.1. Theoretical background: the bridging relation To establish the probabilistic laws connecting Eulerian and Lagrangian statistics, we consider a velocity field u(x, t) of a homogeneous isotropic turbulent flow governed by the incompressible Navier-Stokes equation in the statistically stationary state. A very small (compared to the smallest scales of the flow) neutrally buoyant Lagrangian tracer particle starting from point x 0 at time t 0 is advected by the velocity field according to with the initial condition X(x 0 , t 0 ) = x 0 . As the particle probes the flow over the time lag τ , it encounters velocity fluctuations which can be characterized by the Lagrangian velocity increment u (X(x 0 , t 0 + τ ), t 0 + τ ) − u(x 0 , t 0 ). The statistics of the velocity fluctuations can be conveniently captured in terms of the probability density function (PDF) Here, . . . denotes an ensemble average over all Lagrangian particles and all velocity fields that can also be replaced with a spatio-temporal average for practical purposes. The PDF f L (v; τ ) is a normalized density with respect to the sample-space variable v with a functional dependence on the time lag τ . Note that it does not explicitly depend on the initial time t 0 due to statistical stationarity. Eq. (A.2) essentially contains statistical information of a Lagrangian tracer particle sampling the Eulerian, fixed labframe velocity fluctuations at two points in space and time. In fact, this statement can be made rigorous by considering stationarity have been exploited for the Eulerian increment PDF. As a result, it only depends on the time lag τ and the spatial difference y − x 0 . Also the starting point x 0 has been added and subtracted in the argument of the delta function of the Lagrangian particle position. From (A.6) to (A.7) a change of variables was applied, which finally allows us to identify the conditional Lagrangian dispersion PDF in (A.8). This concludes the derivation of (1) from the main text. The above bridging relation still contains Eulerian two-time information. The goal in the following is to establish a relation between f L (v; τ ) and the Eulerian single-time velocity increment PDF Due to homogeneity, this PDF does not explicitly depend on the spatial coordinate x, and, as before, it does not depend explicitly on the time t due to statistical stationarity. The main challenge here is to relate the Eulerian two-time statistics to the single-time statistics. As outlined in the main text, the temporal effects can be absorbed into the Lagrangian dispersion PDF by introducing an effective dispersion. The physical arguments therefor are given in the main text, here we add some technical details.

Appendix A.2. Bridging relation for long times
We show here that the long-time behavior is captured by the proposed bridging relation. For very long times we obtain Here, we have introduced the standard Lagrangian dispersion PDF f LD (R; τ ) = δ[X(x 0 , t 0 + τ ) − x 0 − R] and the Eulerian two-point-two-time increment PDF f E (v; R, τ ) and assumed that the two PDFs become statistically independent for long times, which is essentially Corrsin's hypothesis and well justified by the physics of the problem. The Eulerian two-point-two-time increment PDF can be related to the single-time increment PDF by applying the random sweeping hypothesis [24,25], i.e. we assume that the temporal decorrelation is dominated by large-scale random advection effects. That means, velocities are carried past a fixed point in space with a random sweeping velocity over a distance R sw (τ ) = τ 0 ds u sw (s), where u sw is the large-scale random-sweeping velocity at the fixed Eulerian position. This allows us to map the velocity back in time according to For simplicity we assume that u sw is an isotropic Gaussian process with zero mean and correlation function u sw (t) · u sw (t ) /3. As a consequence, the random sweeping displacement is also a Gaussian process with ds u sw (s) · u sw (s ) . (A.14) As most of the kinetic energy is contained in the large scales of the flow, we have u 2 sw /3 ≈ u 2 /3 for the single-time variance. For short times, the mean squared displacement grows ballistically, R 2 sw (τ ) ≈ u 2 τ 2 . For long times, significantly larger than the Eulerian integral time scale T , the displacement becomes diffusive, i.e. R 2 sw (τ ) ≈ 2 u 2 T τ . Taking into account these considerations, we can relate the two-time PDF to the single-time PDF via where f sw is the Gaussian random sweeping PDF with variance R 2 sw (τ )/3. In the derivation we have made the assumption of statistical independence of the sweeping velocity and the velocity increment. Combining this result with the long-time bridging relation, we obtain The convolution of the Lagrangian dispersion PDF and the random sweeping PDF yields the effective Lagrangian dispersion PDF, showing that the structure of the bridging relation holds for long times. Both, the random sweeping PDF as well as the standard Lagrangian dispersion PDF are in good approximation Gaussian with diffusively growing variances for large times. Therefore also the effective Lagrangian dispersion PDF is Gaussian with a variance ∼ τ for long times.

Appendix A.3. Bridging relation for short times
To make contact to the Eulerian single-time statistics for short times, the mean Lagrangian displacement X(x 0 , t) − x 0 = X(x 0 , t)|u 0 − x 0 has been introduced in the main text. Because Lagrangian particle transport is dominated by random advection by the large-scale velocity, the velocity is expected to vary only little along the mean Lagrangian trajectory and hence u(X(x 0 , t 0 + τ ), t 0 + τ ) = u(x 0 , t 0 ) + u , where the correction u is negligible for small to moderate times. This allows us to approximately identify the initial velocity of the Lagrangian increment with the velocity along the mean Lagrangian trajectory. As a result the bridging relation becomes for τ not too large. The crucial point is that this expression involves two points at a single time only. The remaining steps to the complete bridging relation are now straightforward: The main difference is that we have exploited the fact that f E (v; R) depends only on the distance magnitude due to isotropy. As a result, we can integrate over the angles from (A.35) to (A.36). In (A.36) we have thereby introduced the effective Lagrangian dispersion PDF of the magnitude R. If a model PDF is introduced for P (R|v, τ ), such as the one discussed in the main text, the additional factor N (τ ) ensures normalization. The corresponding long-time bridging relation can be straightforwardly obtained in the same manner. This concludes the derivation of the bridging relation (5) from the main text.
Appendix A.5. Description of the DNS For this study, an extensive suite of direct numerical simulations of statistically stationary fully developed turbulence in a periodic domain was conducted. A standard MPI-parallel pseudo-spectral scheme is used to solve the vorticity formulation of the Navier-Stokes equations in a periodic cubic box of length 2π (DNS units). Aliasing errors are controlled by using a high-order Fourier smoothing [37]. Time stepping is performed by means of a memory-saving third-order Runge-Kutta method [38]. To maintain the statistically stationary state, a large-scale band-passed Lundgren forcing [39] is applied in the wavenumber range [1.5, 3] (DNS units) with an amplitude of 0.5 (DNS units).
To obtain the Lagrangian data, tracer particles are advected with the flow. Cubic splines are used to interpolate the velocity fields. Time stepping for the particles is performed by means of a fourth-order Adams-Bashforth method (see, e.g., §6.7 in [40]).
As preparation, the Eulerian flow field is advanced in time until a statistically stationary state is obtained. The production runs are launched using this state as an initial condition, and Lagrangian tracers initialized at random locations in the cubic domain (with a uniform distribution). The main simulation parameters obtained for the production runs are summarized in table Appendix A.6.
Appendix A.6. Description of the fitting procedure The bridging relation (5) from the main text requires explicit knowledge of σ(τ ) and µ(τ ). While our theoretical arguments provide estimates for very short and very long times, data from direct numerical simulations is needed to obtain information for the full range of time scales. The fitting procedure to obtain σ(τ ) and µ(τ ) is explained in detail in the following. 1. 2   Table A1.
Main DNS parameters. Direct numerical simulations are run over the time interval [t 0 , t 1 ]; field information is computed on a real space grid of n 3 points and 10 7 tracer trajectories are integrated for each flow. Characteristic parameters are as follows: Taylor-based Reynolds number R λ , root-mean-squared velocity u = (2 dkE(k)/3) 1/2 (given in DNS units, E(k) is the energy spectrum), integral length scale L = π 2u 2 dk k E(k) (given in DNS units), ratios of integral to Kolmogorov scales (η and τ η are computed from the mean kinetic energy dissipation ε and the kinematic viscosity ν, and the integral time scale is computed as T = L/u), and resolution criterion k max η, where k max is the highest wavenumber resolved by our pseudo-spectral code.  Figure A1. Normalization factor N (τ ) required by (5) for the data at different R λ . N (τ ) reaches unity at the integral time T independently of R λ .
From our numerical data, we obtain the PDFs of Lagrangian velocity increment magnitudes f L (v; τ ). Our subsequent fitting is based on structure functions of these PDFs, which serve as a reference; the pth-order structure function is defined as For the fitting, even-order structure functions up to order eight are considered. To evaluate the structure functions from the right-hand side of the bridging relation (5), we furthermore obtain the Eulerian velocity increment magnitude PDFs f E (v; R) directly from velocity field data. Both Lagrangian and Eulerian data are averaged over the full simulation time t 1 − t 0 . σ(τ )/η R λ = 148 R λ = 208 R λ = 259 R λ = 304 R λ = 316 R λ = 429 Figure A2. Results of the fitting procedure (continuous lines) compared with the empirical formula for σ(τ ) (dotted lines) for various Reynolds numbers. Data from different data sets was multiplied with successive powers of 3 for clarity. The data for all the different R λ is matched quite well by the empirical three-power-law formula (7) for σ(τ ) given in the main text.  Figure A3. Comparison of Lagrangian structure functions for R λ = 429: DNS data (filled circles) and results of applying the bridging relation with σ(τ ) and µ(τ ) fitted to minimize (A.38) (continuous lines), as well as using the analytic expressions (7) and (8), respectively, for σ(τ ) and µ(τ ) (dotted lines).
For fixed τ and given initial values of σ and µ, model PDFsf L (v; τ ) are evaluated from (5). Because the bridging relation mixes Eulerian velocity increments for different scales depending on the velocity increment magnitude, the model PDFs need to be explicitly normalized. The normalization factor N (τ ), computed at the end of the fitting procedure, is shown in figure A1. For time lags below one τ η , for which the dependence on the velocity increment is most pronounced, N (τ ) deviates significantly from unity. As expected, this deviation decreases with increasing time lag.
After computing the model PDFs, the corresponding model structure functions we obtain best fits for σ and µ. For the results presented, the "Nelder-Mead" algorithm [41,42] (through the minimize method available in the scipy.optimize [43] package) is used. Repeating this procedure for all samples of τ under consideration, we obtain estimates for σ(τ ) and µ(τ ), which then inspired our analytical parameterizations (7) and (8) presented in the main text. For the analytical µ(τ ) we combine the short-time prediction of τ /2 with an exponential drop-off for larger times to take into account the decreasing correlation with the velocity increment. It has to be pointed out that the precise determination of the drop-off of µ(τ ) becomes technically more difficult at long times where its influence is almost negligible. For σ(τ ), the data indicates three regimes, which we capture analytically by an extension of Batchelor interpolation [44]. We then iterate this procedure taking our analytical parameterizations as improved initial guesses until convergence is reached. In figure A2, a direct comparison between the results of the fitting procedure for σ(τ ) and the corresponding empirical formula (7) for various Reynolds numbers shows good agreement. The resulting structure functions from the DNS data, along with the ones obtained from the fitting procedure as well as from the analytical parameterizations, are shown in figure A3.