At the Batchelor scalar microscale, molecular diffusion rate balances the smooth Kolmogorov strain rate . For this scale is much smaller than the Kolmogorov length scale.
Corrsin scalar microscale 2026-10-14
In the large-diffusivity regime, equate diffusion rate with inertial eddy rate . This gives the Corrsin scalar microscale, which exceeds the Kolmogorov length scale when .
This range lies above both the Kolmogorov length scale and the scalar diffusion cutoff, but below the forcing scale. Molecular viscosity and scalar diffusion have negligible direct influence on transfer at those scales. The Obukhov-Corrsin theory predicts a two-thirds scalar structure function there.
Assume statistically steady, three-dimensional turbulence at large Reynolds number, with local statistical homogeneity and statistical isotropy at scales small compared with the forcing scale . The Kolmogorov 1941 theory additionally assumes a local energy cascade, constant mean energy flux per unit mass , and negligible viscous dissipation within the inertial range. Its dimensional form neglects corrections from internal intermittency.
At separation , let be a typical velocity increment. The eddy turnover time is . Constant transfer rate then gives , hence
Here is the one-dimensional turbulent energy spectrum, obtained from with . The scaling holds only for
The Kolmogorov length scale follows by setting the scale-dependent Reynolds number to one. The associated Kolmogorov microscales give . At larger scales the forcing matters; at smaller scales kinematic viscosity matters.
Let be the one-dimensional turbulent energy spectrum, normalized so that . Large energy-containing eddies receive energy from the forcing. In the inertial range, nonlinear transfer passes this energy towards smaller scales at a mean rate per unit mass. Dimensional analysis of Kolmogorov 1941 theory gives
The dissipation range begins near the inverse of the Kolmogorov length scale, where kinematic viscosity removes the transferred kinetic energy. The large-scale spectrum depends on the forcing; the rising low- branch in this sketch is illustrative, rather than a universal power law.
The tube has neither homogeneous nor isotropic large-scale turbulence: gravity selects a vertical direction, its walls constrain the motion, and its density gradient varies with height. The width bounds the transverse size of the energy-containing eddies. At sufficiently high Reynolds number, much smaller eddies may become approximately locally isotropic and exhibit an inertial range, even though the full tube flow does not satisfy the homogeneous-isotropic assumptions.
To obtain the Reynolds stress, write and , with . Averaging the incompressible Navier-Stokes equations uses and gives the Reynolds-averaged momentum equation
Here the final term is the mean buoyancy force, with the hydrostatic reference removed. Thus averaging the quadratic momentum flux introduces the covariance ; multiplied by , it has the dimensions of a stress. Its divergence transfers momentum between the mean flow and the fluctuations.
An eddy viscosity model represents this Reynolds stress by
The isotropic term can be absorbed into the mean pressure, while the deviatoric term gives a diffusion of mean momentum with turbulent viscosity , where and are representative turbulent speed and mixing length. Unlike kinematic viscosity, this coefficient depends on the flow. The closure assumes a local relation between stress and mean strain; anisotropic, buoyancy-driven large eddies may require a more detailed Reynolds stress model.
For positive initial separation below the Kolmogorov length scale, a phenomenological pair-dispersion model has smooth-flow exponential stretching, inertial-range Richardson pair dispersion, and diffusive large-scale pair dispersion beyond the integral scale of turbulence. A positive Lyapunov exponent is needed for the first regime, scale locality and initial-memory loss for the second, and integrable relative-velocity autocorrelation for the third. The inertial-range crossing time is of order , one outer eddy turnover time.
For scalar diffusivity much smaller than kinematic viscosity, a scalar has structure below the Kolmogorov length scale while the velocity is smooth. The strain rate transfers scalar variance, giving . Integration produces a logarithm plus an additive matching constant.