
[Sponsors] 
March 5, 2013, 04:15 
Wall function's solving procedure.

#1 
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 742
Rep Power: 9 
Hi guys:
In Wall functions. there are lots of functions as follows.What confused me is how to solve this bunch of functions?I know the main idea of wall functions. but how to solve it? Thanks. 

March 5, 2013, 12:16 

#2 
Senior Member
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 196
Rep Power: 6 
Do you mean how to calculate y+ prior to the simulation of the flow, when you create the mesh?
Estimate the shear stress (by calculating the skin friction coefficient with one of the many formulas) then calculate Ut and then you can solve for y+.
__________________
Lefteris 

March 5, 2013, 19:52 

#3  
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 742
Rep Power: 9 
Quote:
Yes, Are you sure about estimating the shear stress? If I have the shear stress,I have the friction velocity,but if I want to know y+,I dont know k, how can I get k? And seems according to this formula,you can caculate y+ only by y and k. Thanks. 

March 6, 2013, 07:30 

#4 
Senior Member
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 196
Rep Power: 6 
Hello Forrest,
(Cμ^0.25)*(k^0.5) is the shear velocity. However, I usually use the formula y+=(ρ*Ut*y)/μ which is the same with the one you use but I estimate the shear velocity through the shear stress and therefore I don't have to use the k.
__________________
Lefteris 

March 6, 2013, 21:44 

#5  
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 742
Rep Power: 9 
Quote:
First I need to caclulate the skin friction coefficient by moody picture and other things, then I get shear stress. according to shear stress,I have shear velocity,and I have yplus. After I have yplus I can get k and U ,epsilon in the P node which can be treated as boundary conditions. but how can I get pressure in node P? I think if I have k,epsilon,U,P,I can integrate the KE model in the main flow and get the solution. 

March 7, 2013, 06:10 

#6 
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,296
Rep Power: 19 
You take all unknown values from the iteration before, or from the initialization. You can not know all of them at the same time, because then you would know the result...
__________________
The skeleton ran out of shampoo in the shower. 

March 22, 2013, 07:33 

#7  
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 742
Rep Power: 9 
Quote:
I attach my procedure, I dont know if its rite...still I am confused about it. To Rodriguez: I dont know where to start my iteration ... 

March 22, 2013, 09:33 

#8 
Senior Member
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 196
Rep Power: 6 
Hello Forrest,
You don't have to make an assumption for Uτ since you can compute it from the formula Uτ=sqrt(τw/ρ), where ρ is the density and τw the shear stress which can be computed from the formula τw=Cf*(0.5ρU^2). In order to estimate Cf, choose a formula from here: http://www.cfdonline.com/Wiki/Skin_...on_coefficient
__________________
Lefteris 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
alphaEqn.H in twoPhaseEulerFoam  cheng1988sjtu  OpenFOAM Bugs  15  May 1, 2016 16:12 
Extrusion with OpenFoam problem No. Iterations 0  Lord Kelvin  OpenFOAM Running, Solving & CFD  8  March 28, 2016 11:08 
Multiple floating objects  CKH  OpenFOAM Running, Solving & CFD  12  March 21, 2016 14:05 
pimpleFoam: turbulence>correct(); is not executed when using residualControl  hfs  OpenFOAM Running, Solving & CFD  3  October 29, 2013 09:35 
problem with Min/max rho  tH3f0rC3  OpenFOAM  7  February 23, 2013 06:37 