At separations well above the integral scale of turbulence, two fluid velocities can decorrelate spatially. If the relative-velocity autocorrelation is stationary and integrable, the displacement identity used in Taylor turbulent dispersion gives the displayed long-time linear mean-square growth. With speed and Lagrangian integral time of order , its relative diffusivity is of order . The characteristic separation grows as , unlike inertial-range Richardson pair dispersion.
Let , assumed finite. The displacement is . By stationarity, its two-time covariance is , with even for a real scalar process. Integrating over the two triangles of the square gives the Taylor turbulent dispersion identity
A sufficient condition for ordinary long-time variance growth is with . The dominated convergence theorem then gives
More generally the necessary asymptotic condition for a finite positive variance-based diffusion coefficient is convergence of to that coefficient. Absolute integrability is sufficient, not necessary. A zero integral does not give nondegenerate ordinary diffusion. Moreover a second-moment calculation alone does not establish a normal distribution or an entire diffusive scaling limit; those require additional probabilistic assumptions.
For the specified positive autocorrelation, direct integration gives, when ,
Thus the power-law velocity-correlation dispersion has the following regimes. If , the Lagrangian integral time is and ordinary variance growth holds, with . At the exact expression is . If ,
For ,
The last two regimes spread faster than ordinary diffusion and have no finite long-time diffusion coefficient. Their slowly decaying velocity memory is responsible. At short times all these correlations instead give the ballistic law .
Take a positive initial separation . In the dissipation range, the velocity field is smooth, so the difference of the two Lagrangian trajectories obeys the linearized separation equation . Chaotic stretching gives a positive Lyapunov exponent of order and a typical separation
This needs a nonzero initial separation: exactly coincident particles remain coincident in a smooth flow.
Once lies in the inertial range, the Kolmogorov 1941 theory gives a relative speed of order . A scale-local typical-separation estimate integrates to
After loss of the entry-scale memory, the statistical Richardson pair dispersion law is . Thus the characteristic separation grows as while . This is a statistical scaling closure, rather than a deterministic equation for every pair.
For , the particle velocities become approximately independent. If their Lagrangian velocity autocorrelations have a finite Lagrangian integral time , the diffusive large-scale pair dispersion has effective relative diffusivity of order . After a further velocity-decorrelation transient,
This is the long-time mechanism of Taylor turbulent dispersion applied to the difference of two decorrelated trajectories. The root mean square separation eventually loses memory of the term.
Reaching at the end of the inertial range takes
Order-one coefficients depend on the dispersion closure. The initial exponential stage can take a long time if is extremely small, but the inertial-range stage is of order one outer eddy turnover time.
Figure 1.
Exponential, Richardson and diffusive regimes of turbulent pair separation
. Details.
The three-regime turbulent pair-separation model in the sketch matches the regimes continuously for illustration. Its coefficients are schematic; the universal claims here are the scaling powers under the stated stretching, locality and decorrelation assumptions.
Since the particle starts at the origin, its Lagrangian trajectory satisfies . Differentiating the mean-square displacement gives
Here we set and assume finite second moments so that averaging and integration commute. Under Lagrangian stationarity write the autocorrelation as . Uniform fluid-tracer sampling in homogeneous incompressible turbulence gives . A second integration yields the exact Taylor turbulent dispersion formula
For times short compared with velocity decorrelation, continuity at zero gives throughout the integral. Thus
The particle retains nearly its initial velocity: this is ballistic transport.
If is integrable, then at long times the displacement derivative approaches , using the Lagrangian integral time. Therefore
This is diffusive transport: successive long-time displacements lose memory of one another. In an isotropic three-dimensional diffusion equation, , so . The shorthand and presumes that is representative of the decorrelation time. With strongly oscillatory correlations, the area-based integral time need not characterize every crossover; the exact integral formula remains the appropriate starting point.
The three questions concern different statistics. Taylor turbulent dispersion follows the absolute displacement of one marked fluid element from its release position, or the spreading of the centre of a cloud. Its input is the velocity autocorrelation along a Lagrangian trajectory. It relates short-time ballistic motion to long-time effective diffusion. This is not the shear-flow mechanism ordinarily called Taylor dispersion.
Richardson pair dispersion instead follows the separation of two fluid elements, or the growth of a cloud relative to its moving centre. Common sweeping by an eddy much larger than the cloud contributes to absolute displacement but largely cancels from relative displacement. The relative diffusivity depends on separation; the classical inertial-range prediction is a mean-square separation regime after memory of the release separation is lost.
Obukhov-Corrsin theory, also called Kolmogorov-Obukhov-Corrsin theory, concerns spatial fluctuations of a passive scalar, their transfer across scales and their destruction by molecular diffusion. Its objects are scalar structure functions and scalar spectra, with an inertial-convective prediction controlled by scalar-variance flux and the turbulent energy supply. Thus Taylor describes one-particle transport, Richardson describes relative spreading, and Obukhov-Corrsin describes scale-by-scale scalar mixing. Following a cloud's location is not the same as predicting its internal concentration fluctuations.