CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

gamma-ReTheta turbulence model for predicting transitional flows

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree23Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 30, 2011, 04:32
Default
  #41
jms
Member
 
José
Join Date: Jan 2011
Posts: 73
Rep Power: 15
jms is on a distinguished road
Airfoil FX-W-343.

geom_FX77W343.png

Comparison for the results with and withou transition model in the boudnary layer.

free_full_OF_cl.png

free_full_OF_cd.png
jms is offline   Reply With Quote

Old   April 30, 2011, 04:33
Default
  #42
jms
Member
 
José
Join Date: Jan 2011
Posts: 73
Rep Power: 15
jms is on a distinguished road
Airfoil FX-77-W-343.

Comparison for the results given by OpenFOAM and ANSYS CFX.

free_OF_CFX_cl.png

free_OF_CFX_cd.png
jms is offline   Reply With Quote

Old   April 30, 2011, 04:35
Default
  #43
jms
Member
 
José
Join Date: Jan 2011
Posts: 73
Rep Power: 15
jms is on a distinguished road
Airfoil FX-77-W-400.

geom_FX77W400.png

Comparison of results of OpenFOAM and ANSYS CFX. Both with transition model in the boundary layer applied.

free_OF_cfx_cl.png

free_OF_cfx_cd.png
jms is offline   Reply With Quote

Old   April 30, 2011, 04:40
Default
  #44
jms
Member
 
José
Join Date: Jan 2011
Posts: 73
Rep Power: 15
jms is on a distinguished road
Airfoil FB-3500-0875

geom_FB3500_0875.png

Comparison of OpenFOAM, ANSYS CFX, XFOIL and Wind Tunnel measurements.

free_OF_cl.png

free_OF_cd.png
jms is offline   Reply With Quote

Old   April 30, 2011, 04:41
Default
  #45
jms
Member
 
José
Join Date: Jan 2011
Posts: 73
Rep Power: 15
jms is on a distinguished road
Airfoil FB-3500-1750

geom_FB3500_1750.png

Comparison of OpenFOAM, ANSYS CFX, XFOIL and Wind Tunnel measurements.

fb3500_1750_cl.png

fb3500_1750_cd.png
jms is offline   Reply With Quote

Old   May 1, 2011, 14:51
Default
  #46
Senior Member
 
Join Date: Nov 2010
Posts: 139
Rep Power: 15
taxalian is on a distinguished road
Send a message via Skype™ to taxalian
hi jsm,
first of all can you comment on the results obtained for NACA0012 airfoil with blunt trailing edge. The results look reasonable in the linear region with OF, CFX versus the literature.
But what happens in the maximum Cl region with CFX results the lift drops dramatically. I guess you did not run that long enough to get a converged solution or their is any other explanation.
On the other why with OF you overpredicted the maximum Cl and OF does not seem to predict the leading edge type stall observed in the literature.
I have the following queries:
--> what bc's you used
--> what turbulence level, is this same for both CFX and OF
--> did you try unsteady in the non linear region

really looking forward to your reply.

taxalian.
taxalian is offline   Reply With Quote

Old   May 5, 2011, 04:30
Default
  #47
Senior Member
 
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 18
FelixL is on a distinguished road
Hello, José,


thanks for sharing your results with us. Especially the NACA0012 looks great! I didn't think the implementation would perform so well. It does seem to have some problems in the non-linear region, though, so probably there still are some improvements left.

A few questions:

- which correlations did you use?
- did you try the suggestion of SULUKSNA et. al. 2009, i.e. keeping lambda = 0 in gammaReThetatSST::ReThetat()?
- do you have any detailed information about the CFX implementation and simulations?


Thank you very much again!


Greetings,
Felix.
FelixL is offline   Reply With Quote

Old   June 1, 2011, 10:10
Default Running the transition model in parallel
  #48
jms
Member
 
José
Join Date: Jan 2011
Posts: 73
Rep Power: 15
jms is on a distinguished road
Dear all,

Has anyone of you tried to run the transition model given by Felix L. in parallel mode? I have tried it but it starts to diverge.
I know it sounds strange but this is how it is. I just one run one case in single mode and it works, then I do the decomposition and run it in parallel the force coefficients of the airfoil I am computing start to diverge.
I have tried to use another compiled model and it works (the k-omegaSST corrected, which only one number is added, no new equation is added as it happens with the transition model).

If any of you has something to say about this, please let me know.

Thanks,


José
jms is offline   Reply With Quote

Old   July 15, 2011, 10:58
Default
  #49
Senior Member
 
aerothermal's Avatar
 
Guilherme da Silva
Join Date: Aug 2010
Location: Sao Paulo - Brazil
Posts: 118
Rep Power: 16
aerothermal is on a distinguished road
Send a message via Skype™ to aerothermal
Hello Felix, Jose and All,

Do you have your test cases available? Both system and constant folder? Mesh file?

I would like to test here and compare with my data simulated during my PhD thesis about aircraft icing. At that time, I did some validation with ERCOFTAC, transonic airfoil, NACA airfoil (two types) and non-isothermal flat plate cases in order to check the results of the transition (prediction and development) models of the boundary layer codes. Therefore, I am interested to see or download any case you may have to compare with my results. Further, I intend to use OpenFoam in aircraft icing cases in order to calculate local Cf and St around it.

You can download my thesis from here:
http://www.teses.usp.br/teses/dispon...-082825/en.php

Regards,

Guilherme

Last edited by aerothermal; July 15, 2011 at 11:01. Reason: header and greetings
aerothermal is offline   Reply With Quote

Old   August 18, 2011, 09:54
Default
  #50
New Member
 
Christian
Join Date: Aug 2011
Posts: 2
Rep Power: 0
chvh is on a distinguished road
Hello,

I need help. I want to compile your model in openfoam 2.01. But i got some error.

Code:
gammaReThetatSST.C: In member function ‘void Foam::incompressible::RASModels::gammaReThetatSST::ReThetat(Foam::volScalarField&) const’:
gammaReThetatSST.C:244: error: no match for ‘operator[]’ in ‘Foam::incompressible::turbulenceModel::nu() const()[cellI]’
gammaReThetatSST.C: In member function ‘virtual void Foam::incompressible::RASModels::gammaReThetatSST::correct()’:
gammaReThetatSST.C:726: error: ‘omega0_’ was not declared in this scope
gammaReThetatSST.C:749: error: ‘k0_’ was not declared in this scope
In file included from gammaReThetatSST.H:109,
                 from gammaReThetatSST.C:26:
/opt/openfoam201/src/turbulenceModels/incompressible/RAS/lnInclude/RASModel.H: In static member function ‘static Foam::autoPtr<Foam::incompressible::RASModel> Foam::incompressible::RASModel::adddictionaryConstructorToTable<RASModelType>::New(const Foam::volVectorField&, const Foam::surfaceScalarField&, Foam::transportModel&, const Foam::word&) [with RASModelType = Foam::incompressible::RASModels::gammaReThetatSST]’:
/opt/openfoam201/src/turbulenceModels/incompressible/RAS/lnInclude/RASModel.H:138:   instantiated from ‘Foam::incompressible::RASModel::adddictionaryConstructorToTable<RASModelType>::adddictionaryConstructorToTable(const Foam::word&) [with RASModelType = Foam::incompressible::RASModels::gammaReThetatSST]’
gammaReThetatSST.C:44:   instantiated from here
/opt/openfoam201/src/turbulenceModels/incompressible/RAS/lnInclude/RASModel.H:126: error: no matching function for call to ‘Foam::incompressible::RASModels::gammaReThetatSST::gammaReThetatSST(const Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&, const Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::transportModel&, const Foam::word&)’
gammaReThetatSST.C:341: note: candidates are: Foam::incompressible::RASModels::gammaReThetatSST::gammaReThetatSST(const Foam::volVectorField&, const Foam::surfaceScalarField&, Foam::transportModel&)
gammaReThetatSST.H:128: note:                 Foam::incompressible::RASModels::gammaReThetatSST::gammaReThetatSST(const Foam::incompressible::RASModels::gammaReThetatSST&)
make: *** [Make/linux64GccDPOpt/gammaReThetatSST.o] Error 1
Could you help me please?
chvh is offline   Reply With Quote

Old   September 6, 2011, 05:09
Default
  #51
Senior Member
 
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 18
FelixL is on a distinguished road
Hello,


sorry about being slow with things lately - I caught the chance of having a cold and thus staying at home today to make the turbulence model compatible with OF 2.0.x. Please bear in mind that the code still is NOT validated!

Changelog:

Code:
    2011-09-06 : - ensured compatibility with OpenFOAM 2.0.x
Greetings,
Felix
Attached Files
File Type: gz gammaReThetatSST.tar.gz (7.2 KB, 144 views)
wiedangel, sail and Sune like this.
FelixL is offline   Reply With Quote

Old   September 12, 2011, 07:39
Default
  #52
New Member
 
Christian
Join Date: Aug 2011
Posts: 2
Rep Power: 0
chvh is on a distinguished road
Hello Felix,

thanks for your work. It was possible to compile so that i'm going to test it in the next days.

Greetings,
Christian
chvh is offline   Reply With Quote

Old   November 10, 2011, 16:49
Default
  #53
Senior Member
 
sail's Avatar
 
Vieri Abolaffio
Join Date: Jul 2010
Location: Always on the move.
Posts: 308
Rep Power: 17
sail is on a distinguished road
Dear Felix,

thanks for the good job and the effort to implement this turbolence model in openFoam. Exactly what i was looking for. I owe you a beer.

I will test and validate it in the next few days and report back.
sail is offline   Reply With Quote

Old   November 16, 2011, 07:26
Default
  #54
Senior Member
 
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 18
FelixL is on a distinguished road
Cheers, mate, a cold beer's always appreciated!

Greetings,
Felix
FelixL is offline   Reply With Quote

Old   February 22, 2012, 14:51
Default
  #55
New Member
 
xiaobao jia
Join Date: Jun 2011
Posts: 7
Rep Power: 15
xiaobao_jia is on a distinguished road
Hello, FelixL

I am using version 1.6-ext and I am a new user. Do you think I can apply this transitional model to komega SST model in version 1.6-ext?
In version 1.6-ext, there are both Komega SST low Rey model and Komega SST model. Which one we need to apply this transitional model?

please help me!


regards

jia
xiaobao_jia is offline   Reply With Quote

Old   November 11, 2012, 08:06
Default
  #56
Member
 
Zifei Yin
Join Date: Sep 2012
Location: Shanghai & Ames
Posts: 33
Rep Power: 14
yzf1215 is on a distinguished road
I have run the code using the T3A testcase, but when I apply the model in airfoil, it blows up. Can somebody give me an testcase of airfoil?
yzf1215 is offline   Reply With Quote

Old   February 18, 2014, 07:08
Default
  #57
Senior Member
 
Artur's Avatar
 
Artur
Join Date: May 2013
Location: Southampton, UK
Posts: 372
Rep Power: 20
Artur will become famous soon enough
Hi Foamers!

First of all, thank you Felix for sharing the code, saved me (and likely many others) lots of effort in implementing it.

I just wanted to post my experience with using the model with OpenFOAM-2.2.2 for future reference if anyone else tries to compile under this version.

List of changes I had to make to make it work (I started off from the most recent version that you posted in this thread):
1. divDevRhoReff method had to be added (simply copied over from standard kOmegaSST)
2. in the correct() method there was a problem with some naming convention at the line:
Code:
    omega_.boundaryField().updateCoeffs();
Closer look at kOmegaSST.C revealed that replacing the line before that as follows solves the issue:
Code:
    //volScalarField G("RASModel::G", nut_*2*S2); // replaced this
    volScalarField G(GName(), nut_*2*S2);
I also took the liberty of updating the test case you posted on page 1 of this thread to match the recent OF conventions of schemes.

Once again thanks a lot for posting this.

Peace,

Artur
Attached Images
File Type: jpg wallShearStress.jpg (21.5 KB, 63 views)
Attached Files
File Type: gz gammaReThetatSST.tar.gz (8.3 KB, 56 views)
File Type: gz T3A.tar.gz (26.3 KB, 85 views)
Artur is offline   Reply With Quote

Old   September 2, 2014, 05:03
Default
  #58
New Member
 
Join Date: Apr 2014
Posts: 1
Rep Power: 0
Guffalo is on a distinguished road
Hi there!
Has anyone of you tried Felix implementation on an airfoil case and had a look on the C_f distribution? Especially transition onset points and length?

I am trying to get a solution with the S809 airfoil right now. It is kind of difficult to get a well converged solution I think. My C_l is not too bad with a 5% error compared to experimental results. But if I look at my pressure distribution the onset point on the upper side has significantly moved towards the leading edge.
I also made some simulations with Ansys CFX for comparison with the exact same mesh and I get good and reasonable results respecting the pressure distribution.

Anyone who faces the same complications?
Guffalo is offline   Reply With Quote

Old   September 11, 2014, 09:21
Default About the Missing Correlations in Gamma-Re_Theta Turbulence Model
  #59
New Member
 
Mazharul Islam
Join Date: Jan 2014
Posts: 6
Rep Power: 12
sbvawt is on a distinguished road
Thanks a lot for sharing your code Artur. We know that there are two proprietary correlations in the original gamma-Re_theta turbulence model proposed by Meter et. al which are unknown to the public and implemented in the commercial code like Fluent. Different researchers (e.g. Sørensen) have attempted to propose their own correlations in the past. May I know which correlations are used in the default code?

I have attempted to use your code for analyzing NACA 4415 and found that though Cd values are closer to the experimental results, but Cl values are worse than kOmegaSST!

Thanks.
sbvawt is offline   Reply With Quote

Old   September 11, 2014, 09:27
Default
  #60
Senior Member
 
Artur's Avatar
 
Artur
Join Date: May 2013
Location: Southampton, UK
Posts: 372
Rep Power: 20
Artur will become famous soon enough
Hi,

I was not the one who created the code, I merely upgraded it to work with OF version 2.2.2. Hopefully the original Author (Felix) will be able to answer your questions.

Sorry,

A
Artur is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Low Reynolds k-epsilon model YJZ ANSYS 1 August 20, 2010 13:57
Centrifugal Pump and Turbulence Model Michiel CFX 12 January 25, 2010 03:20
RSM & Transitional Flows!!! Erika FLUENT 0 March 31, 2006 10:56
Turbulence model Herry Phoenics 1 May 29, 2003 13:19
turbulent separated flows Yin Fidelity CFD 9 February 19, 2003 11:50


All times are GMT -4. The time now is 11:28.