CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Bugs

kOmegaSSTSAS implementation in OF 2.1

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

LinkBack Thread Tools Display Modes
Old   May 13, 2012, 11:48
Default kOmegaSSTSAS implementation in OF 2.1
New Member
Ben Colcord
Join Date: Mar 2012
Location: Atlanta
Posts: 2
Rep Power: 0
RealSoutherner is on a distinguished road

I think there are errors in the implementation of the kOmegaSSTSAS model in OF 2.1, but I'm not certain.

Firstly, the production term G differs from the kOmegaSST model and the basic kOmega model, where I believe they should be consistent. OF 1.6 has a consistent formulation in all 3.

In kOmegaSSTSAS(OF2.1):
S2 = 2*magSqr(symm(fvc::grad(U)));
G = nuSgs*0.5*S2;
nuSgs = a1*k/max(a1*omega, F2()*sqrt(S2));

In kOmegaSST(OF2.1):
S2 = 2*magSqr(symm(fvc::grad(U)));
G = nut*S2;
nut = a1*k/max(a1*omega, F2()*sqrt(S2));

So, G is half the magnitude in the kOmegaSSTSAS model. I can't find anywhere in the literature where this is true, but maybe someone knows. I have tested both formulations on a turbulent channel flow and actually see a slightly better result with the new "wrong" formulation, which makes me question the wrongness.

Secondly, the related term in the omega equation is given as gamma(F1)*0.5*S2. This is how it's shown in Menter et al (Ten Years of Industrial Experience with the SST Turbulence Model, 2003) but in Menter et al (The SST Turbulence Model with Improved Wall Treatment..., 2003) it is given (using OF variables) as gamma(F1)*G/nut. Furthermore, at, the formulation that OF uses is said to be an error in the original paper. This also affects the kOmegaSST model in OF.

Hopefully someone can confirm these errors or knows where in the literature these modifications have been made. Thanks.
RealSoutherner is offline   Reply With Quote

Old   August 17, 2012, 01:31
Senior Member
Join Date: Aug 2010
Location: State College, PA
Posts: 110
Rep Power: 9
Tarak is on a distinguished road

This was a bug in OF and has been rectified through mantis bugtracker with ID 0000575 (you can see the edited version in the git repository). This reduced production term was leading to a spurious reduction of turbulent viscosity.
Tarak is offline   Reply With Quote


Thread Tools
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 On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
No reactions in reactingFoam 2.1 OMN OpenFOAM Running, Solving & CFD 16 April 7, 2015 12:14
a little confused with the implementation of compressibleInterFoam su_junwei OpenFOAM Running, Solving & CFD 2 May 7, 2014 17:33
Implementation issues of fvSchemes / laplacianScheme, in particular gaussLaplacianSch thomek OpenFOAM Programming & Development 0 October 18, 2010 05:10
Request for help so as to debug MAC implementation HectorRedal Main CFD Forum 0 September 13, 2010 16:52
Please explain the implementation of species transport Eqn in reactingFoam kallipygian OpenFOAM Running, Solving & CFD 0 October 13, 2008 07:29

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