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

wall functions.-E value, epsilon calculation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 5, 2009, 05:40
Default wall functions.-E value, epsilon calculation
  #1
New Member
 
Join Date: Mar 2009
Posts: 20
Rep Power: 8
Xabi is on a distinguished road
Hi all!

I am reading through the wall functions code for OF and comparing it with fluent and I noticed some strange things:

1- For the value E instead of using 9.781 as fluent 9 is used.

2- When OF calculates the epsilon for the wall function uses the following formula:
epsilon_[faceCelli] +=Cmu75*pow(k_[faceCelli], 1.5)/(kappa_.value()*RASModel::y_[patchi][facei]);
In Fluent the same formula is used but what is bothering me is the faceCelli term. In fluent epsilon_ is calculated at the point (epsilonp), however, in OF the term epsilon_[faceCelli] seems to refer to the value at the wall.

Thank you in advance,
Xabi
Xabi is offline   Reply With Quote

Old   September 10, 2009, 13:57
Default
  #2
Member
 
Sven Schweikert
Join Date: Jun 2009
Posts: 38
Rep Power: 8
svens is on a distinguished road
Hey Xabi

I'm also interested how OF is dealing with the wall functions.

So far I'm not really deep in it but about the E thing I found out that they overwriting the E=9 to E=9.8 in the respective turbulence Model with this autoCreateEpsilon (-Nut etc.) function.

Cheers,
Sven
svens is offline   Reply With Quote

Old   September 11, 2009, 18:54
Default
  #3
Member
 
Sven Winkler
Join Date: May 2009
Posts: 70
Rep Power: 8
sven is on a distinguished road
Hey Xabi,
I am also interested in OpenFOAM deals with wall functions. In which file did you find the
equation epsilon_[faceCelli] +=Cmu75*pow(k_[faceCelli], 1.5)/(kappa_.value()*RASModel::y_[patchi][facei]); you mention in 2.? Thanks a lot!
sven is offline   Reply With Quote

Old   February 25, 2014, 12:24
Default
  #4
Member
 
Manan
Join Date: Oct 2013
Location: Göteborg
Posts: 37
Rep Power: 3
MaLa is on a distinguished road
Hi

On this note of the epsilonWallFunction, does anybody know the source of the formulation mentioned above by Sven? I wish to know the reasoning behind why epsilon @ the first node near wall is recalculated this way.

If anyone can direct me to the scientific article, which forms the basis of the epsilonWallFunction, I would be grateful. Thanks.
MaLa 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
Modify, G_k,k and epsilon values on wall functions Javier Larrondo FLUENT 0 September 15, 2008 21:34
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug unoder OpenFOAM Installation 11 January 30, 2008 21:30
Wall functions tutlhino OpenFOAM Pre-Processing 0 July 2, 2007 05:04
Warning 097- AB CD-adapco 6 November 15, 2004 05:41
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12


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