CFD Online Discussion Forums - Implementation of dynOneEqEddy LES model

Hi all,

I have doubts about implementation of incompressible::dynOneEqEddy turbulence model.

For example, in the part of calculating ck, symmTensor MM is calculated as

Code:

MM = simpleFilter_(-2.0*delta()*sqrt(KK)*filter_(D))

in 2.2.x original source code.

but I think, MM must be calculated as

Code:

MM = delta()*(filter_(sqrt(k_)*D) - 2*sqrt(KK + filter_(k_))*filter_(D))

written in comment section of header file "dynOneEqEddy.H".

In addition, MM is calculated as the latter in "homogeneousDynOneEqEddy"

Why the former is used?