[shubz]

Hello guys,

Currently, I am trying to implement SA turbulence model with rotation and curvature correction in Openfoam 6.

The link below provided very useful information to do the same

To start with I've got the results with the standard Openfoam 6 SA model.


Now I'm trying with**SA-RCM*, the link below provided some very useful information.

https://hobbyfoam.blogspot.com/2015/...orrection.html


Unfortunately, I cannot directly use the .H and .C files that have been provided because Openfoam 6 has slightly different class implementations to avoid duplication of code.*


Based on this Blog, I have modified them accordingly (attached files)


Sorry, incase my programming terminology is not so great. bear with me.


There are two changes in the .H file*

1) Defining the dimensioned scalars cr1,cr2,cr3.

2) Defining the* f1 function*


In the .C file, I've made the following changes in order

1) In* Protected Member Functions, I've added the fr1 class* --> highlighted as new code*

2) In the Constructors, defined the lookup for cr1 ,cr2 ,cr3

3) In the member functions,* mentioned cr1 ,cr2,cr3 for the read class.

4) In the correct class, I have tried adding fr1,


So with all these changes, the code gets compiled successfully.


However, the simulation crashes when I add the fr1 to the correct class. If I keep the old lines as it is, everything works fine (obviously)*


Also, as you must have noticed, I don't use an if condition for the RCM correction, I want to use it by default but yeah, I would like to include that later on so that I don't have to compile two separate codes.*

I am confident about my setup since with the same setup I have ran the case with the default Openfoam6 SA Model and obtained good results.

Can you please help me out. I can provide you the files required for debugging this.

I don't know why I couldn't attach the files , so here is the Google drive Link where the .H and the .C files that I have modified using the above reference.

https://drive.google.com/open?id=1_L...eInc-klFGIU-76


Kindly help me out.. it's kind of urgent.

Regards

Shubham