// Compute Filter scalar alph = 3.0; //Should be in a dictionary scalarField Lt = sqrt(k_)/(betaStar_*omega_); scalarField lt = alph*Foam::max(Foam::pow(mesh_.V().field(), 1.0/3.0), (mag(U_)*runTime_.deltaT())->internalField()); // Re-calculate viscosity nut_.internalField() = Foam::min(Foam::pow(lt/Lt,4.0/3.0), 1.0)* (a1_*k_/max(a1_*omega_, F2()*sqrt(S2)))->internalField();