code_saturne User's Forum

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

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

Hello,

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

Best Regards
Martin

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