
[Sponsors] 
May 13, 2003, 13:16 
abt wall functions and kepsilon

#1 
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 ! 

May 14, 2003, 04:48 
Re: abt wall functions and kepsilon

#2 
Guest
Posts: n/a

The ke 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 nearwall 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?


May 14, 2003, 05:21 
Re: abt wall functions and kepsilon

#3 
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, gdensity * 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 underrelaxation 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 ! 

May 14, 2003, 05:59 
Re: abt wall functions and kepsilon

#4 
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.


May 14, 2003, 11:32 
Procedure here:

#5 
Guest
Posts: n/a

Please, take a look at section 7.3.4, page 345 in the book of Poinsot and Veynante.


May 14, 2003, 12:12 
Re: Procedure here:

#6 
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. 

May 20, 2003, 14:17 
You should buy it, it is good

#7 
Guest
Posts: n/a

You should buy it, it is good


Thread Tools  
Display Modes  


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 