[Tibo99]

Hi everyone,

I need some help to find the specific location where I can apply a correction factor on the Reynold Stresses equation for the k-epsilon turbulence model.

From the pictures in attachment, the correction factor is apply on the Laplacian term from the implicit part of the Reynolds stresses equation.

I used dOxygen in order to try to located the term, but I don't know from which file and term that I need to start in order to find the path that will lead me to the right location.

Thanks for your help and guidance.

Best Regards,

[agustinvo]

Hi Tibo99,


it is a bit hard to understand your question: are you modifying a k-epsilon model or a Reynolds stress model? As far as I know, if you want to define a specific divDevReff or similiar you can do it directly in your new k-epsilon model. You can take a look into the code of old OpenFOAM versions to have a hint about this.


Regards,
agustinvo

[Tibo99]

Thank you for replying!

Sorry that my post could be a bit confusing. In fact, by talking about the Reynold stresses, that probably bring the confusion.

The reason why I mention the Reynold stresses is that in the following paper:
https://onlinelibrary.wiley.com/doi/....1002/fld.2709,

the authors talk about the Reynold Stress equation (see PDF in attachment for a section of the text) and then apply this specific correction on the Laplacian term for the implicit part of the equation. Also, they use the k-epsilon model.

Since you talk about the term 'divDevReff' and which I found the term in the file 'UEqn.H', I had the intuition that's could be the function I have to follow the path in doxygen in order to find the location where to apply the correction. Unfortunately, I use OF v1906, but from what I've seen in doxygen, it look similar to the new OF v2006 online.

'divDevReff' seems to be a virtual function and I ended up to hit a 'cul de sac' in the file 'IncompressibleTurbulenceModel.C' in doxygen and don't know where to go from there.

Thank again for the discussion.

Reagrds,