
[Sponsors] 
Production limiter only used for k  kOmegaSST  OpenFOAM 2.2.x 

LinkBack  Thread Tools  Search this Thread  Display Modes 
January 13, 2014, 12:07 
Production limiter only used for k  kOmegaSST  OpenFOAM 2.2.x

#1 
New Member
Join Date: Jul 2012
Posts: 1
Rep Power: 0 
I'm trying to understand which version of Menter's komega SST turbulence model OpenFOAM 2.2.2 is using.
Looking at the incompressible komega SST model, the production limiter is only added to the k equation. Code:
// Turbulent kinetic energy equation tmp<fvScalarMatrix> kEqn ( fvm::ddt(k_) + fvm::div(phi_, k_)  fvm::laplacian(DkEff(F1), k_) == min(G, c1_*betaStar_*k_*omega_) // Production limiter  fvm::Sp(betaStar_*omega_, k_) ); Code:
// Turbulent frequency equation tmp<fvScalarMatrix> omegaEqn ( fvm::ddt(omega_) + fvm::div(phi_, omega_)  fvm::laplacian(DomegaEff(F1), omega_) == gamma(F1)*S2 // No limiter  fvm::Sp(beta(F1)*omega_, omega_)  fvm::SuSp ( (F1  scalar(1))*CDkOmega/omega_, omega_ ) ); I found a bug report on the same topic, but I don't think it provides a answer to the question. http://openfoam.org/mantisbt/view.php?id=286 Thanks in advance. 

January 21, 2014, 15:40 

#2 
Member
Eric Robertson
Join Date: Jul 2012
Posts: 95
Rep Power: 13 
...........
Last edited by msuaeronautics; February 16, 2014 at 20:33. 

January 21, 2014, 21:19 

#3 
Member
Eric Robertson
Join Date: Jul 2012
Posts: 95
Rep Power: 13 
The NASA resource states that Menter required production limiters on both equations, but his actual 2003 paper only seems to call for the limiter on k. So, I suppose that the code only does what's been called for.
Last edited by msuaeronautics; February 16, 2014 at 20:33. 

January 22, 2014, 00:50 

#4 
Senior Member
Join Date: Jun 2011
Posts: 163
Rep Power: 13 
I have compared Fluent's and OFs' kwsst model. and found always the OF value of w is higher than Fluent' s w.(between 10 to 100 times).
but other fields (k U p) is very close to Fluent ones. I think this thread is the reason of higher w in OF 

January 22, 2014, 10:04 

#5 
Member
Eric Robertson
Join Date: Jul 2012
Posts: 95
Rep Power: 13 
Can you provide some data showing this?


January 30, 2014, 13:05 

#6 
Senior Member
Join Date: Jun 2011
Posts: 163
Rep Power: 13 
OK
I will provide some examples,as soon as possible 

February 20, 2014, 05:33 

#7 
New Member
Join Date: Apr 2010
Posts: 3
Rep Power: 15 
I'm also interested in why the production limiter is not used in the omega equation, or indeed if it should be.
Does anyone know where the correction referred to on the NASA website comes from? Could it be the following reference (I don't have access to this so cannot check, does anybody else?), or is it just knowledge that Chris Rumsey has? (who is curator of the NASA page). Menter, Florian R. "Review of the shearstress transport turbulence model experience from an industrial perspective." International Journal of Computational Fluid Dynamics 23.4 (2009): 305316. There is another bug report about the same issue and it seems the developer (Sergio) decided there is not a published paper to base the correction on: http://openfoam.org/mantisbt/view.php?id=486 

March 3, 2014, 17:11 

#8 
Member
Eric Robertson
Join Date: Jul 2012
Posts: 95
Rep Power: 13 
I have personally tried using a custom limiter on omega, which appears to keep omega in check, but introduces some nonphysical elements (like hotspots in TKE). Will have to do further testing.


September 21, 2015, 08:55 

#9  
Senior Member
Join Date: Mar 2015
Posts: 250
Rep Power: 11 
Hello all,
the omegaequation has changed with 2.3.x: Code:
// Turbulent frequency equation tmp<fvScalarMatrix> omegaEqn ( fvm::ddt(omega_) + fvm::div(phi_, omega_)  fvm::laplacian(DomegaEff(F1), omega_) == gamma(F1) *min(S2, (c1_/a1_)*betaStar_*omega_*max(a1_*omega_, b1_*F23()*sqrt(S2)))  fvm::Sp(beta(F1)*omega_, omega_)  fvm::SuSp ( (F1  scalar(1))*CDkOmega/omega_, omega_ ) ); Best regards, Kate Quote:


Tags 
komega sst, production limiter, turbulence models 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Confusion with magSqr in turbulence model kOmegaSST  OpenFOAM 2.2.x  A_Pete  OpenFOAM Programming & Development  4  September 18, 2014 04:52 
using production limiter for kw SST model  Madreed  FLUENT  1  September 2, 2013 17:06 
Same SimpleFOAM Case converges with openFOAM 2.1 but diverges with openFOAM 2.0.1  alsdia  OpenFOAM Running, Solving & CFD  3  October 22, 2012 11:25 
Wrong calculation of nut in the kOmegaSST turbulence model  FelixL  OpenFOAM Bugs  27  March 27, 2012 09:02 
kOmegaSST in openfoam 1.6  Gearb0x  OpenFOAM  2  March 3, 2010 06:02 