Reynolds number effects in Lagrangian stochastic models of turbulent dispersion

Abstract

A second-order autoregressive equation is used to model the acceleration of fluid particles in turbulence in order to study the effect of Reynolds number on Lagrangian turbulence statistics. It is shown that this approach provides a good representation of dissipation subrange structure of Lagrangian velocity and acceleration statistics. The parameters of the model, two time scales representing the energy-containing and dissipation scales, are determined by matching the model velocity autocorrelation function to Kolmogorov similarity forms in the inertial subrange and the dissipation subrange. The model is tested against the Lagrangian statistics obtained by Yeung and Pope [J. Fluid Mech. 207, 531 (1989)] from direct numerical simulations of turbulence. Agreement between the model predictions and simulation data for second-order Lagrangian statistics such as the velocity structure function, the acceleration correlation function, and the dispersion of fluid particles is excellent, indicating that the main departures from Kolmogorov’s theory of local isotropy shown by the simulation data are due to low Reynolds number. For Reynolds numbers typical of laboratory experiments and direct numerical simulations of turbulence the root-mean-square dispersion of marked particles is changed from the Langevin equation (i.e., infinite Reynolds number) prediction by up to about 50% at large times. Most of this change can be accounted for by the change in the Lagrangian integral time scale. It is also shown that Reynolds number effects in laboratory dispersion or Lagrangian turbulence measurements can cause significant errors (typically of order 50%) when the value of the Kolmogorov Lagrangian structure function constant C0 is estimated by fitting the predictions of the Langevin equation to these data. A value C0 = 7 is obtained by fitting the new model to the direct simulation data.