Generalised Langevin Model in Correspondence with a Chosen Standard Scalar-Flux Second-Moment Closure

Abstract

In the context of transported joint velocity-scalar probability density function methods, the correspondence between Generalised Langevin Models (GLM) for Lagrangian particle velocity evolution and Eulerian Reynolds-stress turbulence models has been established in the 1990’s by S.B. Pope. It was shown that the GLM representation of a given Reynolds stress model is not unique. It was also shown that a given GLM together with a given mixing model for particle composition evolution implies a differential scalar-flux model. In this paper, we study how extra constraints can be applied on the choice of the GLM coefficients in order to imply a chosen scalar-flux model. This correspondence between GLM-implied and standard scalar-flux models is based on the linear relaxation term and on the mean velocity gradient contributions in the rapid term. In general, GLM-implied models possibly involve more terms (including anisotropy effects in the scalar-flux decay rate and some high-order terms in the rapid-pressure-scrambling term). The proposed form of the GLM supposes a non-constant value for the diffusion coefficient C 0, originally identified as a Kolmogorov constant. Here, the value of C 0 is determined in order to yield the Monin model for linear relaxation of the scalar-flux, and the constant in the rapid-pressure contribution is related to the choice of the parameter β ∗ in the GLM. We finally show how GLM-implied scalar-flux models are in general dependent on the choice of the mixing model and how the proposed GLM can reduce this dependency. These developments are illustrated by results obtained from calculations of the Sydney bluff-body stabilised flame HM1.