[renos]

Hello,

I want to modify turbulent viscosity of the k-omega turbulence model to avoid the overprediction according to the proposed limiter:

nut= min(k/omega,(sqrt(2)/3)*k/(magnitude of the strain tensor))

How I can write this in openfoam?

I have tried with this format, the compile of the turbulence model is without errors, but when I run the solver, then there is a division with zero.

template<class BasicMomentumTransportModel>
void kOmega2006<BasicMomentumTransportModel>::correctNu t
(
const volTensorField& gradU
)
{

this->nut_ = min(k_/omega_, (double(0.81649)*k_)/sqrt(magSqr(symm(gradU)))));
this->nut_.correctBoundaryConditions();
fv:ptions::New(this->mesh_).correct(this->nut_);
}


Anybody knows what's the error?

Kind regards,

Renos

[korenml]

Hello Renos,
did you figured it out?
I would rather to code it in this form

template<class BasicMomentumTransportModel>
void kOmega2006<BasicMomentumTransportModel>::correctNut
(
     const volScalarField& S2
)
{
     this->nut_ = min(k_/omega_, ((sqrt(2)/3)*k_)/sqrt(S2);
     this->nut_.correctBoundaryConditions();
     fv::options::New(this->mesh_).correct(this->nut_);
}

template<class BasicEddyViscosityModel>
void kOmega2006<BaxicMomentumTransportModel>::correctNut
{
     correctNut(2*magSqr(symm(fvc::grad(this->U_))));
}

See the code for kOmegaSST.

Or you can try if the kOmegaSST model is calculating without error, just to check your BC's are set correctly.

Best,
Tomas