CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Initial Turbulence Conditions

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

Like Tree4Likes
  • 2 Post By plucas
  • 2 Post By malaboss

Reply
 
LinkBack Thread Tools Display Modes
Old   April 2, 2013, 22:31
Default Initial Turbulence Conditions
  #1
Member
 
plucas's Avatar
 
Anonymous
Join Date: Dec 2011
Location: Everywhere
Posts: 75
Rep Power: 5
plucas is on a distinguished road
I am doing an automotive analysis using the k-omega SST turbulence model and I have a question on initial conditions.

The two values:
turbulentKE
turbulentOmega

So to calculate these first I calculated Reynolds Number:
Re = U * L / nu = 45 * 4.16 / 1.5e-5 = 12.5 E6

Then initial turbulence %:
I = 0.16 Re ^-1/8 = 0.0208

Now to turbulentKE:
k = 1.5 * (UI)^2 = 1.5 * (45 * 0.0208)^2 = 1.309

Now is the real question. I found 2 equations for epsilon with one being an approximation. After calculating both a get a big difference in numbers and wonder which one I should use to calculate omega.

epsilon ~ (k^1.5) / L = (1.309^1.5) / 4.16 = 0.360

epsilon = (cmu^.75) * (k^1.5) * (l^-1)
cmu = 0.09 (seems to be most commonly used)
l = 0.07 * L = 0.07 * 4.16 = 0.2912
epsilon = (0.09^0.75) * (1.309^1.5) * (0.2912^-1)
epsilon = 0.845

::::
::::
::::

epsilon = 0.360
epsilon = 0.845

So am I doing this right? If so, which epsilon value should I be using?

::::
::::
::::

turbulentOmega = k / epsilon
immortality and JR22 like this.
plucas is offline   Reply With Quote

Old   April 3, 2013, 21:55
Default
  #2
Member
 
plucas's Avatar
 
Anonymous
Join Date: Dec 2011
Location: Everywhere
Posts: 75
Rep Power: 5
plucas is on a distinguished road
Anybody know?
plucas is offline   Reply With Quote

Old   April 4, 2013, 11:12
Default
  #3
Senior Member
 
JR22's Avatar
 
Jose Rey
Join Date: Oct 2012
Posts: 125
Rep Power: 8
JR22 will become famous soon enough
I am also interested in some good input for this question. Although the openFoam forum has more traction, and this question is extremely relevant for openFoam, it might be appropriate for it to be in this forum:
http://www.cfd-online.com/Forums/main/

If you do repost the question there, let me know as I am very interested.
JR22 is offline   Reply With Quote

Old   April 4, 2013, 12:17
Default
  #4
Member
 
plucas's Avatar
 
Anonymous
Join Date: Dec 2011
Location: Everywhere
Posts: 75
Rep Power: 5
plucas is on a distinguished road
I thought about posting there but since this part of the forum (OpenFOAM) gets more foot traffic, I decided to post it here. I just made a new post however in the main.

Initial Turbulence Conditions
plucas is offline   Reply With Quote

Old   April 5, 2013, 17:45
Default
  #5
Senior Member
 
sivakumar selvaraju
Join Date: Mar 2009
Location: Cape Town - South Africa
Posts: 186
Rep Power: 8
sivakumar is on a distinguished road
Send a message via Skype™ to sivakumar
Hi,
I assume that your are doing external aerodynamics, in that case your problem is similar to motorBike case in OF.

According to OF training manual, you need to use the value corresponds to the following formula epsilon = (Cmu^0.75 K^1.5)/ L


Thanks,
Sivakumar

Last edited by sivakumar; April 8, 2013 at 04:55.
sivakumar is offline   Reply With Quote

Old   April 7, 2013, 08:39
Default
  #6
Member
 
Malik
Join Date: Dec 2012
Location: Austin, USA
Posts: 52
Rep Power: 4
malaboss is on a distinguished road
Quote:
Originally Posted by plucas View Post
epsilon ~ (k^1.5) / L = (1.309^1.5) / 4.16 = 0.360
Hi !

Could you tell us where you found this relation ? I've always seen the other one you talked about.


Quote:
Originally Posted by plucas View Post
l = 0.07 * L = 0.07 * 4.16 = 0.2912
For the turbulent length scale value, you choosed the same value for L as the characteristic length used for the Reynolds. If you are doing external aerodynamics, those two lengths could be different.

The characteristic length used for turbulent length scale is a length which depends on your wind tunnel's dimensions.

Could you tell us the meaning of the length L you used for the reynolds and the turbulent length scale ?

I guess it was a dimension of your obstacle. If it is right, then your calculations may not be right : the epsilon values you are calculating are the values of the freestream, before the fluid meets the obstacle : this is a boundary condition. You should have chosen a domain large enough so that this freestream value is not influenced by your obstacle. Then the epsilon value for the free stream should not be influenced by a dimension of your obstacle (the 4.16 value you showed us)


Have a nice day !
plucas and immortality like this.
malaboss is offline   Reply With Quote

Old   April 27, 2013, 21:05
Default
  #7
New Member
 
Tenglubao
Join Date: Feb 2013
Posts: 3
Rep Power: 4
philocfd is on a distinguished road
Quote:
Originally Posted by sivakumar View Post
Hi,
I assume that your are doing external aerodynamics, in that case your problem is similar to motorBike case in OF.

According to OF training manual, you need to use the value corresponds to the following formula epsilon = (Cmu^0.75 K^1.5)/ L


Thanks,
Sivakumar

In the user guide epsilon = (Cmu^0.75 K^1.5)/ l ,l is different from L
philocfd is offline   Reply With Quote

Reply

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
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
calculation stops after few time steps sivakumar OpenFOAM Running, Solving & CFD 7 March 17, 2013 07:37
Compressible Nozzle Flow sebastian OpenFOAM Running, Solving & CFD 13 October 21, 2010 05:25
Error while running rhoPisoFoam.. nileshjrane OpenFOAM Running, Solving & CFD 8 August 26, 2010 12:50
lift and drag on ship superstructures vaina74 OpenFOAM Running, Solving & CFD 3 June 8, 2010 12:30


All times are GMT -4. The time now is 18:44.