[Crigerfra]

Hello,
I am would like to implement the Spalart-Allmaras turbulence model with QCR nonlinear stress model in Openfoam. However, i already wrote a code but when i test my code, the results do not agree with the other data available in literature. I tested the code for the tutorial case pitzDaily. In fact, the case pitzDaily can to be run using other nonlinear stress model already implemented in Openfoam. Can you help me?
This is my code:

void SpalartAllmarasQCR::correctNonlinearStress(const volTensorField& gradU)
{

volSymmTensorField S(symm(gradU));
volTensorField W(skew(gradU));


volScalarField sBar(mag(gradU));



nut_ = nuTilda_*fv1(this->chi());
nut_.correctBoundaryConditions();
fv:Options::New(mesh_).correct(nut_);

nonlinearStress_ = -4*Cr1_/(sBar)*nut_*twoSymm(W&S);



}

[lpz456]

Maybe the sign is reversed, such as twoSymm(S&W)