# Understanding k-omega SST model source code

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 September 8, 2009, 04:45 Understanding k-omega SST model source code #1 New Member   Tomi Honkanen Join Date: Sep 2009 Location: Finland Posts: 1 Rep Power: 0 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.

 September 8, 2009, 07:33 #2 Member   Marco Müller Join Date: Mar 2009 Location: Germany Posts: 94 Rep Power: 9

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post toto Main CFD Forum 1 March 22, 2011 22:57 cbarry OpenFOAM 3 August 18, 2009 10:09 zuozicheng CFX 2 July 12, 2005 10:05 Hafsia Zouhaier Main CFD Forum 2 March 15, 1999 09:23 Heinz Wilkening Main CFD Forum 38 March 5, 1999 12:44

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

 Contact Us - CFD Online - Top