|
[Sponsors] |
March 5, 2013, 03:15 |
Wall function's solving procedure.
|
#1 |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 839
Rep Power: 17 |
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, 11:16 |
|
#2 |
Senior Member
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 337
Rep Power: 15 |
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, 18:52 |
|
#3 | |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 839
Rep Power: 17 |
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, 06:30 |
|
#4 |
Senior Member
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 337
Rep Power: 15 |
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, 20:44 |
|
#5 | |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 839
Rep Power: 17 |
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, 05:10 |
|
#6 |
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26 |
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, 06:33 |
|
#7 | |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 839
Rep Power: 17 |
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, 08:33 |
|
#8 |
Senior Member
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 337
Rep Power: 15 |
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.cfd-online.com/Wiki/Skin_...on_coefficient
__________________
Lefteris |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
problem with Min/max rho | tH3f0rC3 | OpenFOAM | 8 | July 31, 2019 09:48 |
Multiple floating objects | CKH | OpenFOAM Running, Solving & CFD | 14 | February 20, 2019 09:08 |
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 |
pimpleFoam: turbulence->correct(); is not executed when using residualControl | hfs | OpenFOAM Running, Solving & CFD | 3 | October 29, 2013 08:35 |