0-equation turbulence model

[CTrinca]

Hello everyone,

I need to implement a simple 0-equation turbulence model, where the turbulent viscosity is taken to be proportional to the streamwise component of the velocity. In my plans, I wanted to do so by modifying the molecular viscosity in usphyv and running the calculation as if it were laminar. Is this the correct/best way to proceed?
Thanks to anybody answering.
Regards,

Cosimo

[Yvan Fournier]

Hello,

This is probably the simplest way to do it, and probably should work, but for a cleaner approach, adding a model in src/base/tridim.f90 (at the same level as turbke, turrij, ....) might be better.

I'll let the turbulence specialists answer that one in more detail...

Regards,

Yvan

[Martin FERRAND]

Hello,

Your proposition is the best. Do not forget to activate a variable molecular viscosity.

Best Regards
Martin

[CTrinca]

Hello,

Thank you very much for your answers.
I have implemented the model in the way you suggested. It compiles, at least, but I am having convergence problems that I did not have before I introduced it. I am now trying to find out what is actually happening.
Thank you again nonetheless
Regards,

Cosimo