CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Main CFD Forum (https://www.cfd-online.com/Forums/main/)
-   -   When I use the wall functions....! (https://www.cfd-online.com/Forums/main/5629-when-i-use-wall-functions.html)

 maximus January 13, 2003 10:17

When I use the wall functions....!

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 k-e model accompanying with the wall functions in the near wall region. I got the two-layer 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 Log-Law layer and the outer region. In the visous sublayer,it assumes that 而 is zero .However, in the Log-Law 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 log-law 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.

 maximus January 13, 2003 22:56

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

 Michel January 14, 2003 01:12

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 baldwin-lomax there instead and mached with the k-epsolin 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 :)

 Michel January 14, 2003 09:40

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)

 maximus January 14, 2003 11:32

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 so-called log-law 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 log-law 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 B-L model and k-e 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

 maximus January 14, 2003 11:49

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

 maximus January 14, 2003 11:51