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

Understanding k-omega SST model source code

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 8, 2009, 05:45
Default Understanding k-omega SST model source code
  #1
New Member
 
Tomi Honkanen
Join Date: Sep 2009
Location: Finland
Posts: 1
Rep Power: 0
tmhonka is on a distinguished road
Hi,

I'm currently working on my master's thesis and using OpenFOAM 1.5 for aerodynamics calculations (incompressible, subsonic). I was writing down equations used by the K-omega SST turbulence model implemented in OF1.5, when I run into some difficulties. I hope someone could shed some light into this.

The source code of the model clearly states that the implementation is based on this paper:

Menter, F., Esch, T.
"Elements of Industrial Heat Transfer Prediction"
16th Brazilian Congress of Mechanical Engineering (COBEM),
Nov. 2001

I do not have that, but instead I have looked at the following sources:

[1] Menter SST Two-Equation Model from 2003 http://turbmodels.larc.nasa.gov/sst.html
[2] Ten Years of Industrial Experience with the SST Turbulence Model http://cfd.mace.manchester.ac.uk/flo...nter-SST-paper
[3] K-omega SST in OpenFOAM-1.5 http://www.tfd.chalmers.se/~hani/kur...lenceModel.pdf
[4] http://www.cfd-online.com/Wiki/SST_k-omega_model

My questions are:

1) In the source code KOmegaSST.C on line 397 nut is defined as:
nut_ = a1_*k_/max(a1_*omega_, F2()*sqrt(S2));

and on line 344:

S2 = magSqr(symm(fvc::grad(U_)));


However on line 363 in the omega-equation, the term gamma*S^2 is stated as

gamma(F1)*2*S2

Now where does that factor 2 between gamma and S^2 come from?

In Ref [1] it is found that for the 2003 model in fact this formulation would
be correct (S=sqrt(2*S_ij*S_ij)), if this was used then the nut -equation is
inconsistent as the sqrt(2) factor is missing!

2) The G-term defined in line 345 is

nut_*2*S2

As far as I understand it is not defined like this in Refs [1,2,3,4]. How is this
done in the code when G-term is used in the same way as in Refs?

I suppose there is nothing wrong with the code, and just some programming
gimmick is used that I don't understand, but it would be much appreciated if
someone could help me to see how the model was implemented.


tmhonka is offline   Reply With Quote

Old   September 8, 2009, 08:33
Default
  #2
Member
 
Marco Müller
Join Date: Mar 2009
Location: Germany
Posts: 94
Rep Power: 17
marico is on a distinguished road
http://www.cfd-online.com/Forums/ope...komegasst.html
marico 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
where can i find the source code for KH-RT model? toto Main CFD Forum 1 March 22, 2011 22:57
Wall function implementation K Omega SSt cbarry OpenFOAM 3 August 18, 2009 11:09
about Y+ in SST model zuozicheng CFX 2 July 12, 2005 11:05
Use of an hydrodynamic model integrated in the vertical by the PHOENICS code Hafsia Zouhaier Main CFD Forum 2 March 15, 1999 09:23
public CFD Code development Heinz Wilkening Main CFD Forum 38 March 5, 1999 12:44


All times are GMT -4. The time now is 04:47.