
[Sponsors] 
how to calculate the omega at inlet boundary in k omega sst 

LinkBack  Thread Tools  Search this Thread  Display Modes 
September 25, 2014, 11:46 
how to calculate the omega at inlet boundary in k omega sst

#1 
Member
Join Date: Dec 2013
Location: Newcastle
Posts: 54
Rep Power: 10 
Hi guys,
Recently, I am simulating flow past a stationary cylinder by using komegasst. However, I found that there are 3 omega values on inlet from different calculation method/tools. And the values are totally different. 1.omega=k^0.5/(l*Cmu^0.25) where k is the turbulence kinetic energy on inlet, l is the turbulence length scale and Cmu=0.09 with this condition, at re=10000 my omega on inlet boundary is 0.6388 2.Turbulence variables (k, ε, ω) from turbulence intensity (Tu), eddy viscosity ratio (μt/μ), freestream velocity (U∞) and kinematic viscosity (ν) nut/nu = k/ (nu*omega) which from the eddy viscosity=0.1, at this condition, my omega is 60 3. Turbulence variables (k, ε, ω) from turbulence intensity (Tu), lengthscale (TuL) and freestream velocity (U∞) by using the tools on http://www.cfdonline.com/Tools/turbulence.php that result shows my omega is 0.3499 Why are those omega value totally different, and by using different omega value at inlet boundary I got different Drag force coefficient Could you give me some recommendations about which formula I should use for calculating omega? reference value Freestream velocity U=1m/s,Turbulence intensity/level Tu=2%, Turbulence length scale TuL=0.07D=0.07m, Kinematic viscosity=1e04, Diameter of the cylinder=1m, Eddy viscosity ratio=0.1 Many thx, Scabbard 

September 25, 2014, 14:27 

#2 
Senior Member
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 14 
I'm interested too in this problem, same thing for the calculation of the epsilon that you can define based on the eddy viscosity ratio or kinetic energy and turbulence intensity.
check this page too LINK. I can't find the omega file in the sources, that's where we can find the solution. Ok i searched teh sources and the only reference with a law is this and can be found in the omegawallfunction file: Code:
231 // Set omega and G 232 forAll(nutw, faceI) 233 { 234 label cellI = patch.faceCells()[faceI]; 235 236 scalar w = cornerWeights[faceI]; 237 238 scalar omegaVis = 6.0*nuw[faceI]/(beta1_*sqr(y[faceI])); 239 240 scalar omegaLog = sqrt(k[cellI])/(Cmu25*kappa_*y[faceI]); 241 242 omega[cellI] += w*sqrt(sqr(omegaVis) + sqr(omegaLog)); 243 244 G[cellI] += 245 w 246 *(nutw[faceI] + nuw[faceI]) 247 *magGradUw[faceI] 248 *Cmu25*sqrt(k[cellI]) 249 /(kappa_*y[faceI]); 250 } 251 } Last edited by ArathoN; September 25, 2014 at 19:04. 

September 30, 2014, 13:06 

#3 
Senior Member
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 14 
ok i've found the relation in the kOmegaSST.H file from line 251 where it is defining epsilon and there you'll see that:
epsilon_ = betaStar_*k_*omega_ From here you can explicit omega. There is another expression that is valid only for compressible flows (there should be teh rho variable but i didn't see it inn the source file) that is: nut=k/omega. This is valid only away from the wall otherwise you have to apply the relation that a posted in the last post. 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Bubble column  Inlet boundary conditions  oj.bulmer  CFX  8  February 4, 2018 19:31 
Total pressure and mass flow boundary condition at inlet  bscphil  OpenFOAM PreProcessing  3  July 9, 2017 14:39 
GETVAR Error in Multiband Monte Carlo Radiation Simulation with Directional Source  silvan  CFX  3  June 16, 2014 09:49 
Question about heat transfer coefficient setting for CFX  Anna Tian  CFX  1  June 16, 2013 06:28 
An error has occurred in cfx5solve:  volo87  CFX  5  June 14, 2013 17:44 