CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Calculating the Friction Velocity

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

Like Tree1Likes
  • 1 Post By ngj

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 19, 2009, 04:18
Question Calculating the Friction Velocity
  #1
New Member
 
Christopher Barry
Join Date: Aug 2009
Location: Switzerland
Posts: 18
Rep Power: 17
cbarry is on a distinguished road
Hi everyone,

I'm trying to implement new wall functions for the kOmegaSST model in OpenFOAM 1.5.

There is just one major hurlde I keep coming across. Does anyone know how I can calculate the friciton velcoity in OpenFOAM?

Best regards,

Chris
cbarry is offline   Reply With Quote

Old   August 19, 2009, 05:03
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,902
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Chris

It depends on whether or not you want to use low- or high Reynolds version of SST. If you are using low RE, then (nu() + nut()) U.snGrad() computed at the wall yields satisfactory results.
If on the other hand you are having high RE, then an assumption of a log-profile would give reasonable results.

Best regards,

Niels
ngj is online now   Reply With Quote

Old   August 19, 2009, 06:34
Default
  #3
New Member
 
Christopher Barry
Join Date: Aug 2009
Location: Switzerland
Posts: 18
Rep Power: 17
cbarry is on a distinguished road
Thank you very much Niels.
I am trying to implement a universal, non-linear wall function. For y+ > 30, i.e. the log layer, I want the wall function to take the following form:

omega = u*/ (sqrt(beta2) * kappa * y)

where u* = friction velocity

It is thus in effect a high Reynolds approach. How would you recommend I calculate the friction velocity in that case?

Best regards,

Chris
cbarry is offline   Reply With Quote

Old   August 19, 2009, 06:38
Default
  #4
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,902
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Chris

Then the current implementation of SST is what you want. As I understand it, it is doing exactly that.

Best regards,

Niels
ngj is online now   Reply With Quote

Old   August 19, 2009, 07:00
Default
  #5
New Member
 
Christopher Barry
Join Date: Aug 2009
Location: Switzerland
Posts: 18
Rep Power: 17
cbarry is on a distinguished road
Hi Niels,

as far as I can see the current implementation of SST calculates omega in the wall region as follows:

omega= sqrt (k) / (Cmu25 * kappa * y)

When I looked earlier on I couldn't see how the two expressions are equivalent. Are they actually equivalent then? If so, how?

In that case I'll try leaving that expression as it is and adding different functions for the buffer layer and the viscous sublayer.

Thank you for your quick help.

Best regards,

Chris
cbarry is offline   Reply With Quote

Old   August 19, 2009, 08:58
Default
  #6
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,902
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
The following is from [1]:

k = u_f^2 / sqrt(C_mu)

and if you combine this with your equation for omega, you obtain the implemented result.

Best regards,

Niels

[1]: Wall Function Boundary Condition Including Heat Transfer and Compressibility, RH Nichols and CC Nelson, AIAA, vol. 42, no. 6., june 2004.
albet likes this.
ngj is online now   Reply With Quote

Old   August 19, 2009, 09:34
Default
  #7
New Member
 
Christopher Barry
Join Date: Aug 2009
Location: Switzerland
Posts: 18
Rep Power: 17
cbarry is on a distinguished road
Thank you . It all makes sense now
cbarry is offline   Reply With Quote

Old   November 17, 2009, 12:24
Default
  #8
Senior Member
 
Anonymous
Join Date: Mar 2009
Posts: 110
Rep Power: 17
madad2005 is on a distinguished road
Chris,

Did you get any further with your blended wall function implementation?

Regards,

Adriano
madad2005 is offline   Reply With Quote

Reply

Tags
friction velocity, komegasst, shear velocity, wall function

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Velocity in Porous medium : HELP! HELP! HELP! Kali Sanjay Phoenics 0 November 6, 2006 07:10
How to update polyPatchbs localPoints liu OpenFOAM Running, Solving & CFD 6 December 30, 2005 18:27
Variables Definition in CFX Solver 5.6 R P CFX 2 October 26, 2004 03:13
friction velocity kei-tee Main CFD Forum 10 May 25, 2003 20:31
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11


All times are GMT -4. The time now is 16:56.