[gu1]

Hi,

Could someone tell me what these featured lines do?

template<class BasicTurbulenceModel>
void kEpsilon<BasicTurbulenceModel>::correctNut()
{
this->nut_ = Cmu_*sqr(k_)/epsilon_;
this->nut_.correctBoundaryConditions();
fv:: options::New(this->mesh_).correct(this->nut_);

BasicTurbulenceModel::correctNut(); // If I call only correctNut() it gives error while executing the solverFoam.
}

and

void kEpsilon<BasicTurbulenceModel>::correct()
{
if (!this->turbulence_) //Is it necessary?
{
return;
}

fv:: options& fvOptions(fv:: options::New(this->mesh_));

eddyViscosity<RASModel<BasicTurbulenceModel>>::cor rect(); // If I call only correct() it gives error while executing the solverFoam.

I was testing a version I compiled from the kEpsilon model (OF5.0) and both I did not insert (because when I did, it gave error in executing the solver - to my understanding, I was closing the loop) in my code ... I noticed a difference in results, but very small. Unfortunately '' I could not understand '' what they do ... if anyone can teach me.

[anon_q]

I think you can find the answer in this link:
https://pingpong.chalmers.se/public/...o?item=3855255

[gu1]

Hi Linda,

I got to read this file, but he did not answer my questions.

It says that part of the code: eddyViscosity<RASModel<BasicTurbulenceModel>>::cor rect();, "...that it first calls the correct() function of the base class, in case there are some general things to be done.", what general things?

...and about: BasicTurbulenceModel::correctNut();, ...he does not say anything.

[simrego]

Hi!


For the second question: if (!this->turbulence_) //Is it necessary?
Yes, it is. turbulence_ is a switch from the RASModel<BasicTurbulenceModel>. This is the switch in the RASDict. This is usually on. But if you turn it off, you just skip this part of the code (return from the function), since you turned the turbulence off.


eddyViscosity<RASModel<BasicTurbulenceModel>>::cor rect();
Honestly I have no idea what basic things he talking about.
If you check the code:

template<class BasicTurbulenceModel>
 void Foam::eddyViscosity<BasicTurbulenceModel>::correct()
 {
     BasicTurbulenceModel::correct();
 }

This function just calls the correct() function from the turbulence model.
Maybe this is for further developments. So if you have to add something to all turbulence models when you call the correct() function, you can add it in the base class.


I hope i understood the code correctly and I'm right.