CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Wall functions for near wall treatment : k - epsilon model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 8, 2009, 05:49
Default Wall functions for near wall treatment : k - epsilon model
  #1
Member
 
Phoevos
Join Date: Mar 2009
Posts: 84
Rep Power: 7
fivos is on a distinguished road
Hi everyone,

I am building my own CFD code and I want to use a k-e turbulence model. So I would like to ask for some more information from people with more experience, on the subject of the wall functions.
Form what I have already read, the continuity, momentum, k and epsilon equations are solved only for the nodes/cells away from the walls. On the other hand, nodes which are located near the wall (y_plus less than 300) have prescribed values for velocity, K and e and to be more specific :
U_plus= 1/0.41 * ln (y_plus) + 5.1 where U_plus=U_average / U_shear
and y_plus=U_shear*perpendicular_distance_from_wall/kinematic_viscosity
Also K= U_shear ^ 2 / sqrt (0.09) and
epsilon = U_shear ^3 / (0.41* delta )
where delta is an arbitarily set small distance, less than the cell size.

So my question is how do you exactly implement these wall functions and especially the log - law ? I mean : I can understand how to find the shear velocity (which on complex geometries has to be calculated as the directional derivative of the velocity towards the normal of the wall) and consequently calculate the y_plus, K and epsilon values. But what about the velocity? I mean that shear velocity is a scalar value, how can I find the three components of the velocity (u,v,w) for the area in which the log - law is valid ?

Also something else : Has anyone implemented a k-epsilon turbulence model on meshless techniques and especially on SPH?

Thanks in advance.
fivos is offline   Reply With Quote

Old   July 8, 2009, 08:47
Default
  #2
New Member
 
M. Nabi
Join Date: Jun 2009
Posts: 20
Rep Power: 7
mnabi is on a distinguished road
you need to find the parallel velocity to the boundary/solid. you need to find the normal vector to the solid, and then the tangent vector can be found. By mapping the celocity to the tangential vector, you can find the tangential component of velocity. It is the velocity that must be used for the log-law.

U_tan = U - U_norm = U - U.n
U_t=norm2(U_tan)
mnabi is offline   Reply With Quote

Old   July 9, 2009, 03:26
Default
  #3
Senior Member
 
N/A
Join Date: Mar 2009
Posts: 177
Rep Power: 7
harishg is on a distinguished road
The validity depends on your geometry. If you would encounter flow separation, consider low-Re k-e models or k-omega models.
harishg is offline   Reply With Quote

Old   July 13, 2009, 01:23
Default
  #4
Member
 
Phoevos
Join Date: Mar 2009
Posts: 84
Rep Power: 7
fivos is on a distinguished road
Thanks everyone I' ve managed to figure it out.
fivos is offline   Reply With Quote

Old   April 27, 2011, 07:44
Default
  #5
Member
 
jk
Join Date: Jun 2009
Posts: 64
Rep Power: 7
jyothishkumar is on a distinguished road
Dear Mr Fivos,

I have written a code for laminar flow with colocated grid arrangement. Now i am trying to implement k epsilon turbulence model (high Re) in that. I really dont know where i am going wrong in my code. It is implemented in the following way:

Initialising all the variables (utou, tou_w, etc..)

for the near wall grid node following source term is added to the u momentum equation

ro(k).value * (cmeu ^ 0.25) * (kp(k).value ^ 0.5) * delx / uplus(k).value

for the y momentum equation the velocity value (v velocity) is set to zero at the near wall grid node.

for the k equation

source2 = ro(k).value * (cmeu ^ 0.75) * (kp_old(k).value ^ 0.5) * uplus(k).value * 2 * delx

above source term will go into "ap" expression

su_ke = (2 * meut(k).value * sij * delx * dely) - (ro(k).value * ep(k).value * delx * dely) + (tou_w(k).value * u1(k).value * (2 * delx))

above source will go as "Su" in the ap_kp equation

above two source terms will be added to the near wall node equation

for the dissipation equation

ep(k).value = (cmeu ^ 0.75 * (kp(k).value ^ 1.5)) / (0.41 * 0.5 * dely)

found out directly for the epsilon value at the near grid node.

I have a doubt in the source term of my k and epsilon equation especially the production and dissipation terms. I am just using central difference term for finding the Sij terms (vorticity).

Please tell me your suggestion so that i can improve my code to run with the kepsilon model.

many thanks

jyothish
jyothishkumar 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
Wall Treatment of User defined Turbulence Model Yang Chung FLUENT 0 August 23, 2008 11:31
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12
the problem of the wall functions www_sun Phoenics 2 March 13, 2002 19:15
Wall functions Confused Main CFD Forum 1 August 14, 1998 09:31


All times are GMT -4. The time now is 05:36.