Non-Gaussian scalar statistics in homogeneous turbulence

Results are presented of numerical simulations of passive scalar mixing in homogeneous, incompressible turbulent flows. These results are generated via the Linear Eddy Model (LEM) and Direct Numerical Simulation (DNS) of turbulent flows under a variety of different conditions. The nature of mixing and its response to the turbulence field is examined and the single-point probability density function (p.d.f) of the scalar amplitude and the p.d.f.s of the scalar spatial-derivatives are constructed. It is shown that both Gaussian and exponential scalar p.d.f.s emerge depending on the parameters of the simulations and the initial conditions of the scalar field. Aided by the analyses of data, several reasons are identified for the non-Gaussian behaviour of the scalar amplitude. In particular, two mechanisms are identified for causing exponential p.d.f.s: (i) a non-uniform action of advection on the large and the small scalar scales, (ii) the nonlinear interaction of the scalar and the velocity fluctuations at small scales. In the absence of a constant non-zero mean scalar gradient, the behaviour of the scalar p.d.f is very sensitive to the initial conditions. In the presence of this gradient, an exponential p.d.f is not sustained regardless of initial conditions. The numerical results pertaining to the small-scale intermittency (non-Gaussian scalar derivatives) are in accord with laboratory experimental results. The statistics of the scalar derivatives and those of the velocity-scalar fluctuations are also in accord with laboratory measured results.