CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Wall function's solving procedure.

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 5, 2013, 04:15
Default Wall function's solving procedure.
  #1
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17
sharonyue is on a distinguished road
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.
Attached Images
File Type: jpg 1.jpg (8.1 KB, 18 views)
File Type: jpg 2.jpg (4.8 KB, 11 views)
File Type: jpg 3.jpg (3.1 KB, 10 views)
File Type: jpg 4.jpg (5.1 KB, 8 views)
sharonyue is offline   Reply With Quote

Old   March 5, 2013, 12:16
Default
  #2
Senior Member
 
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 337
Rep Power: 15
Aeronautics El. K. is on a distinguished road
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

Aeronautics El. K. is offline   Reply With Quote

Old   March 5, 2013, 19:52
Default
  #3
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17
sharonyue is on a distinguished road
Quote:
Originally Posted by Aeronautics El. K. View Post
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+.
Hi Lefteris,

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.
Attached Images
File Type: jpg 1.jpg (4.0 KB, 5 views)
sharonyue is offline   Reply With Quote

Old   March 6, 2013, 07:30
Default
  #4
Senior Member
 
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 337
Rep Power: 15
Aeronautics El. K. is on a distinguished road
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

Aeronautics El. K. is offline   Reply With Quote

Old   March 6, 2013, 21:44
Default
  #5
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17
sharonyue is on a distinguished road
Quote:
Originally Posted by Aeronautics El. K. View Post
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+.
Okay, Thanks, I got it.

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.
sharonyue is offline   Reply With Quote

Old   March 7, 2013, 06:10
Default
  #6
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
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.
RodriguezFatz is offline   Reply With Quote

Old   March 22, 2013, 07:33
Default
  #7
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17
sharonyue is on a distinguished road
Quote:
Originally Posted by Aeronautics El. K. View Post
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.
Hi,

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 ...
Attached Images
File Type: jpg 2.jpg (50.6 KB, 11 views)
sharonyue is offline   Reply With Quote

Old   March 22, 2013, 09:33
Default
  #8
Senior Member
 
Lefteris
Join Date: Oct 2011
Location: UK
Posts: 337
Rep Power: 15
Aeronautics El. K. is on a distinguished road
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

Aeronautics El. K. is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 22:47.