Code for EARSM by Hellsten and RSM by Manceau, square duct, 2d turbulent channel
 

Hi everybody,

under the link

http://www.istm.kit.edu/english/83.php


you can download some codes and tutorials of turbulence models I implemented in OpenFOAM for different versions. They are implemented for incompressible flow (simpleFoam, pisoFoam,...). I validated them at serveral cases (2d channel, square duct, annular pipe flow, rod bundles).

OpenFOAM 1.7.0
Code: nonlinear algebraic stress model ( EARSM ) by Hellsten and two validationcases: channel flow and square duct. Good agreement of the velocity profile at the 2d channel. The uv Reynolds stress is good (even for a annular pipe flow). (uu, vv, ww) are too low.
Pic: Attachment 23602

OpenFOAM 2.1.0
Code: Reynolds stress model ( RSM ) by Manceau and the same two validation cases. Good agreement of velocity profiles, Reynolds stresses and secondary flows (see the pics)

PICS:
RSM: Reynolds stress profiles of a 2d channel
Attachment 23599

RSM: Square Duct
Attachment 23598
Attachment 23600
Attachment 23601



Post-Processing yplus: To check the y+ values: post-processing tool to calculate them using low Reynolds models.

Good luck ;)

Thomas

Hello Thomas,

thank you very much for providing your files!

I am interested in your square duct simulation using the Reynolds-stress model you implemented. Unfortunately, the link to the case setup (OF211_RSM_Manceau_quadchannel) doesn't work correctly (it's linked to the EARSM case). I wonder if you could fix it or upload your case folder here, since I have some problems with setting the right boundary conditions / fvSchemes.

Thank you in advance!
Alex

Hi,

I have uploaded it now.

Try following link:

https://dl.dropboxusercontent.com/u/...channel.tar.gz

Regards Thomas

Hello Thomas,

thank you very much for your reply. Could you also upload the 3D case of a square duct (OF211_RSM_Manceau_quadchannel) or do you use the same boundary conditions and fvSchemes as for the 2D channel?

Furthermore, I have some questions concerning the 2D case you've just uploaded.
What kind of boundary condition do you use for epsilon? Do you calculate the values for epsilon a priori and impose them as fixedValue?
Is the initialization of the flow field important?
If I get it right, there is no need for any lowReWallfunctions when using your code?

Best wishes,
Alex

Hi,

here is the correct link for the quad-channel.

https://dl.dropboxusercontent.com/u/...channel.tar.gz

It is a low Reynolds-number model with a elliptic factor, no damping functions are necessary... epsilon is calculated by a correction-function within the code. Its put to fixed value, but every time step the boundary values are updated.

I made pre-simulations with the Launder-Sharma-model and used the results as starting conditions, because of the model is very sensitive.
Of course, you have to correct the boundary conditions.

You must read the papers of Manceau to understand the code. (an improved version ....)

Best regards,

Thomas

Hi Thomas,

thank you once again. You helped me a lot!

Alex

Hi Thomas,

I've tried your quadchannel case and it works fine if I run it on a single core. Since I'd like to test it also on a bigger case I'll need to run it in parallel. The problem is that the case cannot be decomposed and it returns me the following error:

Do you have the same problem?

Thanks!
Blaz

Dear Thomas,

Would you please write down exact reference to works of Manceau and Hellsten? I am not sure about title of the papers.

Regards,
Mohammadreza

Expressing the production term 'G'
 

Hi Thomas

Thanks for the code!

Thanks for sharing your code and the results.

Has anybody managed to compile the EARSM code for newer openfoam versions?

I have been struggling with it for some time without any success.

Here are some infos; bit late ;)

the EARSM from Hellsten is the version from his Dissertation:
new two equation turbulence model for aerodynamics applications.

The RSM by Manceau should be the one from the paper:
AN IMPROVED VERSION OF THE ELLIPTIC BLENDING MODEL
APPLICATION TO NON-ROTATING AND ROTATING CHANNEL FLOWS

I'm working with Remi Manceau now, your code is helpful a lot :)
A note:
The current version of the code is not robust.
You need to make it stable.
Split the anisotoropy tensor to deviatic and diagonal and use Sp in deviatic part! check the implementation of Spezial-Sarkar-Gatski RSM as well.

This is my REq:

tmp<fvSymmTensorMatrix> REqn
(
fvm::ddt(R_)
+ fvm::div(phi_, R_)
//+ fvm::SuSp(-fvc::div(phi_), R_)
- fvm::Sp(fvc::div(phi_), R_)
- fvm::laplacian(nu(), R_)
- fvm::laplacian(CMu/sigmak*TEBM*R_,R_) //anisotrope SSG-Modellierung
+ fvm::Sp(anisotrope, R_)
==
P
- fvm::Sp(epsilonterm1, R_)
+isotrope
-epsilonterm2
+phistar
);
////////////////////
If you have problem with it contact me!
Ardalan Javadi

request to re-upload the code
 

Sir Baumann,
none of the links is working in the present time (2019) :confused: so if you are still active, please update the link.
Thank you for your time. :)

code for the RSM by Manceau und the EARSM by Hellsten
 

Hallo, the former links don't work anymore.

Here is the code for the RSM by Manceau und the EARSM by Hellsten.
Its the OF 1.7 and the OpenFOAM 2.11 Version.

Dear Thomas Baumann,
thanks for the codes that you shared. Now I am trying to converting your RSTM Elliptic Blending model in OpenFoam 6: I did only minor syntax changes to your version and then I compiled it successfully. I tried then to reproduce your results for the channel flow at Re_tau=395, but I got very different results with poor accuracy. Could you please share again also the channel test case that you used for the validation?
Thank you very much in advance
Matilde