Structure functions and intermittency
Use velocity increments and high-order moments to distinguish typical fluctuation size from rare intense events, without claiming unresolved exponent predictions.
On this page
The question for today
What you’ll understand
- Define longitudinal increments and distinguish signed from absolute moments.
- Calculate how rare events alter fourth moments without altering variance.
- Explain what intermittency means for scaling claims and statistical convergence.
Two turbulent records can have the same variance and still behave very differently. One may fluctuate moderately almost all the time; another may remain quiet between intense bursts. High-order moments distinguish these possibilities because large events receive rapidly increasing weight. Structure functions apply that idea to velocity differences across a specified separation. They connect statistical scaling with the physical question of how unevenly intense motion and dissipation are distributed.
Measure differences across a stated separation
Choose a separation vector and define the longitudinal increment as the velocity difference projected along it. The absolute pth-order structure function averages the pth power of its magnitude. It depends on separation, direction, flow position, and time unless statistical symmetries eliminate some of those variables. A mean over independent realizations is an ensemble average. Replacing it by a spatial or time average requires a justified sampling assumption; taking many nearby samples does not make them independent.
For even integer p, absolute and signed powers coincide. For odd p they generally do not. In particular, the signed longitudinal third moment can be negative while the absolute third moment cannot. Losing that sign erases information carried by the exact four-fifths relation in its inertial-range regime. Always identify which moment has been plotted before comparing a measured exponent or coefficient with theory. A label saying only “third-order structure function” leaves a consequential ambiguity.
Separate simple similarity from observed intermittency
A simple inertial-range similarity estimate gives , with coefficients carrying the required dimensions. Actual high-order statistics can depart from these exponents, reflecting scale-dependent distribution shapes and concentration of intense events. A useful diagnostic is flatness, the fourth moment divided by the second moment squared. A Gaussian increment has flatness three, but a different constant value alone does not establish scale-dependent turbulent intermittency. Examine how the normalized distribution and its tails change with separation.
For homogeneous isotropic incompressible three-dimensional turbulence, under the high-Reynolds-number inertial-range and suitable stationarity conditions, the signed third-order law is . Its sign is linked to the forward energy flux with this increment convention. The relation does not imply that every third-order absolute moment has that coefficient, nor does it determine all higher-order exponents. Proposed intermittency models can organize observations, but no general exact derivation of every turbulent scaling exponent is established.
Rare events demand more than a long-looking record
A fourth or sixth power can make a few intense events dominate an estimate. A record containing many ordinary fluctuations may therefore give a stable variance while its high-order moment continues changing dramatically. Assess convergence by comparing independent blocks, examining which amplitudes contribute to the integral, and checking sensor saturation and noise. Quietly clipping large values as outliers can remove the very physical events under investigation. Conversely, isolated electronic spikes can imitate a heavy tail and must be distinguished using the measurement chain.
Intermittency also has several contextual meanings. A turbulent patch alternating with a nearly nonturbulent region is not identical to small-scale internal intermittency within an already turbulent region. State the conditioning, separation, geometry, Reynolds number, and measured component. In a finite data range, a fitted slope has uncertainty and may mix forcing, inertial, and dissipation behavior. Reporting that limitation is part of the physical result. A good statistical model explains which features it captures and which rare-event probabilities remain insufficiently constrained.
Work it through
Equal variance, unequal tail weight
Compare two invented dimensionless increment distributions. In A, values +1 and −1 each occur with probability one half. In B, zero occurs with probability three quarters and values +2 and −2 each occur with probability one eighth. Calculate their second moments and flatness.
Both means vanish by symmetry. Distribution A has second moment 1 and fourth moment 1, because squaring or taking the fourth power of either nonzero value gives one. Its flatness is therefore 1.
Distribution B has second moment (1/8)×4+(1/8)×4=1. Its fourth moment is (1/8)×16+(1/8)×16=4. Its flatness is 4 despite the same variance as A.
The difference comes from rare larger amplitudes, not a different mean or second moment. Both signed third moments vanish by symmetry, so neither distribution is offered as a complete inertial-range turbulent increment model obeying the four-fifths law.
Matching variance leaves the high-order statistics undetermined. This calculation isolates tail sensitivity without pretending to reproduce the dynamics or scale dependence of actual turbulence.
The living notebook
Explore the model
Interactive lab · probability
How often will a rare event appear?
Set the success probability p to 0.05 and trials n to 20. Predict the mean event count before reading the exact binomial bars. Then increase n to 30 and compare the chance of zero events with a more common event probability. Each success represents a hypothetical threshold exceedance in an independent sample.
- Mean · np
- 5
- Variance · np(1−p)
- 2.5
- Total probability
- 1
These are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forThe lab displays an exact binomial count distribution, not velocity increments or resampled turbulent data. Real threshold exceedances may be correlated, so the independent-trial model does not establish their true frequency or intermittency exponent.
Pause & reason
Check your understanding
Your turn
Try it yourself
A symmetric dimensionless variable is zero with probability 15/16 and takes +4 and −4 each with probability 1/32. Calculate its second moment, fourth moment, and flatness. Explain why a short record could underestimate the last quantity.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Only the nonzero outcomes contribute. A sample that misses those rare outcomes misses most of the fourth-moment contribution.
Compare with the solution
The second moment is 2×(1/32)×16=1. The fourth moment is 2×(1/32)×256=16, giving flatness 16. Each nonzero event has low probability but a large fourth-power weight. A short or strongly correlated record can contain too few such events for a stable tail estimate; equality of measured variances alone would not demonstrate convergence of the fourth moment.
Keep the idea
Take it with you
- Odd signed moments and odd absolute moments contain different information.
- High-order statistics emphasize rare intense increments that a variance can hide.
- Measured intermittency requires scale dependence, sampling checks, and honest uncertainty about unresolved exponent predictions.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.