
[Sponsors] 
January 13, 2003, 10:17 
When I use the wall functions....!

#1 
Guest
Posts: n/a

Hi: I met a problem when I decided to use the wall functions. A project I'm working with needs to calculate the external airfoil aerodynamics(2D.),and I want to use the stantard ke model accompanying with the wall functions in the near wall region. I got the twolayer and three layer model for wall functions in the Amano's paper("Development of a Turbulence Near Wall Model and its Application to Separated and Reattached Flows" in 1984). However,I don't know how to calculate the turbulence shear stess(而) in the LogLaw layer and the outer region. In the visous sublayer,it assumes that 而 is zero .However, in the LogLaw layer ,而 is accquired by interpolation of 而n and 而w,where 而n is the turbulence shear stress of the north boundary of the first node,and 而w ,of the wall. And here comes the problem.Can I calculate the turbulence shear stress in the loglaw layer through the following relation: 而=Mut*(改Ui/改xj+改Uj/xi)? I wonder if those who have the experience of using wall functions would give me a hand. Thanks.


January 13, 2003, 22:56 
Re: Anybody helps me please!

#2 
Guest
Posts: n/a

Maybe my question is very easy to most of you,but it has puzzled me for a long time beacuse i am new for CFD. Please help me and give me some suggestions. Thanks a lot!!


January 14, 2003, 01:12 
Re: Anybody helps me please!

#3 
Guest
Posts: n/a

I think that u have to iterate there u have the vorticity omega lets say multply it with the laminar coeff of viscosity and get the shear from there calculare the k and epsilon then from the k and the epsilon calculate the turb viscosity and add it to the first equation.Now this time, u multiply the vorticity by the summation of the turbulent and laminar viscosity and get the shear. make a loop for this! I didn't do that when I faced such a problem rather I used the baldwinlomax there instead and mached with the kepsolin in the log layer and it was much better and efficient....however, its your choice. I hope I could help....Good luck...Others might have better ideas...look around


January 14, 2003, 09:40 
Re: Anybody helps me please!

#4 
Guest
Posts: n/a

I thought of it and I thing that if u're not well in the log layer( which should be the case) u shouldnt even consider the shear with the turbulent viscousity coefficient coz the turb vis is negligible in the sublayer and I think u were confused and made me confused at the begining (in my previous reply)


January 14, 2003, 11:32 
Re: Anybody helps me please!

#5 
Guest
Posts: n/a

Hi,Michel: Thanks a lot for your reply and consideration .Maybe I have not expressed my question clealy,so it makes u confused.I will emphasized my question on the two layer model. (1)In Amano's paper I mentioned in the topic message ,the near wall region is divided into two parts:a viscous sublayer region(0<y_plus<11) and an overlap layer(11<y_plus<400).Here the overlap layer is the socalled loglaw layer( in the paper it was called the fully turbulent region).In the visous sublayer,the turbluence shear stress is negligible and the paper assumes that 而=0.However ,in the loglaw layer region where the first point is set,the turbulence shear stress can not be neglected and it is assumed to vary linearly.So the 而 in the log law region are expressed as follows:而=而w+(而n而w)y/yn.It is a linear formulation,but how to get the 而n? Without 而n the formulation is useless and I am puzzled by 而n.I am not sure if 而n can be acquired from the definition of turbulence shear stress.(2)I must get the expression of 而 in all the two regions,because it is necessary in the calculation of the mean generation and destruction rate in the k and e equations for the first point.For 2D conditions,the generation rate of k can be written as P=而(改U/改y+改V/x),and for e, it can be expressed as (C1*e*P/k),then the mean generation rate of k and e can be acquired by integration from the wall to yn,where yn is the north boundary fo the first point. When the mean generation and destruction rate in the k and e equations are acquired,the k,e eqution for the first point can be sovled.(3)To include the mean generation and destruction rate into the k and e equations for the first point,it is said that it can reflect the influence of the molecule viscosity in the near wall region. (4)By the way I am also intrested in the hybrid method with BL model and ke model you have mentioned before. If I can't solve the problem of wall functions well,maybe I will have a try of that hybrid method. (5)I am not sure if I have expressed it clealy this time,Michel. (6)Expecting your response and thanks a lot lot lot. Maximus


January 14, 2003, 11:49 
Re: Anybody helps me please!

#6 
Guest
Posts: n/a

Hi,Michel: Sorry,perhaps I wrote a too long message,so only part of it can be showed here.I do not know if there is a limitation of the message length. I have sent a email to u. Thanks a lot. Maximus


January 14, 2003, 11:51 
Re: Anybody helps me please!

#7 
Guest
Posts: n/a

Hi,Michel: Thanks a lot for your reply and consideration .Maybe I have not expressed my question clealy,so it makes u confused.I will emphasized my question on the two layer model. (1)In Amano's paper I mentioned in the topic,the near wall region is divided into two parts:a viscous sublayer region and an overlap layer.Here the overlap layer is the socalled loglaw layer( in the paper it was called the fully turbulent region).In the visous sublayer,the turbluence shear stress is negligible and the paper assumes that 而=0.However ,in the loglaw layer region where the first point is set,the turbulence shear stress can not be neglected and it is assumed to vary linearly.So the 而 in the log law region are expressed as follows:而=而w+(而n而w)y/yn.It is a linear formulation,but how to get the 而n? Without 而n the formulation is useless and I am puzzled by 而n.I am not sure if 而n can be acquired from the definition of turbulence shear stress.(2)I must get the expression of 而 in all the two regions,because it is necessary in the calculation of the mean generation and destruction rate in the k and e equations for the first point.For 2D conditions,the generation rate of k can be written as P=而(改U/改y+改V/x),and for e, it can be expressed as (C1*e*P/k),then the mean generation rate of k and e can be acquired by integration from the wall to yn,where yn is the north boundary fo the first point. When the mean generation and destruction rate in the k and e equations are acquired,the k,e eqution for the first point can be sovled.(3)To include the mean generation and destruction rate into the k and e equations for the first point,it is said that it can reflect the influence of the molecule viscosity in the near wall region. (4)I am not sure if I have expressed it clealy this time,Michel. (5)By the way I am also intrested in the hybrid method with BL model and ke model you have mentioned before. If I can't solve the problem of wall functions,maybe I will have a try of that hybrid method. (6)Expecting your response and thanks a lot lot lot. Maximus


January 20, 2003, 10:35 
Re: When I use the wall functions....!

#8 
Guest
Posts: n/a

Hi there. Eventhough I've not worked with the Amano model, I used the earlier twolayer model the ChiengLaunder model (Num. Heat Transfer 1980, #3 pp 189207). First of all you should not expect stellar result from these models, the may give some improvements as compared with standard wall functions however at a price of being less stable. I would rather suggest you to use a lowReynolds model. Anyhow in the ChiengLaunder the first node is divided into two parts the laminar and the turbulent part. In the laminar part the turbulent viscosity is zero, while in the turbulent part it varies linearly: tau_turb=tau_wall+(tau_northtau_wall)*y/y_north with tau_wall computed from the loglaw: U\sqrt(k_v)\rho/tau_wall=ln(E*y*\sqrt(k_v)/\nu)/kappa E (=5.0) and kappa (=0.23) is changed compared to standard values. k_v is the turbulent kin energy at the edge of the viscous sublayer. k_v is computed from k/k_v=(y/y_v)**2 and y_v is set from y_v=20\nu/sqrt(k_v), where 20 is the viscous sublayer Reynolds number which is roughly equal to y+=11, (the bufferlayer where the linear and loglaw intersect). ChiengLaunder then perform the integrated production and dissipation term in kequation and set epsilon in the first node, see the paper or contact me. The values for the north face is given from interpolate values from node 1 and node 2. The turbulent shear stress could only be evaluate using mu_t*dU/dy in regions of small variation of dU/dy and mu_t. In the near wall region such a formulation is only approriate for a lowReynolds number model, where the variation is well discretized using several cells. Using wall functions the loglaw is applied instead. The formulation is ok in the loglaw (y+>30).
Best of luck.. Jonas 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Modelling Turbulent Wall Functions  victorwfreire  CFX  11  January 14, 2015 16:48 
LES with wall functions, wall boundary, Re 10000  subba  CFX  4  July 10, 2007 05:04 
Wall Functions  pierre  OpenFOAM Running, Solving & CFD  0  October 1, 2005 13:13 
Wall functions  iterative approach  agg  Main CFD Forum  1  May 13, 2005 12:25 
LES wall functions  barrat  Main CFD Forum  10  September 1, 2004 08:29 