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

abt wall functions and k-epsilon

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 13, 2003, 13:16
Default abt wall functions and k-epsilon
  #1
amv
Guest
 
Posts: n/a
Hi

I am trying to write a code for turbulent flows. I used the standard wall functions and solved a lot of problems, but I am finding it hard to solve the cases where there is abrupt contraction or expansion .e.g. in backward facing step or forward facing step problems, I am finding difficulty of convergence and with only very small underrelaxation factors for k & epsilon, I could get convergence.

Does it mean that I have to go to two layer or three layer approaches (e.g. Amano / chiang) for such cases rather than the standard wall functions.

Thanking you !
  Reply With Quote

Old   May 14, 2003, 04:48
Default Re: abt wall functions and k-epsilon
  #2
mike malin
Guest
 
Posts: n/a
The k-e model plus wall functions should converge relatively easily for this basic case. If you use simple equilibrium wall functions, the values of k and e are fixed at the near-wall point, and so that should not be a problem if you are computing the friction velocity correctly from the log law. The wall friction force should be linearised in the momentum equation to facilitate convergence. If you need heavy relaxation on k and e, it may suggest that you are using inappropriate linearisation practices for the source terms in these equations. Does your code converge if you employ a uniform turbulent viscosity plus wall functions?

  Reply With Quote

Old   May 14, 2003, 05:21
Default Re: abt wall functions and k-epsilon
  #3
amv
Guest
 
Posts: n/a
hi, Thanks for your kind cooperation.

Now the code is behaving much better.

I am using

if y+ > 11.63,

g = (tau_wall ^2)/(density*kappa*u_tau*y_p)

and

epsilon=u_tau^3/(kappa*y_p)

and

if y+ < 11.63 then,

g-density * epsilon = 0 ( source term of tke equation)

epsilon=2 * nu * tke / y_p ^2 (boundary condition for epsilon equation)!

is the above formulation ok ?

I have to reduce under-relaxation factors if i use a finer grid now but for a coarser grid now i am able to go take under relax fact of 0.2 which was only .05 earlier !
  Reply With Quote

Old   May 14, 2003, 05:59
Default Re: abt wall functions and k-epsilon
  #4
mike malin
Guest
 
Posts: n/a
For y+ > 11.63, k=u_tau^2/0.3 with simple equilibrium wall functions. In my experience no convergence problems accrue if the same relations are employed when y+ < 11.63, wherein u_tau is computed for laminar flow. For me a linear relaxation factor of 0.2 or smaller seems rather low for securing convergence. A value of 0.4 or 0.5 ought to be sustainable.
  Reply With Quote

Old   May 14, 2003, 11:32
Default Procedure here:
  #5
RMG
Guest
 
Posts: n/a
Please, take a look at section 7.3.4, page 345 in the book of Poinsot and Veynante.
  Reply With Quote

Old   May 14, 2003, 12:12
Default Re: Procedure here:
  #6
AMV
Guest
 
Posts: n/a
Hi RMG, Can you tell me some details ? I think we don't have any book of these authors

Poinsot and Veynante

in our library.

Thanking you.
  Reply With Quote

Old   May 20, 2003, 14:17
Default You should buy it, it is good
  #7
RMG
Guest
 
Posts: n/a
You should buy it, it is good
  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
Natural convection in a closed domain STILL NEEDING help! Yr0gErG FLUENT 3 June 12, 2013 02:12
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
Wall functions for near wall treatment : k - epsilon model fivos Main CFD Forum 4 April 27, 2011 07:44
Modify, G_k,k and epsilon values on wall functions Javier Larrondo FLUENT 0 September 15, 2008 21:34
abt near wall turbulence wall functions amv Main CFD Forum 0 May 12, 2003 03:38


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